Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260479 3,529,526 +G intergenic (‑40/‑264) secA ← / → h16_A3265 preprotein translocase subunit SecA/Hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604793,529,5241.G100.0% 53.2 / NA 15intergenic (‑38/‑266)secA/h16_A3265preprotein translocase subunit SecA/Hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (8/7);  total (8/7)

GGCGATATTGTTTGATCAGCCGCTCATTGCGGCTGCCGAAGACTTTCTTGAGAAGGCCGGTGATCATCGATTGCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTGCCG  >  AM260479/3529420‑3529651
                                                                                                         |                                                                                                                                
ggcgATATTGTTTGATCAGCCGCTCATTGCGGCTGCCGAAGACTTTCTTGAGAAGGCCGGTGATCATCGATTGCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGTCGGGCg                                                                                                      <  1:588115/134‑1 (MQ=255)
   gATATTGTTTGATCAGCCGCTCATTGCGGCTGCCGAAGACTTTCTTGAGAAGGCCGGTGATCATCGATTGCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAg                                                                                                    >  1:813585/1‑133 (MQ=255)
        tGTTTGATCAGCCGCTCATTGCGGCTGCCGAAGACTTTCTTGAGAAGGCCGGTGATCATCGATTGCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATg                                                                                               >  1:689732/1‑133 (MQ=255)
          tttGATCAGCCGCTCATTGCGGCTGCCGAAGACTTTCTTGAGAAGGCCGGTGATCATCGATTGCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGAt                                                                                            <  1:350069/134‑1 (MQ=255)
                         aTTGCGGCTGCCGAAGACTTTCTTGAGAAGGCCGGTGATCATCGATTGCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCaaa                                                                               <  2:689732/132‑1 (MQ=255)
                          ttGCGGCTGCCGAAGACTTTCTTGAGAAGGCCGGTGATCATCGATTGCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACg                                                                            >  2:540095/1‑134 (MQ=255)
                                                                      ttGCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCTGGc                                <  1:84682/134‑1 (MQ=255)
                                                                           ttCACCTGCGCGGAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCTGGCCCCGt                           >  2:689328/1‑134 (MQ=255)
                                                                                       gAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAg                                                                                                    <  1:657092/49‑1 (MQ=38)
                                                                                       gAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAg                                                                                                    >  2:657092/1‑49 (MQ=38)
                                                                                       gAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAAc                                                                             >  1:1056328/1‑72 (MQ=255)
                                                                                       gAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAAc                                                                             <  2:1056328/72‑1 (MQ=255)
                                                                                                 gcACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTg     >  2:528249/1‑134 (MQ=255)
                                                                                                 gcACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTg     >  2:934071/1‑134 (MQ=255)
                                                                                                    ccccAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTGCCg  <  1:406412/134‑1 (MQ=255)
                                                                                                         |                                                                                                                                
GGCGATATTGTTTGATCAGCCGCTCATTGCGGCTGCCGAAGACTTTCTTGAGAAGGCCGGTGATCATCGATTGCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTGCCG  >  AM260479/3529420‑3529651

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