Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,731,581 G→T 5.1% intergenic (‑424/+19) 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,5810GT 5.1% 1063.0 / 25.2 294intergenic (‑424/+19)rrsG/clpB16S ribosomal RNA of rrnG operon/protein disaggregation chaperone
Reads supporting (aligned to +/- strand):  ref base G (160/117);  new base T (8/9);  total (168/126)
Fisher's exact test for biased strand distribution p-value = 4.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 295 total aligned reads displayed.

GTTGTTCCGTGTCAGTGGTGGCGCATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGCACGTTTCAGAGGAC  >  NC_000913/2731433‑2731730
                                                                                                                                                    |                                                                                                                                                     
gttgttCCGTGTCAGTGGTGGCGCATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGa                                                                                                                                                      >  2:4221118/1‑150 (MQ=255)
 ttggtCCGTGTCAGTGGTGGCGCATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGaa                                                                                                                                                     <  2:837379/150‑1 (MQ=255)
        gtgtCTGTGGTGGCGCATTATAGGGAGTTATTCCGGCCTGACAAGCTAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCg                                                                                                                                              <  1:4301511/150‑1 (MQ=255)
        gtgtCAGTGGTGGCGCATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCg                                                                                                                                              >  1:3570156/1‑150 (MQ=255)
            cAGTGGTGGCGCATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGtttt                                                                                                                                          <  2:3972002/150‑1 (MQ=255)
             aGTGGTGGCGCATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTa                                                                                                                                         <  1:4830776/150‑1 (MQ=255)
              gtggtgGCGCATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTAt                                                                                                                                        <  1:3946451/150‑1 (MQ=255)
              gtggtgGCGCATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACATTTCAGGCAAAACCTCGAAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGATTTAt                                                                                                                                        <  1:4450164/150‑1 (MQ=39)
               tggtggCGCATTTTAGGGAGTTATTCCGGCCTGACAATCTAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGCAAAAAAATCTTAAATATAGTCTTTTCCGTCTAACTTATAGAAAAAAACGAGCCCTTCGGGGCTCGtttttgt                                                                                                                                       <  1:1941058/150‑4 (MQ=255)
               tggtggCGCATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCTTCGGGGCTCGtttttgt                                                                                                                                       <  2:3378469/150‑4 (MQ=39)
               tggtggCGCATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCTTCGGGGCTCGtttttgt                                                                                                                                       <  2:3688119/150‑4 (MQ=39)
               tggtggCGCATTATAGGGAGTTAGTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCTTCGGGGCTCGtttttgt                                                                                                                                       <  2:508843/150‑4 (MQ=255)
                    gcgcATTATAGGGAGTTATTCAGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCTAAGGGCTCGTTTTATCATTTa                                                                                                                                  >  2:1056832/1‑150 (MQ=255)
                     cgcATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTAc                                                                                                                                 <  1:927378/150‑1 (MQ=255)
                      gcATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACt                                                                                                                                >  1:4621983/1‑150 (MQ=255)
                       cATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTg                                                                                                                               >  2:939513/1‑150 (MQ=255)
                       cATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTg                                                                                                                               <  1:1718103/150‑1 (MQ=255)
                         ttATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGa                                                                                                                             >  2:290376/1‑150 (MQ=255)
                            taGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACgg                                                                                                                          >  2:3801297/1‑150 (MQ=255)
                              gggAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCg                                                                                                                        >  1:689849/1‑150 (MQ=255)
                               ggAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGa                                                                                                                       >  2:4233334/1‑150 (MQ=255)
                               ggAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGa                                                                                                                       <  1:3822366/150‑1 (MQ=255)
                                   ttATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAAt                                                                                                                   >  1:2257055/1‑150 (MQ=255)
                                   ttATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAAt                                                                                                                   >  2:485246/1‑150 (MQ=255)
                                     aTTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATcc                                                                                                                 <  2:4107522/150‑1 (MQ=255)
                                      ttCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCg                                                                                                                <  2:4983772/150‑1 (MQ=255)
                                      gtCCGGCCTGACAAGCGAAAAATATAAAAACTTGATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCg                                                                                                                <  1:5360267/149‑1 (MQ=255)
                                       tCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCgg                                                                                                               <  2:2584789/150‑1 (MQ=255)
                                         cGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTc                                                                                                             >  2:1268987/1‑150 (MQ=255)
                                           gCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCtt                                                                                                           >  2:2509462/1‑150 (MQ=255)
                                            ccTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTc                                                                                                          >  2:3811885/1‑150 (MQ=255)
                                             cTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCa                                                                                                         <  2:2333762/150‑1 (MQ=255)
                                                  aaGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAAc                                                                                                    <  2:3906471/150‑1 (MQ=255)
                                                     cGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTc                                                                                                 >  2:1471920/1‑150 (MQ=255)
                                                      gAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTcc                                                                                                >  1:1594177/1‑150 (MQ=255)
                                                       aaaaaTATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCa                                                                                               <  1:2122299/150‑1 (MQ=255)
                                                       aaaaaTATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCa                                                                                               >  1:1962766/1‑150 (MQ=255)
                                                            tataAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGa                                                                                          <  1:451230/150‑1 (MQ=255)
                                                              taAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAAt                                                                                        >  1:2881428/1‑150 (MQ=255)
                                                              taAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAAt                                                                                        >  2:1021568/1‑150 (MQ=255)
                                                                 aaaCTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCAc                                                                                     <  2:3572593/150‑1 (MQ=255)
                                                                     tttATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTa                                                                                 >  1:4055560/1‑150 (MQ=255)
                                                                               gCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCa                                                                       >  2:4296269/1‑150 (MQ=255)
                                                                               gCTCACTTTTAAGGCAAAACATCTTAAATAAATTCTTTTCCGTCTAACTTAAAGACAAAAAATAGCCACGAACGGCTCGTTTTATCATTTACTGGACGGCGAAAATCCGGTCTTCAATAACTACCAGGCGAATCACTTTACCCGGAACCa                                                                       >  2:3512384/1‑150 (MQ=255)
                                                                                 tCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAAt                                                                     >  1:1358604/1‑150 (MQ=255)
                                                                                   aCTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCGGCCGAACACATTACCCCGAACCaaatac                                                                  >  1:3231202/1‑146 (MQ=255)
                                                                                    cTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCa                                                                  >  1:1118008/1‑150 (MQ=255)
                                                                                    cTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCa                                                                  <  1:2791876/150‑1 (MQ=255)
                                                                                     ttttCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCAc                                                                 <  2:4973975/150‑1 (MQ=255)
                                                                                      tttCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCAcc                                                                >  2:3120715/1‑150 (MQ=255)
                                                                                         cAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGa                                                             <  2:2257055/150‑1 (MQ=255)
                                                                                          aGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGAc                                                            >  1:4275542/1‑150 (MQ=255)
                                                                                          aGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGAc                                                            >  1:3493520/1‑150 (MQ=255)
                                                                                            gCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAg                                                          <  2:5392697/150‑1 (MQ=255)
                                                                                              aaaaCATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTa                                                        >  2:3782413/1‑150 (MQ=255)
                                                                                                aaCATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTAtt                                                      <  2:5465007/150‑1 (MQ=255)
                                                                                                      ttAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCCCCAGACAGTATTTGCtgt                                                <  2:3847079/150‑1 (MQ=255)
                                                                                                      ttAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCtgt                                                >  1:2330774/1‑150 (MQ=255)
                                                                                                      ttAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCtgt                                                <  2:1898183/150‑1 (MQ=255)
                                                                                                      ttAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCtgt                                                <  1:1021568/150‑1 (MQ=255)
                                                                                                             taGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATACACCAGACAGTATTTGCTGTGCCAGCg                                         >  2:5380391/1‑150 (MQ=255)
                                                                                                                tCTTTTACGTCTAACTTAAAGAAAAAAACGAGCACCTAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATAAACTTCAAGGCGAATCACTTTACACGGAACCAATTCACAAGACAGTAATTGCTGTGCCAGAGGGt                                      >  2:1190398/1‑150 (MQ=255)
                                                                                                                  ttttCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGttt                                    >  1:4120126/1‑150 (MQ=255)
                                                                                                                  ttttCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGttt                                    >  2:1910729/1‑150 (MQ=255)
                                                                                                                  ttttCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGttt                                    >  2:894179/1‑150 (MQ=255)
                                                                                                                  ttttCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGttt                                    <  2:4682821/150‑1 (MQ=255)
                                                                                                                  ttttCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGGTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGttt                                    >  2:3342930/1‑150 (MQ=255)
                                                                                                                   tttCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGtttt                                   >  1:5150323/1‑150 (MQ=255)
                                                                                                                      ccGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGa                                <  2:1435782/150‑1 (MQ=255)
                                                                                                                      ccGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGGATTTGCTGTGCCAGCGGGTTTTCGa                                >  1:1547264/1‑150 (MQ=255)
                                                                                                                      ccGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATGTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGa                                <  1:735604/150‑1 (MQ=255)
                                                                                                                       cGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGAt                               <  2:925418/150‑1 (MQ=255)
                                                                                                                       cGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGAt                               <  2:2400619/150‑1 (MQ=255)
                                                                                                                          cTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATctg                            <  1:542788/150‑1 (MQ=255)
                                                                                                                             aCTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATctgctg                         >  2:3728053/1‑150 (MQ=255)
                                                                                                                             aCTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATctgctg                         <  2:1362620/150‑1 (MQ=255)
                                                                                                                                   aGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAAt                   <  2:5208102/150‑1 (MQ=255)
                                                                                                                                   aGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAAt                   >  1:3117651/1‑150 (MQ=255)
                                                                                                                                     aCAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTg                 >  2:1844378/1‑150 (MQ=255)
                                                                                                                                     aCAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTg                 <  2:1962766/150‑1 (MQ=255)
                                                                                                                                     aCAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTg                 >  1:538529/1‑150 (MQ=255)
                                                                                                                                     aCAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTg                 >  1:3269027/1‑150 (MQ=255)
                                                                                                                                      cAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGc                >  2:1654970/1‑150 (MQ=255)
                                                                                                                                      cAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGc                >  2:4398600/1‑150 (MQ=255)
                                                                                                                                      cAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGc                <  1:1222913/150‑1 (MQ=255)
                                                                                                                                       taaaaCGAGCCCTTCGGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGCa               >  1:4703767/2‑150 (MQ=255)
                                                                                                                                         aaaCGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGCACg             >  2:678535/1‑150 (MQ=255)
                                                                                                                                          aaCGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGCACGt            >  2:2529600/1‑150 (MQ=255)
                                                                                                                                             gAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGCACGTTTc         >  1:3906156/1‑150 (MQ=255)
                                                                                                                                              aGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGCACGTTTCa        >  1:3800879/1‑150 (MQ=255)
                                                                                                                                                ccccGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGCACGTTTCagt      >  1:3423326/1‑149 (MQ=255)
                                                                                                                                                ccccGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGCACGTTTCaga      <  1:3277806/150‑1 (MQ=255)
                                                                                                                                                ccccGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGCACGTTTCaga      >  1:611975/1‑150 (MQ=255)
                                                                                                                                                ccccGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGCACGTTTCaga      >  2:2829437/1‑150 (MQ=255)
                                                                                                                                                ccccGAAGGGCGCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCCGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGCACGTTTCaga      >  2:4432781/1‑150 (MQ=255)
                                                                                                                                                 cccGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGCACGTTTCagag     >  1:1466176/1‑150 (MQ=255)
                                                                                                                                                  ccGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGCACGTTTCAGAgg    >  2:58040/1‑150 (MQ=255)
                                                                                                                                                   cGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGCACGTTTCAGAGGa   >  1:83737/1‑150 (MQ=255)
                                                                                                                                                   cGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGCACGTTTCAGAGGa   >  2:3637868/1‑150 (MQ=255)
                                                                                                                                                    gAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGCACGTTTCAGAGGAc  <  2:2063277/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GTTGTTCCGTGTCAGTGGTGGCGCATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGCACGTTTCAGAGGAC  >  NC_000913/2731433‑2731730

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: