Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,629,356 C→T A389A (GCG→GCA 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,3560CT98.2% 173.4 / NA 57A389A (GCG→GCAettAenergy‑dependent translational throttle protein EttA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (27/29);  minor base A (0/1);  total (27/30)
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.94e-01

TCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAGAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATGATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGACGCCAGTTTCACCGTTTCACCCAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCACCGTTCGGACCGATGATCCCGACGATCGCTCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGACGATCGCCATAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATA  >  NC_000913/4629117‑4629587
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                       
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)
                                                                      aaagTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGAt                                                                                                                                                        >  2:44225/3‑251 (MQ=255)
                                                                         gTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCGGACCAGAGATCATACGGAACAGGGTCGAttt                                                                                                                                                      >  2:494630/1‑250 (MQ=255)
                                                                                   gcgTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAACGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGt                                                                                                                                           >  2:317541/1‑251 (MQ=255)
                                                                                                    ccccTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTTTTGCCGATCTTCATAATATCCAGCCCGCCGGAATCTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCGGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCcg                                                                                                                          <  1:443089/251‑1 (MQ=255)
                                                                                                      ccTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCcgac                                                                                                                        <  1:274462/251‑1 (MQ=255)
                                                                                                       cTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAACCTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCcgacg                                                                                                                       >  1:231716/1‑251 (MQ=255)
                                                                                                             agttaaagCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACGTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGAGTTACCCGCGCCGTTCGGACCGATGATCCCTACGAttgcc                                                                                                                 <  2:414029/251‑5 (MQ=255)
                                                                                                               ttaaagCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGAttgcccc                                                                                                               >  1:324081/1‑245 (MQ=255)
                                                                                                               ttaaagCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGGTGACGATCTTCATAATATCCAGCACGCCGGCAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGTTCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAAAAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGAttgcccc                                                                                                               >  2:417170/1‑245 (MQ=255)
                                                                                                                 aaagCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCtt                                                                                                             <  1:395388/251‑1 (MQ=255)
                                                                                                                  acgCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGGTGATCCCGCCGATTGCCCCttt                                                                                                            <  2:120414/249‑1 (MQ=255)
                                                                                                                   agctgctaacggatgcgcGGCTTGGCATCACGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCGTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATTCCAGTTTCACCGTTTCACCTAAAGTGATGGGGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCTAGTTACCCGCGCCGTTCGGACCGATGATCCcgacgagtgcccctttc                                                                                                           <  2:102436/237‑12 (MQ=255)
                                                                                                                     cGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGTGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCgg                                                                                                         <  2:42584/251‑1 (MQ=255)
                                                                                                                         cAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATc                                                                                                     <  2:308360/251‑1 (MQ=255)
                                                                                                                             gTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGCTGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTAACCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCg                                                                                                    >  2:76152/1‑248 (MQ=255)
                                                                                                                             gTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCg                                                                                                    >  2:141531/1‑248 (MQ=255)
                                                                                                                             gTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTGCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCAGACGGAACAGGGTCGATTTACCCGCGCCGGTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCgaga                                                                                                 <  1:232791/251‑1 (MQ=255)
                                                                                                                               aGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATc                                                                                                     <  1:160637/245‑1 (MQ=255)
                                                                                                                               aGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATc                                                                                                     >  2:160637/1‑245 (MQ=255)
                                                                                                                                   gcgGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCa                                                                                           <  1:53393/251‑1 (MQ=255)
                                                                                                                                     ggCTTGGCATCTCGGTGTTGCCGATCTTAATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGCTGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAgg                                                                                         >  2:228343/1‑251 (MQ=255)
                                                                                                                                         tGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGTGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGtca                                                                                      <  1:614992/250‑1 (MQ=255)
                                                                                                                                                       tGCCGATCTTCATAATAGCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATcagcaggcg                                                                        <  1:474314/250‑4 (MQ=255)
                                                                                                                                                       tGCCGATCTTCAGAATATCCAGCCCGCCGGAACCTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGAGCGAGAAGCTCAGGTCATCAATcagcaggcg                                                                        <  2:488235/250‑4 (MQ=255)
                                                                                                                                                                tCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTAAGGTCATCAATCAGCAGGCGATca                                                                    <  1:146224/245‑2 (MQ=255)
                                                                                                                                                                     atatCCAGCCCGCCGGCAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGACGATCACCGTAGGAttt                                                         >  1:205850/1‑251 (MQ=255)
                                                                                                                                                                       atCCAGCCCGCCGGAAACTTCTACCCAAACGGATTGGCGGGTATCCATTGAGTCACGGAACTTAACAACCGATGCCAGTTTCACCGTTTCACCTAAATTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGTACAGGGTAGATGTGCCCGCGTCGGTCGGACCGATGATCCCGTCGAGTTCCCCTTTCGGGAACGAGAAGCTCAGATCATCAATCAGCAGGCGTTTACCGTAGGATTTAc                                                       <  2:57901/251‑1 (MQ=255)
                                                                                                                                                                            gCCCTCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGGGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAgg                                                  <  2:25965/251‑1 (MQ=255)
                                                                                                                                                                                     aaaCTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACt                                         <  1:2584/251‑1 (MQ=255)
                                                                                                                                                                                      aaCTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACtt                                        <  2:77788/251‑1 (MQ=255)
                                                                                                                                                                                          tcttcCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCa                                     <  1:143173/250‑1 (MQ=255)
                                                                                                                                                                                                   aCGGTTTTGCGGTTATCCATTGAGGCACGGAACTGATCAACCGAATCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCGGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGGTCGGACCGATGAGCCCGACGAGTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATc                           <  2:230386/251‑1 (MQ=255)
                                                                                                                                                                                                     ggTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGc                         >  1:344944/1‑251 (MQ=255)
                                                                                                                                                                                                             tGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGAcgc                 >  1:440325/1‑251 (MQ=255)
                                                                                                                                                                                                                   ccATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTccc           >  1:164980/1‑250 (MQ=255)
                                                                                                                                                                                                                   ccATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGAAGGTTGCTGACTTCCAGCACTTTATCGCCAAGACGCGGTccc           >  2:540562/1‑250 (MQ=255)
                                                                                                                                                                                                                          gTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGa     >  1:272003/1‑250 (MQ=255)
                                                                                                                                                                                                                             aCGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATa  >  2:488123/1‑250 (MQ=255)
                                                                                                                                                                                                                              cGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCa                                     >  1:93535/1‑214 (MQ=255)
                                                                                                                                                                                                                              cGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCa                                     <  2:93535/214‑1 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                       
TCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAGAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATGATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGACGCCAGTTTCACCGTTTCACCCAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCACCGTTCGGACCGATGATCCCGACGATCGCTCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGACGATCGCCATAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATA  >  NC_000913/4629117‑4629587

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: