Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,421,4430AC40.9% 213.5 / 21.0 115Y40S (TAC→TCC) folXD‑erythro‑7,8‑dihydroneopterin triphosphate 2'‑epimerase and dihydroneopterin aldolase
Reads supporting (aligned to +/- strand):  ref base A (60/8);  new base C (0/47);  total (60/55)
Fisher's exact test for biased strand distribution p-value = 2.67e-24
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 115 total aligned reads displayed.

TCATCCTAAGTAAAACAACGGAGAACCTGCAATGGC‑ACAACCTGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACA‑ACCGTC‑AGGATATTGTTATCAATGTGACGATCCA‑CTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAAC  >  NC_000913/2421294‑2421592
                                                                                                                                                          |                                                                                                                                                     
tCATCCTAAGTAAAACAACGGAGAACCTGCAATGGC‑ACAACCTGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTCAGGAAGAAGAAATTAACA‑ACCGTC‑AGGATATTGTTATCAATGTGACGATCCA‑CTAc                                                                                                                                                      >  2:554218/1‑151 (MQ=255)
    ccTAAGTAAAACAACGGAGAACCTGCAATGGC‑ACAACCTGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACA‑ACCGTC‑AGGATATTGTTATCAATGTGACGATCCA‑CTAcccc                                                                                                                                                   >  2:1923251/1‑150 (MQ=255)
    ccTAAGTAAAACAACGGAGAACCTGCAATGGC‑ACAACCTGCCGCTATTATTCGTACAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACA‑ACCGTC‑AGGATATTGTTATCAATGTGACAATCCA‑CTAcccc                                                                                                                                                   <  2:484203/150‑1 (MQ=255)
      tAAGTAAAACAACGGAGAACCTGCAATGGC‑ACAACCTGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACA‑ACCGTC‑AGGATATTGTTATCAATGTGACGATCCA‑CTACCccgcc                                                                                                                                                >  1:1558574/1‑151 (MQ=255)
       aaGTAAAACAACGGAGAACCTGCAATGGC‑ACAACCTGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACA‑ACCGTC‑AGGATATTGTTATCAATGTGACGATCCA‑CTACCccgc                                                                                                                                                 <  1:258181/149‑1 (MQ=255)
        aGTAAAACAACGGAGAACCTGCAATGGC‑ACAACCTGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACA‑ACCGTC‑AGGATATTGTTATCAATGTGACGATCCA‑CTACCccgccg                                                                                                                                               >  2:198189/1‑150 (MQ=255)
          tAAAACAACGGAGAACCTGCAATGGC‑ACAACCTGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACA‑ACCGTC‑AGGATATTGTTATCAATGTGACGATCCA‑CTACCCCGCCGAt                                                                                                                                             >  2:1555937/1‑150 (MQ=255)
           aaaaCAACGGAGAACCTGCAATGGC‑ACAACCTGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACA‑ACCGTC‑AGGATATTGTTATCAATGTGACGATCCA‑CTACCCCGCCGATa                                                                                                                                            >  2:1875426/1‑150 (MQ=255)
                        aCCTGCAAAGGCAACACCCTCCCGTTATTATTCGTAAAAAGACCT‑TTCTTTTGCGTCCTTTTATCGAAATAAAGAAAAAAAAAATTAACA‑CCCGAA‑AGGATTTTTTTATAAATGTGCCGTCCCC‑CTACCCCGCCGATAAAgcgc                                                                                                                                      <  2:1861412/144‑1 (MQ=255)
                         ccTGAAAGGAC‑ACACCCTCCCTTTATTTTTCTAAAAAAGACCC‑TTCTTTTGCGTCTGTTTATCGAAATAAAGAAAAAAAAAATAAACA‑CCCGTC‑AGGATTTTTTTATAAATGTGACGATCCC‑CTCCCCCGCCGATAAAGCGCGCACTAGCg                                                                                                                             <  1:463357/151‑1 (MQ=255)
                          cTGCAATGGC‑ACAACCTGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACA‑ACCGTC‑AGGATATTGTTATCAATGTGACGATCCA‑CTACCCCGCCGATAAAGCGCGCACTAGCg                                                                                                                             >  2:1505243/1‑150 (MQ=255)
                           tgaaaggcc‑acaCCCTGCGGTTATTATTGGTAAAAAGACCT‑TTGTTTTGCGTCCTTTTATCGAAATAAAGGAA‑AAAAAAATAAAA‑ACCCTCGGGGATATTGTTAAAAATGTGACGACCCA‑CTACCCCGCCGATAAAGCGCGCACTAGCg                                                                                                                             <  1:651772/141‑1 (MQ=255)
                           tGCAATGGC‑ACAACCTGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACA‑ACCGTC‑AGGATATTGTTATCAATGTGACGATCCA‑CTACCCCGCCGATAAAGCGCGCACTAGCGa                                                                                                                            >  2:1588272/1‑150 (MQ=255)
                              aagggC‑ACACCCTGCCGTTTTTATTCGAAAAAAAAACT‑TTCTTTTGCGTACGTTTATCGAAATAAAGGAAAAAAAAATAAACA‑CCCGTC‑AGGTTTTTGTTATAAATGGGCCGTCCCA‑CTCCCCCGCCGATAAAGCGCGCACTAGCGaa                                                                                                                           <  1:1004016/145‑1 (MQ=255)
                               aataacacccccTGCCGTTATTATTCGTAAAAAGAACC‑TTCTTTGGCGTACTTTTATGGAAATAAAGGAAAAAAAAATAACCA‑CCCGTC‑GGGATTTTGTTAACAATGTGACGCTCCA‑CCCCCCCGCCGATAAAGCGCGCACTAGCGAAg                                                                                                                          <  1:436574/139‑1 (MQ=255)
                                gggC‑ACAACCGGCCGTTTTTATTCGTAAAAAAAACC‑TTCGTTTGCGTATGTTTGGCATAATAAAGAAAGAAAAAATTAACA‑CCGGCA‑AGTATATTGTAAACAATGTGACGACCCA‑CTACCCCGCCGATAAAGCGCGCACTAGCGAAGatat                                                                                                                      <  1:454388/150‑1 (MQ=255)
                                   c‑acaACCTGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACA‑ACCGTC‑AGGATATTGTTATCAATGTGACGATCCA‑CTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCa                                                                                                                    >  2:1346426/1‑150 (MQ=255)
                                      caACCTGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACA‑ACCGTC‑AGGATATTGTTATCAATGTGACGATCCA‑CTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACg                                                                                                                 >  2:842304/1‑151 (MQ=255)
                                      caACCTGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACA‑ACCGTC‑AGGATATTGTTATCAATGTGACGATCCA‑CTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACg                                                                                                                 >  2:1880731/1‑151 (MQ=255)
                                       acccTCCCGTTATTATTCGTAAAAAAACCT‑TTCTTTTGCGTACTTTTATCGAAATAAAGAAAAAAAAAATAACCA‑CCCGTA‑GGGATATTGTAATAAATGGGACGACCCC‑CTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAAc                                                                                                                  <  1:1451217/147‑1 (MQ=255)
                                         ccTGCCGTTTTTATTCGTAAAAAGACCC‑TTTGTTTGCGTACTTTTATCGAAATAAAGGAAAAAAAAATAAACA‑ACCGTC‑AGGATTTTGTTATCAATGTGACGATCCA‑CTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACg                                                                                                                 <  1:1620632/148‑1 (MQ=255)
                                           tGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACA‑ACCGTC‑AGGATATTGTTATCAATGTGACGATCCA‑CTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATgcg                                                                                                            >  2:1245899/1‑151 (MQ=255)
                                           tGCCGCTATTATTCGTATAAAAAACC‑TCCTTTGCCGAACGTTTTTCGGAATTAAGGAAGAAGAAATTAACA‑ACGGTA‑AGGATTTTTTTACCATTTTGACG‑TCCACCAACCCACCCGTTAAACCGCGCACAAGTAAAGTTATCAACGATgc                                                                                                             >  2:233665/1‑150 (MQ=255)
                                                   ttattCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACA‑ACCGTC‑AGGATATTGTTATCAATGTGACGATCCA‑CTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAAtt                                                                                                     >  1:666992/1‑150 (MQ=255)
                                                    tattCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACA‑ACCGTC‑AGGATATTGTTATCAATGTGACGATCCA‑CTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTAt                                                                                                   >  1:1596747/1‑151 (MQ=255)
                                                     attCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACA‑ACCGTC‑AGGATATTGTTATCAATGTGACGATCCA‑CTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATc                                                                                                  >  2:1894115/1‑151 (MQ=255)
                                                      ttCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACA‑ACCGTC‑AGGATATTGTTATCAATGTGACGATCCA‑CTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCg                                                                                                 >  1:1871070/1‑151 (MQ=255)
                                                       ttctaaaaaaaaCCTTTCTTTTGCGTACTTTTATCGAAATTAAGAAAAAAAAAATAACCA‑CCCGTA‑AGATTTTTTTTATAAATGGGCCGCCCCA‑CTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTAt                                                                                                   <  1:1821502/143‑1 (MQ=255)
                                                       tCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACA‑ACCGTC‑AGGATATTGTTATCAATGTGACGATCCA‑CTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGc                                                                                                >  2:1338044/1‑151 (MQ=255)
                                                        cataaaaaaaaCC‑TTCTTTTGCGTACTTTTATGGAAATAAAGGAAGAAAAAATAAACA‑CCCGTA‑AGGATATTTTTAAAAATGTGCCGACCCC‑CTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCa                                                                                               <  1:1815106/146‑1 (MQ=255)
                                                           aaaaaaaaCC‑TTCTTTGGCGTACTTTTATCGAAATAAAGAAAAAAAAAATAACCA‑CCCGTC‑AGGATATTTTAATCAAGGTGCCGACCCC‑CTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCg                                                                                            <  1:2070778/149‑1 (MQ=255)
                                                           aaaaaGATCC‑TTCTTGTGCTGATGTTTAGAAGAATAAAGGAAAAAGAAATAAACA‑CCCGAA‑AGGATTTTTTTATCATTGTGAAGACCCA‑CCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCa                                                                                               <  2:2005748/146‑1 (MQ=255)
                                                           aaaaaGACCC‑TTTGTTTGCGTACTTTTATCGAAATTAAGGAAAAAAAAATAAACA‑CCCGTC‑AGGATTTTGTTATCAATGTGACGACCCA‑CTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCAc                                                                                              <  1:1546887/147‑1 (MQ=255)
                                                            taAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACA‑ACCGTC‑AGGATATTGTTATCAATGTGACGATCCA‑CTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGt                                                                                           >  2:2015814/1‑151 (MQ=255)
                                                            aaaaaaatctttCTTTTGTGTACTTTTATGGAAATAAAGAAAAAAGAAATAACCA‑CCCTTA‑AGAATTTTTTAATAAATGTGACATCCCC‑CCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCa                                                                                               <  2:801930/138‑1 (MQ=255)
                                                             aaaaaccc‑TTCTTTTCCTTACTTTTATCGAAATAAAGAAAAAAAAAATAACCA‑CCCTTA‑ATAATTTTTTTATCAATTTGACATCCCC‑CTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTTAATTATCGCACCGt                                                                                           <  2:459517/144‑1 (MQ=255)
                                                             aaaaaaCT‑TTCTTTTGCGTGCTTTAAGGAGAATTAAGAAAGAAGAAATAAACA‑CCCGAA‑AGAATTTTGTTAACAAGGTGACGTTCCC‑CTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTa                                                                                          <  2:364514/151‑1 (MQ=255)
                                                             aaaGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACA‑ACCGTC‑AGGATATTGTTATCAATGTGACGATCCA‑CTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGt                                                                                           >  2:1925892/1‑150 (MQ=255)
                                                              aaGACCC‑TTCGTTTGCGTACGTTTATCGAAATTAAGAAAGAAAAAATAAACA‑CCCGTC‑AGGATATTTTTATAAATGTGCCGACCCC‑CTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTaa                                                                                         <  1:1943810/151‑1 (MQ=255)
                                                              aaGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACA‑ACCGTC‑AGGATATTGTTATCAATGTGACGATCCA‑CTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTa                                                                                          >  1:811664/1‑150 (MQ=255)
                                                                 accc‑TTCTTTTGCGTACTTTTATGGAAATAAAGAAAAAAAAAATAACCA‑CCCGCC‑AGGATATTTTAATAAATGAGACGACCCA‑CCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACg                                                                                       <  2:941721/148‑1 (MQ=255)
                                                                 accc‑TTCTTTGGCGTCCGTTTATCGAAATAAAGAAAAAAAAAATAAACA‑CCCGTC‑AGGATTTTTTTATAAATGTGACGACCCC‑CTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACg                                                                                       <  1:1236670/148‑1 (MQ=255)
                                                                  acctttCTTTTGCGTTCTTTTATCAGAATAAAGGAAAAAAAAATAACCA‑CCCGAC‑AGTATATTTTTAAAATGGTGACGCCCCA‑CCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGa                                                                                      <  2:383173/147‑1 (MQ=255)
                                                                   cc‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACA‑ACCGTC‑AGGATATTGTTATCAATGTGACGATCCA‑CTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGaa                                                                                     >  1:1550767/1‑150 (MQ=255)
                                                                       tCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACA‑ACCGTC‑AGGATATTGTTATCAATGTGACGATCCA‑CTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAAc                                                                                 >  1:269285/1‑151 (MQ=255)
                                                                          tttGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACA‑ACCGTC‑AGGATATTGTTATCAATGTGACGATCCA‑CTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACatt                                                                              >  1:1163787/1‑151 (MQ=255)
                                                                          ttgtcttactttttccGGAATAAAGAAAAAAAAAATTAACA‑CCCGTC‑AGGATATTTTTATCAATGTGACGACCCC‑CACCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACat                                                                               <  1:1701929/135‑1 (MQ=255)
                                                                             gCGTACGTTTATCGAAATAAAGGAAAAAAAAATAAACA‑CCCGTC‑AGGATTTTTTTATCAATGTGACGACCCC‑CTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACattatt                                                                           <  1:1676536/151‑1 (MQ=255)
                                                                              cGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACA‑ACCGTC‑AGGATATTGTTATCAATGTGACGATCCA‑CTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTc                                                                          >  1:239707/1‑151 (MQ=255)
                                                                               gTACGTTTATCGGAATTAAGGAAGAAGAAATTAACA‑ACCGTC‑AGGATATTGTTATCAATGTGACGATCCA‑CTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTc                                                                          >  1:353205/1‑150 (MQ=255)
                                                                               gTACGTTTATCGGAATTAAGGAAGAAGAAATTAACA‑ACCGTC‑AGGATATTGTTATCAATGTGACGATCCA‑CTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTc                                                                          >  2:753085/1‑150 (MQ=255)
                                                                               gTACGTTTATCGGAATTAAGGAAGAAGAAATTAACA‑ACCGTC‑AGGATATTGTTATCAATGTGACGATCCA‑CTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTc                                                                          >  1:661580/1‑150 (MQ=255)
                                                                                ttcttttATGGAAATAAAGGAAGAAAAAATAAACA‑CCCGTC‑AGGATTTTTTTATAAATGTGACAACCCC‑CTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCa                                                                         <  1:251413/146‑1 (MQ=255)
                                                                                 aCGTTTATCGGAATTAAGGAAGAAGAAATTAACA‑ACCGTC‑AGGATATTGTTATCAATGTGACGATCCA‑CTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGc                                                                       >  2:186571/1‑151 (MQ=255)
                                                                                   gTTTATCGGAATTAAGGAAGAAGAAATTAACA‑ACCGTC‑AGGATATTGTTATCAATGTGACGATCCA‑CTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCAt                                                                     >  1:2069928/1‑151 (MQ=255)
                                                                                          ggAATTAAGGAAGAAGAAATTAACA‑ACCGTC‑AGGATATTGTTATCAATGTGACGATCCA‑CTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTag                                                                 >  1:471814/1‑148 (MQ=255)
                                                                                            aaTAAAGAAAAAAAAAATAAACA‑CCCGTC‑AGAATTTTGTTATCAATGGGACGATCCA‑CTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTagag                                                               <  1:218258/148‑1 (MQ=255)
                                                                                            aaTAAAGAAAAAAAAAATAAACA‑CCCGTC‑AGAATATTGTTATCAATGTGACGATCCA‑CTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAAt                                                            <  1:1480967/151‑1 (MQ=255)
                                                                                              ttAAGGAAGAAGAAATTAACA‑ACCGTC‑AGGATATTGTTATCAATGTGACGATCCA‑CTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATa                                                           >  1:1185160/1‑150 (MQ=255)
                                                                                              taaaGAAAAAAAAAATAAACA‑CCCGTG‑AGAATATTGTTATAAATGTGACGATCCC‑CTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATa                                                           <  1:898613/148‑1 (MQ=255)
                                                                                               tAAGGAAGAAGAAATTAACA‑ACCGTC‑AGGATATTGTTATCAATGTGACGATCCA‑CTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAAc                                                         >  1:465387/1‑151 (MQ=255)
                                                                                                aaGAAAAAAAAAATAACCA‑CCCGCA‑AGGATTTTTTTATCAAGGTGACGACCCC‑CACCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAgaga                                                              <  2:433263/145‑1 (MQ=255)
                                                                                                 aGGAAGAAGAAATTAACA‑ACCGTC‑AGGATATTGTTATCAATGTGACGATCCA‑CTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCg                                                       >  2:595413/1‑151 (MQ=255)
                                                                                                      gaagaaATTAACA‑ACCGTC‑AGGATATTGTTATCAATGTGACGATCCA‑CTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTtct                                                  >  2:1745830/1‑151 (MQ=255)
                                                                                                       aaaaaaTAAACA‑CCCTTC‑AGAATTTTTTTATAAATGTGACACCCCC‑CACCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTtctc                                                 <  2:1568425/148‑1 (MQ=255)
                                                                                                         aaaaaTAACACACCGCC‑AGAATTTTGTTATAAATGTGCCGACCCC‑CTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTctct                                                <  1:1933363/150‑1 (MQ=255)
                                                                                                            aTTAAAA‑ACCGTC‑AGGCTATTGTTATCCAAGTGACGATCCA‑CTACCCCGCCGTTAAAGCGCGCACTAGCGAAGATATCAACGATGAGCGGAATTATCTCACCGTAACGAAAAACCTTCTTCAGAATGAAGAGAATAACCGTCTCTCttt                                              >  2:1077217/1‑149 (MQ=255)
                                                                                                              tAACA‑ACCGTC‑AGGATATTGTTATCAATGTGACGATCCA‑CTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTg                                          >  1:1282851/1‑151 (MQ=255)
                                                                                                              tAACA‑ACCGTC‑AGGATATTGTTATCAATGTGACGATCCA‑CTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTg                                          >  1:1555130/1‑151 (MQ=255)
                                                                                                                   acccGTA‑AGAATTTTTTTATAAATGTGCCAACCCC‑CCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGaaa                                      <  1:777787/149‑1 (MQ=255)
                                                                                                                   acccGTA‑AGAATTTTGTTACAATTGTGACGACCCC‑CTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGaaa                                      <  1:325546/149‑1 (MQ=255)
                                                                                                                    acCGTC‑AGGATATTGTTATCAATGTGACGATCCA‑CTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGaaa                                      >  2:542799/1‑150 (MQ=255)
                                                                                                                     ccGTC‑AGGATTTTGTTATAAATGTGCCCCCCCC‑CCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGaaa                                      <  2:449169/149‑1 (MQ=255)
                                                                                                                     ccGTA‑AGAATATTTTAATAAATGTAACGACCCC‑CTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGaaaa                                     <  2:622919/150‑1 (MQ=255)
                                                                                                                       gTC‑AGGATATTGTTATCAATGTGACGATCCA‑CTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAAtt                                  >  2:2027918/1‑151 (MQ=255)
                                                                                                                           aGGATTTTTTTATAAATGGGACGACCCA‑CTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTaa                                <  1:1460601/150‑1 (MQ=255)
                                                                                                                              atatTGTTATCAATGTGACGATCCA‑CTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTc                             >  1:1535604/1‑150 (MQ=255)
                                                                                                                               tttttttaaaaaTGTGACGACCCA‑CTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTc                             <  1:1947902/139‑1 (MQ=255)
                                                                                                                               ttttGTTATCGGTGGGAAGACCCA‑CACCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTGTCTCTTTGCt                                           <  2:906500/133‑1 (MQ=255)
                                                                                                                               tatTGTTATCAATGTGACGATCCA‑CTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTc                             >  2:1546746/1‑149 (MQ=255)
                                                                                                                               tatTGTTATCAATGTGACGATCCA‑CTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTc                             >  2:311980/1‑149 (MQ=255)
                                                                                                                                 ttGTTATCAATGTGACGATCCA‑CTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAgg                          >  1:1040544/1‑150 (MQ=255)
                                                                                                                                      ataaaggtgACGACCCA‑CCACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGAtgtg                     <  2:1338529/144‑1 (MQ=255)
                                                                                                                                      ataaaggtgACAACCCC‑CTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGAtgtg                     <  2:1385464/144‑1 (MQ=255)
                                                                                                                                         atggtgACGACCCC‑CTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTc                  <  1:1137879/147‑1 (MQ=255)
                                                                                                                                         aaTGTGCCGACCCA‑CTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGAtgtg                     <  1:2181964/147‑1 (MQ=255)
                                                                                                                                         aaTGTGACGATCCA‑CTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCg                 >  2:858410/1‑151 (MQ=255)
                                                                                                                                         aaTGTGACGATCCA‑CTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCg                 >  2:693936/1‑151 (MQ=255)
                                                                                                                                         aaTGTGACGACCCC‑CTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTc                  <  1:1293967/150‑1 (MQ=255)
                                                                                                                                         aaTGTGACGACCCA‑CTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCg                 <  1:1479383/151‑1 (MQ=255)
                                                                                                                                          aagtgccaaTCCC‑CTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATGACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGAtgtg                     <  2:960356/138‑1 (MQ=255)
                                                                                                                                            gtgACGATCCA‑CTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGat               >  2:1956841/1‑150 (MQ=255)
                                                                                                                                            gtgACGACCCA‑CTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGat               <  1:1567237/150‑1 (MQ=255)
                                                                                                                                             ggacgacCCA‑CTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGat               <  1:1620477/148‑1 (MQ=255)
                                                                                                                                               acaccccA‑CTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGatat             <  1:1100922/144‑1 (MQ=255)
                                                                                                                                                 gATCCA‑CTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGc          >  1:2026167/1‑150 (MQ=255)
                                                                                                                                                    ccA‑CTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGt      >  2:1623590/1‑151 (MQ=255)
                                                                                                                                                     cA‑CTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTg     >  1:1450062/1‑151 (MQ=255)
                                                                                                                                                        cTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGa    <  2:1667746/150‑1 (MQ=255)
                                                                                                                                                         tACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAAc  >  1:630827/1‑151 (MQ=255)
                                                                                                                                                          |                                                                                                                                                     
TCATCCTAAGTAAAACAACGGAGAACCTGCAATGGC‑ACAACCTGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACA‑ACCGTC‑AGGATATTGTTATCAATGTGACGATCCA‑CTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAAC  >  NC_000913/2421294‑2421592

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