Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,281,524 C→T 5.0% pseudogene (454/555 nt) agaA → pseudogene, N‑acetylgalactosamine‑6‑phosphate deacetylase fragment; putative N‑acetylgalactosamine‑6‑phosphate deacetylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,281,5240CT 5.0% 258.0 / 2.2 80pseudogene (454/555 nt)agaApseudogene, N‑acetylgalactosamine‑6‑phosphate deacetylase fragment; putative N‑acetylgalactosamine‑6‑phosphate deacetylase
Reads supporting (aligned to +/- strand):  ref base C (34/41);  new base T (2/2);  total (37/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.17e-01

CGTACCGCGTCTGGTGGGCTGGCGGGCAGTACGCTGTCTGTTGATGCGGCAGTGCGCAATATGGTCGAGTTGACGGGCGTAACGCCTGCGGAAGCCATCCATATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCT  >  NC_000913/3281377‑3281672
                                                                                                                                                   |                                                                                                                                                    
cGTACCGCGTCTGGTGGGCTGGCGGGCAGTACGCTGTCTGTTGATGCGGCAGTGCGCAATATGGTCGAGTTGACGGGCGTAACGCCTGCGGAAGCCATCCATATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATTGTGTTc                                                                                                                                                      <  6:87233/148‑1 (MQ=255)
     cgcgTCTGGTGGGCTGGCGGGCAGTACGCTGTCTGTTGATGCGGCAGTGCGCAATATGGTCGAGTTGACGGGCGTAACGCCTGCGGAAGCCATCCATATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGAt                                                                                                                                                <  4:318181/149‑1 (MQ=255)
      gcgTCTGGTGGGCTGGCGGGCAGTACGCTGTCTGTTGATGCGGCAGTGCGCAATATGGTCGAGTTGACGGGCGTAACGCCTGCGGAAGCCATCCATATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATc                                                                                                                                               <  1:195448/149‑1 (MQ=255)
             gttggCTGGCGGTCAGTTCGCTGGATGTTGCTGCGGCAGTGCGAAATATTGTCGAGTTGACGGTCGTAACGGCTGCGGAAGCGTTCCATATTGCGTCTCTGCATCCGTTGCGTATGCTGGTTGTTGATTGTGTTCTGGGATCGCTTaaa                                                                                                                                        <  6:187525/146‑1 (MQ=255)
              tggGCTGGCGGGCAGTACGCTGTCTGTTGATGCGGCAGTGCGCAATATGGTCGAGTTGACGGGCGTAACGCCTGCGGAAGCCATCCATATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAAc                                                                                                                                       >  3:115387/1‑149 (MQ=255)
                 gCTGGCGGGCAGTACGCTGTCTGTTGATGCGGCAGTGCGCAATATGGTCGAGTTGACGGGCGTAACGCCTGCGGAAGCCATCCATATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCgg                                                                                                                                    <  7:308438/149‑1 (MQ=255)
                 gCTGGCGGGCAGTACGCTGTCTGTTGATGCGGCAGTGCGCAATATGGTCGAGTTGACGGGCGTAACGCCTGCGGAAGCCATCCATATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCgg                                                                                                                                    <  8:242920/149‑1 (MQ=255)
                  cTGGCGGTCAGTTCGCTTGATGTTGATGCGGCAGTGCGCAATATGGTCGCGTTGCTGGGCGTAACGCCTGCGGAAGCTCTCCATATGGCGTAGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCggg                                                                                                                                   <  8:397893/149‑1 (MQ=255)
                   tGGCGGGCAGTACGCTGTCTGTTGATGCGGCAGTGCGCAATATGGTCGAGTTGACGGGCGTAACGCCTGCGGAAGCCATCCATATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTc                                                                                                                                                      >  1:122614/1‑129 (MQ=255)
                   tGGCGGGCAGTACGCTGTCTGTTGATGCGGCAGTGCGCAATATGGTCGAGTTGACGGGCGTAACGCCTGCGGAAGCCATCCATATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTc                                                                                                                                                      <  2:122614/129‑1 (MQ=255)
                           aGTACGCTGTCTGTTGATGCGGCAGTGCGCAATATGGTCGAGTTGACGGGCGTAACGCCTGCGGAAGCCATCCATATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACgcgc                                                                                                                          <  8:53475/149‑1 (MQ=255)
                             tACGCTGTCTGTTGATGCGGCAGTGCGCAATATGGTCGAGTTGACGGGCGTAACGCCTGCGGAAGCCATCCATATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACgcgc                                                                                                                          >  7:229622/1‑147 (MQ=255)
                              aCGCTGGCTGTTGATGCGGCAGTGCGCAATATGGTCGAGTTGACGGGCGTAACGCCTGCGGAAGCCATCCATATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCa                                                                                                                        <  1:237005/148‑1 (MQ=255)
                                gCTGTCTGTTGATGCGGCAGTGCGCAATATGGTCGAGTTGACGGGCGTAACGCCTGCGGAAGCCATCCATATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTTTGGGATCGCTTAAACCGGGCAAACGCGCCAGcg                                                                                                                     >  8:44682/1‑149 (MQ=255)
                                gCTGTCTGTTGATGCGGCAGTGCGCAATATGGTCGAGTTGACGGGCGTAACGCCTGCGGAAGCCATCCATATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTTTGGGATCGCTTAAACCGGGCAAACGCGCCAGcg                                                                                                                     >  8:51210/1‑149 (MQ=255)
                                 ctgtctgtTGATGCGGCAGTGCGCAATATGGTCGAGTTGACGGGCGTAACGACTGCGGAAGCCATCCATATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGcgt                                                                                                                    <  2:178736/149‑1 (MQ=255)
                                       gtTGATGCGGCAGTGCGCAGTATGGTCGAGTTGACGGGCGTAACGCCTGCGGAAGCCATCCATATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTgc                                                                                                              >  8:78981/1‑149 (MQ=255)
                                       gtTGATGCGGCAGTGCGCAATATGGTCGAGTTGACGGGCGTAACGCCTGCGGAAGCCATCCATATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTgc                                                                                                              >  5:252620/1‑149 (MQ=255)
                                        ttGATGCGGCAGTGCGCAATATGGTCGAGTTGACGGGCGTAACGCCTGCGGAAGCCATCCATATGGCGTCGCTGCATCCGGCGCGATTGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGtt                                                                                                                <  8:146406/146‑1 (MQ=255)
                                        ttGATGCGGCAGTGCGCAATATGGTCGAGTTGACGGGCGTAACGCCTGCGGAAGCCATCCATATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTgcg                                                                                                             >  4:338261/1‑149 (MQ=255)
                                            tGCGGCAGTGCGCAATATGGTCGAGTTGACGGGCGTAACGCCTGCGGAAGCCATCCATATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTg                                                                                                          >  3:49045/1‑148 (MQ=255)
                                            tGCGGCAGTGCGCAATATGGTCGAGTTGACGGGCGTAACGCCTGCGGAAGCCATCCATATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGAGCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTgg                                                                                                         <  5:224660/149‑1 (MQ=255)
                                             gCGGCAGTGCGCAATATGGTCGAGTTGACGGGCGTAACGCCTGCGGAAGCCATCCATATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGa                                                                                                        >  7:32350/1‑149 (MQ=255)
                                               ggCAGTGCGCAATATGGTCGAGTTGACGGGCGTAACGCCTGCGGAAGCCATCCATATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGtt                                                                                                                >  7:185940/1‑139 (MQ=255)
                                               ggCAGTGCGCAATATGGTCGAGTTGACGGGCGTAACGCCTGCGGAAGCCATCCATATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGtt                                                                                                                <  8:185940/139‑1 (MQ=255)
                                                  aGTGCGCAATATGGTCGAGTTGACGGGCGTAACGCCTGCGGAAGCCATCCATATGGCGTCGCGGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTGCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCg                                                                                                   <  6:96501/149‑1 (MQ=255)
                                                    tGCGCAATATGGTCGAGTTGACGGGCGTAACGCCTGCGGAAGCCATCCATATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCggg                                                                                                 >  5:215762/1‑149 (MQ=255)
                                                    tGCGCAATATGGTCGAGTTGACGGGCGTAACGCCTGCGGAAGCCATCCATATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGGTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCggg                                                                                                 <  2:278050/149‑1 (MQ=255)
                                                      cgcAATATGGTCGAGTTGACGGGCGTAACGCCTGCGGAAGCCATCCATATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGc                                                                                                >  2:36734/1‑148 (MQ=255)
                                                      cgcAATATGGTCGAGTTGACGGGCGTAACGCCTGCGGAAGCCATCCATATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGc                                                                                                >  8:146774/1‑148 (MQ=255)
                                                         aaTATGGTCGAGTTGACGGGCGTAACGCCTGAGGAAGCCATCCATATTGCGTAGATGCATCCGTCGCGAATGCTGGGTGTTGATGGTGTTCTGTGATCGCGTAAACCGGGCAAACTAGACAGCGTCGTTGCGCTGGATAGAGGGCTaa                                                                                             <  2:99479/148‑2 (MQ=255)
                                                          atatGGTCGAGTTGACGGGCGTAACGCCTGCGGAAGCCATCCATATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACAt                                                                                           <  6:215762/149‑1 (MQ=255)
                                                             tGGTCGAGTTGACGGGCGTAACGCCTGCGGAAGCCATCCATATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTTTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACAtgt                                                                                         <  7:44682/148‑1 (MQ=255)
                                                             tGGTCGAGTTGACGGGCGTAACGCCTGCGGAAGCCATCCATATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTTTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACAtgt                                                                                         <  7:51210/148‑1 (MQ=255)
                                                             tGGTCGAGTTGACGGGCGTAACGCCTGCGGAAGCCATCCATATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACAtgt                                                                                         <  4:49045/148‑1 (MQ=255)
                                                                  gAGGTGACGGGCGTAACGCCTGCGGAAGCCATCCATATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGcaac                                                                                    <  1:157249/148‑1 (MQ=255)
                                                                   aGTTGACGGGCGTAACGCCTGCGGAAGCCATCCATATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGcaaca                                                                                   >  7:400568/1‑148 (MQ=255)
                                                                     ttGACGGGCGTAACGCCTGCGGAAGCCATCCATATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAAt                                                                                <  6:310762/149‑1 (MQ=255)
                                                                          gggCGTAACGCCTGCGGAAGCCATCCATATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGa                                                                           >  3:98748/1‑149 (MQ=255)
                                                                              gTAACGCCTGCGGAAGCCATCCATATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCa                                                                       <  4:98748/149‑1 (MQ=255)
                                                                               tAACGCCTGAGGAAGCTCTGCATATCTGGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGGTCTGGGATCGCTTACACCGGGCAAACGCGCCAGCGGCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAg                                                                      <  4:46003/149‑1 (MQ=255)
                                                                                 aCGCCTGCGGAAGCCATCCATATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCa                                                                       <  6:192828/146‑1 (MQ=255)
                                                                                 aCGCCTGCGGAAGCCATCCATATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCa                                                                       >  5:192828/1‑146 (MQ=255)
                                                                                 aCGCCTGCGGAAGCCATCCATATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAgg                                                                     <  8:49977/148‑1 (MQ=255)
                                                                                 aCGCCTGCGGAAGCCATCCATATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAgg                                                                     <  6:279753/148‑1 (MQ=255)
                                                                                  cGCCTGCGGATGCCATCCATATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGt                                                                   <  2:211019/149‑1 (MQ=255)
                                                                                       gCGGAAGCCATCCATATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAAt                                                               >  2:131058/1‑148 (MQ=255)
                                                                                        cGGAAGCCATCCATATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTa                                                             <  8:400568/149‑1 (MQ=255)
                                                                                          gAAGCCATCCATATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGc                                                                                                                                             <  5:334235/67‑1 (MQ=255)
                                                                                          gAAGCCATCCATATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGc                                                                                                                                             >  6:334235/1‑67 (MQ=255)
                                                                                            aGCCATCCATATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCtt                                                         <  2:354467/149‑1 (MQ=255)
                                                                                               catccatATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGa                                                                           <  4:62409/128‑1 (MQ=255)
                                                                                                atccatATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGtt                                                     >  5:298356/1‑149 (MQ=255)
                                                                                                atccatATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGtt                                                     >  7:40778/1‑149 (MQ=255)
                                                                                                 tccatATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGttt                                                    >  3:175066/1‑149 (MQ=255)
                                                                                                        ggCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGt                                             >  6:60441/1‑149 (MQ=255)
                                                                                                                   aTCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGa                                                 >  6:186746/1‑134 (MQ=255)
                                                                                                                   aTCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGa                                                 <  5:186746/134‑1 (MQ=255)
                                                                                                                   aTCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTa                                  >  8:387603/1‑149 (MQ=255)
                                                                                                                   aTCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTa                                  <  7:387603/149‑1 (MQ=255)
                                                                                                                      cGGCGCGAATGCTGGGTGTTTATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTAAATTTGCAACAAATCTGGATTCTGGGGCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTAtt                                <  8:229622/148‑1 (MQ=255)
                                                                                                                      cGGCGCGAATGCTGGGTGTTGATGGTGTTCTTGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTAtt                                <  1:28653/148‑1 (MQ=255)
                                                                                                                       ggCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTgg                              >  6:77874/1‑149 (MQ=255)
                                                                                                                       ggCGCGAATGCTGGGTGTTGATGGTGTGCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTgg                              <  2:279522/149‑1 (MQ=255)
                                                                                                                         cgcgAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTGCGTTTTGATAGTTTGCTCCTTTATTGGGc                            <  1:36734/149‑1 (MQ=255)
                                                                                                                          gcgAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGcc                           >  8:2748/1‑149 (MQ=255)
                                                                                                                          gcgAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGcc                           >  6:367392/1‑149 (MQ=255)
                                                                                                                          gcgAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTCGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGcc                           >  1:323460/1‑149 (MQ=255)
                                                                                                                          gcgAATGCTGGGAGTACATGGTCTTCACGGATCGCTTAAACCGCGCAAACGCGCCAGCGTCGTTGCCCTCGATAGCGGCCTCCATCTGCAACAAATCTGGATTCAGGGTCAATTACCTTCCGTATGATAGTTTGCTCCTTTATTgcgc                            >  8:125325/1‑148 (MQ=255)
                                                                                                                                  tGGGTGTTGATGGTGTTATGGTATCGATTAAACCGGGCAAACGCGCCAGCTTCTTTGCGCTGGATAGATGGCTACATGTTCAACAAATATGGATTAAGTGTCAATTAGCTTCGTTATGATAGTTTGATCCTTTATTGGGCCTTCaagtc                   >  2:207469/1‑145 (MQ=255)
                                                                                                                                    ggTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCTTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTcc                  <  4:115387/148‑1 (MQ=255)
                                                                                                                                      tgtTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTccccc               <  3:338261/149‑1 (MQ=255)
                                                                                                                                        ttGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCAc                      >  7:323268/1‑140 (MQ=255)
                                                                                                                                        ttGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCAc                      <  8:323268/140‑1 (MQ=255)
                                                                                                                                        ttGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCg              >  5:128975/1‑148 (MQ=255)
                                                                                                                                          gATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTAtt                                <  2:72237/128‑1 (MQ=255)
                                                                                                                                          gATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTAtt                                >  1:72237/1‑128 (MQ=255)
                                                                                                                                          gATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGt             >  6:115415/1‑147 (MQ=255)
                                                                                                                                           aTGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCt                                                                              <  5:103742/81‑1 (MQ=255)
                                                                                                                                           aTGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAg          >  7:183267/1‑149 (MQ=255)
                                                                                                                                           aTGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGCATAGCGGGCTACATGTGCAACAAATCt                                                                              >  6:103742/1‑81 (MQ=255)
                                                                                                                                                 ttCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCcttt    <  7:2748/149‑1 (MQ=255)
                                                                                                                                                 ttCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCctt     >  4:361061/1‑148 (MQ=255)
                                                                                                                                                   cTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCctttct  >  6:33190/1‑149 (MQ=255)
                                                                                                                                                   cTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCctttct  >  7:302568/1‑149 (MQ=255)
                                                                                                                                                   cTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCctttc   <  5:367392/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CGTACCGCGTCTGGTGGGCTGGCGGGCAGTACGCTGTCTGTTGATGCGGCAGTGCGCAATATGGTCGAGTTGACGGGCGTAACGCCTGCGGAAGCCATCCATATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCT  >  NC_000913/3281377‑3281672

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑