Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,731,582 A→C 5.3% intergenic (‑425/+18) rrsG ← / ← clpB 16S ribosomal RNA of rrnG operon/protein disaggregation chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,731,5820AC 5.3% 824.1 / 21.8 238intergenic (‑425/+18)rrsG/clpB16S ribosomal RNA of rrnG operon/protein disaggregation chaperone
Reads supporting (aligned to +/- strand):  ref base A (121/101);  new base C (6/8);  total (129/109)
Fisher's exact test for biased strand distribution p-value = 4.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 241 total aligned reads displayed.

GTTGTTCCGTGTCAGTGGTGGCGCATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGCACGTTTCAGAGGA  >  NC_000913/2731433‑2731729
                                                                                                                                                     |                                                                                                                                                   
gttgttCCGTGTCAGTGGTGGCGCATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGa                                                                                                                                                     >  1:915751/1‑150 (MQ=255)
gttgttCCGTGTCAGTGGTGGCGCATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGa                                                                                                                                                     >  1:3029164/1‑150 (MQ=255)
      ccGTGTCAGTGGTGGCGCATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCt                                                                                                                                               >  2:1120944/1‑150 (MQ=255)
      ccGTGTCAGTGGTGGCGCATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCt                                                                                                                                               >  1:2651533/1‑150 (MQ=255)
         tgtCAGTGGTGGCGCATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCTTCGGGGCTCGt                                                                                                                                            <  2:1227522/150‑1 (MQ=255)
           tCAGTGGTGGCGCATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCTTCGGGGCTCGttt                                                                                                                                          <  2:1281251/150‑1 (MQ=255)
            cAGTGGTGGCGCATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGtttt                                                                                                                                         >  1:417448/1‑150 (MQ=255)
               gggtggCGCATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTGATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCTTCGGGGCTCGtttttgt                                                                                                                                      <  1:442708/149‑4 (MQ=39)
                    gcgcATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTa                                                                                                                                 <  2:773014/150‑1 (MQ=255)
                     cgcATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTAc                                                                                                                                <  2:894595/150‑1 (MQ=255)
                     cgcATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTAc                                                                                                                                >  1:1604217/1‑150 (MQ=255)
                     cgcATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTAc                                                                                                                                <  2:1864031/150‑1 (MQ=255)
                      gcATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACt                                                                                                                               >  2:1986297/1‑150 (MQ=255)
                            taGGGAGTTAGTCCGGCCGGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACgg                                                                                                                         <  2:556927/150‑1 (MQ=255)
                               ggAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAATCGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGa                                                                                                                      >  2:2960302/1‑150 (MQ=255)
                               ggAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGa                                                                                                                      >  2:1225757/1‑150 (MQ=255)
                               ggAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTGTTATCATTTACTGGACGGCgc                                                                                                                      >  1:3253719/1‑149 (MQ=255)
                                     aTTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTAATGGACGGCGACAATcc                                                                                                                <  2:1710757/150‑1 (MQ=255)
                                      ttCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATCGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCg                                                                                                               <  1:2477664/150‑1 (MQ=255)
                                      ttCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCg                                                                                                               >  1:1477496/1‑150 (MQ=255)
                                      ttCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCg                                                                                                               <  2:2399128/150‑1 (MQ=255)
                                        ccGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGt                                                                                                             <  1:2241899/150‑1 (MQ=255)
                                        ccGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGt                                                                                                             <  1:95555/150‑1 (MQ=255)
                                         cGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTc                                                                                                            >  1:1822088/1‑150 (MQ=255)
                                         cGGCCTGACAAGCGAAAAATATAAAAACTTTATCGGTCGCTCACTTTTCAGGCAAAAAATCTTAAATATAGTCTTTTCCGTCTATCTTATAGACAAAAACGAGCCCCGAAGGTCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTc                                                                                                            <  1:366051/150‑1 (MQ=255)
                                             cTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCa                                                                                                        <  1:169658/150‑1 (MQ=255)
                                              tGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCAt                                                                                                       <  1:1986968/150‑1 (MQ=255)
                                              tGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCAt                                                                                                       <  1:1986771/150‑1 (MQ=255)
                                                aCAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTa                                                                                                     >  1:1348646/1‑150 (MQ=255)
                                                   aGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACt                                                                                                  <  1:755824/150‑1 (MQ=255)
                                                    gCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGCAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACtt                                                                                                 >  1:1573145/1‑150 (MQ=255)
                                                    gCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTAATGGACGGCGACAATCCGGTCTTCATTAACtt                                                                                                 >  1:3094828/1‑150 (MQ=255)
                                                    gCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTAATGGACGGCGACAATCAGGTCTTCATTAACtt                                                                                                 >  1:1726364/1‑150 (MQ=255)
                                                     cGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAGTCCGGTCTTCATTAACTTc                                                                                                >  1:2306425/1‑150 (MQ=255)
                                                     cGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTc                                                                                                >  1:478021/1‑150 (MQ=255)
                                                       aaaaaTATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCa                                                                                              >  1:3695891/1‑150 (MQ=255)
                                                              taAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAAt                                                                                       >  2:4057956/1‑150 (MQ=255)
                                                               aaaaaCTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAATTATAGTCTTTTCCGTCTAACTTATAGACAATATCTTGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATc                                                                                      >  2:2702195/1‑150 (MQ=255)
                                                               aaaaaCTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATc                                                                                      >  2:4023362/1‑150 (MQ=255)
                                                               aaaaaCTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATc                                                                                      >  1:1859283/1‑150 (MQ=255)
                                                               aaaaaCTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATc                                                                                      >  1:2320502/1‑150 (MQ=255)
                                                                aaaaCTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCa                                                                                     <  2:3805314/150‑1 (MQ=255)
                                                                 aaaCTTTCTCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCAc                                                                                    <  1:4098581/150‑1 (MQ=255)
                                                                  aaCTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACt                                                                                   <  1:1850641/150‑1 (MQ=255)
                                                                  aaCTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACt                                                                                   <  2:447244/150‑1 (MQ=255)
                                                                            ttCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGaa                                                                         <  2:2443665/150‑1 (MQ=255)
                                                                            ttCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGaa                                                                         <  2:1920535/150‑1 (MQ=255)
                                                                            ttCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGaa                                                                         <  2:3170306/150‑1 (MQ=255)
                                                                             tCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAAc                                                                        <  1:2877283/150‑1 (MQ=255)
                                                                                     ttttCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCAc                                                                >  1:1175790/1‑150 (MQ=255)
                                                                                       ttCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCa                                                              <  2:1800624/150‑1 (MQ=255)
                                                                                         cAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGa                                                            <  2:4035102/150‑1 (MQ=255)
                                                                                           ggCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACa                                                          >  2:3948111/1‑150 (MQ=255)
                                                                                             cAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGt                                                        >  2:3780733/1‑150 (MQ=255)
                                                                                              aaaaCATCTGAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTa                                                       <  2:41417/150‑1 (MQ=255)
                                                                                               aaaCATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGTGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTAACAGGCGAATCACTTTACACGGAACCAAATCACCACACAGTAt                                                      >  2:2529387/1‑150 (MQ=255)
                                                                                               aaaCATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTAt                                                      >  2:4082788/1‑150 (MQ=255)
                                                                                               aaaCATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTAt                                                      <  1:3270083/150‑1 (MQ=255)
                                                                                               aaaCATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTAt                                                      >  2:2643726/1‑150 (MQ=255)
                                                                                                 aCATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTAttt                                                    <  1:1761301/150‑1 (MQ=255)
                                                                                                      ttAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCtgt                                               >  1:1728415/1‑150 (MQ=255)
                                                                                                      ttAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCtgt                                               <  1:3556397/150‑1 (MQ=255)
                                                                                                      ttAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCtgt                                               <  1:809363/150‑1 (MQ=255)
                                                                                                      ttAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCtgt                                               <  1:4093780/150‑1 (MQ=255)
                                                                                                      ttAAATATAGTCTTGTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCtgt                                               <  1:2333771/150‑1 (MQ=255)
                                                                                                       tAAATATAGTCTTTTCCGTCTACCTGATAGACAAAAACGAGCCCCGCTGGGCTCGTTTTATCATGTACTGGAGGGCGACAATACGGTCTTCATTAACTTCCTGGCGAATCACTTTACCCGGAACCAATTCACCAGTCTGTATTTGCtgtg                                              <  2:587726/150‑1 (MQ=255)
                                                                                                        aaaTATAGTCTTGTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGc                                             <  1:3547454/150‑1 (MQ=255)
                                                                                                         aaTATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGcc                                            >  2:345793/1‑150 (MQ=255)
                                                                                                         aaTATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGcc                                            >  2:3890038/1‑150 (MQ=255)
                                                                                                           tataGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAg                                          <  1:3939528/150‑1 (MQ=255)
                                                                                                             taGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAAACCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCg                                        >  1:2709611/1‑150 (MQ=255)
                                                                                                               gtctggtCCGTCTACCTTATAGACAAAAACGAGCCCCGACGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCggg                                      <  2:3293588/144‑1 (MQ=255)
                                                                                                               gTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCggg                                      >  1:1070767/1‑150 (MQ=255)
                                                                                                                tCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGt                                     >  1:1885501/1‑150 (MQ=255)
                                                                                                                  ttgaccGTCTAACGTATAGACAAAAACGAGCCCCGCCGGGCTCGTTTTATCATTTACTGGCCGGCGACAATCCGGTCTTCATTAACTTCCAGGCGCATCACTTTACCCGGAACCCATTCACCAGACAGTATTTGCTGTGCCAGCGGGttt                                   <  2:2741422/146‑1 (MQ=255)
                                                                                                                   tttCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGtttt                                  >  1:2471280/1‑150 (MQ=255)
                                                                                                                   tttCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGtttt                                  >  1:773459/1‑150 (MQ=255)
                                                                                                                    ttCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTc                                 >  2:802882/1‑150 (MQ=255)
                                                                                                                      ccGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGa                               <  1:3675254/150‑1 (MQ=255)
                                                                                                                       cGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGAt                              <  2:3161760/150‑1 (MQ=255)
                                                                                                                       cGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGAt                              <  1:573949/150‑1 (MQ=255)
                                                                                                                        gTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATc                             <  1:3837196/150‑1 (MQ=255)
                                                                                                                               ttATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATctgctgct                      >  1:2182761/1‑150 (MQ=255)
                                                                                                                               gtataGACAAACACGAGCCCCGACGGGCTCGTTTTATCATTGACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACGCTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATctgctgct                      <  2:2919114/149‑1 (MQ=255)
                                                                                                                                   aGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAAt                  <  2:435379/150‑1 (MQ=255)
                                                                                                                                      cAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGc               >  2:3637371/1‑150 (MQ=255)
                                                                                                                                       taaaaCGAGCCCTTCGGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGCa              >  2:1378985/2‑150 (MQ=255)
                                                                                                                                       taaaaCGAGCCCTTCGGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGCa              >  2:3592142/2‑150 (MQ=255)
                                                                                                                                       aaaaaCGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGCa              <  2:2211139/150‑1 (MQ=255)
                                                                                                                                       aaaaaCGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGATGTGCCAGCGGGTTTTCGAACTGCTGCTGAATTGCa              >  1:4084286/1‑150 (MQ=255)
                                                                                                                                        aaaaCGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGCAc             >  2:667558/1‑150 (MQ=255)
                                                                                                                                           aCGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGCACGtt          >  1:1023044/1‑150 (MQ=255)
                                                                                                                                             gAGCCCCGAAGGGCTCGTTTTAACATTTACTTGACGGCGACAATCCGGTCTTCAATAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCAGCAGCTGAATTGCACGTTTc        >  1:3466736/1‑150 (MQ=255)
                                                                                                                                              aGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGCACGTTTCa       <  1:1539390/150‑1 (MQ=255)
                                                                                                                                               gCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGCACGTTTCag      >  2:3393212/1‑150 (MQ=255)
                                                                                                                                                ccccGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGCACGTTTCaga     >  2:158174/1‑150 (MQ=255)
                                                                                                                                                ccccGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGCACGTTTCaga     <  1:1131230/150‑1 (MQ=255)
                                                                                                                                                 cccGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGCACGTTTCagag    >  2:2239891/1‑150 (MQ=255)
                                                                                                                                                 cccGAAGGGCTCGTTTTATCAAATACATGACGGCGACAATCCCGTCTTAAATAACTTCCACGCGAAACACTTTACCACCAAACAAATCACCACACAGTAATTGCTGTTCCAGCGGGATTTCTATATTCTGCTGAAATGCACGTTTCatag    >  1:1659651/1‑150 (MQ=255)
                                                                                                                                                   cGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGCACGTTTCAGAGGa  <  2:2266742/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GTTGTTCCGTGTCAGTGGTGGCGCATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGCACGTTTCAGAGGA  >  NC_000913/2731433‑2731729

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: