Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,189,651 G→A 12.1% R55C (CGC→TGC)  tatC‑I ← Sec‑independent protein translocase protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,189,6510GA12.1% 84.0 / 3.3 33R55C (CGC→TGC) tatC‑ISec‑independent protein translocase protein
Reads supporting (aligned to +/- strand):  ref base G (15/14);  new base A (2/2);  total (17/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01

CCAGAAAACCCCAGGCCTGATGCAGCAGCAGCGGCATGGCCAGGAACAGGGCGCACATCGCGGTCAGCTTGAACGGTGCCAGAAAGGGCGAGGTGACGCTGGTGGCGATCATGCTCGCGCCTTCCGGCAAAAAGCGCCGTAGCGGCTCGGATATGAGCGTGTACAGTGTCTGGGCAAAGGGGAACAGGCCGGCGAACACCAATGCCACCAGCGCCAGGCAACGCACCAGGCGTTTGCGCAGGTCGCGCAGGTGTTCGGTCAGGGGCATGCTCGCG  >  NC_002947/1189516‑1189790
                                                                                                                                       |                                                                                                                                           
ccAGAAAACCCCAGGCCTGATGCAGCAGCAGCGGCATGGCCAGGAACAGGGCGCACATCGCGGTCAGCTTGAACGGTGCCAGAAAGGGCGAGGTGACGCTGGTGGCGATCATGCTCGCGCCTTCCGGCAAAAAGCGCCGTAGCGGCTCg                                                                                                                                <  3:54076/149‑1 (MQ=255)
   gAAAACCCCAGGCCTGATGCAGCAGCAGCGGCATGGCCAGGAACAGGGCGCACATCGCGGTCAGCTTGAACGGTGCCAGAAAGGGCGAGGTGACGCTGGTGGCGATCATGCTCGCGCCTTCCGGCAAAAAGCGCCGTAGCGGCTCGGat                                                                                                                             <  1:62394/149‑1 (MQ=255)
           cAGGCCTGATGCAGCAGCAGCGGCATGGCCAGGAACAGGGCGCACATCGCGGTCAGCTTGAACGGTGCCAGAAAGGGCGAGGTGACGCTGGTGGCGATCATGCTCGCGCCTTCCAGCAAAAAGCGCCGTAGCGGCTCGGATATGAGCgt                                                                                                                     >  2:64056/1‑149 (MQ=255)
                 tGATGCAGCAGCAGCGGCATGGCCAGGAACNGGGCGCACATCGCGGTCAGCTTGAACGGTGCCAGAAAGGGCGAGGTGACGCTGGTGGCGATCATGCTCGCGCCTTCCGGCAAAAAGCGCCGTAGCGGCTCGGATATGAGCGTGTACAg                                                                                                               >  1:23764/1‑149 (MQ=255)
                   aTGCAGCAGCAGCGGCATGGCCAGGAACAGGGCGCACATCGCGGTCAGCTTGAACGGTGCCAGAAAGGGCGAGGTGACGCTGGTGGCGATCATGCTCGCGCCTTCCGGCAAAAAGCGCCGTAGCGGCTCGGATATGAGCGTGTACAgt                                                                                                              <  1:101377/148‑1 (MQ=255)
                       agcagcagcGGCATGGCCAGGAACAGGGCGCACATCGCGGTCAGCTTGAACGGTGCCAGAAAGGGCGAGGTGACGCTGGTGGCGATCATGCTCGCGCCTTCCGGCAAAAAGCGCCGTAGCGGCTCGGATATGAGCGTGTACAGTGTCTg                                                                                                         >  4:171836/1‑149 (MQ=255)
                       agcagcagcGGCATGGCCAGGAACAGGGCGCACATCGCGGTCAGCTTGAACGGTGCCAGAAAGGGCGAGGTGACGCTGGTGGCGATCATGCTCGCGCCTTCCGCCAAAAAGCGCCGTAGCGGCTCGCATATGAGCGTGTACAGTGTCTg                                                                                                         >  4:8988/1‑149 (MQ=255)
                       agcagcagcGGCATGGCCAGGAACAGGGCGCACATCGCGGTCAGCTTGAACGGTGCCAGAAAGGGCGAGGTGACGCTGGTGGCGATCATGCTCGCGCCTTCCAGCAAAAAGCGCCGTAGCGGCTCGGATATGAGCGTGTACAGTGTCTg                                                                                                         <  1:64056/149‑1 (MQ=255)
                         cagcagcGGCATGGCCAGGAACAGGGCGCACATCGCGGTCAGCTTGAACGGTGCCAGAAAGGGCGAGGTGACGCTGGTGGCGATCATGCTCGCGCCTTCCGGCAAAAAGCGCCGTAGCGGCTCGGATATGAGCGTGTACAGTGTCTggg                                                                                                       <  1:165152/149‑1 (MQ=255)
                                ggCATGGCCAGGAACAGGGCGCACATCGCGGTCAGCTTGAACGGTGCCAGAAAGGGCGAGGTGACGCGGGTGGCGATCATGCTCGCGCCTTCCGGCAAAAAGCGCCGTAGCGGCTCGGATATGAGCGTGTACAGTGTCTGGGCAAAggg                                                                                                <  3:114896/149‑1 (MQ=255)
                                        cAGGAACAGGGCGCACATCGCGGTCAGCTTGAACGGTGCCAGAAAGGGCGAGGTGACGCTGGTGGCGATCATGCTCGCGCCTTCCGGCAAAAAGCGCCGTAGCGGCTCGGATATGAGCGTGTACAGTGTCTGGGCAAAGGGGAACAGGc                                                                                        >  8:103710/1‑149 (MQ=255)
                                            aaCAGGGCGCACATCGCGGTCAGCTTGAACGGTGCCAGAAAGGGCGAGGTGACGCTGGTGGCGATCATGCTCGCGCCTTCCGGCAAAAAGCGCCGTAGCGGCTCGGATATGAGCGTGAACAGTGTCTGGGCAAAGGGGAACAGGCCgg                                                                                     <  1:22286/148‑1 (MQ=255)
                                                 ggCGCACATCGCGGTCAGCTTGAACGGTGCCAGAAAGGGCGAGGTGACGCTGGTGGCGATCATGCTCGCGCCTTCCGGCAAAAAGCACCGTAGCGGCTCGGATATGAGCGTGTACAGTGTCTggg                                                                                                       <  4:101679/125‑1 (MQ=255)
                                                 ggCGCACATCGCGGTCAGCTTGAACGGTGCCAGAAAGGGCGAGGTGACGCTGGTGGCGATCATGCTCGCGCCTTCCGGCAAAAAGCACCGTAGCGGCTCGGATATGAGCGTGTACAGTGTCTggg                                                                                                       <  4:97695/125‑1 (MQ=255)
                                                 ggCGCACATCGCGGTCAGCTTGAACGGTGCCAGAAAGGGCGAGGTGACGCTGGTGGCGATCATGCTCGCGCCTTCCGGCAAAAAGCACCGTAGCGGCTCGGATATGAGCGTGTACAGTGTCTggg                                                                                                       >  3:97695/1‑125 (MQ=255)
                                                 ggCGCACATCGCGGTCAGCTTGAACGGTGCCAGAAAGGGCGAGGTGACGCTGGTGGCGATCATGCTCGCGCCTTCCGGCAAAAAGCACCGTAGCGGCTCGGATATGAGCGTGTACAGTGTCTggg                                                                                                       >  3:101679/1‑125 (MQ=255)
                                                  gcgcACATCGCGGTCAGCTTGAACGGTGCCAGAAAGGGCGAGGTGACGCTGGTGGCGATCATGCTCGCGCCTTCCGGCAAAAAGCGCCGTAGCGGCTCGGATATGAGCGTGTACAGTGTCTGGGCAAAGGGGAACAGGCCGGCGAacac                                                                              >  4:208650/1‑149 (MQ=255)
                                                    gcACATCGCGGTCAGCTTGAACGGTGCCAGAAAGGGCGAGGTGACGCTGGTGGCGATCATGCTCGCGCCTTCCGGCAAAAAGCGCCGTAGCGGCTCGGATATGAGCGTGTACAGTGTCTGGGCAAAGGGGAACAGGCCGGCGAACAcc                                                                             >  2:23705/1‑148 (MQ=255)
                                                        aTCGCGGTCAGCTTGAACGGTGCCAGAAAGGGCGAGGTGACGCTGGTGGCGATCATGCTCGCGCCTTCCGGCAAAAAGCGCCGTAGCGGCTCGGATATGAGCGTGTACAGTGTCTGGGCAAAGGGGAACAGGCCGGCGAACACCAATGc                                                                        >  8:210508/1‑149 (MQ=255)
                                                           gcgGTCAGCTTGAACGGTGCCAGAAAGGGCGAGGTGACGCTGGTGGCGATCATGCTCGCGCCTTCCGGCAAAAAGCGCCGTAGCGGCTCGGATATGAGCGTGTACAGTGTCTGGGCAAAGGGGAACAGGCCGGCGAACACCAATGcca                                                                      <  7:103710/148‑1 (MQ=255)
                                                                   cTTGAACGGTGCCAGAAAGGGCGAGGTGACGCTGGTGGCGATCATGCTCGCGCCTTCCGGCAAAAAGCGCCGTAGCGGCTCGGATATGAGCGTGTACAGTGTCTGGGCAAAGGGGAACAGGCCGGCGAACACCAATGCCACCAgcgc                                                               >  5:161895/1‑147 (MQ=255)
                                                                                aGAAAGGGCGAGGTGACGCTGGTGGCGATCATGCTCGCGCCTTCCGGCAAAAAGCGCCGTAGCGGCTCGGATATGAGCGTGTACAGTGTCTGGGCAAAGGGGAACAGGCCGGCGAACACCAATGCCACCAGCGCCAGGCAACGCACCAg                                                <  7:210508/149‑1 (MQ=255)
                                                                                 gAAAGGGCGAGGTGACGCTGGTGGCGATCATGCTCGCGCCTTCCGGCAAAAAGCGCCGTAGCGGCTCGGATATGAGCGTGTACAGTGTCTGGGCAAAGGGGAACAGGCCGGCGAACACCAATGCCACCAGCGCCAGGCAACGCACCAgg                                               >  1:49621/1‑149 (MQ=255)
                                                                                                       ggCGATCATGCTCGCGCCTTCCGGCAAAAAGCGCCGTAGCGGCTCGGATATGAGCGTGTACAGTGTCTGGGCAAAGGGGAACAGGCCGGCGAACACCAATGCCACCAGCGCCAGGCAACGCACCAGGCGTTTGCGCAGGTCGCGCAGgt                         <  3:8988/149‑1 (MQ=255)
                                                                                                           aTCATGCTCGCGCCTTCCGGCAAAAAGCGCCGTAGCGGCTCGGATATGAGCGTGTACAGTGTCTGGGCAAAGGGGAACAGGCCGGCGAACACCAATGCCACCAGCGCCAGGCAACGCACCAGGCGTTTGCGCAGGTCGCGCAGGTGTTc                     <  1:23705/149‑1 (MQ=255)
                                                                                                              aTGCTCGCGCCTTCCGGCAAAAAGCGCCGTAGCGGCTCGGATATGAGCGTGTACAGTGTCTGGGCAAAGGGGAACAGGCCGGCGAACACCAATGCCACCAGCGCCAGGCAACGCACCAGGCGTTTGCGCAGGTCGCGCAGGTGTTCGGt                  <  3:171836/149‑1 (MQ=255)
                                                                                                                  tCGCGCCTTCCGGCAAAAAGCGCCGTAGCGGCTCGGATATGAGCGTGTACAGTGTCTGGGCAAAGGGGAACAGGCCGGCGAaca                                                                               >  8:53879/1‑84 (MQ=255)
                                                                                                                  tCGCGCCTTCCGGCAAAAAGCGCCGTAGCGGCTCGGATATGAGCGTGTACAGTGTCTGGGCAAAGGGGAACAGGCCGGCGAaca                                                                               <  7:53879/84‑1 (MQ=255)
                                                                                                                         ttCCGGCAAAAAGCGCCGTAGCGGCTCGGATATGAGCGTGTACAGTGTCTGGGCAAAGGGGAACAGGCCGGCGAACACCAATGCCACCAGCGCCAGGCAACGCACCAGGCGTTTGCGCAGGTCGCGCAGGTGTTCGGTCAGGGGCATGc       >  4:188602/1‑149 (MQ=255)
                                                                                                                             ggCAAAAAGCGCCGTAGCGGCTCGGATATGAGCGTGTACAGTGTCTGGGCAAAGGGGAACAGGCCGGCGACCACCAATGCCACCAGCGCCAGGCAACGCACCAGGCGTTTGCGCAGGTCGCGCAGGTGTTCGGTCAGGGGCATGCTcgc   >  7:91967/1‑149 (MQ=255)
                                                                                                                             ggCAAAAAGCGCCGTAGCGGCTCGGATATGAGCGTGTACAGTGTCTGGGCAAAGGGGAACAGGCCGGCGAACACCAATGCCACCAGCGCCAGGCAACGCACCAGGCGTTTGCGCAGGTCGCGCAGGTGTTCGGTCAGGGGCATGCTcgc   <  6:161895/149‑1 (MQ=255)
                                                                                                                              gCAAAAAGCGCCGTAGCGGCTCGGATATGAGCGTGTACAGTGTCTGGGCAAAGGGGAACAGGCCGGCGAACACCAATGCCACCAGCGCCAGGCAACGCACCAGGCGTTTGCGCAGGTCGCGCAGGTGTTCGGTCAGGGGCATGCTcgcg  >  6:124320/1‑149 (MQ=255)
                                                                                                                              gCAAAAAGCGCCGTAGCGGCTCGGATATGAGCGTGTACAGTGTCTGGGCAAAGGGGAACAGGCCGGCGAACACCAATGCCACCAGCGCCAGGCAACGCACCAGGCGTTTGCGCAGGTCGCGCAGGTGTTCGGTCAGGGGCATGCTcgcg  >  4:112680/1‑149 (MQ=255)
                                                                                                                                       |                                                                                                                                           
CCAGAAAACCCCAGGCCTGATGCAGCAGCAGCGGCATGGCCAGGAACAGGGCGCACATCGCGGTCAGCTTGAACGGTGCCAGAAAGGGCGAGGTGACGCTGGTGGCGATCATGCTCGCGCCTTCCGGCAAAAAGCGCCGTAGCGGCTCGGATATGAGCGTGTACAGTGTCTGGGCAAAGGGGAACAGGCCGGCGAACACCAATGCCACCAGCGCCAGGCAACGCACCAGGCGTTTGCGCAGGTCGCGCAGGTGTTCGGTCAGGGGCATGCTCGCG  >  NC_002947/1189516‑1189790

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