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% 89.5 / NA 24intergenic (‑53/‑222)h16_A1136/grpEconserved hypothetical protein/molecular chaperone GrpE
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base T (12/12);  total (12/12)

CCGGAAAGGCAGATAAAGACGGTCATGACCGCCAATGCGGCCACAAGGCAGAGTGTCACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGCCTGGCCTTTGTAACCGTGTGTTTCCGCGTCACGCCAGCGCCACTCT  >  AM260479/1239924‑1240172
                                                                                                                                  |                                                                                                                       
ccGGAAAGGCAGATAAAGACGGTCATGACCGCCAATGCGGCCACAAGGCAGAGTGTCACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCAtttt                                                                                                                     <  1:430138/134‑1 (MQ=255)
    aaaGGCAGATAAAGACGGTCATGACCGCCAATGCGGCCACAAGGCAGAGTGTCACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCAtttttgt                                                                                                                  >  2:335329/1‑130 (MQ=255)
          aGATAAAGACGGTCATGACCGCCAATGCGGCCACAAGGCAGAGTGTCACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCAttttt                                                                                                            <  1:335329/134‑1 (MQ=255)
               aaGACGGTCATGACCGCCAATGCGGCCACAAGGCAGAGTGTCACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTttatta                                                                                                        <  2:190530/133‑1 (MQ=255)
               aaGACGGTCATGACCGCCAATGCGGCCACAAGGCAGAGTGTCACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTttatta                                                                                                        <  1:468136/133‑1 (MQ=255)
               aaGACGGTCATGACCGCCAATGCGGCCACAAGGCAGAGTGTCACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTttatta                                                                                                        <  1:862888/133‑1 (MQ=255)
                                           caAGGCAGAGTGTCACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATgcagc                                                                            >  1:130579/1‑133 (MQ=255)
                                                         aCGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTttat                                                                                                          >  1:223996/1‑89 (MQ=255)
                                                         aCGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTttat                                                                                                          <  2:223996/89‑1 (MQ=255)
                                                                     ttATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCAttt                                                                                                              >  2:83264/1‑73 (MQ=255)
                                                                     ttATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCAttt                                                                                                              <  1:83264/73‑1 (MQ=255)
                                                                                  aCCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAAt                                                                    >  1:519131/1‑102 (MQ=255)
                                                                                  aCCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAAt                                                                    <  2:519131/102‑1 (MQ=255)
                                                                                  aCCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCata                                                   >  1:673134/1‑119 (MQ=255)
                                                                                  aCCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCata                                                   <  2:673134/119‑1 (MQ=255)
                                                                                      tCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGCCTGGCCTTTGTAACCg                                >  1:302355/1‑134 (MQ=255)
                                                                                       ccGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGCCTGGCCTTTGTAACCgt                               >  2:663970/1‑134 (MQ=255)
                                                                                                  acaAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCAttt                                                                                                              <  1:824852/44‑1 (MQ=38)
                                                                                                  acaAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCAttt                                                                                                              >  2:824852/1‑44 (MQ=38)
                                                                                                              cATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATgcagc                                                                            <  2:116039/66‑1 (MQ=255)
                                                                                                              cATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATgcagc                                                                            >  1:116039/1‑66 (MQ=255)
                                                                                                                    cacaACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCagcaaag                                                                        <  2:217046/64‑1 (MQ=255)
                                                                                                                    cacaACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCagcaaag                                                                        >  1:217046/1‑64 (MQ=255)
                                                                                                                    cacaACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGCCTGGCCTTTGTAACCGTGTGTTTCCGCGTCACGCCAGCGCCActct  >  1:812441/1‑134 (MQ=255)
                                                                                                                         cccccGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCg                                                              <  1:466874/69‑1 (MQ=255)
                                                                                                                         cccccGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCg                                                              >  2:466874/1‑69 (MQ=255)
                                                                                                                                  |                                                                                                                       
CCGGAAAGGCAGATAAAGACGGTCATGACCGCCAATGCGGCCACAAGGCAGAGTGTCACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGCCTGGCCTTTGTAACCGTGTGTTTCCGCGTCACGCCAGCGCCACTCT  >  AM260479/1239924‑1240172

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