Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,755,086 G→C 8.3% G471A (GGG→GCG) ‡ 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,0860GC 8.3% 152.5 / 4.1 48G471A (GGG→GCG) ‡czcA‑IIcation efflux system protein
Reads supporting (aligned to +/- strand):  ref base G (21/23);  new base C (2/2);  total (23/25)
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.73e-01

CCCATGCCCAGCAGCGCCACGGTCGCATGCTGACCCGCGGCGAACGGTTCCATGAGGTGTTTGCCGCTGCCCGCGAAGCACGCCGGCCGCTGATCTATGGGCAGTTGATCATCATGGTGGTGTACCTGCCGATCTTTGCCCTGACCGGGGTCGAGGGCAAGATGTTCCACCCTATGGCCTTTACCGTGGTCATGGCCCTGTTGGGCGCGATGATCCTTTCGGTCACGTTCGTTCCGGCAGCCATTGCCCTGTTCGTCACCGGCAAGGTCAAGGAAGAGGAAGGCCTGGTGATGC  >  NC_002947/2754939‑2755232
                                                                                                                                                    |                                                                                                                                                   
cccATGCCCAGCAGCGCCACGGTCGCATGCTGACCCGCGGCGAACGGTTCCATGAGGTGTTTGCCGCTGCCCGCGAAGCACGCCGGCCGCTGATCTATGGGCAGTTGATCATCATGGTGGTGTACCTGCCGATCTTTGCCCTGACCggg                                                                                                                                                    >  1:79607/1‑149 (MQ=255)
 ccATGCCCAGCAGCGCCACGGTCGCATGCTGACCCGCGGCGAACGGTTCCATGAGGTGTTTGCCGCTGCCCGCGAAGCACGCCGGCCGCTGATCTATGGGCAGTTGATCATCATGGTGGTGTACCTGCCGATCTTTGCCCTGACCgggg                                                                                                                                                   >  4:207849/1‑149 (MQ=255)
     gCCCAGCAGCGCCACGGTCGCATGCTGACCCGCGGCGAACGGTTCCATGAGGTGTTTGCCGCTGCCCGCGAAGCACGCCGGCCGCTGATCTATGGGCAGTTGATCATCATGGTGGTGTACCTGCCGATCTTTGCCCTGACCGGGGTc                                                                                                                                                 >  3:128730/1‑147 (MQ=255)
            agcGCCACGGTCGCATGCTGACCCGCGGCGAACGGTTCCATGAGGTGTTTGCCGCTGCCCGCGAAGCACGCCGGCCGCTGATCTATGGGCAGTTGATCATCATGGTGGTGTACCTGCCGATCTTTGCCCTCGACCCCGGTCAGGGCaaa                                                                                                                                       <  5:266920/149‑2 (MQ=25)
               gcCACGGTCGCATGCTGACCCGCGGCGAACGGTTCCATGAGGTGTTTGCCGCTGCCCGCGAAGCACGCCGGCCGCTGATCTATGGGCAGTTGATCATCATGGTGGTGTACCTGCCGATCTTTGCCCTCGACCCCGGTCAGGGCaaagat                                                                                                                                    <  5:50481/149‑5 (MQ=25)
               gcCACGGTCGCATGCTGACCCGCGGCGAACGGTTCCATGAGGTGTTTGCCGCTGCCCGCGAAGCACGCCGGCCGCTGATCTATGGGCAGTTGATCATCATGGTGGTGTACCTGCCGATCTTTGCCCTCGACCCCGGTCAGGGCa                                                                                                                                         >  6:50481/1‑144 (MQ=25)
                    ggTCGCATGCTGACCCGCGGCGAACGGTTCCATGAGGTGTTTGCCGCTGCCCGCGAAGCACGCCGGCCGCTGATCTATGGGCAGTTGATCATCATGGTGGTGTACCTGCCGATCTTTGCCCTGACCGGGGTTGAGGGCAAGATGTTCCa                                                                                                                               <  8:156250/149‑1 (MQ=255)
                    ggTCGCATGCTGACCCGCGGCGAACGGTTCCATGAGGTGTTTGCCGCTGCCCGCGAAGCACGCCGGCCGCTGATCTATGGGCAGTTGATCATCATGGTGGTGTACCTGCCGATCTTTGCCCTGACCGGGGTCGAGGGCAAGATGTTCCa                                                                                                                               >  2:34723/1‑149 (MQ=255)
                     gTCGCATGCTGACCCGCGGCGAACGGTTCCATGAGGTGTTTGCCGCTGCCCGCGAAGCACGCCGGCCGCTGATCTATGGGCAGTTGATCATCATGGTGGTGTACCTGCCGATCTTTGCCCTGACCGGGGTCGAGGGCAAGATGTTCCAc                                                                                                                              <  2:269938/149‑1 (MQ=255)
                        gCATGCTGACCCGCGGCGAACGGTTCCATGAGGTGTTTGCCGCTGCCCGCGAAGCACGCCGGCCGCTGATCTATGGGCAGTTGATCATCATGGTGGTGTACCTGCCGATCTTTGCCCTGACCGGGGTCGAGGGCAAGATGTTCCAccc                                                                                                                            >  4:83028/1‑148 (MQ=255)
                         cATGCTGACCCGCGGCGAACGGTTCCATGAGGTGTTTGCCGCTGCCCGCGAAGCACGCCGGCCGCTGATCTATGGGCAGTTGATCATCATGGTGGTGTACCTGCCGATCTTTGCCCTGACCGGGGTCGAGGGCAAGATGTTCCACCCTa                                                                                                                          >  3:183738/1‑149 (MQ=255)
                              tGACCCGCGGCGAACGGTTCCATGAGGTGTTTGCCGCTGCCCGCGAAGCACGCCGGCCGCTGATCTATGGGCAGTTGATCATCATGGTGGTGTACCTGCCGATCTTTGCCCTGACCGGGGTCGAGGGCAAGATGTTCCACCCTATGGcc                                                                                                                     >  5:12000/1‑149 (MQ=255)
                                aCCCGCGGCGAACGGTTCCATGAGGTGTTTGCCGCTGCCCGCGAAGCACGCCGGCCGCTGATCTATGGGCAGTTGATCATCATGGTGGTGTACCTGCCGATCTTTGCCCTGACCGGGGTCGAGGGCAAGATGTTCCACCCTATGGCCtt                                                                                                                   >  4:55599/1‑149 (MQ=255)
                                  ccGCGGCGAACGGTTCCATGAGGTGTTTGCCGCTGCCCGCGAAGCACGCCGGCCGCTGATCTATGGGCAGTTGATCATCATGGTGGTGTACCTGCCGATCTTTGCCCTGACCGGGGTCGAGGGCAAGATGTTCCACCCTATGGCCTTTa                                                                                                                 >  8:182117/1‑149 (MQ=255)
                                     cggcgAACGGTTCCATGAGGTGTTTGCCGCTGCCCGCGAAGCACGCCGGCCGCTGATCTATGGGCAGTTGATCATCATGGTGGTGTACCTGCCGATCTTTGCCCTGACCGGGGTCGAGGGCAAGATGTTCCACCCTATGGCCTTTACCg                                                                                                              <  3:55599/149‑1 (MQ=255)
                                         gAACGGTTCCATGAGGTGTTTGCCGCTGCCCGCGAAGCACGCCGGCCGCTGATCTATGGGCAGTTGATCATCATGGTGGTGTACCTGCCGATCTTTGCCCTGACCGGGGTCGAggg                                                                                                                                            <  5:49877/116‑1 (MQ=255)
                                         gAACGGTTCCATGAGGTGTTTGCCGCTGCCCGCGAAGCACGCCGGCCGCTGATCTATGGGCAGTTGATCATCATGGTGGTGTACCTGCCGATCTTTGCCCTGACCGGGGTCGAggg                                                                                                                                            >  6:49877/1‑116 (MQ=255)
                                          aaCGGTTCCATGAGGTGTTTGCCGCTGCCCGCGAAGCACGCCGGCCGCTGATCTATGGGCAGTTGATCATCATGGTGGTGTACCTGCCGATCTTTGCCCTGACCGGGGTCGAGGGCAAGATGTTCCACCCTATGGCCTTTACCGTGGTc                                                                                                         <  3:207849/149‑1 (MQ=255)
                                            cGGTTCCATGAGGTGTTTGCCGCTGCCCGCGAAGCACGCCGGCCGCTGATCTATGGGCAGTTGATCATCATGGTGGTGTACCTGCCGATCTTTGCCCTGACCGGGGTCGAGGGCAAGATGTTCCACCCTATGGCCTTTACCGTGGTCAt                                                                                                       >  4:92981/1‑149 (MQ=255)
                                              gTTCCATGAGGTGTTTGCCGCTGCCCGCGAAGCACGCCGGCCGCTGATCTATGGGCAGTTGATCATCATGGTGGTGTACCTGCCGATCTTTGCCCTGACCGGGGTCGAGGGCAAGATGTTCCACCCTATGGCCTTTACCGTGGTCATgg                                                                                                     <  7:182117/149‑1 (MQ=255)
                                               ttCCATGAGGTGTTTGCCGCTGCCCGCGAAGCACGCCGGCCGCTGATCTATGGGCAGTTGATCATCATGGTGGTGTACCTGCCGATCTTTGCCCTGACCGGGGTCGAGGGCAAGATg                                                                                                                                    <  7:50144/117‑1 (MQ=255)
                                               ttCCATGAGGTGTTTGCCGCTGCCCGCGAAGCACGCCGGCCGCTGATCTATGGGCAGTTGATCATCATGGTGGTGTACCTGCCGATCTTTGCCCTGACCGGGGTCGAGGGCAAGATg                                                                                                                                    >  8:50144/1‑117 (MQ=255)
                                               ttCCATGAGGTGTTTGCCGCTGCCCGCGAAGCACGCCGGCCGCTGATCTATGGGCAGTTGATCATCATGGTGGTGTACCTGCCGATCTTTGCCCTGACCGGGGTCGAGGGCAAGATGTTCCACCCTATGGCCTTTACCGTGGTCATGGc                                                                                                    >  8:76899/1‑149 (MQ=255)
                                                   aTGAGGTGTTTGCCGCTGCCCGCGAAGCACGCCGGCCGCTGATCTATGGGCAGTTGATCATCATGGTGGTGTACCTGCCGATCTTTGCCCTGACCGGGGTCGAGGGCAAGATGTTCCACCCTATGGCCTTTACCGTGGTCATGGCCCTg                                                                                                <  5:261288/149‑1 (MQ=255)
                                                      aGGTGTTTGCCGCTGCCCGCGAAGCACGCCGGCCGCTGATCTATGGGCAGTTGATCATCATGGTGGTGTACCTGCCGATCTTTGCCCTGACCGGGGTCGAGGGCAAGATGTTCCACCCTATGGCCTTTACCGTGGTCATGGCCCTGTTg                                                                                             <  2:30887/149‑1 (MQ=255)
                                                           tttGCCGCTGCCCGCGAAGCACGCCGGCCGCTGATCTATGGGCAGTTGATCATCATGGTGGTGTACCTGCCGATCTTTGCCCTGACCGGGGTCGAGGGCAAGATGTTCCACCCTATGGCCTTTACCGTGGTCATGGCCCTGTTGGgcgc                                                                                        >  3:301847/1‑149 (MQ=255)
                                                           tttGCCGCTGCCCGCGAAGCACGCCGGCCGCTGATCTATGGGCAGTTGATCATCATGGTGGTGTACCTGCCGATCTTTGCCCTGACCGGGGTCGAGGGCAAGATGTTCCACCCTATGGCCTTTACCGTGGTCATGGCCCTGTTGGgcgc                                                                                        <  7:76899/149‑1 (MQ=255)
                                                               ccGCTGCCCGCGAAGCACGCCGGCCGCTGATCTATGGGCAGTTGATCATCATGGTGGTGTACCTGCCGATCTTTGCCCTGACCGGGGTCGAGGGCAAGATGTTCCACCCTATGGCCTTTACCGTGGTCATGGCCCTGTTGGGCGCgatg                                                                                    <  5:168379/149‑1 (MQ=255)
                                                                cGCTGCCCGCGAAGCACGCCGGCCGCTGATCTATGGGCAGTTGATCATCATGGTGGTGTACCTGCCGATCTTTGCCCTGACCGGGGTCGAGGGCAAGATGTTCCACCCTATGGCCTTTACCGTGGTCATGGCCCTGTTGGGCGCgatga                                                                                   <  3:83028/149‑1 (MQ=255)
                                                                  cTGCCCGCGAAGCACGCCGGCCGCTGATCTATGGGCAGTTGATCATCATGGTGGTGTACCTGCCGATCTTTGCCCTGACCGGGGTCGAGGGCAAGATGTTCCACCCTATGGCCTTTACCGTGGTCATGGCCCTGTTGGGCGCGATGATc                                                                                 <  4:183738/149‑1 (MQ=255)
                                                                  cTGCCCGCGAAGCACGCCGGCCGCTGATCTATGGGCAGTTGATCATCATGGTGGTGTACCTGCCGATCTTTGCCCTGACCGGGGTCGAGGGCAAGATGTTCCACCCTATGGCCTTTACCGTGGTCATGGCCCTGTTGGGCGCGATGATc                                                                                 <  3:92981/149‑1 (MQ=255)
                                                                        gcgAAGCACGCCGGCCGCTGATCTATGGGCAGTTGATCATCATGGTGGTGTACCTGCCGATCTTTGCCCTGACCGGGGTCGAGGGCAAGATGTTCCACCCTATGGCCTTTACCGTGGTCATGGCCCTGTTGGGCGCGATGATCCTTTCg                                                                           <  2:79607/149‑1 (MQ=255)
                                                                         cgAAGCACGCCGGCCGCTGATCTATGGGCAGTTGATCATCATGGTGGTGTACCTGCCGATCTTTGCCCTGACCGGGGTCGAGGGCAAGATGTTCCACCCTATGGCCTTTACCGTGGTCATGGCCCTGTTGGGCGCGATGATCCTTTCgg                                                                          >  2:40187/1‑149 (MQ=255)
                                                                         cgAAGCACGCCGACCGCTGATCTATGGGCAGTTGATCATCATGGTGGTGTACCTGCCGATCTTTGCCCTGACCGGGGTCGAGGGCAAGATGTTCCACCCTATGGCCTTTACCGTGGTCATGGCCCTGTTGGGCGCGATGATCCTTTCgg                                                                          >  1:175547/1‑149 (MQ=255)
                                                                                                 agggCAGTTGATCATCATGGTCGACTACCAGCCGAACTTTGACCTGACCGGGGTCAAGGCCAAGATGTTCCACACTATCGCCTTTACCGTGGTCATGGCCCAGTTGGGCTCTAAGATCCTTACGGTCACTTTCATACCGGCAGCCAtt                                                   >  8:122201/2‑148 (MQ=255)
                                                                                                      aGTTGATCATCATGGTGGTGTACCTGCCGATCTTTGCCCTGACCGGGGTCGAGGGCAAGATGTTCCACCCTATGGCCTTTACCGTGGTCATGGCCCTGTTGGGCGCGATGATCCTTTCGGTCACGTTCGTTCCGGCAGCCATTGCCCTg                                             >  1:190769/1‑149 (MQ=255)
                                                                                                           atcatcATGGTGGTGTACCTGCCGATCTTTGCCCTGACCGGGGTCGAGGGCAAGATGTTCCACCCTATGGCCTTTACCGTGGTCATGGCCCTGTTGGGCGCGATGATCCTTTCGGTCACgttcgttc                                                              <  8:178626/127‑1 (MQ=255)
                                                                                                           atcatcATGGTGGTGTACCTGCCGATCTTTGCCCTGACCGGGGTCGAGGGCAAGATGTTCCACCCTATGGCCTTTACCGTGGTCATGGCCCTGTTGGGCGCGATGATCCTTTCGGTCACgttcgttc                                                              >  7:178626/1‑127 (MQ=255)
                                                                                                              atcatGGTGGTGTACCTGCCGATCTTTGCCCTGACCGGGGTCGAGGGCAAGATGTTCCACCCTATGGCCTTTACCGTGGTCATGGCCCTGTTGGGCGCGATGATCCTTTCGGTCACGTTCGTTCCGGCAGCCATTGCCCTGTTCGTCAc                                     <  2:190769/149‑1 (MQ=255)
                                                                                                               tcatGGTGGTGTACCTGCCGATCTTTGCCCTGACCGGGGTCGAGGGCAAGATGTTCCACCCTATGGCCTTTACCGTGGTCATGGCCCTGTTGGGCGCGATGATCCTTTCGGTCACGTTCGTTCCGGCAGCCATTGCCCTGTTCGTCAcc                                    >  1:180873/1‑149 (MQ=255)
                                                                                                                    gtggttTACCTGCCTATCTTTGCCCTGACCGGGGTCGAGGGCAAGATGTTCCACCCTATGGCCTTTACCGTGGTCATGGCCCTGTTGGGCGCGATGATCCTTTCGGTCACGTTCGTTCCGGCAGCCATTGCCCTGTTCGTCACCGGCaa                               <  7:122201/149‑1 (MQ=255)
                                                                                                                      ggtgTACCTGCCGATCTTTGCCCTGACCGGGGTCGAGGGCAAGATGTTCCACCCTATGGCCTTTACCGTGGTCATGGCCCTGTTGGGCGCGATGATCCTTTCGGTCACGTTCGTTCCGGCAGCCATTGCCCTGTTCGTCACCGGCAAgg                             <  2:180873/149‑1 (MQ=255)
                                                                                                                          tACCTGCCGATCTTTGCCCTGACCGGGGTCGAGGGCAAGATGTTCCACCCTATGGCCTTTACCGTGGTCATGGCCCTGTTGGGCGCGATGATCCTTTCGGTCACGTTCGTTCCGGCAGCCATTGCCCTGTTCGTCACCGGCAAGGTCaa                         <  1:34723/149‑1 (MQ=255)
                                                                                                                             cTGCCGATCTTTGCCCTGACCGGGGTCGAGGGCAAGATGTTCCACCCTATGGCCTTTACCGTGGTCATGGCCCTGTTGGGCGCGATGATCCTTTCGGTCACGTTCGTTCCGGCAGCCATTGCCCTGTTCGTCACCGGCAAGGTCAagg                       >  6:260126/1‑148 (MQ=255)
                                                                                                                                 cGATCTTTGCCCTGACCGGGGTCGAGGGCAAGATGTTCCACCCTATGGCCTTTACCGTGGTCATGGCCCTGTTGGGCGCGATGATCCTTTCGGTCACGTTCGTTCCGGCAGCCATTGcc                                                >  2:272133/1‑119 (MQ=255)
                                                                                                                                    atcttGCCCTCGACCCCGGTCAGGGCAAAGATGTTCCACCCTATGGCCTTTACCGTGGTCATGGCCCTGTTGGGCGCGATGATCCTTTCGGTCACGTTCGTTCCGGCAGCCATTGCCCTGTTCGTCACCGGCAAGGTCAaggaagagg                 >  4:260392/4‑148 (MQ=255)
                                                                                                                                     cTTTGCCCTGACCGGGGTCGAGGGCAAGATGTTCCACCCTATGGCCTTTACCGTGGTCATGGCCCTGTTGGGCGCGATGATCCTTTCGGTCACGTTCGTTCCGGCAGCCATTGccc                                               >  6:297675/1‑116 (MQ=255)
                                                                                                                                     cTTTGCCCTGACCGGGGTCGAGGGCAAGATGTTCCACCCTATGGCCTTTACCGTGGTCATGGCCCTGTTGGGCGCGATGATCCTTTCGGTCACGTTCGTTCCGGCAGCCATTGccc                                               <  5:297675/116‑1 (MQ=255)
                                                                                                                                        tGCCCTGACCGGGGTCGAGGGCAAGATGTTCCACCCTATGGCCTTTACCGTGGTCATGGCCCTGTTGGGCGCGATGATCCTTTCGGTCACGTTCGTTCCGGCAGCCATTGCCCTGTTCGTCACCGGCAAGGTCAAGGAAGAGGAAgg             <  6:12000/147‑1 (MQ=255)
                                                                                                                                                  cGGGGTCGAGGGCAAGATGTTCCACCCTATGGCCTTTACCGTGGTCATGGCCCTGTTGGGCGCGATGATCCTTTCGGTCACGTTCGTTCCGGCAGCCATTGCCCTGTTCGTCACCGGCAAGGTCAAGGAAGAGGAAGGCCTGGTGATGc  <  4:303651/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CCCATGCCCAGCAGCGCCACGGTCGCATGCTGACCCGCGGCGAACGGTTCCATGAGGTGTTTGCCGCTGCCCGCGAAGCACGCCGGCCGCTGATCTATGGGCAGTTGATCATCATGGTGGTGTACCTGCCGATCTTTGCCCTGACCGGGGTCGAGGGCAAGATGTTCCACCCTATGGCCTTTACCGTGGTCATGGCCCTGTTGGGCGCGATGATCCTTTCGGTCACGTTCGTTCCGGCAGCCATTGCCCTGTTCGTCACCGGCAAGGTCAAGGAAGAGGAAGGCCTGGTGATGC  >  NC_002947/2754939‑2755232

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: