Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,229,1900TA22.0% 410.3 / 53.0 224intergenic (+125/+248)yohP/dusCputative membrane protein YohP/tRNA‑dihydrouridine(16) synthase
Reads supporting (aligned to +/- strand):  ref base T (34/136);  new base A (48/0);  total (86/138)
Fisher's exact test for biased strand distribution p-value = 2.47e-26
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 226 total aligned reads displayed.

TGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCT‑GGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCG  >  NC_000913/2229041‑2229338
                                                                                                                                                      |                                                                                                                                                      
tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAAt                                                                                                                                                        <  1:1436885/150‑1 (MQ=255)
 ggCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATc                                                                                                                                                       <  1:1391176/150‑1 (MQ=255)
    gTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATcgcc                                                                                                                                                    <  1:2002652/150‑1 (MQ=255)
               gatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTTGCGCAAGGAGCGC‑GCCAGGGGCGGCCCAAATCCGccccccccac                                                                                                                                         >  2:1105524/1‑148 (MQ=255)
                  atTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGCGGGGGGGGGCCAAAC‑CCGccccccccccccc                                                                                                                                      >  1:401024/1‑145 (MQ=255)
                   ttttAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCCAAAGCCCCCGCCCCCAACAcc                                                                                                                                     >  2:657367/1‑150 (MQ=38)
                   ttttAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAAACGCCGCCGCCCCCAGCAcc                                                                                                                                     >  2:1067667/1‑150 (MQ=255)
                    tttAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAAAACCCGCCCCCCCCAGCAccc                                                                                                                                    >  2:2012767/1‑150 (MQ=38)
                        aaaTTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGCGGGGGGGGCCCAAAGCCCGCCCCCCCCGCAccccccg                                                                                                                                >  2:1259240/1‑148 (MQ=255)
                        aaaTTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGAGGGGGGGGCCCCACGCCCCCCCCCCCAcccccccccg                                                                                                                                >  2:784232/1‑148 (MQ=255)
                         aattaattAATGTCATCAGGTCCTAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCCGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCggg                                                                                                                               <  1:657606/150‑1 (MQ=255)
                         aattaattAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCggg                                                                                                                               <  2:1713096/150‑1 (MQ=255)
                           ttaattaaTGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGCCGGGGGGGGCCAA‑CGCCGCCCGCCCCCCCACCCCCcgggg                                                                                                                             >  2:241303/1‑148 (MQ=255)
                              attaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCT‑g                                                                                                                         <  2:1839693/150‑1 (MQ=255)
                                taatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCT‑GGc                                                                                                                       <  1:528076/150‑1 (MQ=255)
                                  atGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCCAACGCCCCCCCCCCCCGCACCCCCCGGCTCA‑GGCGa                                                                                                                     >  2:543861/1‑150 (MQ=37)
                                   tGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCT‑GGCGaa                                                                                                                    <  2:313788/150‑1 (MQ=255)
                                    gTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCT‑GGCGAAc                                                                                                                   <  2:1037443/150‑1 (MQ=255)
                                     tcatcaGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCT‑GGCGAACa                                                                                                                  <  2:736180/150‑1 (MQ=255)
                                      catcaGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGCGGGGGGGGCCCAAACCCGCCCCCCCCAGCACCCCCCGGCGCT‑GGGGAaaaa                                                                                                                 >  2:1227753/1‑150 (MQ=255)
                                          aGGTCCTAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCT‑GGCGAACAAAATc                                                                                                             <  2:1457500/150‑1 (MQ=255)
                                          aGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCCTATAATCACTGGCGCAAGGAGCGC‑GGGGGGGGGGGGCCAACCCCCCCCCCCCCAGCACCCCCCCGCGCG‑GGCGCACacaaaa                                                                                                             >  1:615194/1‑148 (MQ=255)
                                          aGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGGGGGGGGGGGCCAAAGCCGCCCCCCCCCGCCCCCCCCGGCGCG‑GGGGCACaaaaaa                                                                                                             >  2:154337/1‑148 (MQ=255)
                                              cctaaaaTAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCTTGTGCATTGCTTCATATAATCCCTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCT‑GGCGAACAAAATcgccg                                                                                                         <  1:1473687/147‑1 (MQ=255)
                                               cGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGCGGGGGGGGCCCAACGCCCCCCCCCCCACCACCCCCCCGCGCA‑GGCGAACaaaaaaacccc                                                                                                        >  1:1567910/1‑143 (MQ=255)
                                               cGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGAGGGGGGGGGCCAAAGCCCCCCCCCCCCGCCCCCCCCCGCGCA‑GGGGAACaaaaaaacccc                                                                                                        >  2:691865/1‑143 (MQ=255)
                                                gAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGGGGGGGGGGCCCAAACCCGCCCCCCCCCGCCCCCCCCGGGGCTGGGGGAACaaaaaaccccg                                                                                                        >  2:165245/1‑143 (MQ=25)
                                                gAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCCAACGCCCCCCCCCCCCCCAACCCCCGGCGCT‑GGGGAACAAAAAagccgcg                                                                                                       >  1:2132987/1‑149 (MQ=25)
                                                 aaaaTAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCCATCGCCGCCCCCCCCCGCACCCCCCGGCTCT‑GGGGAACCAAAAcgccgccg                                                                                                      >  2:1902765/1‑148 (MQ=25)
                                                  aaaTAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCGGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCCGGGGGGGGGCAATACCCCCCGCCCCCAGCCCCCCCCGGGGCA‑TGCGAACAAAAAAGCCGCGgc                                                                                                     >  2:448600/1‑150 (MQ=255)
                                                       aCGATAATATTTCAGTCTCTCATCCTGTTGCGCTCCTTTCATGTGCATTGCTACATATAACCCCTGGCGCACGGAGCGC‑GCAGGGGGCGGCCAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCT‑GGCGAACAAAATCGCCGCTGCGCGGt                                                                                                <  1:995903/150‑1 (MQ=17)
                                                                     gTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCT‑GGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATc                                                                                  <  2:197680/150‑1 (MQ=255)
                                                                     gTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCT‑GGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATc                                                                                  <  2:933948/150‑1 (MQ=255)
                                                                      tctctcATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCCGGGGGGGGCCAAACGCCGCCCCCCCCAGCACCCCCCGGCTCT‑GGCGACCAAAAAAGCCGCTGCGCGGTGCCCTCGGCGTAtac                                                                                 >  2:1580742/1‑148 (MQ=255)
                                                                        tctcATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCT‑GGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCt                                                                               <  1:308946/150‑1 (MQ=255)
                                                                          tcATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCT‑GGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTa                                                                             <  1:1262044/150‑1 (MQ=255)
                                                                           cATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCT‑GGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTAc                                                                            <  2:1782037/150‑1 (MQ=255)
                                                                           cATCCTGTTGCGCTCCTGTCATGTGCATTGATTAATATAATCACTGGCGCAAGGAGCGC‑GCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCT‑GGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTAc                                                                            <  1:830229/150‑1 (MQ=255)
                                                                            aTCCTGTTGCGCTCATGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCT‑GGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACg                                                                           <  2:627285/150‑1 (MQ=255)
                                                                                tGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGGGGGGGGGGCCCAAAGCCCCCCCCCCCAGCACCCCCCGGGGCT‑GGGGGACAAAAAAGCCCCGGCGCGGGGGCCGCGGCGTATACCATaacgcga                                                                       >  2:371219/1‑144 (MQ=12)
                                                                                      gcTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGCGGGGGGGGGCCAA‑CCCCGCCCCCCCCCGCACCCCCCGGGGCT‑GGGGGACAAAATAGCCCCGGGGCGGTGGCCACGGCGTATACCATAAGGGTAACCggg                                                                 >  2:481364/1‑150 (MQ=255)
                                                                                       cTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGGGGGGGGGGCCCATCGCCCCCGGCCCCAGCACCCCCCGGGGCT‑GGCGAGCACAATAGCCCCGGGGGGGTGGCCTCGGCGGTTCCCCTACGGCT‑CCCggggg                                                               >  2:1200175/1‑148 (MQ=255)
                                                                                       cTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCT‑GGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTc                                                               <  1:1221260/150‑1 (MQ=255)
                                                                                           tGTCATGTGCATTGCTTCATATAATCCCTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCT‑GGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGgc                                                           <  1:467380/150‑1 (MQ=255)
                                                                                           tGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCT‑GGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGgc                                                           <  2:1941751/150‑1 (MQ=255)
                                                                                               aTGTGCATTGCTTCATATACTCACTGGCGCACGTGGGGC‑GAAGGGGGCGGCAAATCGCCGCCGCCCCCTGCACCCCCGGGCTCT‑GGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGa                                                       <  1:970820/150‑1 (MQ=14)
                                                                                                tgtgCATTGCTTCCTATCATCCCTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCCCCGCCGCCCCCTGCACCCCCGGGCTCT‑GGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAg                                                      <  2:1404901/150‑1 (MQ=255)
                                                                                                  tgCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCT‑GGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGt                                                    <  1:1972131/150‑1 (MQ=255)
                                                                                                  tgCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCT‑GGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGt                                                    <  1:1600151/150‑1 (MQ=255)
                                                                                                  tgCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCT‑GGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGt                                                    <  1:44750/150‑1 (MQ=255)
                                                                                                   gCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCT‑GGCGAACAAACTCGCCGCTGCGCTGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTa                                                   <  1:426100/150‑1 (MQ=255)
                                                                                                      ttGCTTCATATAATCACTGGCGCATGGAGCGC‑GCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCT‑GGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACa                                                <  1:657238/150‑1 (MQ=255)
                                                                                                      ttGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGCGGGGGCGGGCAAAACCCGCCGCCCCCAGCACCCCCGGGCTCT‑TGCGAACAAAAACGCCGCGGCGCGGGGCCCCCGGCGTATACCATAAGGCT‑ACCCCGTCGGGCGCGCGGGAACa                                                >  2:1159800/1‑150 (MQ=255)
                                                                                                      ttGCTTAATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCT‑GGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACa                                                <  1:228457/150‑1 (MQ=255)
                                                                                                       tGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCT‑GGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACAt                                               <  1:811131/150‑1 (MQ=255)
                                                                                                        gCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCT‑GGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACATc                                              <  2:805244/150‑1 (MQ=255)
                                                                                                        gCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCT‑GGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACATc                                              <  2:970155/150‑1 (MQ=255)
                                                                                                         cTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCACTCGCCGCCGCCCCCTGCACCCCCGGGCTCT‑GGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTCCGGCT‑ACCGGGTCGGGCGCGAGGTAACATcc                                             <  2:1337786/150‑1 (MQ=255)
                                                                                                         cTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCT‑GGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACATcc                                             <  2:218626/150‑1 (MQ=40)
                                                                                                            cATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCT‑GGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACATCCCTg                                          <  2:466841/150‑1 (MQ=38)
                                                                                                             atatAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCT‑GGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACATCCCTGt                                         <  1:185512/150‑1 (MQ=38)
                                                                                                             atatAAACACTGGCGAATGGAGGGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCT‑GGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACATCCCTGt                                         <  1:759829/150‑1 (MQ=255)
                                                                                                                 aaTCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAAACGCCGCCGCCCCCAGCACCCCCCGGCTCT‑TGCGAACAAAAACGCCGCGGCGCGGGGCCCCCCGCGTATACCATAAGGCA‑ACCCGGGCGGGCGCGGGGGTACCTCCCAGTaaaa                                     >  2:1464456/1‑150 (MQ=2)
                                                                                                                  aTCACTGGCGCAAGGAGCGC‑GGAGGGGGGGGGCCAACGCCGCCGCCCCCGGGCCCCCCCGGCGCT‑GGCGCACAAAAAAGCCCCGGGGCGGTGGCCCCGGCGTATCCCATAAGGCG‑ACCCGGGGGGGGGCGGGGGAAAAACCCCGGaaaaa                                    >  2:1749573/1‑149 (MQ=255)
                                                                                                                    cACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCT‑GGCAAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAcgc                                  <  2:1224502/150‑1 (MQ=17)
                                                                                                                    cACGGGCGCATGGAGGGC‑GCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCT‑GGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACATCCCTGTAAACcgc                                  <  1:1831347/150‑1 (MQ=11)
                                                                                                                      cTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCT‑GGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACgcgc                                <  1:101534/150‑1 (MQ=32)
                                                                                                                      cTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCT‑GGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACgcgc                                <  2:1349546/150‑1 (MQ=32)
                                                                                                                      cTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCT‑GGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACgcgc                                <  1:101431/150‑1 (MQ=32)
                                                                                                                      cTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCT‑GGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACgcgc                                <  1:970432/150‑1 (MQ=32)
                                                                                                                      cTGGCGAAAGGAGCGC‑GCAGGGGGCGGCAAATCGCCGCCGCCCCCTGCACCCCCGGGCTCT‑GGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACgcgc                                <  1:1689594/150‑1 (MQ=17)
                                                                                                                       tGGGGCATGGAGCGC‑GCTGGGGGCGGCCATGCGCCGCCGCCCCCTGCACCCCCGGGCTCT‑GGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGcc                               <  1:152469/150‑1 (MQ=11)
                                                                                                                             aaGGAGCGC‑GCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCT‑GGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGc                         <  1:1113431/150‑1 (MQ=17)
                                                                                                                              aGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCT‑GGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGcc                        <  1:1257317/150‑1 (MQ=32)
                                                                                                                                  gcgc‑gcAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCT‑GGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCcaca                    <  2:1884538/150‑1 (MQ=32)
                                                                                                                                  gcgc‑gcAGGGGGCGGCCAATCACCGCCGCCCCCTGCACCCCCGGGCTCT‑GGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCcaca                    <  1:1917606/150‑1 (MQ=17)
                                                                                                                                       gccgggggCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCT‑GGCGAACAAAATCGCCGCGGCGCGGTGCCCACGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACATACCAGTAAAAAGCGCCCTCAGCCCACAACCa                >  2:1048687/4‑150 (MQ=11)
                                                                                                                                        cAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCT‑GGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAt               <  2:1477331/150‑1 (MQ=32)
                                                                                                                                         aGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCT‑GGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAtg              >  2:403472/1‑150 (MQ=32)
                                                                                                                                          gggggCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCT‑GGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAtgt             <  2:264034/150‑1 (MQ=32)
                                                                                                                                          gggggCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCT‑GGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAtgt             >  1:574606/1‑150 (MQ=32)
                                                                                                                                          gggggCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCT‑GGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAtgt             <  1:2085141/150‑1 (MQ=32)
                                                                                                                                          gggggCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCT‑GGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAtgt             <  1:1030294/150‑1 (MQ=32)
                                                                                                                                          gggggCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCT‑GGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAtgt             <  2:1078155/150‑1 (MQ=32)
                                                                                                                                           ggggCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCT‑GGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAtgtg            >  2:1605155/1‑150 (MQ=32)
                                                                                                                                           ggggCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCT‑GGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAtgtg            <  2:1249203/150‑1 (MQ=32)
                                                                                                                                            gggCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCT‑GGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTgg           <  2:500856/150‑1 (MQ=32)
                                                                                                                                             ggcggcCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCT‑GGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTggg          <  1:1685520/150‑1 (MQ=32)
                                                                                                                                             ggcggcCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCT‑GGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTggg          <  2:1420019/150‑1 (MQ=32)
                                                                                                                                             ggcggcCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCT‑GGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACACCCATGTggg          <  1:2012767/150‑1 (MQ=17)
                                                                                                                                              gcggcCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCT‑GGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCATACGGCA‑ACCGGGTCGGGCGCGAGGTAACAACCCTGTAAAACGCGCCCTCAGCCCACAACCAAGAGGGc         >  1:957365/1‑150 (MQ=11)
                                                                                                                                                  ccAATCGCCGCCGCCCCCTGCACCCCCGGGCTCT‑GGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATGGATGTGGGCTgag     <  1:1198372/150‑3 (MQ=12)
                                                                                                                                                  ccAATCGCCGCCGCCCCCTGCACCCCCGGGCTCT‑GGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGcc     >  2:1890727/1‑150 (MQ=32)
                                                                                                                                                  ccAATCGCCGCCGCCCCCTGCACCCCCGGGCTCT‑GGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGcc     <  2:1299365/150‑1 (MQ=32)
                                                                                                                                                   cAATCGCCGCCGCCCCCTGCACCCCCGGGCTCT‑GGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGccc    >  2:872251/1‑150 (MQ=32)
                                                                                                                                                   cAATCGCCGCCGCCCCCTGCACCCCCGGGCTCT‑GGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGccc    >  2:954463/1‑150 (MQ=32)
                                                                                                                                                    aaTCGCCGCCGCCCCCTGCACCCCCGGGCTCT‑GGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGcc                        >  1:1821819/1‑129 (MQ=32)
                                                                                                                                                    aaTCGCCGCCGCCCCCTGCACCCCCGGGCTCT‑GGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGcc                        <  2:1821819/129‑1 (MQ=32)
                                                                                                                                                    aaTCGCCGCCGCCCCCTGCACCCCCGGGCTCT‑GGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGcccc   >  2:823034/1‑150 (MQ=32)
                                                                                                                                                    aaTCGCCGCCGCCCCCTGCACCCCCGGGCTCT‑GGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGcccc   >  2:425423/1‑150 (MQ=32)
                                                                                                                                                      |                                                                                                                                                      
TGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCT‑GGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCG  >  NC_000913/2229041‑2229338

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑