Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,240,425 A→G K56E (AAA→GAA) 
F549F (TTT→TTC
ymgM →
cvrA ←
protein YmgM
putative K(+):H(+) antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,240,4250AG95.0% 61.4 / NA 20K56E (AAA→GAA) 
F549F (TTT→TTC
ymgM
cvrA
protein YmgM
putative K(+):H(+) antiporter
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (12/7);  minor base T (0/1);  total (12/8)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAACCGGCGATAGATGCTTCATTAAATTTAAGATTCAGCTTCCTCTTCAGCCACCCGTACGCCAATCTTCAACACTTCATTATCTTCTTTCTCGGCCACCGTCCAGATCATCCCGGCAAACTCTACCTGGTCACCGACAACCGGTGCCGCGCCTAACAACTGCTGGACAATTTCACCCAGCGTTTGCTGCTTATCACGATATTCTCGCCCGTCTTCCAGGCCATATATCAGCGCAACATCAGCA  >  NC_000913/1240308‑1240551
                                                                                                                     |                                                                                                                              
cAACCGGCGATAGATGCTTCATTAAATTTAAGATTCAGCTTCCTCTTCAGCCACCCGTACGCCAATCTTCAACACTTCATTATCTTCTTTCTCGGCCACCGTCCAGATCATCCCGGCGAACTCTACCTGGTCa                                                                                                                 <  2:268074/133‑1 (MQ=255)
        gATAGATGCTTCATTCAATTTACGATTCAGCTTCCTCTTCAGCCACCCGTACGCCAATCTTCAACACTTCATTATCTTCTTTCTCGGCCACCGTCCAGATCATCCCGGCGAACTCTACCTGGTCACCGACAAc                                                                                                         <  1:176068/133‑1 (MQ=255)
               gCTTCATTAAATTTAAGATTCAGCTTCCTCTTCAGCCACCCGTACGCCAATCTTCAACACTTCATTATCTTCTTTCTCGGCCACCGTCCAGATCATCCCGGCGAACTCTACCTGGTCACCGACAACCGGTGcc                                                                                                  >  1:296060/1‑133 (MQ=255)
                   cATTAAATTTAAGATTCAGCTTCCTCTTCAGCCACCCGTACGCCAATCTTCAACACTTCATTATCTTCTTTCTCGGCCACCGTCCAGATCATCCCGGCGAACTCTACCTGGTCACCGACAACCGGTGCCgcgc                                                                                              >  2:445489/1‑133 (MQ=255)
                   cATTAAATTTAAGATTCAGCTCCTTCTTCAGCCCCCCGTCCGCCAATCTTCAACTCTTCATTATCTTCCTTCTCGGCCGCCGTACAGATCTTCCCCGCGGACTCTTCCTGGTCGCCGACAACCGGTGCCgcgc                                                                                              <  2:478774/133‑1 (MQ=255)
                          tttAAGATTCAGCTTCCTCTTCAGCCACCCGTACGCCAATCTTCAACACTTCATTATCTTCTTTCTCGGCCACCGTCCAGATCATCCCGGCGAACTCTACCTGGTCACCGACAACCGGTGCCGCGCCTaacaa                                                                                       <  2:68064/133‑1 (MQ=255)
                                aTTCAGCTTCTTCTTGAGCCACCCATACTCCACTCTTCAACTCTTCATTATCTTCTTTCTCGGCCACCGTCCAGATCATCCCGGCGAACTCTACCTGGTCACCGACAACCGGTGCCGCGCCTAACAActgctg                                                                                 <  2:131659/133‑1 (MQ=255)
                                         ccTCTTCAGCCACCCGTACGCCAATCTTCAACACTTCATTATCTTCTTTCTCGGCCACCGTCCAGATCATCCCGGCGAACTCTACCTGGTCACCGACAACCGATGCCGCGACTAACACCTGCTGGACAATTTc                                                                        >  1:180111/1‑133 (MQ=255)
                                               cAGCCACCCGTACGCCAATCTTCAACACTTCATTATCTTCTTTCTCGGCCACCGTCCAGATCATCCCGGCGAACTCTACCTGGTCACCGACAACCGGTGCCGCGCCTAACAACTGCTGGACAATTTCACCCAg                                                                  >  1:418280/1‑133 (MQ=255)
                                                     cccGTACGCCAATCTTCAACACTTCATTATCTTCTTTCTCGGCCACCGTCCAGATCATCCCGGCGAACTCTACCTGGTCACCGACAACCGGTGCCGCGCCTAACAACTGCTGGACAATTTCACCCAGCGTTtg                                                            >  1:48389/1‑133 (MQ=255)
                                                     cccGTACGCCAATCTGCAACACTTTATTGTCTGCTTTCTCGTCCACCGTCCAGATCATCCCGGCGAACTCTACCTGGTCACCGACAACCGGTGCCGCGCCTAACAACTGCTGGACAATTTCACCCAGCGTTtg                                                            <  2:169349/133‑1 (MQ=255)
                                                                aTCTTCAACACTTCATTATCTTCTTTCTCGGCCACCGTCCAGATCATCCCGGCGAACTCTACCTGGTCACCGACAACCGGTGCCGCGCCTAACAACTGCTGGACAATTTCACCCAGCGTTTGCTGCTTATCAc                                                 >  2:84382/1‑133 (MQ=255)
                                                                aTCTTCAACACTTCATTATCTTCTTTCTCGGCCACCGTCCAGATCATCCCGGCGAACTCTACCTGGTCACCGACAACCGGTGCCGCGCCTAACAACTGCTGGACAAATGCACCCAGCTTTTGCCGCTTATCAc                                                 >  2:237011/1‑133 (MQ=255)
                                                                         acTTCATTATCTTCTTTCTCGGCCACCGTCCAGATCATCCCGGCGAACTCTACCTGGTCACCGACAACCGGTGCCGCGCCTAACAACTGCTGGACAATTTCACCCAGCGTTTGCTGCTTATCACGATATtctc                                        >  1:484414/1‑133 (MQ=255)
                                                                             cATTATCTTCTTTCTCGGCCACCGTCCAGATCATCCCGGCGAACTCTACCTGGTCACCGACAACCGGTGCCGCGCCTAACAACTGCTGGACAATTTCACCCAGCGTTTGCTGCTTATCACGATATTCTCGccc                                    >  2:55703/1‑133 (MQ=255)
                                                                                     tctttctCGGCCACCTTCCAGATCGTCCGGGCTAACTCTACCTGGTCACCGACAACCGGTGCCGCGCCTAACAACTGCTGGACAATTTCACCCAGCGTTTGCTGCTTATCACGATATTCTCGCCCGTCTTCCa                            <  2:180111/133‑1 (MQ=255)
                                                                                         tctcGGCCACCGTCCAGATCATCCCGGCGAACTCTACCTGGTCACCGACAACCGGTGCCGCGCCTAACAACTGCTGGACAATTTCACCCAGCGTTTGCTGCTTATCACGATATTCTCGCCCGTCTTCCAGGcc                        >  1:512071/1‑133 (MQ=255)
                                                                                         tctcGGCCACCGTCCAGATCATCCCGGCGAACTCTACCTGGTCACCGACAACCGGTGCCGCGCCTAACAACTGCTGGACAATTTCACCCAGCGTTTGCTGCTCATCACGAGATTCTCGCCCGTCGTCCAGGcc                        >  1:223049/1‑133 (MQ=255)
                                                                                                       cAGATCATCCCGGCGTACCCTAACTGGTCACCGACAACCGGTGCCGCGCCTAACAACTGCTGGACAATTTCACCCAGCGTTTGCTGCTTATCACGATATTCTCGCCCGTCTTCCAGGCCATATATCAGCGCaa          <  1:237011/133‑1 (MQ=255)
                                                                                                       cAGATCATCCCGGCGAACTCTACCTGGTCACCGACAACCGGTGCCGCGCCTAACAACTGCTGGACAATTTCACCCAGCGTTTGCTGCTTATCACGATATTCTCGCCCGTCTTCCAGGCCATATATCAGCGCaa          <  1:84382/133‑1 (MQ=255)
                                                                                                               cccGGCGAACTCTACCTGGTCACCGACAACCGGTGCCGCGCCTAACAACTGCTGGACAATTTCACCCAGCGTTTGCTGCTTATCACGATATTCTCGCCCGTCTTCCAGGCCATATATCAGCGCAACATCAGCa  >  1:116642/1‑133 (MQ=255)
                                                                                                                     |                                                                                                                              
CAACCGGCGATAGATGCTTCATTAAATTTAAGATTCAGCTTCCTCTTCAGCCACCCGTACGCCAATCTTCAACACTTCATTATCTTCTTTCTCGGCCACCGTCCAGATCATCCCGGCAAACTCTACCTGGTCACCGACAACCGGTGCCGCGCCTAACAACTGCTGGACAATTTCACCCAGCGTTTGCTGCTTATCACGATATTCTCGCCCGTCTTCCAGGCCATATATCAGCGCAACATCAGCA  >  NC_000913/1240308‑1240551

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