Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260479 1,240,057 +T intergenic (‑57/‑218) h16_A1136 ← / → grpE conserved hypothetical protein/molecular chaperone GrpE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604791,240,0531.T100.0% 134.5 / NA 36intergenic (‑53/‑222)h16_A1136/grpEconserved hypothetical protein/molecular chaperone GrpE
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base T (18/18);  total (18/18)

GGCAGATAAAGACGGTCATGACCGCCAATGCGGCCACAAGGCAGAGTGTCACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGCCTGGCCTTTGTAACCGTGTGTTTCCGCGTCACGCCAGCGCCACTCT  >  AM260479/1239931‑1240172
                                                                                                                           |                                                                                                                       
ggCAGATAAAGACGGTCATGACCGCCAATGCGGCCACAAGGCAGAGTGTCACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCAtt                                                                                                               <  1:1132049/134‑1 (MQ=255)
        aaGACGGTCATGACCGCCAATGCGGCCACAAGGCAGAGTGTCACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTttatta                                                                                                        <  1:576453/133‑1 (MQ=255)
        aaGACGGTCATGACCGCCAATGCGGCCACAAGGCAGAGTGTCACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTAc                                                                                                       >  1:523205/1‑134 (MQ=255)
                   gACCGCCAATGCGGCCACAAGGCAGAGTGTCACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCt                                                                                              >  1:564862/1‑132 (MQ=255)
                           aTGCGGCCACAAGGCAGAGTGTCACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGc                                                                                    >  1:355431/1‑134 (MQ=255)
                             gCGGCCACAAGGCAGAGTGTCACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGAtt                                                                                      >  1:1057852/1‑130 (MQ=255)
                             gCGGCCACAAGGCAGAGTGTCACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGAtt                                                                                      <  2:1057852/130‑1 (MQ=255)
                                gCCACAAGGCAGAGTGTCACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATgc                                                                               <  2:355431/134‑1 (MQ=255)
                                                  aCGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGt                                                             <  2:564862/134‑1 (MQ=255)
                                                     aCTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGa                                                          >  2:880343/1‑134 (MQ=255)
                                                       tGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCAtttt                                                                                                                     <  1:1012600/72‑1 (MQ=255)
                                                       tGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCAtttt                                                                                                                     >  2:1012600/1‑72 (MQ=255)
                                                            gCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTg                                                                                         >  1:540753/1‑96 (MQ=255)
                                                            gCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTg                                                                                         <  2:540753/96‑1 (MQ=255)
                                                            gCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCg                                                              >  1:459563/1‑123 (MQ=255)
                                                            gCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCg                                                              <  2:459563/123‑1 (MQ=255)
                                                              ttATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGAt                                                         >  1:587971/1‑126 (MQ=255)
                                                              ttATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGAt                                                         <  2:587971/126‑1 (MQ=255)
                                                                aTGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATg                                                 <  2:523205/132‑1 (MQ=255)
                                                                 tGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGcc                                               <  1:880343/133‑1 (MQ=255)
                                                                  gCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGCCTg                                             >  1:474492/1‑134 (MQ=255)
                                                                           aCCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGcc                                               >  2:253080/1‑123 (MQ=255)
                                                                           aCCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGcc                                               <  1:253080/123‑1 (MQ=255)
                                                                                    aaGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGCCTGGCCTTTGTAACCgtgtgt                           >  1:631854/1‑134 (MQ=255)
                                                                                           acaAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATgc                                                                               <  2:663903/75‑1 (MQ=255)
                                                                                           acaAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATgc                                                                               >  1:663903/1‑75 (MQ=255)
                                                                                             aaaGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCTATTGCATATGCCTGGCCTTTGTAACCGTGTGTTTCCGCGt                   >  2:125384/1‑133 (MQ=255)
                                                                                                     aCCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGCCTGGCCTTTGTAACCGTGTGTTTCCGCGTCACGCCAgc          <  2:631854/134‑1 (MQ=255)
                                                                                                      ccATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGCCTGGCCTTTGTAACCGTGTGTTTCCGCGTCACGCCAgcg         <  1:429945/134‑1 (MQ=255)
                                                                                                       cATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATgc                                                                               >  2:575171/1‑63 (MQ=255)
                                                                                                       cATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATgc                                                                               <  1:575171/63‑1 (MQ=255)
                                                                                                       cATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATgc                                                                               >  2:348548/1‑63 (MQ=255)
                                                                                                       cATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATgc                                                                               <  1:348548/63‑1 (MQ=255)
                                                                                                             cacaACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCg                                                              <  1:781937/74‑1 (MQ=255)
                                                                                                             cacaACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCg                                                              <  1:187073/74‑1 (MQ=255)
                                                                                                             cacaACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCg                                                              >  2:187073/1‑74 (MQ=255)
                                                                                                             cacaACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCg                                                              >  2:781937/1‑74 (MQ=255)
                                                                                                             cacaACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCATCGATTGCATATGCCTGGCCTTTGTAACCGTGTGTTTCCGCGTCACGCCAGCGCCActct  >  2:459659/1‑134 (MQ=255)
                                                                                                                           |                                                                                                                       
GGCAGATAAAGACGGTCATGACCGCCAATGCGGCCACAAGGCAGAGTGTCACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGCCTGGCCTTTGTAACCGTGTGTTTCCGCGTCACGCCAGCGCCACTCT  >  AM260479/1239931‑1240172

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