Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913370,8760CG40.9% 120.9 / 31.4 66S200R (AGC→AGG) mhpB2,3‑dihydroxyphenylpropionate 1,2‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base C (33/6);  new base G (0/27);  total (33/33)
Fisher's exact test for biased strand distribution p-value = 9.04e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.28e-01

TACCCGCATGTTGGGTGAAGCCATTGGACGTTTCACCAGCACTCTCAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTT‑CCCAT‑CAGCCGCCGGT‑GCCC‑GAA‑CTGGCGAAA‑GCCGATGCCCAT‑ATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACG  >  NC_000913/370726‑371022
                                                                                                                                                             |                                                                                                                                                  
tACCCGCATGTTGGGTGAAGCCATTGGACGTTTCACCAGCACTCTCAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTT‑CCCAT‑CAGCCGCCGGT‑GCCC‑GAA‑CTGGCGAAA‑GCCGATGCCCAT‑ATGCGCGACCGTCTGTTGGGGAGc                                                                                                                                                    >  2:519012/1‑151 (MQ=255)
    cGCATGTTGGGTGAAGCCATTGGACGTTTCACCAGCACTCTCAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTT‑CCCAT‑CAGCCGCCGGT‑GCCC‑GAA‑CTGGCGAAA‑GCCGATGCCCAT‑ATGCGCGACCGTCTGTTGGGGAGCggg                                                                                                                                                 >  1:527084/1‑150 (MQ=255)
        tGTTGGGTGAAGCCATTGGACGTTTCACCAGCACTCTCAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTT‑CCCAT‑CAGCCGCCGGT‑GCCC‑GAA‑CTGGCGAAA‑GCCGATGCCCAT‑ATGCGCGACCGTCTGTTGGGGAGCGGGaa                                                                                                                                               >  2:868593/1‑148 (MQ=255)
           tGGGTGAAGCCATTGGACGTTTCACCAGCACTCTCAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTT‑CCCAT‑CAGCCGCCGGT‑GCCC‑GAA‑CTGGCGAAA‑GCCGATGCCCAT‑ATGCGCGACCGTCTGTTGGGGAGCGGGAAAGAtt                                                                                                                                          <  2:127024/150‑1 (MQ=255)
           tGGGTGAAGCCATTGGACGTTTCACCAGCACTCTCAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTT‑CCCAT‑CAGCCGCCGGT‑GCCC‑GAA‑CTGGCGAAA‑GCCGATGCCCAT‑ATGCGCGACCGTCTGTTGGGGAGCGGGAAAGAtt                                                                                                                                          <  1:954702/150‑1 (MQ=255)
                                tCACCAGCACTCTCAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTT‑CCCAT‑CAGCCGCCGGT‑GCCC‑GAA‑CTGGCGAAA‑GCCGATGCCCAT‑ATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAAtt                                                                                                                    >  2:1436416/1‑151 (MQ=255)
                                tCACCAGCACTCTCAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTT‑CCCAT‑CAGCCGCCGGT‑GCCC‑GAA‑CTGGCGAAA‑GCCGATGCCCAT‑ATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAAtt                                                                                                                    >  2:741644/1‑151 (MQ=255)
                                 caccacCCCTCAAAAAAACCGGGTGCTGTTCC‑GGGTTTCGGGGGGTTTTCCCCAT‑CGGCCGCGGGG‑CCCC‑GAG‑CGGGGAAAA‑GCCGATCCCCAT‑ATGCGCGCCCGTTTGTGGGGGAGGGGGAAATATTTCCCCGCCAGTGAGCGCGAATTg                                                                                                                   <  1:1344212/151‑1 (MQ=255)
                                  aCCAGCACTCTCAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTT‑CCCAT‑CAGCCGCCGGT‑GCCC‑GAA‑CTGGCGAAA‑GCCGATGCCCAT‑ATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGc                                                                                                                  >  2:346977/1‑151 (MQ=255)
                                  aCCAGCACTCTCAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTT‑CCCAT‑CAGCCGCCGGT‑GCCC‑GAA‑CTGGCGAAA‑GCCGATGCCCAT‑ATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGc                                                                                                                  >  1:966721/1‑151 (MQ=255)
                                   ccacccctctAAATAACCGGGTGTTTTTCCTGGTTTCGGGTGGGTTTT‑CCCAT‑CGCCGGCCGGT‑GCCC‑GAA‑CGGGAGAAA‑GCCGATCCCAAT‑CTGCCCGACCTTTTTTGGGGGGGGGGAAAAGATTTCCCCGCCAGTGAGCGCGAATTg                                                                                                                   <  1:1476136/143‑1 (MQ=255)
                                   ccAGCACTCTCAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTT‑CCCAT‑CAGCCGCCGGT‑GCCC‑GAA‑CTGGCGAAA‑GCCGATGCCCAT‑ATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCg                                                                                                                 >  1:1402311/1‑151 (MQ=255)
                                     aGCACTCTCAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTT‑CCCAT‑CAGCCGCCGGT‑GCCC‑GAA‑CTGGCGAAA‑GCCGATGCCCAT‑ATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTc                                                                                                               >  1:540840/1‑151 (MQ=255)
                                           caaaacaaccgcgTGTTGTGGCGGGGTTCCGG‑GGGCTTT‑CCCCTGCGGCGGCGGGC‑GCCC‑GAA‑GGGGAAAGG‑GCGGAGCCCCAT‑AGGCGCGCCCTTTGGTGGGGGGGGGGGAAATATTGTCCCGCCAGTGAGCGCGAATTGCGTCAGCaa                                                                                                          <  2:1202074/141‑1 (MQ=255)
                                            tcAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTT‑CCCAT‑CAGCCGCCGGT‑GCCC‑GAA‑CTGGCGAAA‑GCCGATGCCCAT‑ATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACg                                                                                                        >  2:1017974/1‑151 (MQ=255)
                                            tcAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTT‑CCCAT‑CAGCCGCCGGT‑GCCC‑GAA‑CTGGCGAAA‑GCCGATGCCCAT‑ATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACg                                                                                                        >  1:119453/1‑151 (MQ=255)
                                              aaTAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTT‑CCCAT‑CAGCCGCCGGT‑GCCC‑GAA‑CTGGCGAAA‑GCCGATGCCCAT‑ATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACggg                                                                                                      >  1:859743/1‑151 (MQ=255)
                                                 aaaCGCGTGCTGTTCCTGGGTTCCGGTGGGCTTT‑CCCAT‑CAGCCGCCGGT‑GCCC‑GAA‑CTGGCGAAA‑GCCGATGCCCAT‑ATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTg                                                                                                    >  2:1390003/1‑150 (MQ=255)
                                                   aCGCGTGCTGTTCCTGGGTTCCGGTGGGCTTT‑CCCAT‑CAGCCGCCGGT‑GCCC‑GAA‑CTGGCGAAA‑GCCGATGCCCAT‑ATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGAtt                                                                                                 >  2:1224844/1‑151 (MQ=255)
                                                    cctctgcgtttCCTGGTTTGCGGGGGTCTTC‑CCCTT‑CCGCCGCC‑GT‑TCCC‑GAACCTGGAG‑AA‑GCCCCGGCCCATAACGCGCGCCCGTTTGTGGGGGGGGGGAAAATATTTTCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGAt                                                                                                  <  1:541900/140‑1 (MQ=255)
                                                          ttGTTCCTGGGTCCGGGTGGGTTTC‑CCAAC‑CAGCCGCGGGG‑GCCC‑GAA‑CGGGCGAAA‑GCCGACGCCCAT‑AGGCGCGCCCGTTTGTGGGGGGGCGGAAAAATTTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGAtt                                                                                                 <  1:310194/143‑1 (MQ=255)
                                                               ccTGGGTTCCGGTGGGCTTT‑CCCAT‑CAGCCGCCGGT‑GCCC‑GAA‑CTGGCGAAA‑GCCGATGCCCAT‑ATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTg                                                                                       >  1:881945/1‑149 (MQ=255)
                                                                cTGGGTTCCGGTGGGCTTT‑CCCAT‑CAGCCGCCGGT‑GCCC‑GAA‑CTGGCGAAA‑GCCGATGCCCAT‑ATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTg                                                                                       >  1:25197/1‑148 (MQ=255)
                                                                cTGGGCTCCGGGTGGCCTT‑CCCCA‑CCGCCCCGGGAGGCCC‑GAG‑CTGGAGAAA‑GCTGATGCCCTT‑GGGCCCGACCGTTTGTGGGGGGGGGAAAAATTTTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTg                                                                                       <  2:944735/149‑1 (MQ=255)
                                                                  gggTCCGGGGGGGTTTC‑CCCAC‑CGGCCGCCGGG‑GCCG‑GAC‑GTGGGAAAG‑GCGGGGGCCGAT‑CGGCGCGCCGGGTTGTGGGGGGGGGGAAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTgag                                                                                     <  2:173811/148‑1 (MQ=255)
                                                                   ggTTCCGGTGGGCTTT‑CCCAT‑CAGCCGCCGGT‑GCCC‑GAA‑CTGGCGAAA‑GCCGATGCCCAT‑ATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGt                                                                                 >  2:717683/1‑151 (MQ=255)
                                                                      tCCGGTGGGCTTT‑CCCAT‑CAGCCGCCGGT‑GCCC‑GAA‑CTGGCGAAA‑GCCGATGCCCAT‑ATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTttg                                                                              >  2:324665/1‑151 (MQ=255)
                                                                        cGGTGGGCTTT‑CCCAT‑CAGCCGCCGGT‑GCCC‑GAA‑CTGGCGAAA‑GCCGATGCCCAT‑ATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTttgtt                                                                            >  1:1301322/1‑151 (MQ=255)
                                                                             ggCTTC‑CCCAG‑CGGCCGCGGGT‑CCCC‑AAC‑CTGGCAAAA‑GCCGATGCCCTT‑GTGCGCGACCGTGTTTGGGGGGGGGGAAAAGATTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGa                                                                       <  2:632778/151‑1 (MQ=255)
                                                                                ttt‑CCCAT‑CAGCCGCCGGT‑GCCC‑GAA‑CTGGCGAAA‑GCCGATGCCCAT‑ATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATc                                                                     >  1:193590/1‑150 (MQ=255)
                                                                                ttt‑CCCAT‑CAGCCGCCGGT‑GCCC‑GAA‑CTGGCGAAA‑GCCGATGCCCAT‑ATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATc                                                                     >  1:881797/1‑150 (MQ=255)
                                                                                   ccccAT‑CGCCCGCGGGT‑CCCA‑AAA‑GTGGAAAAC‑GCCGATGCCAAT‑TTGCGCGCCCGTTTGTGGGGGGGGGGAAAATATTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATc                                                                     <  1:1217422/147‑1 (MQ=255)
                                                                                    cccAT‑CAGCCGCCGGT‑GCCC‑GAA‑CTGGCGAAA‑GCCGATGCCCAT‑ATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCag                                                                   >  1:1211945/1‑149 (MQ=255)
                                                                                     ccAT‑CAGCCGCCGGT‑GCCC‑GAA‑CTGGCGAAA‑GCCGATGCCCAT‑ATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCag                                                                   >  2:866915/1‑148 (MQ=255)
                                                                                      cccaccgcccgcgGG‑GCCCGGAA‑CGGGCAAAA‑GCCGATCCCAAT‑ATGCGCGCCCGTTTGTGGGGGGGGGGAAAAGTTTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCaga                                                                  <  1:1139037/139‑1 (MQ=255)
                                                                                         ccAGCCGCCGGT‑GCCC‑AAC‑CTGGCAAAA‑CCCGACCCCCAT‑ATGCCCGACCGTTTGTGGGGGAGGGGAAAAATTTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAgaga                                                                <  1:1055662/147‑1 (MQ=255)
                                                                                             ccgccgGT‑GCCC‑GAA‑CTGGAGAAA‑GCCGATGCCAAT‑ATGCGCGACCGTTTGTGGGGGGGGGGGAAAGTTTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCt                                                           <  1:487623/149‑1 (MQ=255)
                                                                                             ccgccgGG‑GCCA‑G‑A‑CGGGGGAAAGGCCGATCCCCAT‑ATGCGCGCCCGTCTGTGGGGGGGGGGAAAAGATTCCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGc                                                         <  1:129236/151‑1 (MQ=255)
                                                                                                 cgGT‑GCCC‑GAA‑CTGGCGAAA‑GCCGATGCCCAT‑ATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATcc                                                     >  2:978431/1‑151 (MQ=255)
                                                                                                         c‑GAA‑CTGGCGAAA‑GCCGATGCCCAT‑ATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAAc                                              >  1:832206/1‑151 (MQ=255)
                                                                                                            aa‑CTGGCGAAA‑GCCGATGCCCAT‑ATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAAccc                                            >  2:311932/1‑151 (MQ=255)
                                                                                                               ggggagAAA‑GCCGACCCCCAT‑GGGCGCGACGTTTTGTGGGGGAGCGGGAAAGATTCCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGa                                          <  2:369900/146‑1 (MQ=255)
                                                                                                                 ggaaaaa‑GCCGATCCCCAT‑ATGCGCGCCCGTTTGTGGGGGGGGGGAAAAGTTTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGAt                                         <  1:868609/146‑1 (MQ=255)
                                                                                                                     aaa‑GCCGATGCCCAT‑ATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTggg                                    >  1:261316/1‑151 (MQ=255)
                                                                                                                     aaa‑GCCGATCCCCAT‑GGGCGCGCCCGTCTGTGGGGGGGCGGGAAAGTTTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGAttt                                       <  1:1386329/148‑1 (MQ=255)
                                                                                                                      aaaccggtgccccAT‑AGGCGCCCCCGTTTGTGGGGGAGGGGGAAAGTTTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTg                                      <  1:358987/139‑1 (MQ=255)
                                                                                                                       acgcggATGCCCAT‑ATGGGGGACGGTTTGTGGGGGGGGGGGAAAGTTTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTg                                      <  1:710051/143‑1 (MQ=255)
                                                                                                                          ccgttccctat‑agGCGCGACCGTTGGTGGGGGGGGGGGAAAGTTTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATaa                                <  2:980281/142‑1 (MQ=255)
                                                                                                                          ccgacccccAT‑ACGCGCGCCGGTTTGGGGGGGAGGGGAAAATTTTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATaa                                <  2:649518/145‑1 (MQ=255)
                                                                                                                          ccGTTGCCGAT‑ATGCGCGCCGTTCTGTGGGGGTGCGGAAAAGTTTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTg                                      <  1:205912/145‑1 (MQ=255)
                                                                                                                          ccGATGCCCAT‑ATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATa                                 >  2:91503/1‑150 (MQ=255)
                                                                                                                           ccgttccccagttGCGCGCCCGTTTGGGGGGGAGGGAAAAAGTTTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATa                                 <  2:643726/138‑1 (MQ=255)
                                                                                                                            gATGCCCAT‑ATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAAcc                              >  2:719699/1‑151 (MQ=255)
                                                                                                                              atcccaatgtGCGCGACCGTTTGTGGGGGGGGGGAAAATTTTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCa                             <  1:467762/142‑1 (MQ=255)
                                                                                                                                cccAT‑AGGCGCGCCCGTTTGTGGGGGGGGGGAAAAGTTTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCa                             <  1:1032094/148‑1 (MQ=255)
                                                                                                                                   at‑atGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCAt                       >  1:1120703/1‑151 (MQ=255)
                                                                                                                                   at‑agGCGCGCCCGTCTGTGGGGGGGGGGAAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCAt                       <  1:1143929/151‑1 (MQ=255)
                                                                                                                                        gcgcgcCGGTGTGTGGGGGGGGGGAAAAGTTTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCa                        <  2:1489677/146‑1 (MQ=255)
                                                                                                                                         cggacccGTGTGTGGGGGGGGGGAAAAGTTTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACtt                  <  1:596370/146‑1 (MQ=255)
                                                                                                                                          gcgccggtaggtgggggAGGGGGAAATATTCACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGa                     <  1:1015299/134‑1 (MQ=255)
                                                                                                                                           cgcccGTTTGTGGGGGGGGGGAAAATTTTCCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATg                      <  1:810250/142‑1 (MQ=255)
                                                                                                                                                 tCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTgg            >  1:1212963/1‑149 (MQ=255)
                                                                                                                                                 tCTGTGGGGGAGGGGAAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGc               <  1:1305664/146‑1 (MQ=255)
                                                                                                                                                  cTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGc         >  1:458743/1‑151 (MQ=255)
                                                                                                                                                   tGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGc         >  2:379041/1‑150 (MQ=255)
                                                                                                                                                    gtgggggggggGAAAAGTTTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGc         <  1:130635/146‑1 (MQ=255)
                                                                                                                                                        gggAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGAc   >  1:890927/1‑151 (MQ=255)
                                                                                                                                                         ggAGGGGGAAAGATTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACg  <  1:231709/151‑1 (MQ=255)
                                                                                                                                                             |                                                                                                                                                  
TACCCGCATGTTGGGTGAAGCCATTGGACGTTTCACCAGCACTCTCAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTT‑CCCAT‑CAGCCGCCGGT‑GCCC‑GAA‑CTGGCGAAA‑GCCGATGCCCAT‑ATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACG  >  NC_000913/370726‑371022

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