Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,629,701 T→C E274E (GAA→GAG ettA ← energy‑dependent translational throttle protein EttA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,629,7010TC97.1% 219.1 / ‑4.3 70E274E (GAA→GAGettAenergy‑dependent translational throttle protein EttA
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (36/32);  total (37/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGATGATCCCGACGATCGCTCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGACGATCGCCATAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGTTCTTCAAAGCGCGCCAGACGTGCTTTACCTTTCGACTGACGGCCTTTAGTACCTTGACGTACCCATTCCAGCTCTTTCTCAATCGACTTACGACGCGCCGCTTCTTGTGAAGCTTCCTGCGCCAGGCGCTGATCTTTCTGCTCCAGCCAGGAGGAGTAGTTACCTTCCCACGGAATACCTTCACCGCGGTCAAGTTCGAGGATCCAGCCTGCAACGTTATCGAGGAAGTAACGGTCGTGGGTAATCGCCACAACGGTGCCTTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCG  >  NC_000913/4629456‑4629938
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                             
ccGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGCTCTTCAAAGCGTGCCAGACGTGCTTTACCTTTCGACTGGCGGCCTTTAGTTCCCTGACGTACCCACTCCAg                                                                                                                                                                                                                                          <  1:477414/251‑1 (MQ=255)
        cccGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGCTCTTCAAAGCGTGCCAGACGTGCTTTACCTTTCGACTGGCGGCCTTTAGTTCCCTGACGTACCCACTCCAGCTCTTtc                                                                                                                                                                                                                                   <  1:195578/250‑1 (MQ=255)
         ccGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGCTCTTCAAAGCGTGCCAGACGTGCTTTACCTTTCGACTGGCGGCCTTTAGTTCCCTGACGTACCCACTCCAGCTCTTtctc                                                                                                                                                                                                                                 >  2:22147/1‑251 (MQ=255)
            acgaTTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGCTCTTCAAAGCGTGCCAGACGTGCTTTACCTTTCGACTGGCGGCCTTTAGTTCCCTGACGTACCCACTCCAGCTCTTtctcgat                                                                                                                                                                                                                              >  2:276628/1‑248 (MQ=255)
                   ccccTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGCTCTTCAAAGCGTGCCAGACGTGCTTTACCTTTCGACTGGCGTCCTTTAGTTCCCTGACGTACCCACTCCAGCTCTTTCTCGATCGACTTa                                                                                                                                                                                                                       <  1:430716/249‑1 (MQ=255)
                            ggATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGCTCTTCAAAGCGTGCCAGACGTGCTTTACCTTTCGACTGGCGGCCTTTAGTTCCCTGACGTACCCACTCCAGCTCTTTCTCGATCGACTTACGACGcgcc                                                                                                                                                                                                              <  1:325891/251‑1 (MQ=255)
                             gATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGCTCTTCAAAGCGTGCCAGACGTGCTTTACCTTTCGACTGGCGGCCTTTAGTTCCCTGACGTACCCACTCCAGCTCTTTCTCGATCGACTTACGACGcgccg                                                                                                                                                                                                             >  1:441517/1‑251 (MQ=255)
                                  agaAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAACCAGTTCGTTGGTTTCGTTACGTTTCTGCTATTCGGTGCTGTTCAGCTCTTCAAAGCGTGCCAGACGTGCTTTACCTTTCGACTGGCGGCCTTTAGTTCCCTGTCGTCCCCACTCCAGCTCTTTCTCGATCGACTTACGCCGCGCCGcttct                                                                                                                                                                                                        >  2:344192/1‑251 (MQ=255)
                                    aaGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGCTCTTCAAAGCGTGCCAGACGTGCTTTACCTTTCGACTGGCGGCTTTTAGTTCCCTGACGTACCCACTCCAGCTCTTTCTCGATCGACTTACGACGCGCCGCTTCTtg                                                                                                                                                                                                      <  1:85895/251‑1 (MQ=255)
                                         cAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCTCTTTATCGCCCAGACGCGGCCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGTTGTTGTTCAGCTCTTCAAAGCTTGCCAGACGTGCTTTACCTTTCGACTGGCGGCCTTTAGTTCCCTGACGTACCCACTCCAGCTCTTTCTCGCTCGACCTACGACGCGCCGCTTCTTGTGaa                                                                                                                                                                                                  <  1:71717/250‑1 (MQ=255)
                                         cAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGCTCTTCAAAGCGTGCCAGACGTGCTTTACCTTTCGACTGGCGGCCTTTAGTTCCCTGACGTACCCACTCCAGCTCTTTCTCGATCGACTTACGACGCGCCGCTTCTTGTGaa                                                                                                                                                                                                  <  1:337463/250‑1 (MQ=255)
                                          aGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGCTCTTCAAAGCGTGCCAGACGTGCTTTACCTTTCGACTGGCGGCCTTTAGTTCCCTGACGTACCCACTCCAGCTCTTTCTCGATCGACTTACGACGCGCCGCTTCTTGTGAAGc                                                                                                                                                                                                <  1:172466/251‑1 (MQ=255)
                                              catcaATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGCTCTTCAAAGCGTGCCAGACGTGCTTTACCTTTCGACTGGCGGCCTTTAGTTCCCTGACGTACCCACTCCAGCTCTTTCTCGATCGACTTACGACGCGCCGCTTCTTGTGAAGCTTcc                                                                                                                                                                                            <  1:476232/251‑1 (MQ=255)
                                               atcaatcaGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGCTCTTCAAAGCGTGCCAGACGTGCTTTACCTTTCGACTGGCGGCCTTTAGTTCCCTGACGTACCCACTCCAGCTCTTTCTCGATCGACTTACGACGCGCCGCTTCTTGTGAAGCTTCCt                                                                                                                                                                                           >  1:432635/1‑251 (MQ=255)
                                                   atcaGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGATCTTCAAAGCGTGCCAGACGTGCTTTACCTTTCGACTGGCGGCCTTTAGTACCCTGACGTACCCACTCCAGCTCTTTCTCGATCGACTTACGCCGCGCCGCTTCTTGTGAAGCTTCCTgcgc                                                                                                                                                                                       >  1:172895/1‑251 (MQ=255)
                                                         aggcGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGCTCTTCAAAGCGTGCCAGACGTGCTTTACCTTTCGACTGGCGGCCTTTAGTTCCCTGACGTACCCACTCCAGCTCTTTCTCGATCGACTTACGACGCGCCGCTTCTTGTGAAGCTTCCTGCGCCAGgcg                                                                                                                                                                                 <  1:24907/248‑1 (MQ=255)
                                                            cGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGCTCTTCAAAGCGTGCCAGACGTGCTTTACCTTTCGACTGGCGGCCTTTAGTTCCCTGACGTACCCACTCCAGCTCTTTCTCGATCGACTTACGACGCGCCGCTTCTTGTGAAGCTTCCTGCGCCAGGCGCTg                                                                                                                                                                              >  1:128603/1‑251 (MQ=255)
                                                             gatcaccgtAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGCTCTTCAAAGCGTGCCAGACGTGCTTTACCTTTCGACTGGCGGCCTTTAGTTCCCTGACGTACCCACTCCAGCTCTTTCTCGATCGACTTACGACGCGCCGCTTCTTGTGAAGCTTCCTgcg                                                                                                                                                                                        >  1:187751/9‑240 (MQ=255)
                                                             gatcaccgtAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGCTCTGCAAAGCGTGCCAGACGTGCTTTACCTTTCGACTGGCGGCCTTTAGTTCCCTGACGTACCCACTCCAGCTCTTTCTCGATCGACTTACGACGCGCCGCTTCTTGTGAAGCTTCCTgcg                                                                                                                                                                                        <  2:187751/232‑1 (MQ=255)
                                                             gATCGCCATAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGATCTTCAAAGCGCGCCAGACGTGCTTTACCTTTCGACTGACGGCCTTTAGTACATTGACGTACCCATTCCAGCTCTTTCTCAATCGACTTACGACGCGACGCTTCTTGTGAAGCTTCCTGCGCCAGTCGCTg                                                                                                                                                                              >  1:105327/1‑250 (MQ=255)
                                                                caccgtAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGCTCTTCAAAGCGTGCCAGACGTGCTTTACCTTTCGACTGGCGGCCTTTAGTTCCCTGACGTACCCACTCCAGCTCTTTCTCGATCGACTTACGACGCGCCGCTTCTTGTGAAGCTTCCTGCGCCAGGCGCTGATCt                                                                                                                                                                          <  1:343404/246‑1 (MQ=255)
                                                                    ttcttctttcctcAGGTTGCTTATTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGCTCTTCAAAGCGTGCCAGACGTGCTTTACCTTTCGACTGGCGGCCTTTAGTTCCCTGACGTACCCACTCCAGCTCTTTCTCGATCGACTTACGACGCGCCGCTTCTTGTGAAGCTTCCTGCGCCAGGCGCTGATCTTTCt                                                                                                                                                                      >  1:203944/13‑251 (MQ=255)
                                                                     tAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGCTCTTCAAAGCGTGCCAGACGTGCTTTACCTTTCGACTGGCGGCCTTTAGTTCCCTGACGTACCCACTCCAGCTCTTTCTCGATCGACTTACGACGCGCCGCTTCTTGTGAAGCTTCCTGCGCCAgg                                                                                                                                                                                   >  2:432409/1‑237 (MQ=255)
                                                                     tAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGCTCTTCAAAGCGTGCCAGACGTGCTTTACCTTTCGACTGGCGGCCTTTAGTTCCCTGACGTACCCACTCCAGCTCTTTCTCGATCGACTTACGACGCGCCGCTTCTTGTGAAGCTTCCTGCGCCAgg                                                                                                                                                                                   <  1:432409/237‑1 (MQ=255)
                                                                         aTTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGCTCTTCAAAGCGTGCCAGACGTGCTTTACCTTTCGACTGGCGGCCTTTAGTTCCCTGACGTACCCACTCCAGCTCTTTCTCGATCGACTTACGACGCGCCGCTTCTTGTGAAGCTTCCTGCGCCAGGCGCTGATCTTTCTGCTcc                                                                                                                                                                 >  1:448758/1‑251 (MQ=255)
                                                                         aTTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGCTCTTCAAAGCGTGCCAGACGTGCTTTACCTTTCGACTGGCGGCCTTTAGTTCCCTGACGTACCCACTCCAGCTCTTTCTCGATCGACTTACGACGCGCCGCTTCTTGTGAAGCTTCCTGCGCCAGGCGCTGATCTTTCTGCTcc                                                                                                                                                                 <  1:22147/251‑1 (MQ=255)
                                                                            tACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGCTCTTCAAAGCGTGCCAGACGTGCTTTACCTTTCGACTGGCGGCCTTTAGTTCCCTGACGTACCCACTCCAGCTCTTTCTCGATCGACTTACGACGCGCCGCTTCTTGTGAAGCTTCCTGCGCCAGGCGCTGATCTTTCTGCTccagc                                                                                                                                                              >  1:18557/1‑251 (MQ=255)
                                                                                    ttGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGCTCTTCAAAGCGTGCCAGACGTGCTTTACCTTTCGACTGGCGGCCTTTAGTTCCCTGACGTACCCACTCCAGCTCTTTCTCGATCGACTTAAGACGCGCCGCTTCTTGTGAAGCTTCCTGCGCCAGGCGCTGATCTTTCTGCTCCAGCCcggaggg                                                                                                                                                      >  2:335304/1‑250 (MQ=255)
                                                                                        tGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGCTCTTCAAAGCGTGCCAGACGTGCTTTACCTTTCGACTGGCGGCCTTTAGTTCCCTGACGTACCCACTCCAGCTCTTTCTCGATCGACTTACGACGAGCCGCTTCTTGTGAAGCTTCCTGCGCCAGGCGCTGATCTTTCTGCTCCAGCCAGGAGGagt                                                                                                                                                    >  1:254897/1‑249 (MQ=255)
                                                                                         gtcTTCCAGCATTTTAACGCCTAGACGCGGTCCAGGTGAAATAACCAGTCCTTTGGGTTGGTTGCTTTTATGATAGTCGGTGCTGTTCAGCTCTTCAAAGCGTGCCAGACGTTCGTTACCTTTCGACTGGCGGTCTTTAGTTCCCTGACGTACCCACTCCAGCTCTTTCTCGAACGACTTACGACTCGCCGCGTCGTGTGAAGCTTCCTGCGCCAGGCGCTGATCTTTCTGCTCCAGCCAGGAGGagta                                                                                                                                                   <  2:416924/247‑1 (MQ=255)
                                                                                           cTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGCTCTTCAAAGCGTGCCAGACGTGCTTTACCTTTCGACTGGCGGCCTTTAGTTCCCTGACGTACCCACTCCAGCTCTTTCTCGATCGACTTACGACGCGCCGCTTCTTGTGAAGCTTCCTGCGCCAGGCGCTGATCTTTCTGCTCCAGCCAGGAGGAGTAGttg                                                                                                                                               >  1:71536/1‑250 (MQ=255)
                                                                                              ccAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGCTCTTCAAAGCGTGCCAGACGTGCTTTACCTTTCGACTGGCGGCCTTTAGTTCCCTGACGTACCCACTCCAGCTCTTTCTCGATCGACTTACGACGCGCCGCTTCTTGTGAAGCTTCCTGCGCCAGGCGCTGATCTTTCTGCTCCAGCCAGGAGGAGTAGtt                                                                                                                                                >  2:158646/1‑247 (MQ=255)
                                                                                              ccAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGCTCTTCAAAGCGTGCCAGACGTGCTTTACCTTTCGACTGGCGGCCTTTAGTTCCCTGACGTACCCACTCCAGCTCTTTCTCGATCGACTTACGACGCGCCGCTTCTTGTGAAGCTTCCTGCGCCAGGCGCTGATCTTTCTGCTCCAGCCAGGAGGAGTAGtt                                                                                                                                                <  1:158646/247‑1 (MQ=255)
                                                                                                 gCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGCTCTTCAAAGCGTGCCAGACGTGCTTTACCTTTCGACTGGCGGCCTTTAGTTCCCTGACGTACCCACTCCAGCTCTTTCTCGATCGACTTACGACGCGCCGCTTCTTGTGAAGCTTCCTGCGCCAGGCGCTGATCTTTCTGCTCCAGCCAGGAGGAGTAGTTGCCTTcc                                                                                                                                         <  1:481546/251‑1 (MQ=255)
                                                                                                         ccGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGCTCTTCAAAGCGTGCCAGACGTGCTTTACCTTTCGACTGGCGGCCTTTAGTTCCCTGACGTACCCACTCCAGCTCTTTCTCGATCGACTTACGACGCGCCGCTTCTTGTGAAGCTTCCTGCGCCAGGCGCTGATCTTTCTGCTCCAGCCAGGAGGAGTAGTTGCCTTCCCACGGAAt                                                                                                                                 <  2:188214/250‑1 (MQ=255)
                                                                                                          cGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGCTCTTCAAAGCGTGCCAGACGTGCTTTACCTTTCGACTGGCGGCCTTTAGTTCCCTGACGTACCCACTCCCGCTCTTTCTCGCTCGACTTACGACGCGCCGCTTCTTGTGAAGCTTCCTGCGCCAGGCGCTGATCTTTCTGCTCCAGCCACGCGGAGTAGTTGCCTTCCCACGGacac                                                                                                                                >  2:326584/1‑248 (MQ=255)
                                                                                                            cccAGACGCGGACCAGGCGCAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGCTCTTCAAAGCGTGCCAGACGTGCTTTACCTTTCGAATGGCGGCCTTTAGTTCCCTGACGTACCCACTCCAGATCTTTATCGATCGACTTCCGACGCGCCGCTTCTTGTGCAGCTTCCTGCGCCCGGCGCTGATCTTTCTGCTCCAGCCAGGAGGGGTAGTTGCCTTCCCACGGCATAcc                                                                                                                              >  2:439484/1‑251 (MQ=255)
                                                                                                                   gcgGGCCAGGTTGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGCTCTTCAAAGCGTGCCAGACGTGCTTTACCTTTCGACTGGCGGCCTTTAGTTCCCTGACGTACCCACTCCAGCTCTTTCTCGATCGACTTACGACGCGCCGCTTCTTGTGAAGCTTCCTGCGCCAGGCGCTGATCTTTCTGCTCCAGCCAGGAGGAGTAGTTGCCTTCCCACGGAATACCTTCACcg                                                                                                                       <  2:58464/251‑1 (MQ=255)
                                                                                                                            gtgtaaTAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGTTCTTCAAAGCGCGCCAGACGTGCTTTACCTTTCGACTGACGGCCTTTAGTACCTTGACGTACCCATCCCAGCTCTTTCTCAATCGACATACGACGCGCCGCTCCTTGTGAAGCTTCCGGCGCCAGGCGCTGATCTTTCTGCTCCAGCCAGGAGGAGTAGTTACCTTCCCACGGAATACCTTCACCGCGGTCAAGt                                                                                                              <  1:34904/247‑1 (MQ=255)
                                                                                                                                aaTAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGCTCTTCAAAGCGTGCCAGACGTGCTTTACCTTTCGACTGGCGGCCTTTAGTTCCCTGACGTACCCACTCCAGCTCTTTCTCGATCGACTTACGACGCGCCGCTTCTTGTGAAGCTTCCTGCGCCAGGCGCTGATCTTTCTGCTCCAGCCAGGAGGAGTAGTTGCCTTCCCACGGAATACCTTCACCGCGGTCAAGTTCg                                                                                                           >  1:240962/1‑250 (MQ=255)
                                                                                                                                 aTAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGCTCTTCAAAGCGTGCCAGACGTGCTTTACCTTTCGACTGGCGGCCTTTAGTTCCCTGACGTACCCACTCCAGCTCTTTCTCGATCGACTTACGACGCGCCGCTTCTTGTGAAGCTTCCTGCGCCAGGCGCTGATCTTTCTGCTCCAGCCAGGAGGTGTAGTTGCCTTCCCACGGAATACCTTCACCGCGGTCAAGTTCg                                                                                                           >  1:195277/1‑249 (MQ=255)
                                                                                                                                      cAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGCTCTTCAAAGCGTGCCAGACGTGCTTTACCTTTCGACTGGCGGCCTTTAGTTCCCTGACGTACCCACTCCAGCTCTTTCTCGATCGACTTACGACGCGCCGCTTCTTGTGAAGCTTCCTGCGCCAGGCGCTGATCTTTCTGCTCCAGCCAGGAGGAGTAGTTGCCTTCCCACGGAATACCTTCACCGCGGTCAAGTTCGAGGATcc                                                                                                    >  2:421311/1‑251 (MQ=255)
                                                                                                                                      cAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGCTCTTCAAAGCGTGCCAGACGTGCTTTACCTTTCGACTGGCGGCCTTTAGTTCCCTGACGTACCCACTCCAGCTCTTTCTCGATCGACTTACGACGCGCCGCTTCTTGTGAAGCTTCCTGCGCCAGGCGCTGATCTTTCTGCTCCAGCCAGGAGGAGTAGTTGCCTTCCCACGGAATACCTTCACCGCGGTCAAGTTCGAGGATcc                                                                                                    >  1:47545/1‑251 (MQ=255)
                                                                                                                                       aGTTCTTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGCTCTTCAAAGCGTGCCAGACGTGCTTTACCTTTCGACTGGCGGCCTTTAGTTCCCTGACGTACCCACTCCAGCTCTTTCTCGATCGACTTACGACGCGCCGCTTCTTGTGAAGCTTCCTGCGCCAGGCGCTGATCTTTCTGCTCCAGCCAGGAGGAGTAGTTGCCTTCCCACGGAATACCTTCACCGCGGTCAAGTTCGAGGATCCa                                                                                                   <  2:448758/251‑1 (MQ=255)
                                                                                                                                            gTTGGTTTCGTTAAGTTTCTGATATTCGGGGCTGGTCAGCCCTTCAAAGCGTGCCAGACGTGCTTTACCTTGCGACTGGCGGCCTTTAGTCCCCTGACGTACCCACACCAGCTCTTTCGCGATCGACTTACGACGCGCCGCTTCTTGTGAAGCTTCCTGCGCCAGGCGCTGATCTTTCTGCTCCAGCCAGGAGGAGTAGTTGCCTTCCCACGGAATACCTTCACCGCGGTCAAGTTCGAGGATCCAGCCTg                                                                                              <  2:203944/251‑1 (MQ=255)
                                                                                                                                                      ttACGTTTCTGATATTCGGTGCTGTTCAGCTCTTCAAAGCGTGCCAGACGTGCTTTACCTTTCGACTGGCGGCCTTTAGTTCCCTGACGTACCCACTCCAGCTCTTTCTCGATCGACTTACGACGCGCCGCTTCTTGTGAAGCTTCCTGCGCCAGGCGCTGATCTTTCTGCTCCAGCCAGGAGGAGTAGTTGCCTTCCCACGGAATACCTTCACCGCGGTCAAGTTCGAGGATCCAGCCTGCAACGTTATc                                                                                    >  1:102935/1‑251 (MQ=255)
                                                                                                                                                        aCGTTTCTGATATTCGGTGCTGTTCAGCTCTTCAAAGCGTGCCAGACGTGCTTTACCTTTCGACTGGCGGCCTTTAGTTCCCTGACGTACCCACTCCAGCTCTTTCTCGATCGACTTACGACGCGCCGCTTCTTGTGAAGCTTCCTGCGCCAGGCGCTGATCTTTCTGCTCCAGCCAGGAGGAGTAGTTGCCTTCCCACGGAATACCTTCACCGCGGTCAAGTTCGAGGATCCAGCCTGCAACGTTATCGa                                                                                  <  1:176877/251‑1 (MQ=255)
                                                                                                                                                          gTTTCTGATATTCGGTGCTGTTCAGCTCTTCAAAGCGTGCCAGACGTGCTTTACCTTTCGACTGGCGGCCTTTAGTTCCCTGACGTACCCACTCCAGCTCTTTCTCGATCGACTTACGACGCGCCGCTTCTTGTGAAGCTTCCTGCGCCAGGCGCTGATCTTTCTGCTCCAGCCAGGAGGAGTAGTTGCCTTCCCACGGAATACCTTCACCGCGGTCAAGTTCGAGGATCCAGCCTGCAACGTTATCGAgg                                                                                <  1:24316/251‑1 (MQ=255)
                                                                                                                                                                   atTCGGTGCTGTTCAGCTCTTCAAAGCGTGCCAGACGTGCTTTACCTTTCGACTGGCGGCCTTTAGTTCCCTGACGTACCCACTCCAGCTCTTTCTCGATCGACTTACGACGCGCCGCTTCTTGTGAAGCTTCCTGCGCCAGGCGCTGATCTTTCTGCTCCAGCCAGGAGGAGTAGTTGCCTTCCCACGGAATACCTTCACCGCGGTCAAGTTCGAGGATCCAGCCTGCAACGTTATCGAGGAAGTAACg                                                                        <  2:128603/250‑1 (MQ=255)
                                                                                                                                                                      cGGTGCTGTTCCGCTCTTCAACGCGTGCCAGACGTGCTTTCCCTTTCGACTGGCGGCCTTTAGTTCCCTGACGTACCCACTCCAGCTCTTTCTCGATCGACTTACGACGCGCCGCTTCTTGTGAAGCTTCCTGCGCCAGGCGCTGATCTTTCTGCTCCAGCCAGGAGGAGTAGTTGCCTTCCCACGGAATACCTTCACCGCGGTCAAGTTCGAGGATCCAGCCTGCAACGTTATCGAGGAAGTAACGGTCg                                                                    <  2:282156/251‑1 (MQ=255)
                                                                                                                                                                            tGTTCAGCTCTTCAAAGCGTGCCAGACGTGCTTTACCTTTCGACTGGCGGCCTTTAGTTCCCTGACGTACCCACTCCAGCTCTTTCTCGATCGAATTACGACGCGCCGCTTCTTGTGAAGCTTCCTTCGCCAGGCGCTGATCTTTCTGCTCCAGCCAGGAGGAGTAGTTGCCTTCCCACGGAATACCTTCACCGAGGTCAAGTTCGAGGCTCCAGCCTGCAACGTTATCGAGGTAGTAACGGTCGgggga                                                               >  2:151591/1‑249 (MQ=255)
                                                                                                                                                                             gTTCAGCTCTTCAAAGCGTGCCAGACGTGCTTTACCTTTCGACTGGCGGCCTTTAGTTCCCTGACGTACCCACTCCAGCTCTTTCTCGATCGACTTACGACGCGCCGCTTCTTGTGAAGCTTCCTGCGCCAGGCGCTGATCTTTCTGCTCCAGCCAGGAGGAGTAGTTGCCTTCCCACGGAATACCTTCACCGCGGTCAAGTTCGAGGATCCAGCCTGCAACGTTATCGAGGAAGTAACGGTCGTGGGTa                                                              <  1:405914/250‑1 (MQ=255)
                                                                                                                                                                                        cAAAGCGTGCCAGACGTGCTTTACCTTTCGACTGGCGGCCTTTAGTTCCCTGACGTACCCACTCCAGCTCTTTCTCGATCGACTTACGACGCGCCGCTTCTTGTGAAGCTTCCTGCGCCAGGCGCTGATCTTTCTGCTCCAGCCAGGAGGAGTCGTTGCCTTCCCACGGAATACCTTCACCGCGGTCAAGTTCGAGGATCCAGCCTGCAACGTTATCGAGGAAGTAACGGTCGTGGGTaaacgacaccc                                                    >  2:255313/1‑240 (MQ=255)
                                                                                                                                                                                        cAAAGCGTGCCAGACGTGCTTTACCTTTCGACTGGCGGCCTTTAGTTCCCTGACGTACCCACTCCAGCTCTTTCTCGATCGACTTACGACGCGCCGCTTCTTGTGAAGCTTCCTGCGCCAGGCGCTGATCTTTCTGCTCCAGCCAGGAGGAGTAGTTGCCTTCCCACGGAATACCTTCACCGCGGTCAAGTTCGAGGATCCAGCCTGCAACGTTATCGAGGAAGTAACGGTCGTGGGTAATCGCCACCACg                                                  >  2:38430/1‑251 (MQ=255)
                                                                                                                                                                                          aaGCGTGCCAGACGTGCTTTACCTTTCGACTGGCGGCCTTTAGTTCCCTGACGTACCCACTCCAGCTCTTTCTCGATCGACTTACGACGCGCCGCTTCTTGTGAAGCTTCCTGCGCCAGGCGCTGATCTTTCTGCTCCAGCCAGGAGGAGTAGTTGCCTTCCCACGGAATACCTTCACCGCGGTCAAGTTCGAGGATCCAGCCTGCAACGTTATCGAGGAAGTAACGGTCGTGGGTAATCGCCACCACGGt                                                >  2:373860/1‑251 (MQ=255)
                                                                                                                                                                                            gcgtgcCAGACGTGCTTTACCTTTCGACTGGCGGCCTTTAGTTCCCTGACGTACCCACTCCAGCGCGTTCTCGATCGACTTACGACGCGCCGCTTCTTGTGAAGCTTCCTGCGCCAGGCGCTGATCTTTCTGCTCCAGCCAGGAGGAGTAGTTGCCTTCCCACGGAATACCTTCACCGCGGTCAAGTTCGAGGATCCAGCCTGCAACGTTATCGAGGAAGTAACGGTCGTGGGTAATCGCCACCACGGTGc                                              <  1:464129/247‑1 (MQ=255)
                                                                                                                                                                                              gtgcCAGACGTGCTTTACCTTTCGACTGGCGGCCTTTAGTTCCCTGACGTACCCACTCCAGCTCTTTCTCGATCGACTTACGACGCGCCGCTTCTTGTGAAGCTTCCTGCGCCAGGCGCTGATCTTTCTGCTCCAGCCAGGAGGAGTAGTTGCCTTCCCACGGAATACCTTCACCGCGGTCAAGTTCGAGGATCCAGCCTGCAACGTTATCGAGGAAGTAACGGTCGTGGGTAATCGCCACCACGGTGccc                                            >  1:213356/3‑250 (MQ=255)
                                                                                                                                                                                                 ccAGACGTGTTTTCCCTTTCGACGGGCGGCCTTTATTTCCCTGACGTACCCACTCCAGCTCTTTCTCGTTCGACTTACGACGCGCCGCTTCTTGTGAAGCTTCCTGCGCCAGGCGCTGATCTTTCTGCTCCAGCCAGGAGGAGTAGTTGCCTTCCCACGGAATACCTTCACCGCGGTCAAGTTCGAGGATCCAGCCTGCAACGTTATCGAGGAAGTAACGGTCGTGGGTAATCGCCACCACGGTGCCCTCg                                         <  2:295865/251‑1 (MQ=255)
                                                                                                                                                                                                       gTGCTTTACCTTTCGACTGGCTGCCTTTAGTTCCCTGACGTACCCACTCCAGCTCTTTCTCTATCGACTTACGACGCGCCGCTTCTTGTGAAGCTTTCTGCGCCAGGCGCTGATCTTTCTGCTCCAGCCAGGAGGAGTAGTTGCCTTCCCACGGAATACCTTCACCGCGGTCAAGTTCGAGGATCCAGCCTGCAACGTTATCGAGGAAGTAACGGTCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTCg                                   <  1:219398/251‑1 (MQ=255)
                                                                                                                                                                                                       gTGCTTTACCTTTCGACTGGCGGCCTTTAGTTCCCTGACGTACCCACTCCAGCTCTTTCTCGATCGACTTACGACGCGCCGCTTCTTGTGAAGCTTCCTGCGCCAGGCGCTGATCTTTCTGCTCCAGCCAGGAGGAGTAGTTGCCTTCCCACGGAATACCTTCACCGCGGTCAAGTTCGAGGATCCAGCCTGCAACGTTATCGAGGAAGTAACGGGCGTGGGGAATCGCCACCACGGTGCCCTCGAAGTCg                                   >  1:360194/1‑251 (MQ=255)
                                                                                                                                                                                                           tttACCTTTCGACTGGCGGCCTTTAGTTCCCTGACGTACCCACTCCAGCTCTTTCTCGATCGACTTACGACGCGCCGCTTCTTGTGAAGCTTCCTGCGCCAGGCGCTGATCTTTCTGCTCCAGCCAGGAGGAGTAGTTGCCTTCCCACGGAATACCTTCACCGCGGTCAAGTTCGAGGATCCAGCCTGCAACGTTATCGAGGAAGTAACGGTCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTCGTGc                                >  1:83221/1‑250 (MQ=255)
                                                                                                                                                                                                            ttACCTTTCGACTGGCGGCCTTTAGTTCCCTGACGTACCCACTCCAGCTCTTTCTCGATCGACTTACGACGCGCCGCTTCTTGTGAAGCTTCCTGCGCCAGGCGCTGATCTTTCTGCTCCAGCCAGGAGGAGTAGTTGCCTTCCCACGGAATACCTTCACCGCGGTCAAGTTCGAGGATCCAGCCTGCAACGTTATCGAGGAAGTAACGGTCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTCGTGc                                <  1:421311/249‑1 (MQ=255)
                                                                                                                                                                                                                    cGACTGGCGGCCTTTAGTTCCCTGACGTACCCACTCCAGCTCTTTCTCGATCGACTTACGACGCGCCGCTTCTTGTGAAGCTTCCTGCGCCAGGCGCTGATCTTTCTGCTCCAGCCAGGAGGAGTAGTTGCCTTCCCACGGAATACCTTCACCGCGGTCAAGTTCGAGGATCCAGCCTGCAACGTTATCGAGGAAGTAACGGTCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTCGTGCAGGAAGCGtt                      <  1:125364/251‑1 (MQ=255)
                                                                                                                                                                                                                        tggcGGCCTTTAGTTCCCTGACGTACCCACTCCAGCTCTTTCTCGATCGACTTACGACGCGCCGCTTATTGTGAAGCTTCCTGCGCCAGGCGCTGCTCTTTCTGCTCCAGCCAGGAGGAGTAGTTGCCTTCCCACGGAATACCTTCACCGCGGTCAAGTTCGAGGATCTCGCCTGCAACGTTATCGAGGAAGTAACGGTCGTAGGTAATCGCCACCATGTTGCCATCGAAGGCGTgacgagacggttac                    >  2:400184/4‑236 (MQ=255)
                                                                                                                                                                                                                         ggcGGCCTTTAGTTCCCTGACGTACCCACTCCAGCTCTTTCTCGATCGACTTACGACGCGCCGCTTCTTGTGAAGCTTCCTGCGCCAGGCGCTGATCTTTCTGCTCCAGCCAGGAGGCGTAGTTGCCTTCCCACGGAATACCTTCACCGCGGTCAAGTTCGAGGTTCCAGCCTGCACCGTTATCGGGGAAGTACCGGTCGTGGGTAATCGCCCCCCCGGGGCCCTCGACGTCTGGCAGGAAGCGTTcacgc                 >  2:35395/3‑247 (MQ=255)
                                                                                                                                                                                                                             gCCTTTAGTTCCCTGACGTACCCACTCCAGCTCTTTCTCGATCGACTTACGACGCGCCGCTTCTTGTGAAGCTTCCTGCGCCAGGCGCTGATCTTTCTGCTCCAGCCAGGAGGAGTAGTTGCCTTCCCACGGAATACCTTCACCGCGGTCAAGTTCGAGGATCCAGCCTGCAACGTTATCGAGGAAGTAACGGTCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAgg             >  1:262514/1‑251 (MQ=255)
                                                                                                                                                                                                                                       ccctGACGTACCCACTCCAGCTCTTTCTCGATCGACTTACGACGCGCCGCTTCTTGTGAAGCTTCCTGCGCCAGGCGCTGATCTTTCTGCTCCAGCCAGGAGGAGTAGTTGCCTTCCCACGGAATACCTTCACCGCGGTCAAGTTCGAGGATCCAGCCTGCAACGTTATCGAGGAAGTAACGGTCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTc   >  1:421439/4‑251 (MQ=255)
                                                                                                                                                                                                                                        cctGACGTACCCACTCCAGCTCTTTCTCGATCGACTTACGGCGCGCCGCTTCTTGTGAAGCTTCCTGCGCCAGGCGCTGATCTTTCTGCTCCAGCCAGGAGGAGTAGTTGCCTTCCCACGGAATACCTTCACCGCGGTCAAGTTCGAGGATCCAGCCTGCAACGTTATCGAGGAAGTAACGGTCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCg  >  1:262862/3‑251 (MQ=255)
                                                                                                                                                                                                                                        cctGACGTACCCACTCCAGCTCTTTCTCGATCGACTTACGACGCGCCGCTTCTTGTGAAGCTTCCTGCGCCAGGCGCTGATCTTTCTGCTCCAGCCAGGAGGAGTAGTTGCCTTCCCACGGAATACCTTCACCGCGGTCAAGTTCGAGGATCCAGCCTGCAACGTTATCGAGGAAGTAACGGTCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCg  <  1:169027/249‑1 (MQ=255)
                                                                                                                                                                                                                                        cctGACGTACCCACTACAGCTCTTTCTCGATCGACTTACGACGCGCCGCTTCTTGTGAAGCTTCCTGCGCCAGGCGCTGATCTTTCTGCTCCAGCCAGGAGGAGTAGTTGCCTTCCCACGGAATACCTTCACCGCGGTCAAGTTCGAGGATCCAGCCTGCAACGTTATCGAGGAAGTAACGGTCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCg  >  1:262863/3‑251 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                             
CCGATGATCCCGACGATCGCTCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGACGATCGCCATAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGTTCTTCAAAGCGCGCCAGACGTGCTTTACCTTTCGACTGACGGCCTTTAGTACCTTGACGTACCCATTCCAGCTCTTTCTCAATCGACTTACGACGCGCCGCTTCTTGTGAAGCTTCCTGCGCCAGGCGCTGATCTTTCTGCTCCAGCCAGGAGGAGTAGTTACCTTCCCACGGAATACCTTCACCGCGGTCAAGTTCGAGGATCCAGCCTGCAACGTTATCGAGGAAGTAACGGTCGTGGGTAATCGCCACAACGGTGCCTTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCG  >  NC_000913/4629456‑4629938

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: