Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,629,188 G→A L445L (CTC→CTT 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,1880GA98.6% 216.4 / NA 70L445L (CTC→CTTettAenergy‑dependent translational throttle protein EttA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (31/38);  minor base C (1/0);  total (32/38)
Fisher's exact test for biased strand distribution p-value = 4.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

ACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAATCCAGAATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAGAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATGATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGACGCCAGTTTCACCGTTTCACCCAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTC  >  NC_000913/4628948‑4629434
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                      
aCCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCa                                                                                                                                                                                                                                               <  1:501408/250‑1 (MQ=255)
aCCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCa                                                                                                                                                                                                                                               >  2:203503/1‑250 (MQ=255)
    tCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAAcg                                                                                                                                                                                                                                            <  1:433219/249‑1 (MQ=255)
     cGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAAcgcg                                                                                                                                                                                                                                          <  1:557084/250‑1 (MQ=255)
     cGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGt                                                                                                                                                                                                                                         >  1:120414/1‑251 (MQ=255)
     cGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGAGGGGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAAcgcg                                                                                                                                                                                                                                          <  2:554936/250‑1 (MQ=255)
      gaagaaCTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGtt                                                                                                                                                                                                                                        >  2:436596/1‑251 (MQ=255)
          aaCTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTAcc                                                                                                                                                                                                                                    >  1:409890/1‑251 (MQ=255)
              cAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCt                                                                                                                                                                                                                                  <  2:484227/249‑1 (MQ=255)
                   ttACCTTTATCCTGGGAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGAAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCacc                                                                                                                                                                                                                           >  2:62545/1‑249 (MQ=255)
                   ttACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCa                                                                                                                                                                                                                             <  2:236384/249‑1 (MQ=255)
                             ccTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTttaa                                                                                                                                                                                                                   >  1:308360/1‑249 (MQ=255)
                              cTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGt                                                                                                                                                                                                                                         <  1:216120/226‑1 (MQ=255)
                              cTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGt                                                                                                                                                                                                                                         >  2:216120/1‑226 (MQ=255)
                                 gtagtCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTttaaagtta                                                                                                                                                                                                              >  2:201138/5‑250 (MQ=255)
                                   agtCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAACGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTttaaagttaaa                                                                                                                                                                                                            <  1:298647/248‑1 (MQ=255)
                                       caggaTGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGc                                                                                                                                                                                                       >  1:399878/5‑251 (MQ=255)
                                        aggaggtgCGGGGCGAGACGGTTGAGGAACTAACGGGCGTGCGAGATCACCATCGCACAGCCAGGGAACACCCGTAGGGCGCTTTCCAGCGCGCGCAGGGTTGCGATATCCAGGTCGTTGGTTGGTTCGTCGAGTAGCAGCATGTTGTCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGATCCCCTCCGGAAAGTTCACCAACGCGTTGACCCGGATCAACCCCTGTAAATTTAAAGCGGct                                                                                                                                                                                                      <  2:54991/246‑2 (MQ=255)
                                          gaTGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCaa                                                                                                                                                                                                    <  1:259165/250‑1 (MQ=255)
                                                 gTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCCGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGTTATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCCGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCCACCCCTTTAAGGTTAAAGCGGCCAACGTAGgc                                                                                                                                                                                             >  2:627079/1‑251 (MQ=255)
                                                  tGGCGATACGGTCGAGGAACCAACGGTCGTGCGCGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGgcg                                                                                                                                                                                            >  2:143173/1‑251 (MQ=255)
                                                       aTACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCt                                                                                                                                                                                       >  2:368924/1‑251 (MQ=255)
                                                        tACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCtt                                                                                                                                                                                      >  1:262269/1‑251 (MQ=255)
                                                         aCGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCtt                                                                                                                                                                                      <  1:62907/250‑1 (MQ=255)
                                                         aCGGGCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCtt                                                                                                                                                                                      <  1:17750/250‑1 (MQ=255)
                                                              cGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCAtc                                                                                                                                                                                <  1:247283/251‑1 (MQ=255)
                                                                aGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCACCGCGCGCAGGGTTTCGATATCCACGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCACCGCAGGCGCGGCTTGGCAtctc                                                                                                                                                                              >  2:614992/1‑251 (MQ=255)
                                                                 ggAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCg                                                                                                                                                                             <  1:345073/251‑1 (MQ=255)
                                                                    aCCAACGGTCGGGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGgtg                                                                                                                                                                          <  2:409890/251‑1 (MQ=255)
                                                                       atggggggTGCGAGATAACCATCGCACAGCCCGGGAACTGCAGCAGTGCGTTTCCCAGCGCGCGCAGGGTTTCGATATCCATGTCGTTGGTTGGGTCGTCCAGCAGGAGGATGCTGCCGCCAACCTGCAGCAGCTTCGCCGGAGGCAGACGACCGCGTGCACCACCGGAAAGTGCACGAACGCGTTGACGCTGATCAAGCCCTTTAAAGTTATATCGTGCAACGTAGGCGTGGCGTGGCATCGCGGTGTTg                                                                                                                                                                       <  2:298096/244‑1 (MQ=255)
                                                                            gcGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGAt                                                                                                                                                                  <  1:520315/250‑1 (MQ=255)
                                                                             cGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCTGTGTTGCCGAt                                                                                                                                                                  <  2:38127/250‑1 (MQ=255)
                                                                                 cGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGGCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGGTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCa                                                                                                                                                             <  1:58016/251‑1 (MQ=255)
                                                                                       aaCCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAatat                                                                                                                                                       >  2:53393/1‑251 (MQ=255)
                                                                                       aaCCATCGAACAGCCCGGGAACTCCAGCAGCGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGAAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCCACGCGTTTACCCTGATCACCCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAatat                                                                                                                                                       >  2:292089/1‑251 (MQ=255)
                                                                                             cgaacaGCGCGGTAACTCCTGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATTCTGGGCGTTGGGTGGGTCGGCGAGCAGCAGCATGTTGCCGCCAACCTGCAGGAGCCTCGCCAGATGCAGACGCCCGCGTGCACCACCGGAAAGTTCATCAACGCGGTTATCCGGATCAACCTCTTTAAAGTTAAAGCGGCCAACGTATGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGcc                                                                                                                                                 <  1:592389/248‑1 (MQ=255)
                                                                                             cGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGcc                                                                                                                                                 <  1:592420/251‑1 (MQ=255)
                                                                                                  aGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCccgccg                                                                                                                                            <  1:215973/251‑1 (MQ=255)
                                                                                                             tCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCACAGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCAAGCCCGCCGGAAACttcttc                                                                                                                                 >  2:2584/1‑251 (MQ=255)
                                                                                                              ccAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACttcttc                                                                                                                                 <  1:353722/250‑1 (MQ=255)
                                                                                                              ccAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGGCGTTGGTTTGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCTGTGTTGCCGATCTTCATATTATCCAGCCCGCCGGAAACTTCTTcc                                                                                                                                <  2:559128/251‑1 (MQ=255)
                                                                                                                   agGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTGCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTGCATAATATCCAGCCCGCCGGACTCTTCTTCCCAAAc                                                                                                                           <  1:350501/251‑1 (MQ=255)
                                                                                                                    gggCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACg                                                                                                                          >  1:40771/1‑251 (MQ=255)
                                                                                                                        gTTTTCCAGAGCGCGCAGGGTTTCGATATACAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACgggtt                                                                                                                      >  2:404662/1‑248 (MQ=255)
                                                                                                                                gcgcgcgcAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGttt                                                                                                              >  2:274462/1‑250 (MQ=255)
                                                                                                                                       cAGGGTTTCGATATCCAGGTCGATGGTTGGTTCGACGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTg                                                                                                       <  2:194317/251‑1 (MQ=255)
                                                                                                                                        aGGGTTTCGCTAGCCAGGGCGTTGGTGGGCTCGTCGAGCAGCAGCATGTTGCCGCCAACCTTCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCCACGCGTTTACCCTGAGCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCTGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCGTCCCAAACGGTTTTGCTGTTATCCATTGa                                                                                                      <  2:444829/251‑1 (MQ=255)
                                                                                                                                           gTTTCGATATCCAGGTCGTTGGTTGGTTCGTCTAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGGTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGTAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTc                                                                                                   >  2:146224/1‑251 (MQ=255)
                                                                                                                                             ttCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCa                                                                                                  <  1:368924/250‑1 (MQ=255)
                                                                                                                                              tCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCAc                                                                                                 <  1:68144/250‑1 (MQ=255)
                                                                                                                                                 atatCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGa                                                                                              >  1:135580/1‑250 (MQ=255)
                                                                                                                                                    tCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGGAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGATGTTATCCATTGAGCCACggtaatg                                                                                          >  2:344292/1‑246 (MQ=255)
                                                                                                                                                            gttgggtggtgctgcgaccagcagcaTGTTGCCGCCAACCGGCAGCAGTTTCGCCAGATGCAGACGACCGCGCTTACCACCGGAAAGGTCACCAACGCGTTTCCCCTGATCAACCCCTGTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATAGCCAGCCCGCCGTAAACTACTGCCCAAACGGGTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCg                                                                                  <  2:518260/233‑1 (MQ=255)
                                                                                                                                                                         tcgAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACAAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCCCGGAACTGATCAACCGATGCCAgtttcac                                                                     >  1:39338/1‑251 (MQ=255)
                                                                                                                                                                          cgAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAggttcacc                                                                    <  1:513705/251‑1 (MQ=255)
                                                                                                                                                                           gAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAgtttcaccg                                                                   >  1:42584/1‑251 (MQ=255)
                                                                                                                                                                              cagcagcaTGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGCCGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACGTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAgtttcaccgttt                                                                <  1:547331/251‑1 (MQ=255)
                                                                                                                                                                                 cagcaTGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCAATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAgtttcaccgtttcac                                                             <  1:203503/251‑1 (MQ=255)
                                                                                                                                                                                 cagcaTGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAACCGTCCAACGTAGGCGCGGCTTTGCATCTCGGTGTTGCCTATCTTCATAATATCCAGCCCGCCGGACACTTCTTCCCACACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGTTCAACCGATGCCAgtttcaccgtttcac                                                             >  1:203855/1‑251 (MQ=255)
                                                                                                                                                                                  agcaTGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCAACTCGGTGTTGCCGATCTTCATAATAACCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAgtttcaccgtttcacc                                                            <  1:627079/251‑1 (MQ=255)
                                                                                                                                                                                   gcaTGTTGCCGCCAACCGGCAGCAGCGTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCGTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAgtttcaccgtttcacct                                                           <  1:248854/251‑2 (MQ=255)
                                                                                                                                                                                                   cTGCAGCAGCTTCGCCAGCTGCAGACGACCGAGCTAACCACCGGACCGTTCACCAACGCGTTTACCCTGCTCACCCCCTTTAAAGATAAAGCGGACAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGGTCTTCATCATATCCACCCAGCCGGCAACTTCTTCCCAAACGGTCTTGCTGTTCTCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGc                                           >  2:571838/1‑251 (MQ=255)
                                                                                                                                                                                                   cTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAACGATCACCAACGCGTATACCCTGATCAACCAATTTACAGATAAAGCGGACCACGTAAGCGCGGATTGGAATCTCGGAGTTGCCGATCTTCATAATATCCAGCACGCCGGAAACTTCTTCCCAAACGTTTTTGCTGTTATCCATTGAGTCACGGACCTGATCAACCGATGCCAGTTTCACCGTTTCACCTACAGTGTTGGTGCCGc                                           >  1:302338/1‑251 (MQ=255)
                                                                                                                                                                                                          agcTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCgg                                    >  2:199122/1‑251 (MQ=255)
                                                                                                                                                                                                                         aGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGAcc                        >  1:77788/1‑248 (MQ=255)
                                                                                                                                                                                                                             gacCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATc                 <  1:201138/251‑1 (MQ=255)
                                                                                                                                                                                                                                      cgccaccGGAAAGTTCACCACCGCGTTTACCCTGATCAACTCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCTGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCATCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCTATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTTCCGCTGTCCGGATGTTCCTGAACAGAGATCATACGGaa         <  1:10176/248‑1 (MQ=255)
                                                                                                                                                                                                                                      caccacCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCGCGGTGTGGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTTATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGaa         <  1:62545/250‑1 (MQ=255)
                                                                                                                                                                                                                                       accaccGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAAc        >  1:496365/1‑250 (MQ=255)
                                                                                                                                                                                                                                            cGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCAAGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTACAGTGCTGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCCTACGGAACAGGGTc  >  2:443089/1‑251 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                      
ACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAATCCAGAATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAGAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATGATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGACGCCAGTTTCACCGTTTCACCCAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTC  >  NC_000913/4628948‑4629434

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