Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 454,873 C→G 7.6% intergenic (+84/‑260) bolA → / → tig stationary‑phase morphogene, transcriptional repressor for mreB; also regulator for dacA, dacC, and ampC/peptidyl‑prolyl cis/trans isomerase (trigger factor)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913454,8730CG 7.6% 1399.3 / 15.1 440intergenic (+84/‑260)bolA/tigstationary‑phase morphogene, transcriptional repressor for mreB; also regulator for dacA, dacC, and ampC/peptidyl‑prolyl cis/trans isomerase (trigger factor)
Reads supporting (aligned to +/- strand):  ref base C (182/213);  new base G (35/2);  total (225/215)
Fisher's exact test for biased strand distribution p-value = 2.72e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 452 total aligned reads displayed.

TGGGAAGGGTTGCAGGACACCGTCTTTGCCTCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAACCGTTCTCGACTCATAAAAGTGATGCCGCTATAATGCCGCGTCCTATTTGAATGCTTTCGGGATGATTCT  >  NC_000913/454724‑455022
                                                                                                                                                     |                                                                                                                                                     
tGGGAAGGGTTGCAGGACACCGTCTTTGCCTCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTc                                                                                                                                                       <  2:5648266/150‑1 (MQ=255)
ggggACGGGTTGCCGGACACCGTCTTTGCCTCTCCGCCCGGGCGTGGCGCCGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTc                                                                                                                                                       <  1:6999941/149‑1 (MQ=255)
  ggAAGGGTTGCAGGACACCGTCTTTGCCTCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTg                                                                                                                                                     >  1:481401/1‑150 (MQ=255)
  ggAAGGGTTGCAGGACACCGTCTTTGCCTCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTg                                                                                                                                                     >  2:4374172/1‑150 (MQ=255)
  ggAAGGGTTGCAGGACACCGTCTTTGCCTCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTGTTGTTTTGTCTg                                                                                                                                                     >  2:102215/1‑150 (MQ=255)
   gAAGGGTTGCAGGACACCGTCTTTGCCTCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTgg                                                                                                                                                    >  1:7162293/1‑150 (MQ=255)
   gAAGGGTTGCAGGACACCGTCTTTGCCTCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTgg                                                                                                                                                    >  2:354014/1‑150 (MQ=255)
        gTTGCAGGACACCGTCTTTGCCTCTCCTCCCTGTCGTGGAGCAGGAAGCATCTCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGAtttt                                                                                                                                               <  2:834338/150‑1 (MQ=255)
         ttGCAGGACACCGTCTTTGCCTCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTg                                                                                                                                              <  2:4654818/150‑1 (MQ=255)
         ttGCAGGACACCGTCTTTGCCTCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTGTTGTTTTGTCTGGATTTTg                                                                                                                                              >  1:4931775/1‑150 (MQ=255)
              ggACACCGTCTTTGCCTCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGGTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTTTTTTGTGTGGAtttttgggtt                                                                                                                                         >  1:2142642/1‑144 (MQ=255)
               gACACCGTCTTTGCCTCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACTCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCttt                                                                                                                                        >  1:3507508/1‑150 (MQ=255)
               gACACCGTCTTTGCCTCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTTTGTGGATTTTGCGCttt                                                                                                                                        >  1:2752516/1‑150 (MQ=255)
               gACACCGTCTTTGCCTCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGGGGTGTTTTGTCTGGAGTTTGCGCttt                                                                                                                                        >  2:3906254/1‑150 (MQ=255)
                 cacCGTCTTTGCCTCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCttttt                                                                                                                                      <  2:2685950/150‑1 (MQ=255)
                 cacCGTCTTTGCCTCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCttttt                                                                                                                                      <  2:6825346/150‑1 (MQ=255)
                  acCGTCTTTGCCTCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGGTTTTGCGCTTTTTg                                                                                                                                     >  1:1400062/1‑150 (MQ=255)
                  acCGTCTTTGCCTCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGGATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGGTTTGTCTGGATTTTTCGCTTTTTg                                                                                                                                     >  2:3277236/1‑150 (MQ=255)
                   ccGTCTTTGCCTCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGc                                                                                                                                    <  1:1747928/150‑1 (MQ=255)
                     gTCTTTGCCTCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGccc                                                                                                                                  <  2:1059589/150‑1 (MQ=255)
                         ttGCCTCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCa                                                                                                                              <  2:3378345/150‑1 (MQ=255)
                            ccTCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTTTCTGGATTTTGCGCTTTTTTACCAGCAtta                                                                                                                           >  2:1878734/1‑149 (MQ=255)
                            ccTCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTGTTGTTTTGTCTGGATTTTGCGCTATTTTCCCAGCAATc                                                                                                                           >  2:1267082/1‑150 (MQ=255)
                             ctctCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTTTTTTTTCTGGATTTTGCGCTTTTTGCCCAGCATTCa                                                                                                                          >  1:469068/1‑150 (MQ=255)
                             ctctCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCa                                                                                                                          >  2:3941202/1‑150 (MQ=255)
                             ctctCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGGTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCa                                                                                                                          >  1:2349530/1‑150 (MQ=255)
                              tctcCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGGAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAg                                                                                                                         >  1:6416464/1‑150 (MQ=255)
                              tctcCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAg                                                                                                                         <  2:649080/150‑1 (MQ=255)
                               ctcctcCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGa                                                                                                                        >  2:1152776/1‑150 (MQ=255)
                               ctcctcCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTGTTGTTTTGTGTGGATTTTTCGCTTTTTGCCCAGCAATAAGa                                                                                                                        >  1:4866377/1‑150 (MQ=255)
                                tcctccCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGAc                                                                                                                       >  1:1139612/1‑150 (MQ=255)
                                tcctccCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGGGGTGGTTTGTCTGGATTTTGCGCTTTTTGCCCCGCATTCCGAc                                                                                                                       >  2:6237068/1‑150 (MQ=255)
                                   tccCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGaa                                                                                                                    <  1:4187511/150‑1 (MQ=255)
                                   tccCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGaa                                                                                                                    >  1:4021922/1‑150 (MQ=255)
                                    cccTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTGGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGaaa                                                                                                                   <  2:60714/150‑1 (MQ=255)
                                     ccTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTACGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGaaaa                                                                                                                  <  2:2114112/150‑1 (MQ=255)
                                       tGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTTTTTTGTCTGGTTTTTTCGCTTTTTTCACAGCATTAAGACGaaaaat                                                                                                                >  1:1288990/1‑148 (MQ=255)
                                       tGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTTAGACGaaaaat                                                                                                                >  1:1576323/1‑148 (MQ=255)
                                         tCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCCGACGAAAAtttc                                                                                                              >  1:5600018/1‑148 (MQ=255)
                                         tCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGc                                                                                                              <  1:6149154/150‑1 (MQ=255)
                                          cGTGGAGCTGGATGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGcc                                                                                                             <  1:5504470/150‑1 (MQ=255)
                                           gTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGccc                                                                                                            >  2:1371351/1‑150 (MQ=255)
                                            tGGAGCAGTAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCg                                                                                                           <  1:4355804/150‑1 (MQ=255)
                                            tGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTTCCCg                                                                                                           >  1:3289055/1‑150 (MQ=255)
                                            tGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTGTTGTTTTGTGTGGATTTTGCGCTTTTTGCCCAGCATTCAGGCGGAAATTGCCCg                                                                                                           >  1:681613/1‑150 (MQ=255)
                                             ggAGCAGGAAGCATCGCTTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCgg                                                                                                          <  1:2554787/150‑1 (MQ=255)
                                                gCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGaa                                                                                                       <  2:5388667/150‑1 (MQ=255)
                                                 cAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTCTGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCCGCATTCAGACGAAAATTGCCCGGGAAt                                                                                                      >  1:3918741/1‑150 (MQ=255)
                                                 cAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAAt                                                                                                      <  2:5996067/150‑1 (MQ=255)
                                                  aGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAAtt                                                                                                     <  1:5298771/150‑1 (MQ=255)
                                                  aGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTGTTGTTTTGTCTGGATTTTGCGCTTTTTTCCCAGCATTCAGACGAAAATTGCCCGGGAAtt                                                                                                     >  1:3756626/1‑150 (MQ=255)
                                                    gAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTTTGTGGTTTTTGCGCTTTTTGGCCAGCATTCAGACGAAAATTGCCCGGGAATtgt                                                                                                   >  2:205283/1‑150 (MQ=255)
                                                          tCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGaaaaa                                                                                             >  2:6337920/1‑150 (MQ=255)
                                                           cgcgTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGaaaaaa                                                                                            <  2:6728203/150‑1 (MQ=255)
                                                               tAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTGGGGTTTGTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGGCCGGGGATTGTTTAAAAATaag                                                                                        >  1:5975154/1‑148 (MQ=255)
                                                                  aaaCGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGa                                                                                     <  2:6588570/150‑1 (MQ=255)
                                                                    aCGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATTAAAAACGGCCTGCGGGCCGTTTTGTTTTTTGTGGATTTTGCGCTTTTTGGCCACCATTAAGACGAAAATTGCCCGGGAATTTTTAAAAAATACGCTACa                                                                                   >  2:3079634/1‑150 (MQ=255)
                                                                    aCGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTTCCCAGCATTAAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACa                                                                                   >  2:1988422/1‑150 (MQ=255)
                                                                    aCGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACa                                                                                   >  2:2795547/1‑150 (MQ=255)
                                                                    aCGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGCAATTGTGAAAAAATACGCGACa                                                                                   >  2:5173705/1‑150 (MQ=255)
                                                                         tttGCAACGGTGGGCGCTTTTCCTGTTTGTTTCTAAATTTTTTATGAAAAAGGGCGGGGGGGCTTTTTTTTTTTGTCTGTTTTTTGCGCTTTTTGCCCAGCTTTCAGACTAAATTTGCCGGGTATTTTTGAAAAAATACGCGACAgcgcg                                                                              <  1:652502/150‑1 (MQ=255)
                                                                           tGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCa                                                                            <  2:481401/150‑1 (MQ=255)
                                                                              aaCTGTCGGCGCTTTTCCAGTATTTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTATTGATTTGTCTTGATTTATCGCTTTTTGCCCACCATTAAGACAAAAAAAGCAAGGGAAATTTGAAAAAATACACGAAAGAGCCCAATa                                                                         >  2:6436435/1‑150 (MQ=255)
                                                                                cTGTGGGCGATTTTCCTGTATGTTGATAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTTTTTTGTCTGTATTTTGCGCTTTTTGCCCAGCTTTCAGACTAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAAc                                                                       <  1:1651514/150‑1 (MQ=255)
                                                                                cTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAAc                                                                       <  2:1714391/150‑1 (MQ=255)
                                                                                 tGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAAcc                                                                      <  2:2509323/150‑1 (MQ=255)
                                                                                  gtgggCGCTTTTCCCGTTTGTTGCTAATGATTTTTTGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAACCg                                                                     <  1:5799912/147‑1 (MQ=255)
                                                                                  gTCGGCGCTTTTCCAGTATTTTGATAATTATTTTATGAAAAACGGCCTGGGGGCCTTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGATTTGTGAAAAAATACGCGACAGCGCGCAATAACCg                                                                     <  1:540397/150‑1 (MQ=255)
                                                                                  gTCGGCGCTTTTCCAGTATGTTGCTAATTATTTTATTAAAAACGGCCTGGGGGCGTTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACTAAATTTGCCCGGGATTTGTTAAAAAATACGCGACAGCGCGCAATAACCg                                                                     <  1:669306/150‑1 (MQ=255)
                                                                                  gTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAACCg                                                                     <  2:6914419/150‑1 (MQ=255)
                                                                                     ggCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAACCGTtc                                                                  >  2:3750705/1‑150 (MQ=255)
                                                                                      gcgcTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAACCGTtcg                                                                 >  2:260737/1‑149 (MQ=255)
                                                                                         cTTTTCCATTCTTTTGCTATAGTTTTTATGAAAAACGGCCTGGGGGCTTTTTTTTTTTGTCTTTATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAACCGTTCTCGa                                                              <  1:1878734/150‑1 (MQ=255)
                                                                                         cTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTTTCTGGATTTTGCGCTTTTTTCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAAGAACCGtgcgcgc                                                              >  2:5571336/1‑144 (MQ=255)
                                                                                         cTTTTCCAGTATGTTGATAAAGATTTTATTAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAACCGTTCTCGa                                                              <  1:4542591/150‑1 (MQ=255)
                                                                                          ttttCCCGTAGGTTGCTACAGATTTTATGAAAAACGGCCTGCGGGCCGTTGTGTTTTGTCTGGATGTTGCGCTTTTTGCCCAGCATTCGGACGACGATGGCCCGGGAGTTGTGTGTATATTCGCGTCAGGGCGCGATGGCCGGTCGCGAc                                                             >  2:4439598/1‑150 (MQ=255)
                                                                                          ttttCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAACCGTTCTCGAc                                                             >  2:1071313/1‑150 (MQ=255)
                                                                                           tttCCAGTATGTTGCTAAAGATTTTATGAAAATCGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAACCGTTCTCGACt                                                            <  1:1008606/150‑1 (MQ=255)
                                                                                             tCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAACCGTTCTCGACTCa                                                          >  2:4218091/1‑150 (MQ=255)
                                                                                               cAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAACCGTTCTCGACTCATa                                                        >  2:2422919/1‑150 (MQ=255)
                                                                                                aGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAACCGTTCTCGACTCATaa                                                       <  2:3148941/150‑1 (MQ=255)
                                                                                                  tATGTTGCAAAAGATTTTATGAAAAACGGCCTGCGGGCCGATTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAACCGTTCTCGACTCATaaaa                                                     <  1:3866793/150‑1 (MQ=255)
                                                                                                       tGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAACCGTTCTCGACTCATAAAAGTGAt                                                <  1:769306/150‑1 (MQ=255)
                                                                                                         cTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAACCGTTCTCGACTCATAAAAGTGATGc                                              <  2:1476188/150‑1 (MQ=255)
                                                                                                          tAAATATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAACCGTTCTCGACTCATAAAAGTGATGcc                                             <  2:5423960/150‑1 (MQ=255)
                                                                                                           aaaGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGACTTGTGAAAAAATACGCGACAGCGCGCAATAACCGTTCTCGACTCATAAAAGTGATGCCg                                            >  2:7004350/1‑150 (MQ=255)
                                                                                                                          aaCGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAACCGTTCTCGACTCATAAAAGTGATGCCGCTATAATGCCGCGTc                             <  1:6049003/150‑1 (MQ=255)
                                                                                                                                cTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAACCGTTCTCGACTCATAAAAGTGATGCCGCTATAATGCCGCGTCCTAttt                       <  1:2985324/150‑1 (MQ=255)
                                                                                                                                   cGGGCCGTTTTGTTTTGTCTGGATTTTGCGATTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGCATTGTGAAAAAATACGCGACAGCGCGCAATAACCGTTCTCGACTCATAAAAGTGATGCCGCTATAATGCCGCGTCCTATTTGaa                    <  1:2657390/150‑1 (MQ=255)
                                                                                                                                    gggCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAACCGTTCTCGACTCATAAAAGTGATGCCGCTATAATGCCGCGTCCTATTTGAAt                   <  2:507831/150‑1 (MQ=255)
                                                                                                                                    gggCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAACCGTTCTCGACTCATAAAAGTGATGCCGCTATAATGCCGCGTCCTATTTGAAt                   <  2:4078901/150‑1 (MQ=255)
                                                                                                                                          ttttgttttgTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAACCGTTCTCGACTCATAAAAGTGATGCCGCTATAATGCCGCGTCCTATTTGAATGCTTTc             >  2:2927709/1‑150 (MQ=255)
                                                                                                                                           tttgttttgtCTGGTTTTTGCGCTTTTTGCCCAGCTTTCAGTCGAAAATTGCCCGTGAATTGTGAAAAAATACGCGACAGCGCGCAATTACCGTTCTCGACTCATAAAAGTGTTGCCGCTATTATGCCGCGTCCTATTTGAATGCTTTCg            <  1:2950811/150‑1 (MQ=255)
                                                                                                                                            ttgttttgtCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAACCGTTCTCGACTCATAAAAGTGATGCCGCTATAATGCCGCGTCCTATTTGAATGCTTTCgg           <  1:1910715/150‑1 (MQ=255)
                                                                                                                                              gttttgtCTGTATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAACCGTTCTCGACTCATAAAAGTGATGCCGCTATAATGCCGCGTCCTATTTGAATGCTTTCGGga         <  1:1260764/150‑1 (MQ=255)
                                                                                                                                              gttttgtCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAACCGTTCTCGACTCATAAAAGTGATGCCGCTATAATGCCGCGTCCTATTTGAATGCTTTCGGga         <  2:4163350/150‑1 (MQ=255)
                                                                                                                                              gttttgtCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAACCGTTCTCGACTCATAAAAGTGATGCCGCTATAATGCCGCGTCCTATTTGAATGCTTTCGGga         <  1:1964146/150‑1 (MQ=255)
                                                                                                                                               ttttgtCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAATTTGCCCGGGAATGTTGAAAAAATACGCGACAGCGCGCAATAACCGTTCTCGACTCATAAAAGTGATGCCGCTATAATGCCGCGTCCTATTTGAATGCTTTCGGgat        <  2:6099004/150‑1 (MQ=255)
                                                                                                                                               ttttgtCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAACCGTTCTCGACTCATAAAAGTGATGCCGCTATAATGCCGCGTCCTATTTGAATGCTTTCGGgat        >  2:2307861/1‑150 (MQ=255)
                                                                                                                                               ttttgtCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAACCGTTCTCGACTCATAAAAGTGATGCCGCTATAATGCCGCGTCCTATTTGAATGCTTTCGGgat        >  2:1806023/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGGGAAGGGTTGCAGGACACCGTCTTTGCCTCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAACCGTTCTCGACTCATAAAAGTGATGCCGCTATAATGCCGCGTCCTATTTGAATGCTTTCGGGATGATTCT  >  NC_000913/454724‑455022

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: ‑