Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260479 3,529,607 +C intergenic (‑121/‑183) 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,6061.C100.0% 148.0 / NA 38intergenic (‑120/‑184)secA/h16_A3265preprotein translocase subunit SecA/Hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (20/18);  total (20/18)

GTGATCATCGATTGCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTGCCGGATGTCGGCCACGCGTGTGGTCCGCAGATTGCCCGCGTGATTGCCCGTGCGCCCATTGGCCACAGTTTGGGGCGGCCGC  >  AM260479/3529479‑3529730
                                                                                                                                 |                                                                                                                            
gTGATCATCGATTGCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCAcggccgg                                                                                                                          >  1:394623/1‑130 (MQ=255)
   atcatcGATTGCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGGATGTGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAAcg                                                                                                                       <  1:88814/134‑1 (MQ=255)
       tcGATTGCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAACGCGcc                                                                                                                   <  2:394623/134‑1 (MQ=255)
         gATTGCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCt                                                                                                                  <  2:115513/133‑1 (MQ=255)
          aTTGCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAAcg                                                                                                                       >  1:1182441/1‑127 (MQ=255)
          aTTGCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAAcg                                                                                                                       <  2:1182441/127‑1 (MQ=255)
           ttGCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCTGGc                                                                                                               <  2:1175909/134‑1 (MQ=255)
           ttGCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCTGGc                                                                                                               >  2:291171/1‑134 (MQ=255)
               tttCACCTGCGCGGAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCTGGCCCCg                                                                                                           >  2:162226/1‑134 (MQ=255)
                       gcgcGGAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCTGGCCCCGTCCCGc                                                                                                     >  2:239857/1‑132 (MQ=255)
                            gAACCGGGCGGCATCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAggg                                                                                                >  2:198662/1‑132 (MQ=255)
                            gAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAAcg                                                                                                                       >  1:84854/1‑109 (MQ=255)
                            gAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAAcg                                                                                                                       <  2:84854/109‑1 (MQ=255)
                                      gcACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTg                                                                                    >  2:1122316/1‑134 (MQ=255)
                                                     gggCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAAcg                                                                                                                       >  2:378057/1‑84 (MQ=255)
                                                     gggCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAAcg                                                                                                                       <  1:378057/84‑1 (MQ=255)
                                                      ggCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTGCCGGATGTCGGCCAcg                                                                    <  1:239857/134‑1 (MQ=255)
                                                       gctgctCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTGCCGGATGTCGGCCAcgc                                                                   <  2:762962/134‑1 (MQ=255)
                                                          gctCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAAcg                                                                                                                       <  2:186679/79‑1 (MQ=255)
                                                          gctCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAAcg                                                                                                                       >  1:186679/1‑79 (MQ=255)
                                                               ggcggcCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTGCCGGATGTCGGCCACGCGTGTGGTc                                                           >  1:1135229/1‑134 (MQ=255)
                                                                       ggCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTGCCGGATGTCGGCCACGCGTGTGGTCCGCAGAtt                                                   <  2:1135229/134‑1 (MQ=255)
                                                                                 ggATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTGCCGGATGTCGGCCACGCGTGTGGTCCGCAGATTGCCcgcg                                            >  2:851005/1‑131 (MQ=255)
                                                                                 ggATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTGCCGGATGTCGGCCACGCGTGTGGTCCGCAGATTGCCcgcg                                            <  1:851005/131‑1 (MQ=255)
                                                                                      gATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTGCCGGATGTCGGCCACGCGTGTGGTCCGCAGATTGCCCGCGTGATTGcc                                    >  1:288592/1‑134 (MQ=255)
                                                                                                aaaCGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTGCCGGATGt                                                                            >  1:110615/1‑84 (MQ=255)
                                                                                                aaaCGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTGCCGGATGt                                                                            <  2:110615/84‑1 (MQ=255)
                                                                                                aaaCGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTGCCGGATGTCGGCCAcgc                                                                   <  1:316195/93‑1 (MQ=255)
                                                                                                aaaCGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTGCCGGATGTCGGCCAcgc                                                                   >  2:316195/1‑93 (MQ=255)
                                                                                                aaaCGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTGCCGGATGTCGGCCACGCgtgt                                                               >  1:406138/1‑97 (MQ=255)
                                                                                                aaaCGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTGCCGGATGTCGGCCACGCgtgt                                                               <  2:406138/97‑1 (MQ=255)
                                                                                                                       ggggTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGc                                                                                       >  1:124547/1‑50 (MQ=255)
                                                                                                                       ggggTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGc                                                                                       <  2:124547/50‑1 (MQ=255)
                                                                                                                       ggggTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTGCCGGATGTCGGCCACGCGTGTGGTCCGCAGATTGCCCGCGTGATTGCCCGTGCGCCCAtt                        >  2:220575/1‑113 (MQ=255)
                                                                                                                       ggggTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTGCCGGATGTCGGCCACGCGTGTGGTCCGCAGATTGCCCGCGTGATTGCCCGTGCGCCCAtt                        <  1:220575/113‑1 (MQ=255)
                                                                                                                       ggggTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTGCCGGATGTCGGCCACGCGTGTGGTCCGCAGATTGCCCGCGTGATTGCCCGTGCGCCCATTGGCCACAGTTTGGGGCGgccg   <  1:733549/134‑1 (MQ=255)
                                                                                                                        gggTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTGCCGGATGTCGGCCACGCGTGTGGTCCGCAGATTGCCCGCGTGATTGCCCGTGCGCCCATTGGCCACAGTTTGGGGCGgccgc  >  2:158792/1‑134 (MQ=255)
                                                                                                                        gggTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTGCCGGATGTCGGCCACGCGTGTGGTCCGCAGATTGCCCGCGTGATTGCCCGTGCGCCCATTGGCCACAGTTTGGGGCGgccgc  >  2:868472/1‑134 (MQ=255)
                                                                                                                        gggTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTGCCGGATGTCGGCCACGCGTGTGGTCCGCAGATTGCCCGCGTGATTGCCCGTGCGCCCATTGGCCACAGTTTGGGGCGgcccc  >  2:1175820/1‑132 (MQ=255)
                                                                                                                                 |                                                                                                                            
GTGATCATCGATTGCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTGCCGGATGTCGGCCACGCGTGTGGTCCGCAGATTGCCCGCGTGATTGCCCGTGCGCCCATTGGCCACAGTTTGGGGCGGCCGC  >  AM260479/3529479‑3529730

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: