Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 267,472 C→G 5.4% P178P (CCG→CCC yafZ ← CP4‑6 prophage; conserved protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913267,4720CG 5.4% 2556.1 / 22.3 748P178P (CCG→CCCyafZCP4‑6 prophage; conserved protein
Reads supporting (aligned to +/- strand):  ref base C (403/292);  new base G (2/51);  total (405/343)
Fisher's exact test for biased strand distribution p-value = 2.94e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 757 total aligned reads displayed.

GGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCACCTCATACGCGCCTTCAATCACCTGACTCACCACGTCCCCCTTGTGTGGCACCCGCACCTCGCCAAACGACTCACCGCAGAC  >  NC_000913/267324‑267621
                                                                                                                                                    |                                                                                                                                                     
     tcCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGt                                                                                                                                                 >  2:1102538/1‑150 (MQ=255)
      ccTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTg                                                                                                                                                <  1:462611/150‑1 (MQ=255)
         gAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCa                                                                                                                                             >  1:5159149/1‑150 (MQ=255)
         gAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCa                                                                                                                                             >  2:1316556/1‑150 (MQ=255)
          aaTACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCaa                                                                                                                                            >  1:2055270/1‑150 (MQ=255)
            tACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACa                                                                                                                                          <  1:10057335/150‑1 (MQ=255)
                cTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCaa                                                                                                                                      <  1:8797790/150‑1 (MQ=255)
                 tGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAAc                                                                                                                                     <  2:6896777/150‑1 (MQ=255)
                 tGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAAc                                                                                                                                     >  1:8224710/1‑150 (MQ=255)
                 tGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAAc                                                                                                                                     >  1:7281578/1‑150 (MQ=255)
                   gTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGa                                                                                                                                   <  2:4172249/150‑1 (MQ=255)
                     aCGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACt                                                                                                                                 >  2:6120137/1‑150 (MQ=255)
                     aCGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACt                                                                                                                                 <  2:2729482/150‑1 (MQ=255)
                        tGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGAGTGCa                                                                                                                              >  2:6659175/1‑150 (MQ=255)
                          gTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATg                                                                                                                            >  2:4838517/1‑150 (MQ=255)
                            ccACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGc                                                                                                                          >  1:6095619/1‑150 (MQ=255)
                            ccACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGc                                                                                                                          >  1:3646212/1‑150 (MQ=255)
                               caGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAACGCGTTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCACCAGCAACGACTGCATGGCATc                                                                                                                       <  2:10742391/150‑1 (MQ=255)
                                 ggTCATTGCTCTCATCCAGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGACGGGGGTGGCAACAGCAACGACTGCATCGCATccc                                                                                                                     >  1:7096032/1‑150 (MQ=255)
                                         ctctCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTctcc                                                                                                             <  2:5972211/150‑1 (MQ=255)
                                          tctcATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTctcct                                                                                                            >  2:10050026/1‑150 (MQ=255)
                                          tctcATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTctcct                                                                                                            >  1:5010090/1‑150 (MQ=255)
                                          tctcATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTctcct                                                                                                            >  2:7483146/1‑150 (MQ=255)
                                          tctcATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTctcct                                                                                                            >  1:2907487/1‑150 (MQ=255)
                                          tctaATCCTGCCAGCGTCTAGGGGATAGGATCTGCGATTCAGTCACCGGCAGGTGGTCTTCACCAAAGCGGTATGTGAGTGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTctcct                                                                                                            >  2:3614577/1‑150 (MQ=21)
                                           ctcATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCCCCGGCTGGTGGTCTTCACCAAAGCGGTCTGTGAGGGCGGCTTTTGCCAGGGCCTGCTGTGCCGGGGGTGGCACCAGCCCCGACTGCATGGCATCCCGTTTctccgc                                                                                                           >  2:9782137/1‑148 (MQ=255)
                                           ctcATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTctcctc                                                                                                           <  2:4584175/150‑1 (MQ=255)
                                                ccTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCGTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCACCGACTGCATGGCATCCCGTGTCTCCTCCAccc                                                                                                      >  2:3873616/1‑150 (MQ=255)
                                                 ctgcctgcGGGGTGGGTCGAGGATCTGCTATTCAGCCCCGGGCTGGTGGTCTACACCAACGCGGTTTGTGCGGGGGGCTTTTGCCTGTGCCTGCTGGGCGGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCg                                                                                                     <  1:3237519/150‑1 (MQ=21)
                                                 cTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCg                                                                                                     <  2:7952189/150‑1 (MQ=255)
                                                  tGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCGGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCgg                                                                                                    <  2:3101450/150‑1 (MQ=255)
                                                  tGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCgg                                                                                                    >  1:6965818/1‑150 (MQ=255)
                                                   gCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTCTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGt                                                                                                   <  2:1948942/150‑1 (MQ=255)
                                                       gcggcgAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCaaa                                                                                               >  2:3176067/1‑150 (MQ=255)
                                                       gcggcgAGGGGAGAGCATCTGCCATTCAGTCACCGGCTGGTGGTCTACACCAAACCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCaaa                                                                                               <  2:5190411/150‑1 (MQ=255)
                                                           cgAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATc                                                                                           >  1:11078957/1‑150 (MQ=255)
                                                            gAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATcc                                                                                          >  1:23416/1‑150 (MQ=255)
                                                            gAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATcc                                                                                          >  2:330731/1‑150 (MQ=255)
                                                                ggAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAg                                                                                      >  2:4445107/1‑150 (MQ=255)
                                                                 gagaGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCGGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTGCTCCACCCGGTCATAAATCCCCAGc                                                                                     >  2:9008470/1‑150 (MQ=255)
                                                                 gagaGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGc                                                                                     >  2:5433714/1‑150 (MQ=255)
                                                                     ggATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCACCt                                                                                 >  1:7460164/1‑150 (MQ=255)
                                                                      tttCTGCTATTCTGTCCCGGGGTGGTGGTCTTCACCAACGGGGTTTGTGTGGGGGGCTTTTGCCTGTGCCTGCTGTGCGGGGGGTGGCAAAAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCACCTc                                                                                <  1:4284948/148‑1 (MQ=21)
                                                                       aTCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCACCTCa                                                                               <  2:4291057/150‑1 (MQ=255)
                                                                       aTCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCACCTCa                                                                               >  2:1541943/1‑150 (MQ=255)
                                                                          tGCGATTCAGTCACCGGCTGGTGGTCTTCACCAACGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCACCTCATAc                                                                            <  2:10431724/150‑1 (MQ=255)
                                                                          tGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCACCTCATAc                                                                            >  2:7828980/1‑150 (MQ=255)
                                                                          tGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCACCTCATAc                                                                            >  2:2571294/1‑150 (MQ=255)
                                                                              aTTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTTTGTGTGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCACCTCATACgcgc                                                                        <  2:5046132/150‑1 (MQ=255)
                                                                              aTTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCACCTCATACgcgc                                                                        <  2:8393634/150‑1 (MQ=255)
                                                                                    tCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCACCTCATACGCGCCTTCaa                                                                  >  2:5319456/1‑150 (MQ=255)
                                                                                    tCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCACCTCATACGCGCCTTCaa                                                                  >  2:1618530/1‑150 (MQ=255)
                                                                                    tCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCACCTCATACGCGCCTTCaa                                                                  >  2:5740997/1‑150 (MQ=255)
                                                                                    tCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCACCTCATACGCGCCTTCaa                                                                  >  1:2999207/1‑150 (MQ=255)
                                                                                    tCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCACCTCATACGCGCCTTCaa                                                                  >  2:9711665/1‑150 (MQ=255)
                                                                                     cACCGGCTGGTGGTCTTCACCAAAGGGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCACCTCATACGCGCCTTCAAt                                                                 <  2:2246066/150‑1 (MQ=255)
                                                                                         ggCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCACCTCATACGCGCCTTCAATCAcc                                                             <  1:6657014/150‑1 (MQ=255)
                                                                                          gCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCACCTCATACGCGCCTTCAATCACCt                                                            >  1:3888426/1‑150 (MQ=255)
                                                                                          gCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCACCTCATACGCGCCTTCAATCACCt                                                            <  1:9901282/150‑1 (MQ=255)
                                                                                            tggtggTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCACCTCATACGCGCCTTCAATCACCTGa                                                          <  2:9443672/150‑1 (MQ=255)
                                                                                                 gtCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCACCTCATACGCGCCTTCAATCACCTGACTcac                                                     >  2:1228909/1‑150 (MQ=255)
                                                                                                  tCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCACCTCATACGCGCCTTCAATCACCTGACTcacc                                                    >  2:9898797/1‑150 (MQ=255)
                                                                                                   cTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCACCTCATACGCGCCTTCAATCACCTGACTcacca                                                   >  1:1417566/1‑150 (MQ=255)
                                                                                                    ttCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCACCTCATACGCACCTTCAATCACCTGACTcaccac                                                  >  1:5649253/1‑150 (MQ=255)
                                                                                                     tCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCACCTCATACGCGCCTTCAATCACCTGACTCACCACg                                                 >  1:6480344/1‑150 (MQ=255)
                                                                                                     tCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCACCTCATACGCGCCTTCAATCACCTGACTCACCACg                                                 >  2:8125848/1‑150 (MQ=255)
                                                                                                       aCCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCACCTCATACGCGCCTTCAATCACCTGACTCACCACGTc                                               >  1:1239004/1‑150 (MQ=255)
                                                                                                        ccAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCTCCTCATACGCGCCTTCAATCACCTGACTCACCACGTcc                                              >  1:10829733/1‑150 (MQ=255)
                                                                                                        ccAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCACCTCATACGCGCCTTCAATCACCTGACTCACCACGTcc                                              >  1:8573091/1‑150 (MQ=255)
                                                                                                        caaaaGGGGTTTGTGTGGGGGGCTTTTGCCTGTGCGTGGTGTGCGGGGGGTGGAAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCACCTCATACGCGCCTTCAATCACCTGACTCACCACGTcc                                              <  1:4827991/148‑1 (MQ=21)
                                                                                                         cAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCACCTCATACGCGCCTTCAATCACCTGACTCACCACGTccc                                             >  1:2506038/1‑150 (MQ=255)
                                                                                                          aacgCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCGGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCACCTCATACGCGCCTTCAATCACCTGACTCACCACGTcccc                                            <  1:9751774/147‑1 (MQ=255)
                                                                                                           acgCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCACCTCATACGCGCCTTCAATCACCTGACTCACCACGTccccc                                           <  2:2365814/148‑1 (MQ=255)
                                                                                                           aaGCGGTTTGTGAGGGCGGATTTTGCCAGTGCGTGCTGTGCGGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCACCTCATACGCGCCTTCAATCACCTGACTCACCACGTccccc                                           <  2:10315137/150‑1 (MQ=255)
                                                                                                           aaGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCACCTCATACGCGCCTTCAATCACCTGACTCACCACGTccccc                                           <  2:1794132/150‑1 (MQ=255)
                                                                                                                gTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCACCTCATACGCGCCTTCAATCACCTGACTCACCACGTCCCCCTagtg                                      >  1:9097500/1‑150 (MQ=255)
                                                                                                                 tATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCACCTCATACGCGCCTTCAATCACCTGACTCACCACGTCCCCCTtgtgt                                     <  2:295023/150‑1 (MQ=255)
                                                                                                                 tATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCACCTCATACGCGCCTTCAATCACCTGACTCACCACGTCCCCCTtgtgt                                     <  2:7142474/150‑1 (MQ=255)
                                                                                                                 tATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTACATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCACCTCATACGCGCCTTCAATCACCTGACTCACCACGTCCCCCTtgtgt                                     <  1:4405551/150‑1 (MQ=255)
                                                                                                                  aTGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCACCTCATACGCGCCTTCAATCACCTGACTCACCACGTCCCCCTtgtgtg                                    <  1:2772726/150‑1 (MQ=255)
                                                                                                                  aTGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCACCTCATACGCGCCTTCAATCACCTGACTCACCACGTCCCCCTtgtgtg                                    <  2:9275203/150‑1 (MQ=255)
                                                                                                                   tgtgAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTtctc                                                                                                              >  1:5304650/1‑75 (MQ=255)
                                                                                                                     tgAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCACCTCATACGCGCCTTCAATCACCTGACTCACCACGTCCCCCTTGTGTGGCa                                 <  2:6395892/150‑1 (MQ=255)
                                                                                                                      gAGGGCGGCTTTTGCCTGTGCCTGCTGTGCGGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCACCTCATACGCGCCTTCAATCACCTGACTCACCACGTCCCCCTTGTGTGGCAc                                <  2:4713150/150‑1 (MQ=255)
                                                                                                                         ggcggcTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAAAGACAGCATGGCATCCCGTTTCTCCACCACCCGGTCAAAAATCCCCAGCACCTCATCCGCGCCTTCAATAACCAGACACACCACGTCCCCCAAGTGTGGCAcccc                             >  1:10583137/1‑149 (MQ=21)
                                                                                                                          gcggcTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCACCTCATACGCGCCTTCAATCACCTGACTCACCACGTCCCCCTTGTGTGGCACCCGc                            >  1:9335705/1‑150 (MQ=255)
                                                                                                                           cggcTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCACCTCATACGCGCCTTCAATCACCTGACTCACCACGTCCCCCTTGTGTGGCACCCGCa                           <  1:6924076/150‑1 (MQ=255)
                                                                                                                               ttttGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCACCTCATACGCGCCTTCAATCACCTGACTCACCACGTCCCCCTTGTGTGGCACCCGCACCTc                       >  1:10246087/1‑150 (MQ=255)
                                                                                                                               ttttGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCACCTCATACGCGCCTTCAATCACCTGACTCACCACGTCCCCCTTGTGTGGCACCCGCACCTc                       >  1:1522445/1‑150 (MQ=255)
                                                                                                                                tttGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCACCTCATACGCGCCTTCAATCACCTGACTCACCACGTCCCCCTTGTGTGGCACCCGCACCTCg                      <  1:9595576/150‑1 (MQ=255)
                                                                                                                                    ccAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCACCTCATACGCGCCTTCAATCACCTGACTCACCACGTCCCCCTTGTGTGGCACCCGCACCTCGCCaa                  >  2:7217403/1‑150 (MQ=255)
                                                                                                                                    ccAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCACCTCATACGCGCCTTCAATCACCTGACTCACCACGTCCCCCTTGTGTGGCACCCGCACCTCGCCaa                  >  1:7785482/1‑150 (MQ=255)
                                                                                                                                    ccAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCACCTCATACGCGCCTTCAATCACCTGACTCACCACGTCCCCCTTGTGTGGCACCCGCACCTCGCCaa                  >  2:8711469/1‑150 (MQ=255)
                                                                                                                                     cAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCACCTCATACGCGCCTTCAATCACCTGACTCACCACGTCCCCCTTGTGTGGCACCCGCACCTCGCCaaa                 >  1:2854418/1‑150 (MQ=255)
                                                                                                                                      aGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCACCTCATACGCGCCTTCAATCACCTGACTCACCACGTCCCCCTTGTGTGGCACCCGCACCTCGCCAAAc                <  1:1818159/150‑1 (MQ=255)
                                                                                                                                       gTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCACCTCATACGCGCCTTCAATCACCTGACTCACCACGTCCCCCTTGTGTGGCACCCGCACCTCGCCAAACg               >  1:550165/1‑150 (MQ=255)
                                                                                                                                          cgtgctgTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCCCCCGGTCAAAAATCCCCAGCACCTCATACGCGCCTTCAATCACCTGACTCACCACGTCCCCCTTGTGTGGCACCCGCACCTCGCCAAACGACt            <  2:5593170/148‑1 (MQ=255)
                                                                                                                                          ccTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCACCTCATACGCGCCTTCAATCACCTGACTCACCACGTCCCCCTTGTGTGGCACCCGCACCTCGCCAAACGACt            >  2:1059947/1‑150 (MQ=255)
                                                                                                                                             gctgTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCACCTCATACGCGCCTTCAATCACCTGACTCACCACGTCCCCCTTGTGTGGCACCCGCACCTCGCCAAACGACTCAc         >  1:6330749/1‑150 (MQ=255)
                                                                                                                                                 tgCGGGGGGGGGAACCCGCAACGACTGCAGGGCATCCGGTTTCTCCTCCACCCGGAAAAAAAGCCCCAGCACCTCATCCGCGCCTTCAATCCCCTGCCTCACCACGGCCCCCTTGTGTGGCACCCGCACCTCGCCACACGACTCACCGCa     <  1:3409622/150‑1 (MQ=18)
                                                                                                                                                    |                                                                                                                                                     
GGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCACCTCATACGCGCCTTCAATCACCTGACTCACCACGTCCCCCTTGTGTGGCACCCGCACCTCGCCAAACGACTCACCGCAGAC  >  NC_000913/267324‑267621

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: