Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,629,515 A→G R336R (CGT→CGC 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,5150AG98.6% 213.2 / ‑4.8 69R336R (CGT→CGCettAenergy‑dependent translational throttle protein EttA
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (35/33);  total (36/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 = 9.40e-01

ATCTTCATGATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGACGCCAGTTTCACCGTTTCACCCAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCACCGTTCGGACCGATGATCCCGACGATCGCTCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGACGATCGCCATAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGTTCTTCAAAGCGCGCCAGACGTGCTTTACCTTTCGACTGACGGCCTTTAGTACCTTGACGTACCCATTCCAGCTCTTTCTCAATCGACTTACGACGCGCCGCTTCTTGTGAAGCTTCCTGCGC  >  NC_000913/4629273‑4629757
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                  
accTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCGCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATcacc                                                                                                                                                                                                                                            <  1:365662/249‑4 (MQ=255)
aTCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATcacc                                                                                                                                                                                                                                            >  1:188214/1‑248 (MQ=255)
        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)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                  
ATCTTCATGATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGACGCCAGTTTCACCGTTTCACCCAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCACCGTTCGGACCGATGATCCCGACGATCGCTCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGACGATCGCCATAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGTTCTTCAAAGCGCGCCAGACGTGCTTTACCTTTCGACTGACGGCCTTTAGTACCTTGACGTACCCATTCCAGCTCTTTCTCAATCGACTTACGACGCGCCGCTTCTTGTGAAGCTTCCTGCGC  >  NC_000913/4629273‑4629757

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: