Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,755,080:1 +C 10.5% coding (1406/3147 nt) czcA‑II → cation efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,755,0801.C10.5% 138.2 / 9.4 38coding (1406/3147 nt)czcA‑IIcation efflux system protein
Reads supporting (aligned to +/- strand):  ref base . (19/15);  new base C (2/2);  total (21/17)
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.00e-01

CCCATGCCCAGCAGCGCCACGGTCGCATGCTGACCCGCGGCGAACGGTTCCATGAGGTGTTTGCCGCTGCCCGCGAAGCACGCCGGCCGCTGATCTATGGGCAGTTGATCATCATGGTGGTGTACCTGCCGATCTTTGCCCTGACCGGGGTCGAGGGCAAGATGTTCCACCCTATGGCCTTTACCGTGGTCATGGCCCTGTTGGGCGCGATGATCCTTTCGGTCACGTTCGTTCCGGCAGCCATTGCCCTGTTCGTCACCGGCAAGGTCAAGGAAGAGGAAGGCCTGGT  >  NC_002947/2754939‑2755227
                                                                                                                                              |                                                                                                                                                    
cccATGCCCAGCAGCGCCACGGTCGCATGCTGACCCGCGGCGAACGGTTCCATGAGGTGTTTGCCGCTGCCCGCGAAGCACGCCGGCCGCTGATCTATGGGCAGTTGATCATCATGGTGGTGTACCTGCCGATCTTTGCCCTGACCggg                                                                                                                                               <  5:47940/149‑1 (MQ=255)
                 cACGGTCGCATGCTGACCCGCGGCGAACGGTTCCATGAGGTGTTTGCCGCTGCCCGCGAAGCACGCCGGCCGCTGATCTATGGGCAGTTGATCATCATGGTGGTGTACCTGCCGATCTTTGCCCTCGACCCCGGTCAGGGCaaagatcg                                                                                                                              <  1:142005/149‑7 (MQ=25)
                           tGCTGACCCGCGGCGAACGGTTCCATGAGGTGTTTGCCGCTGCCCGCGAAGCACGCCGGCCGCTGATCTATGGGCAGTTGATCATCATGGTGGTGTACCTGCCGATCTTTGCCCTGACCGGGGTCGAGGGCAAGATGTTCCACCCTATg                                                                                                                   >  3:134825/1‑149 (MQ=255)
                            gCTGACCCGCGGCGAACGGTTCCATGAGGTGTTTGCCGCTGCCCGCGAAGCACGCCGGCCGCTGATCTATGGGCAGTTGATCATCATGGTGGTGTACCTGCCGATCTTTGCCCTGACCGGGGTCGAGGGCAAg                                                                                                                                  >  5:222553/1‑133 (MQ=255)
                            gCTGACCCGCGGCGAACGGTTCCATGAGGTGTTTGCCGCTGCCCGCGAAGCACGCCGGCCGCTGATCTATGGGCAGTTGATCATCATGGTGGTGTACCTGCCGATCTTTGCCCTGACCGGGGTCGAGGGCAAg                                                                                                                                  <  6:222553/133‑1 (MQ=255)
                              tGACCCGCGGCGAACGGTTCCATGAGGTGTTTGCCGCTGCCCGCGAAGCACGCCGGCCGCTGATCTATGGGCAGTTGATCATCATGGTGGTGTACCTGCCGATCTTTGCCCTGACCGGGGTCGAGGGCAAGATGTTCCACCCTATGGcc                                                                                                                >  5:172677/1‑149 (MQ=255)
                              tGACCCGCGGCGAACGGTTCCATGAGGTGTTTGCCGCTGCCCGCGAAGCACGCCGGCCGCTGATCTATGGGCAGTTGATCATCATGGTGGTGTACCTGCCGATCTTTGCCCTGACCGGGGTCGAGGGCAAGATGTTCCACCCTATGGcc                                                                                                                >  7:132790/1‑149 (MQ=255)
                                   cgcgGCGAACGGTTCCATGAGGTGTTTGCCGCTGCCCGCGAAGCACGCCGGCCGCTGATCTATGGGCAGTTGATCATCATGGTGGTGTACCTGCCGATCTTTGCCCTGACCGGGGTCGAGGGCAAGATGTTCCACCCTATGGCCTTTAc                                                                                                           >  1:240259/1‑149 (MQ=255)
                                       gcgAACGGTTCCATGAGGTGTTTGCCGCTGCCCGCGAAGCACGCCGGCCGCTGATCTATGGGCAGTTGATCATCATGGTGGTGTACCTGCCGATCTTTGCCCTGACCGGGGTCGAGGGCAAGATGTTCCACCCTATGGCCTTTACCGTg                                                                                                       <  6:165688/149‑1 (MQ=255)
                                                     gAGGTGTTTGCCGCTGCCCGCGAAGCACGCCGGCCGCTGATCTATGGGCAGTTGATCATCATGGTGGTGTACCTGCCGATCTTTGCCCTGACCGGGGTCGAGGGCAAGATGTTCCACCCTATGGCCGTTACCGTGGTCATGGCCCTGtt                                                                                         >  5:241546/1‑149 (MQ=255)
                                                     gAGGTGTTTGCCGCTGCCCGCGAAGCACGCCGGCCGCTGATCTATGGGCAGTTGATCATCATGGTGGTGTACCTGCCGATCTTTGCCCTGACCGGGGTCGAGGGCAAGATGTTCCACCCTATGGACTTTACCGTGGTCATGGCCCTGtt                                                                                         >  8:247208/1‑149 (MQ=255)
                                                      aGGTGTTTGCCGCTGCCCGCGAAGCACGCCGGCCGCTGATCTATGGGCAGTTGATCATCATGGTGGTGTACCTGCCGATCTTTGCCCTGACCGGGGTCGAGGGCAAGAt                                                                                                                                >  6:202100/1‑109 (MQ=255)
                                                      aGGTGTTTGCCGCTGCCCGCGAAGCACGCCGGCCGCTGATCTATGGGCAGTTGATCATCATGGTGGTGTACCTGCCGATCTTTGCCCTGACCGGGGTCGAGGGCAAGAt                                                                                                                                <  5:202100/109‑1 (MQ=255)
                                                         tgtTTGCCGCTGCCCGCGAAGCACGCCGGCCGCTGATCTATGGGCAGTTGATCATCATGGTGGTGTACCTGCCGATCTTTGCCCTGACCGGGGTCGAGGGCAAGATGTTCCACCCTATGGCCTTTACCGTGGTCATGGCCCTGTTGGgc                                                                                     <  6:172677/149‑1 (MQ=255)
                                                         tgtTTGCCGCTGCCCGCGAAGCACGCCGGCCGCTGATCTATGGGCAGTTGATCATCATGGTGGTGTACCTGCCGATCTTTGCCCTGACCGGGGTCGAGGGCAAGATGTTCCACCCTATGGCCTTTACCGTGGTCATGGCCCTGTTGGgc                                                                                     <  8:132790/149‑1 (MQ=255)
                                                                cgatgaacgcgAAGCACGCCATCCGCTGATATAACTTCAGTAGATCATCATGGAGGTCTACCTGCAGATATTTGAACTGCCCGGGGTCTAGGGAAAGATGTTACAACATATGGCCTTTACAGTCATCATGGCCCTGTTGGGCGCgagga                                                                              >  4:241971/8‑149 (MQ=255)
                                                                 gCTGCCCGCGAAGCACGCCGGCCGCTGATCTATGGGCAGTTGATCATCATGGTGGTGTACCTGCCGATCTTTGCCCTGACCGGGGTCGAGGGCAAGATGTTCCACCCTATGGCCTTTACCGTGGTCATGGCCCTGTTGGGCGCgatgat                                                                             <  4:134825/149‑1 (MQ=255)
                                                                      ccGCGAAGCACGCCGGCCGCTGATCTATGGGCAGTTGATCATCATGGTGGTGTACCTGCCGATCTTTGCCCTCGACCCCGGTCAGGGCa                                                                                                                                     <  4:92982/89‑1 (MQ=18)
                                                                      ccGCGAAGCACGCCGGCCGCTGATCTATGGGCAGTTGATCATCATGGTGGTGTACCTGCCGATCTTTGCCCTCGACCCCGGTCAGGGCa                                                                                                                                     >  3:92982/1‑89 (MQ=18)
                                                                      ccGCGAAGCACGCCGGCCGCTGATCTATGGGCAGTTGATCATCATGGTGGTGTACCTGCCGATCTTTGCCCTGACCGGGGTCGAGGGCAACGATGTTCCACCCTATGGCCTTTACCGTGGTCATgg                                                                                                <  1:176779/126‑1 (MQ=255)
                                                                      ccGCGAAGCACGCCGGCCGCTGATCTATGGGCAGTTGATCATCATGGTGGTGTACCTGCCGATCTTTGCCCTGACCGGGGTCGAGGGCAAGATGTTCCACCCTATGGCCTTTACCGTGGTCATgg                                                                                                >  2:176779/1‑125 (MQ=255)
                                                                      ccGCGAAGCACGCCGGCCGCTGATCTATGGGCAGTTGATCATCATGGTGGTGTACCTGCCGATCTTTGCCCTGACCGGGGTCGAGGGCAAGATGTTCCACCCTATGGCCTTTACCGTGGTCATGGCCCTGTTGGGCGCGATGATCCttt                                                                        >  5:189971/1‑149 (MQ=255)
                                                                       cgcgAAGCACGCCGGCCGCTGATCTATGGGCAGTTGATCATCATGGTGGTGTACCTGCCGATCTTTGCCCTGACCGGGGTCGAGGGCAAGATGTTCCACCCTATGGCCTTTACCGTGGTCATGGCCCTGTTGGGCGCGATGATCCTTTc                                                                       >  4:212836/1‑149 (MQ=255)
                                                                          gAAGCACGCCGGCCGCTGATCTATGGGCAGTTGATCATCATGGTGGTGTACCTGCCGATCTTTGCCCTGAc                                                                                                                                                   >  4:69361/1‑71 (MQ=255)
                                                                          gAAGCACGCCGGCCGCTGATCTATGGGCAGTTGATCATCATGGTGGTGTACCTGCCGATCTTTGCCCTGAc                                                                                                                                                   <  3:69361/71‑1 (MQ=255)
                                                                           aaGCACGCCGGCCGCTGATCTATGGGCAGTTGATCATCATGGTGGTGTACCTGCCGATCTTTGCCCTGACCGGGGTCGAGGGCAAGATGTTCCACCCTATGGCCTTTACCGTGGTCATGGCCCTGTTGGGCGCGATGATCCTTTCGGTc                                                                   <  6:259289/149‑1 (MQ=255)
                                                                                                   ggCAGTTGATCATCATGGTGGTGTACCTGCCGATCTTTGCCCTGACCGGGGTCGAGGGCAAGATGTTCCACCCTATGGCCTTTACCGTGGTCATGGCCCTGTTGGGCGCGATGATCCTTTCGGTCACGTTCGTTCCGGCAGCCATTGcc                                           <  2:240259/149‑1 (MQ=255)
                                                                                                   ggCAGTTGATCATCATGGTGGTGTACCTGCCGATCTTTGCCCTGACCGGGGTCGAGGGCAAGATGTTCCACCCTATGGCCTTTACCGTGGTCATGGCACTGATGGGCGCTATGATACTTTCGGTAACGTTCGTTCAGGCAGACATTGcc                                           >  4:199679/1‑149 (MQ=21)
                                                                                                          gATCATCATGGTGGTGTACCTGCCGATCTTTGCCCTGACCGGGGTCGAGGGCAAGATGTTCCACCCTATGGCCTTTACCGTGGTCATGGCCCTGTTGGGCGCGATGATCCTTTCGGTCACGTTCGTTCCGGCAGCCATTGCCCTGTTc                                     <  7:247208/148‑1 (MQ=255)
                                                                                                            tcatcaTGGTGGTGTACCTGCCGATCTTTGCCCTGACCGGGGTCGAGGGCAAGATGTTCCACCCTATGGCCTTTACCGTGGTCATGGCCCTGTTGGGCGCGATGATCCTTTCGGTCACGTTCGTTCCGGCAGCCATTGCCCTGTTCGt                                   >  4:229349/1‑148 (MQ=255)
                                                                                                                     gggtgTCCCTGCCGATCTGTGCCCTGACCGGGGTCGAGGGCAAGATGTTCCACCCTATGGCCTTTACCGTGGTCATGGCCCTGTTGGGCGCGATGATCCTTTCGGTCACGTTCGTTCCGGCAGCCATTGCCCTGTTCGTCACCGGCAAg                         <  6:189971/148‑1 (MQ=255)
                                                                                                                            ccTGCCGATCTTTGCCCTGACCGGGGTCGAGGGCAAGATGTTCCACCCTATGGCCTTTACCGTGGTCATGGCCCTGTTGGGCGCGATGATCCTTTCGGTCACGTTCGTTCCGGCAGCCAtt                                              >  4:3448/1‑121 (MQ=255)
                                                                                                                            ccTGCCGATCTTTGCCCTGACCGGGGTCGAGGGCAAGATGTTCCACCCTATGGCCTTTACCGTGGTCATGGCCCTGTTGGGCGCGATGATCCTTTCGGTCACGTTCGTTCCGGCAGCCAtt                                              <  3:3448/121‑1 (MQ=255)
                                                                                                                            ccTGCCGATCTTTGCCCTGACCGGGGTCGAGGGCAAGATGTTCCACCCTATGGCCTTTACCGTGGTCATGGCCCTGTTGGGCGCGATGATCCTTTCGGTCACGTTCGTTCCGGCAGCCATTGCCCTGTTCGTCACCGGCAAGGTCAag                   <  3:212836/148‑1 (MQ=255)
                                                                                                                                 cGATCTTTGCCCTGACCGGGGTCGAGGGCAAGATGTTCCACCCTATGGCCTTTACCGTGGTCATGGCCCTGGTGGGCGCGATGATCCTTTCGGTCACGTTCGTTCCGGCAGCCATTGACCTGTTCGTCACCGGCAAGGTCAaggaagag             >  7:59914/1‑149 (MQ=255)
                                                                                                                                  gATCTTTGCCCTGACCGGGGTCGAGGGCAAGATGTTCCACCCTATGGCCTTTACCGTGGTCATGGCCCTGTTGGGCGCGATGATCCTTTCGGTCACGTTCGTTCCGGCAGCCATTGCCCTGTTCGTCACCGGCAAGGTCAaggaagag             >  5:109508/1‑148 (MQ=255)
                                                                                                                                   aTCTTTGCCCTGACCGGGGTCGAGGGCAAGATGTTCCACCCTATGGCCTTTACCGTGGTCATGGCCCTGTTGGGCGCGATGATCCTTTCGGTCACGTTCGTTCCGGCAGCCATTGCCCTGTTCGTCACCGGCAAGGTCAaggaagagg            >  3:180352/1‑148 (MQ=255)
                                                                                                                                     tcttGCCCTCGACCCCGGTCAGGGCAAGATGTTCCACCCTATGGCCTTTACCGTGGTCATGGCCCTGTTGGGCGCGATGATCCTTTCGGTCACGTTCGTTCCGGCAGCCATTGCCCTGTTCGTCACCGGCACGGTCAaggaagaggaag         >  1:122073/3‑149 (MQ=255)
                                                                                                                                             tGACCGGGGTCGAGGGCAAGATGTTCCACCCTATGGCCTTTACCGTGGTCATGGCCCTGTTGGGCGCGATGATCCTTTCGGTCACGTTCGTTCCGGCAGCCATTGCCCTGTTCGTCACCGGCAAGGTCAAGGAAGAGGAAGGCCTGGt  >  1:100176/1‑148 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
CCCATGCCCAGCAGCGCCACGGTCGCATGCTGACCCGCGGCGAACGGTTCCATGAGGTGTTTGCCGCTGCCCGCGAAGCACGCCGGCCGCTGATCTATGGGCAGTTGATCATCATGGTGGTGTACCTGCCGATCTTTGCCCTGACCGGGGTCGAGGGCAAGATGTTCCACCCTATGGCCTTTACCGTGGTCATGGCCCTGTTGGGCGCGATGATCCTTTCGGTCACGTTCGTTCCGGCAGCCATTGCCCTGTTCGTCACCGGCAAGGTCAAGGAAGAGGAAGGCCTGGT  >  NC_002947/2754939‑2755227

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: