Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,629,521 G→A G334G (GGC→GGT 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,5210GA97.1% 202.8 / ‑6.2 70G334G (GGC→GGTettAenergy‑dependent translational throttle protein EttA
Reads supporting (aligned to +/- strand):  ref base G (1/0);  major base A (34/33);  minor base T (2/0);  total (37/33)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGACGCCAGTTTCACCGTTTCACCCAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCACCGTTCGGACCGATGATCCCGACGATCGCTCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGACGATCGCCATAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGTTCTTCAAAGCGCGCCAGACGTGCTTTACCTTTCGACTGACGGCCTTTAGTACCTTGACGTACCCATTCCAGCTCTTTCTCAATCGACTTACGACGCGCCGCTTCTTGTGAAGCTTCCTGCGCCAGGCGCTG  >  NC_000913/4629281‑4629766
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                     
aattcccAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTGTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGAtt                                                                                                                                                                                                                                             <  2:287072/246‑1 (MQ=255)
  tatCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGACGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGGCGATTTACCCGCGCCGTTCGGACCGATGATCCCTACGATTTCCCCTTTCGTGGTCGGGAAGCGCCGGTCGTCAATACGCAGGCGATCTCCGTAGGtttt                                                                                                                                                                                                                                            >  2:219398/1‑250 (MQ=255)
      cAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGCGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCAACGTAGGATTTACgac                                                                                                                                                                                                                                       >  2:93500/1‑249 (MQ=255)
       accccGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAg                                                                                                                                                                                                                                      <  2:343219/249‑1 (MQ=255)
        gCCCGCCGGAAACGTCTTCCCAAACGGATTGGCTGTTATCCATTGAGTCACGGAACGGCTCAACCGATGCCAGTGTCACCGTTTCACCGAAAGTGAGGGTGCCGCTGGCCGGCTGTTCCTGTCCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTGCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCGCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAgg                                                                                                                                                                                                                                     <  2:332707/251‑1 (MQ=255)
           cgccgGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTg                                                                                                                                                                                                                                  <  1:424237/251‑1 (MQ=255)
            gccgGACACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGCTCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGc                                                                                                                                                                                                                                 >  2:24907/1‑251 (MQ=255)
                 aaaCTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGACGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTgtct                                                                                                                                                                                                                            >  2:125364/1‑248 (MQ=255)
                 aaaCTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACt                                                                                                                                                                                                                            <  1:388173/251‑1 (MQ=255)
                   aCTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCACCCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGACCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTc                                                                                                                                                                                                                          >  2:24316/1‑251 (MQ=255)
                         tcCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAg                                                                                                                                                                                                                       <  1:444191/248‑1 (MQ=255)
                         tcCCAAACGGTTTTGCTGTTAGCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCGCCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGc                                                                                                                                                                                                                      <  2:113038/249‑1 (MQ=255)
                          cccAAACGGTTGTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACt                                                                                                                                                                                                                   <  1:469342/251‑1 (MQ=255)
                           ccAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACtt                                                                                                                                                                                                                  >  1:178472/1‑251 (MQ=255)
                              aaCGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTAt                                                                                                                                                                                                               >  2:476232/1‑251 (MQ=255)
                              aaCGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATAGAGAAGCTCAGGTCATCAATCAGCAGGCGATCAACGTCGGATTTACGCAGGTTGCTGACTTCCAGCACtttct                                                                                                                                                                                                               >  1:319025/1‑249 (MQ=255)
                                cGGTTTTGCTGTTATCCAGTGAGTCACGGAACGGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCg                                                                                                                                                                                                             <  2:395025/251‑1 (MQ=255)
                                      tGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGCCCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGCTTGCCCCTTTCGGGATCGCGAAGCTCAGGTCATCAATCAGCAGGCGCTCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTCTCGCCCAgc                                                                                                                                                                                                       >  2:172466/1‑250 (MQ=255)
                                                    gAGTCCCGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGGCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCGCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAtgtg                                                                                                                                                                                          <  2:339060/250‑1 (MQ=255)
                                                          cGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGAGCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATTGCCAAGACGCGGTCCAGGTGGAATaa                                                                                                                                                                                    >  1:55808/1‑250 (MQ=255)
                                                           ggAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAACGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGACCAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGCTCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAAc                                                                                                                                                                                  >  2:437661/1‑251 (MQ=255)
                                                             aaCTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAAc                                                                                                                                                                                  >  1:58464/1‑249 (MQ=255)
                                                                 gATCAACCGATGCCAGTTTCCCCGTTCCCCCCAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTCCCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGtt                                                                                                                                                                              <  2:167771/249‑1 (MQ=255)
                                                                    cAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGTTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGtt                                                                                                                                                                          <  2:105289/250‑1 (MQ=255)
                                                                      aCCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGt                                                                                                                                                                       >  1:282156/1‑251 (MQ=255)
                                                                       ccGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGtt                                                                                                                                                                      <  1:480315/251‑1 (MQ=255)
                                                                        cgatgCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGTCCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGttt                                                                                                                                                                     <  1:107935/247‑1 (MQ=255)
                                                                                     accgtttcaccTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGAGCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTg                                                                                                                                                        <  1:355752/251‑1 (MQ=255)
                                                                                       cgtttcaccTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGat                                                                                                                                                      <  1:122637/251‑1 (MQ=255)
                                                                                        gtttcaccTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGACTAAACAGTTCGTTGGTTTCGTTACGTTTCTGat                                                                                                                                                      >  2:300696/1‑250 (MQ=255)
                                                                                        gtttcaccCAAGGTGATGGTGCCGCTGTCCGTCTGTTCCTAACCAGAGATCATACGGAACAGGGTCGATTTACCCGCACCGTTCAGACCGATGAACCCGACGATCGCTCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGACGATCTCCATAGGATTTACGCAGGTTGCTGCCTTCCAGCACTTTATCGACCAGACGCGGTCCAGGCGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGat                                                                                                                                                      >  1:229818/1‑250 (MQ=255)
                                                                                                   aGTGATGGTGCCGATGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTAGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGAAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGAGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTACTGt                                                                                                                                          >  1:173250/1‑251 (MQ=255)
                                                                                                                        ctgtacatGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGCTCGAGAAGCTCAAGCCATCAATCAGCAGGCGATCACCGGAGGATTTACGCAGGATGCTGTCTTCCAGCACTTTATCGCCCAGACGCGGTCCAAGAGGAATAAACAGTTCGTTGGTTTCGTTACGGTCCTGATATTCGGTGCTGTTCAGctcttcctagcttgccc                                                                                                                     >  2:115098/8‑240 (MQ=255)
                                                                                                                          gTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGCTCTTCAAAGCGTGCCa                                                                                                                     <  1:300696/249‑1 (MQ=255)
                                                                                                                                 aCCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGACCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGCTCTTCAAAGCGTGCCAGACGTGCtt                                                                                                            >  2:169027/1‑251 (MQ=255)
                                                                                                                                   cAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGCTCTTCAAAGCGTGCCAGACGTGCttt                                                                                                           >  1:295865/1‑250 (MQ=255)
                                                                                                                                   cAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGCTCTTCAAAGCGTGCCAGACGTGCttt                                                                                                           >  2:195578/1‑250 (MQ=255)
                                                                                                                                          cATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGCTCTTCAAAGCGTGCCAGACGTGCTTTACCTTTCg                                                                                                   >  2:176877/1‑251 (MQ=255)
                                                                                                                                              cGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGCTGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGACTAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGCTCTTCAAAGCGTGCCAGACGTGCTTTACCTTTCTACTg                                                                                               >  2:167989/1‑251 (MQ=255)
                                                                                                                                              cGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGCTCTTCAAAGCGTGCCAGACGTGCTTTACCTTTCGACTg                                                                                               >  2:405914/1‑251 (MQ=255)
                                                                                                                                               ggatcAGGGTCGATTTTCGCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCATGTTGCTTACTTCCAGCTCTTTTTCGCCCAGACGCGGTTCAGTTGGAATAAACAGTTCGTTGGTTTCTTTACGTTTCTGATATTCGTTGCTGTTCAGCTCTTCAAAGCGTGCCAGACGTGCTTTACCTTTCGACTgg                                                                                              <  1:449982/247‑2 (MQ=255)
                                                                                                                                                 aaCAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCCGGTGGAATAAACAGTTCGTTGGTTTCGTTa                                                                                                                                                                >  2:145865/1‑183 (MQ=255)
                                                                                                                                                 aaCAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTa                                                                                                                                                                <  1:145865/183‑1 (MQ=255)
                                                                                                                                                 aaCAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGTTATTCGGTGCTGTTCAGCTCTTCAAAGCGTGCCAGACGGGCTTTCCCTTTCGACTggcg                                                                                            >  2:123459/1‑248 (MQ=255)
                                                                                                                                                    aGGGTCGATTTACCCGCGCCGTGCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCGCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGCTCTTCAAAGCGTGCCAGACGTGCTTTACCTTTCGACTGGCGGcc                                                                                         <  2:260529/251‑1 (MQ=255)
                                                                                                                                                      ggTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGCTCTTCAAAGCGTGCCAGACGTGCTTTACCTTTCGACTGGCGGCCtt                                                                                       >  1:279155/1‑251 (MQ=255)
                                                                                                                                                               aCCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGCTCTTCAAAGCGTGCCAGATGTGCTTTACCTTTCGACTGGCGGCCTTTAGttccct                                                                              >  1:347215/1‑246 (MQ=255)
                                                                                                                                                                   gggccGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGCTCTTCAAAGCGTGCCAGACGTGCTTTACCTTTCGACTGGCGGCCTTTAGTTCCCTGACg                                                                          <  2:234097/248‑1 (MQ=255)
                                                                                                                                                                          tCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGCTCTTCAAAGCGTGCCAGACGTGCTTTACCTTTCGa                                                                                                  >  2:104835/1‑220 (MQ=255)
                                                                                                                                                                          tCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGCTCTTCAAAGCGTGCCAGACGTGCTTTACCTTTCGa                                                                                                  <  1:104835/220‑1 (MQ=255)
                                                                                                                                                                               ccGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGCTCTTCAAAGCGTGCCAGACGTGCTTTACCTTTCGACTGGCGGCCTTTAGTTCCCTGACGTACCCACTCCAg                                                              <  1:477414/251‑1 (MQ=255)
                                                                                                                                                                                   tgatCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGCTCTTCAAAgcgtg                                                                                                                        >  2:170215/1‑187 (MQ=255)
                                                                                                                                                                                   tgatCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGCTCTTCAAAgcgtg                                                                                                                        <  1:170215/189‑3 (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)
                                                                                                                                                                                                                                            gATCGCCATAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGATCTTCAAAGCGCGCCAGACGTGCTTTACCTTTCGACTGACGGCCTTTAGTACATTGACGTACCCATTCCAGCTCTTTCTCAATCGACTTACGACGCGACGCTTCTTGTGAAGCTTCCTGCGCCAGTCGCTg  >  1:105327/1‑250 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                     
GATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGACGCCAGTTTCACCGTTTCACCCAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCACCGTTCGGACCGATGATCCCGACGATCGCTCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGACGATCGCCATAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGTTCTTCAAAGCGCGCCAGACGTGCTTTACCTTTCGACTGACGGCCTTTAGTACCTTGACGTACCCATTCCAGCTCTTTCTCAATCGACTTACGACGCGCCGCTTCTTGTGAAGCTTCCTGCGCCAGGCGCTG  >  NC_000913/4629281‑4629766

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