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,1880GA96.0% 231.5 / NA 76L445L (CTC→CTTettAenergy‑dependent translational throttle protein EttA
Reads supporting (aligned to +/- strand):  ref base G (0/2);  major base A (36/35);  minor base C (3/0);  total (39/37)
Fisher's exact test for biased strand distribution p-value = 2.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAATCCAGAATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAGAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATGATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGACGCCAGTTTCACCGTTTCACCCAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTC  >  NC_000913/4628946‑4629434
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                      
ttACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCAcc                                                                                                                                                                                                                                                >  2:428778/1‑251 (MQ=255)
   ccgtCGAAGAACACAATTTTACCTGCACCCTGGTAGTCCAGGATGTGCGTGGCGCTTCGGTCGAGGAACCAACGGTCGTGTGAGAGAACCATCGCACAGCCCGGTAACGCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCATCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGTTCACCCCCGGAAAGTTCACCAAc                                                                                                                                                                                                                                             <  2:67165/248‑1 (MQ=255)
   ccTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAAc                                                                                                                                                                                                                                             <  1:235967/251‑1 (MQ=255)
   ccTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAAc                                                                                                                                                                                                                                             >  1:388400/1‑251 (MQ=255)
      tCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAAcgc                                                                                                                                                                                                                                           <  2:309921/250‑1 (MQ=255)
            aaCTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTAcc                                                                                                                                                                                                                                    <  2:391172/251‑1 (MQ=255)
                 aaCTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGAt                                                                                                                                                                                                                               <  2:178866/251‑1 (MQ=255)
                     ttACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAAc                                                                                                                                                                                                                           >  2:81027/1‑251 (MQ=255)
                     ttACCTTCAACCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCCGCAGCTTCGCCAGATGCAGACGACCGCGCTCCCCACCGGAACGTTCACCAACGCGTTTCCCCTGATCAAc                                                                                                                                                                                                                           >  2:197030/1‑251 (MQ=255)
                          ttCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCCTGAGACATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTCTTCAAGAGCGCGCAGGGTTTCGATATCCAGGACGTTGGTTGGTTCGTCGAGCAGCAGCATGATGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGCCCGCGCTCAACACCGGACAGTTCACCAACGCGTTTACCCTGATCAACCCCtt                                                                                                                                                                                                                      >  1:214858/1‑251 (MQ=255)
                           tCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCt                                                                                                                                                                                                                       <  1:206899/249‑1 (MQ=255)
                            cATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCttt                                                                                                                                                                                                                     >  1:470675/1‑250 (MQ=255)
                            cATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCttt                                                                                                                                                                                                                     <  1:407893/250‑1 (MQ=255)
                              tCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGtt                                                                                                                                                                                                                                                    >  1:127493/1‑217 (MQ=255)
                                     agtCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTttaaagttaaag                                                                                                                                                                                                           <  1:428778/249‑1 (MQ=255)
                                     agtCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTttaaagttaa                                                                                                                                                                                                             >  1:472339/3‑249 (MQ=255)
                                      gtCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGc                                                                                                                                                                                                          >  2:195174/2‑251 (MQ=255)
                                        ccAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCgg                                                                                                                                                                                                        >  1:414594/1‑251 (MQ=255)
                                          aggatgggCGTGGCGATACGGTCGAGGAACCAACGGGCGTGCGAGATAACCATCACACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGGTTCGATATCCAGGTCGGTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACAGGAAAGTGCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGct                                                                                                                                                                                                      <  2:381907/244‑2 (MQ=255)
                                                 gCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAg                                                                                                                                                                                               <  1:373100/251‑1 (MQ=255)
                                                           aCGGTCTAGGAACCCACGGTCGTGCGCGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCTCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACACACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCACCGTCGGCGCGGCTTg                                                                                                                                                                                     >  2:107935/1‑251 (MQ=255)
                                                            cGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTgg                                                                                                                                                                                    >  1:234097/1‑251 (MQ=255)
                                                                cGAGGAACCCACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCAtc                                                                                                                                                                                >  2:444191/1‑251 (MQ=255)
                                                                 gAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCAtct                                                                                                                                                                               >  1:343219/1‑251 (MQ=255)
                                                                     aaCCAACGGTCGTGCGAGATAACCTTCCCACCGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCGAGGTCGGTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGTCCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCTACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGgt                                                                                                                                                                           <  2:178681/251‑1 (MQ=255)
                                                                         aaCGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTGGGCATCTCGGTGTTg                                                                                                                                                                       <  2:30400/251‑1 (MQ=255)
                                                                         aaCGGTCGGGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCATGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCCCGGTGTTg                                                                                                                                                                       >  2:369934/1‑251 (MQ=255)
                                                                          aCGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTAGGCATCTCGGTGTTGc                                                                                                                                                                      >  2:464129/1‑251 (MQ=255)
                                                                           cgttCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGcc                                                                                                                                                                     <  1:230044/248‑1 (MQ=255)
                                                                           cGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGcc                                                                                                                                                                     <  1:348484/251‑1 (MQ=255)
                                                                           cGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGcc                                                                                                                                                                     <  1:197030/251‑1 (MQ=255)
                                                                           cGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGcc                                                                                                                                                                     <  1:310750/251‑1 (MQ=255)
                                                                                gggCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCt                                                                                                                                                                <  2:404432/249‑1 (MQ=255)
                                                                                 tGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGCTCtt                                                                                                                                                               >  1:287072/1‑251 (MQ=255)
                                                                                       aTAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCAAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGGGTTGCCGGTCTTCAtaat                                                                                                                                                         >  2:337463/1‑248 (MQ=255)
                                                                                       aTAACAATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCAtaat                                                                                                                                                         >  1:395025/1‑248 (MQ=255)
                                                                                           ccATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATcc                                                                                                                                                     <  1:438666/251‑1 (MQ=255)
                                                                                           ccATCGCACAGCACGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCACCGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATCATATcc                                                                                                                                                     >  2:469342/1‑251 (MQ=255)
                                                                                                            aaCTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAActtc                                                                                                                                    >  1:339060/1‑251 (MQ=255)
                                                                                                                    cagGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCaa                                                                                                                             >  1:418323/1‑250 (MQ=255)
                                                                                                                    cagGGCGGTTTCCAGCGCGTGCAGTGTTTCGATATTCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCCCCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATAGCCAGCCCGCCGGAAACTTCTTCCCaaa                                                                                                                            <  1:295549/251‑1 (MQ=255)
                                                                                                                     agGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAAc                                                                                                                           <  1:195174/251‑1 (MQ=255)
                                                                                                                               ccAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCt                                                                                                                  <  1:81027/250‑1 (MQ=255)
                                                                                                                                                tCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGACAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTcccg                                                                                                >  2:343404/1‑248 (MQ=255)
                                                                                                                                                   atagccAGGACGTTGGTTGGTTCGGCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGaa                                                                                             <  2:224243/247‑1 (MQ=255)
                                                                                                                                                         aGGTCGTTGGTTGGGTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATc                                                                                       <  2:418323/251‑1 (MQ=255)
                                                                                                                                                          ggTCGTGGGTGGGTCCGTCGAGCAGCAGCAGGTTTCCGCCAACCTGCATCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCa                                                                                      <  2:472339/251‑1 (MQ=255)
                                                                                                                                                           gTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCa                                                                                      >  1:272911/1‑250 (MQ=255)
                                                                                                                                                             cGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAGAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATGATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAAcc                                                                                   <  1:349583/251‑1 (MQ=255)
                                                                                                                                                                          gtcgAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGCACTGATCAACCGATGCCAgtttc                                                                       >  2:137508/1‑250 (MQ=255)
                                                                                                                                                                            cgAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGCCGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCCTAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTAt                                                                                                             >  2:60749/1‑210 (MQ=255)
                                                                                                                                                                            cgAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTAt                                                                                                             <  1:60749/210‑1 (MQ=255)
                                                                                                                                                                            cgAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAgtttcac                                                                     <  1:229929/250‑1 (MQ=255)
                                                                                                                                                                             gAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGt                                                                                                    <  2:422896/218‑1 (MQ=255)
                                                                                                                                                                             gAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGt                                                                                                    >  1:422896/1‑218 (MQ=255)
                                                                                                                                                                                  gcagcaTGTTGCCGCCAACCTGCAGCAGCTGCGCCAGATGCAGACGACCGCGCGCACTACCGGAAAGTTCACCAACGCTTTTCCCCTGATCAACCCCGTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTGGCCGATCTTCATAATATCCAGCCCGCCGGAAACTGCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAggttcaccgtttca                                                              <  2:347729/251‑1 (MQ=255)
                                                                                                                                                                                      caTGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAgtttcaccgtttcacct                                                           >  1:260529/1‑249 (MQ=255)
                                                                                                                                                                                                 cAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCtt                                                                                                                                                                                                                      >  2:301134/1‑84 (MQ=255)
                                                                                                                                                                                                 cAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCtt                                                                                                                                                                                                                      <  1:301134/84‑1 (MQ=255)
                                                                                                                                                                                                  aaCCTGCCGCAACTTCGCCGGATGCAGAAGAGCGCGCTCACCACCGGACCGTTCACAAACGCGGTTACCCTGATCAACCCCTTTGAAGGTAAAGCGGCAAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATATTCATAATATCCAGCCCTCCGGAAACTTCTTCCCACACGGGTTTACTGTTATCCGTTGAGTCACGTAACTGAGCAACAGATGCCCGTTTCCCCGTTTCACCTAAAGTGGTGGTGc                                              >  2:71717/1‑251 (MQ=255)
                                                                                                                                                                                                    ccagcagcagcaTCGCCAGATACAGACGACCGCACTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGACCAACGTAGGCGCGGCTTGGAATCTTGGTGTTGCCGATCTTCATAATATCCAGACCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCTATTGAGTAACGGAACTGATCAACCGATGCAAGTTTCACCGTTTCACCGAAATGGATGGTGCCg                                            >  2:449982/4‑251 (MQ=255)
                                                                                                                                                                                                        cagcagCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCCCGGAACTGATCAACCGCTGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGt                                        >  2:325891/1‑251 (MQ=255)
                                                                                                                                                                                                        cagcagCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGACCGCCGGAAACTTCTTACCAAACGGTTTTGATGTTATCAATTGAGGCACGGAACTTATCAACCAATGCCAGTTTCACCGTTTCACATAAAGTGATGGTGCCGCTgc                                        >  2:430716/1‑250 (MQ=255)
                                                                                                                                                                                                           cagcTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCg                                     >  1:341211/1‑251 (MQ=255)
                                                                                                                                                                                                                tCGCCAGATGCAGACGACCGCGCTCACAACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGAGCCGCTGTCCGGCTGt                                >  1:417478/1‑251 (MQ=255)
                                                                                                                                                                                                                    cagagtgatacgaaCGCGCTCACCACCGGAGAGTTCACCAACGCCTGTACCCCTATCAACCCTAGTAAAGTTGAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATGATATCCAGCCCGCCGGAAACTTCCTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGACCAACCGACGCCAGTGTCACCGTTTCAACCAAAGTGATGGTGCCGCTGTCCGGCTGTTCCt                            <  2:229571/242‑1 (MQ=255)
                                                                                                                                                                                                                        tGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGAGATGGTGCAGCTGTCCGGCTGTTCCTGAcc                        >  1:416924/1‑251 (MQ=255)
                                                                                                                                                                                                                           aGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAACGTTAAAGCGGcc                                                                                                                                                                                                      >  2:204707/1‑74 (MQ=255)
                                                                                                                                                                                                                           aGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGcc                                                                                                                                                                                                      <  1:204707/74‑1 (MQ=255)
                                                                                                                                                                                                                              cgacCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGAt                  <  1:426072/251‑1 (MQ=255)
                                                                                                                                                                                                                                acCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCgg                                                                                                                                           >  2:190907/1‑128 (MQ=255)
                                                                                                                                                                                                                                acCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCgg                                                                                                                                           <  1:190907/128‑1 (MQ=255)
                                                                                                                                                                                                                                acCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCa                <  2:388400/251‑1 (MQ=255)
                                                                                                                                                                                                                                     gcTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACgg           <  2:272911/251‑1 (MQ=255)
                                                                                                                                                                                                                                          ccaccGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAAc        >  1:105289/1‑249 (MQ=255)
                                                                                                                                                                                                                                              cGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTc  <  2:414594/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                      
TTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAATCCAGAATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAGAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATGATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGACGCCAGTTTCACCGTTTCACCCAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTC  >  NC_000913/4628946‑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: