Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 289,257 T→A 5.1% intergenic (‑95/+44) yagI ← / ← argF CP4‑6 prophage; putative DNA‑binding transcriptional regulator/ornithine carbamoyltransferase 2, chain F; CP4‑6 prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913289,2570TA 5.1% 481.3 / 11.2 137intergenic (‑95/+44)yagI/argFCP4‑6 prophage; putative DNA‑binding transcriptional regulator/ornithine carbamoyltransferase 2, chain F; CP4‑6 prophage
Reads supporting (aligned to +/- strand):  ref base T (105/24);  new base A (5/2);  total (110/27)
Fisher's exact test for biased strand distribution p-value = 6.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

Only 100 of 144 total aligned reads displayed.

CGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAACACCTCGTCCGTCACC  >  NC_000913/289108‑289406
                                                                                                                                                     |                                                                                                                                                     
 gTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCtgtg                                                                                                                                                      <  1:1709375/150‑1 (MQ=255)
 gTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCtgtg                                                                                                                                                      <  2:140339/150‑1 (MQ=255)
  ttGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTgg                                                                                                                                                     <  2:3200249/150‑1 (MQ=255)
    gAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGga                                                                                                                                                   <  1:2233110/150‑1 (MQ=255)
      aCAGGTCGAGGATCTGCAACTCACGTTAAACAGACTGAATAATCGTCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGgaga                                                                                                                                                 <  1:95373/150‑1 (MQ=255)
      aCAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGgaga                                                                                                                                                 <  2:2811260/150‑1 (MQ=255)
      aCAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGgaga                                                                                                                                                 <  1:95360/150‑1 (MQ=255)
       cAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGagag                                                                                                                                                <  1:1351131/150‑1 (MQ=255)
        aGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAgg                                                                                                                                               <  2:3976302/150‑1 (MQ=255)
        aGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAgg                                                                                                                                               <  1:2031971/150‑1 (MQ=255)
        aGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCACAGGGAGAgg                                                                                                                                               <  2:1511251/150‑1 (MQ=255)
        aGGTCGAGGATCTGCAACGCACGTTCAACAGACCGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAATACGTACCTTTTTCGGGATGAAAAATCAATGAACTGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAgg                                                                                                                                               <  1:3852886/150‑1 (MQ=255)
        aGGTCGAGGATCTGCAACGCACGTACAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAgg                                                                                                                                               <  1:3352791/150‑1 (MQ=255)
           tCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGcc                                                                                                                                            >  1:1201670/1‑150 (MQ=255)
           tCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGcc                                                                                                                                            >  1:4121585/1‑150 (MQ=255)
           tCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCATTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGcg                                                                                                                                            <  2:572888/150‑2 (MQ=255)
            cGAGGCTCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAACCGGTCCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCg                                                                                                                                           <  2:105793/150‑1 (MQ=255)
             gAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCgg                                                                                                                                          >  1:3997503/1‑150 (MQ=255)
             gAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCgg                                                                                                                                          >  1:1478993/1‑150 (MQ=255)
             gAGGATCTGCAACGCACGTTCAACAGACGGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCgg                                                                                                                                          <  2:3886134/150‑1 (MQ=255)
                gATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGt                                                                                                                                       >  1:2051268/1‑150 (MQ=255)
                   cTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAg                                                                                                                                    <  2:4075805/150‑1 (MQ=255)
                   cTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAg                                                                                                                                    <  2:4075774/150‑1 (MQ=255)
                   cTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAg                                                                                                                                    <  2:2286212/150‑1 (MQ=255)
                    tGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAgg                                                                                                                                   >  2:2958836/1‑150 (MQ=255)
                     gCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAggg                                                                                                                                  >  2:2462777/1‑150 (MQ=255)
                     gCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAggg                                                                                                                                  >  1:243726/1‑150 (MQ=255)
                     gCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCAGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAggg                                                                                                                                  >  1:318883/1‑150 (MQ=255)
                          gCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCa                                                                                                                             >  1:452594/1‑150 (MQ=255)
                          gCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCa                                                                                                                             >  2:1469196/1‑150 (MQ=255)
                          gCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCa                                                                                                                             >  1:3405695/1‑150 (MQ=255)
                               ttCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAAcgcg                                                                                                                        >  1:3365129/1‑150 (MQ=255)
                                  aaCAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAg                                                                                                                     >  1:2925626/1‑150 (MQ=255)
                                  aaCAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAg                                                                                                                     >  1:2623351/1‑150 (MQ=255)
                                   actgACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAACACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCGGTGGGAGAGGGCCGGGGTGAGGGCCCCAACGCGCAGc                                                                                                                    >  2:1196704/4‑150 (MQ=255)
                                   aCAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGc                                                                                                                    >  1:2009359/1‑150 (MQ=255)
                                    cAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCg                                                                                                                   >  1:3820657/1‑150 (MQ=255)
                                    cAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCg                                                                                                                   >  1:796805/1‑150 (MQ=255)
                                     aGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCgg                                                                                                                  >  2:3269678/1‑150 (MQ=255)
                                     aGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCgg                                                                                                                  >  2:3269645/1‑150 (MQ=255)
                                     aGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCgg                                                                                                                  >  2:1441741/1‑150 (MQ=255)
                                       aCTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGAGGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCCGCGGAc                                                                                                                >  2:4086335/1‑150 (MQ=255)
                                         tGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCCCAACGCGCAGCGGAccc                                                                                                              >  2:2484630/1‑150 (MQ=255)
                                          gAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCa                                                                                                             >  2:3608736/1‑150 (MQ=255)
                                           aataatCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCaa                                                                                                            >  2:4102572/1‑150 (MQ=255)
                                                tCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACt                                                                                                       >  2:1309759/1‑150 (MQ=255)
                                                tCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACt                                                                                                       >  1:549350/1‑150 (MQ=255)
                                                 cGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTc                                                                                                      >  1:2865653/1‑150 (MQ=255)
                                                 cGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTc                                                                                                      >  1:2865472/1‑150 (MQ=255)
                                                     aTAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCa                                                                                                  >  1:662183/1‑150 (MQ=255)
                                                     aTAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCa                                                                                                  >  1:662251/1‑150 (MQ=255)
                                                              gCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGcc                                                                                         >  1:220674/1‑150 (MQ=255)
                                                                aTGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCcat                                                                                       >  2:218429/1‑150 (MQ=255)
                                                                aTGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCcat                                                                                       >  2:218196/1‑150 (MQ=255)
                                                                aTGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCcat                                                                                       >  1:3865639/1‑150 (MQ=255)
                                                                   tCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCcatcat                                                                                    >  1:4142800/1‑150 (MQ=255)
                                                                     cACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCAtcatca                                                                                  >  1:3632943/1‑150 (MQ=255)
                                                                     cACGCTGGAATTAAGGCGAAAACGTACCTTTTTAGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCatcatcc                                                                                  >  1:378106/1‑149 (MQ=255)
                                                                          tGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGACATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGCCCCAATCACTCCCCAAGCGGTGCCAGCCTCCCCGcc                                                                             >  1:2299960/1‑150 (MQ=38)
                                                                          tGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGcc                                                                             >  1:898022/1‑150 (MQ=255)
                                                                          tGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGcc                                                                             >  2:1384284/1‑150 (MQ=255)
                                                                                    gCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTAt                                                                   >  2:320346/1‑150 (MQ=255)
                                                                                      gAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGc                                                                 >  1:745937/1‑150 (MQ=255)
                                                                                      gAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGc                                                                 >  2:3070013/1‑150 (MQ=255)
                                                                                            gTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCgg                                                           >  1:2518341/1‑150 (MQ=255)
                                                                                              aCCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGtt                                                         >  2:246268/1‑150 (MQ=255)
                                                                                                 tttttCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTc                                                      >  1:3923154/1‑150 (MQ=255)
                                                                                                 tttttCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGGTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTc                                                      >  2:2815255/1‑150 (MQ=255)
                                                                                                     tCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGcc                                                  >  1:213062/1‑150 (MQ=255)
                                                                                                     tCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGcc                                                  >  1:1876235/1‑150 (MQ=255)
                                                                                                        ggATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCCTCCGGTTTTCCGCCTgg                                               >  1:3499510/1‑150 (MQ=255)
                                                                                                        ggATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTgg                                               >  1:3901540/1‑150 (MQ=255)
                                                                                                          aTGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTc                                             >  2:4040722/1‑150 (MQ=255)
                                                                                                           tGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCg                                            >  2:371424/1‑150 (MQ=255)
                                                                                                            gAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGa                                           >  1:25268/1‑150 (MQ=255)
                                                                                                            gAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGa                                           >  2:3965246/1‑150 (MQ=255)
                                                                                                             aaaaaTCAATTAATTGGCGCCGGTGTTCGCCCTCTCCCTGTGGGAGCGGGCCGGGGTGCGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGaa                                          <  2:1384365/150‑1 (MQ=255)
                                                                                                             aaaaaTCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGCGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGaa                                          >  2:457909/1‑150 (MQ=255)
                                                                                                             aaaaaTCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGaa                                          >  1:3089969/1‑150 (MQ=255)
                                                                                                              aaaaTCAATGACCTGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGAGAGGGCACCAACGCGCAGCGGACCCCCTCACTCCCCAAGCGTTGCCATCATCACCGCCGTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAac                                         >  2:2073388/1‑150 (MQ=38)
                                                                                                                    aaTGAAATGGAGCCGGTGTTCTCCCTCTCCCTTTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCCTCCTCACCGCCTTAATCGTATGCATCCGGTTGTCCGCCTGGTCGAACACGATg                                   >  2:2731391/1‑150 (MQ=255)
                                                                                                                     aTGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGc                                  >  1:2677235/1‑150 (MQ=255)
                                                                                                                      tGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCt                                 >  1:2498963/1‑150 (MQ=255)
                                                                                                                               gCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGa                        >  1:998972/1‑150 (MQ=255)
                                                                                                                                  ggTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGGTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTc                     >  1:3649342/1‑150 (MQ=255)
                                                                                                                                  ggTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACACCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTc                     >  1:3649405/1‑150 (MQ=255)
                                                                                                                                        ggccctctcccaCAGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAacac               >  2:1749680/3‑150 (MQ=255)
                                                                                                                                         tccctctccctGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAACAcc              >  2:210261/1‑150 (MQ=255)
                                                                                                                                          ccctctccctGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAACACCt             >  2:3845552/1‑150 (MQ=255)
                                                                                                                                           cctctccctGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAACACCTc            >  2:1367150/1‑150 (MQ=255)
                                                                                                                                           cctctccctGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAACACCTc            >  2:22819/1‑150 (MQ=255)
                                                                                                                                           cctctcccaCAGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAACACCTc            >  2:1695950/1‑150 (MQ=255)
                                                                                                                                           cctctcccaCAGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAACACCTc            >  1:25641/1‑150 (MQ=255)
                                                                                                                                           cctctcccaCAGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAACACCTc            >  1:25626/1‑150 (MQ=255)
                                                                                                                                           cctccccgtGTGGGTGTGGGCGGGGGTGTGGGCCCCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAACACCTc            <  1:1683066/150‑1 (MQ=255)
                                                                                                                                               tccctGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAACACCTcgtcc        >  1:1264298/1‑150 (MQ=255)
                                                                                                                                                     tgGGGGCGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAACACCTCGTCCGTCAcc  <  1:1024695/150‑1 (MQ=255)
                                                                                                                                                     tgGGAGCGGCCCGGGGTGAGCCCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAAACGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAACACCTCGTCCGTCAcc  <  2:2065488/150‑1 (MQ=255)
                                                                                                                                                     tgGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAAACGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAACACCTCGACCGTCAcc  >  2:3145910/1‑150 (MQ=255)
                                                                                                                                                     tgGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAAACACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCTCCGCCGACTCAAACACCTCGTCCGTCAcc  >  2:262534/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAACACCTCGTCCGTCACC  >  NC_000913/289108‑289406

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: