Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,534,893 G→A 8.9% A297V (GCA→GTA)  mrpD ← K(+)/H(+) antiporter subunit D

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,534,8930GA 8.9% 125.7 / 3.0 45A297V (GCA→GTA) mrpDK(+)/H(+) antiporter subunit D
Reads supporting (aligned to +/- strand):  ref base G (22/19);  new base A (2/2);  total (24/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.04e-01

AAGGGTCGAGTTGACCAGGTAGAACAGCGCAGCGCCGGTCAGGGCCGCCTGGGCAAAACCCACCGCACCCAGCAGCGTGCCGGCCGAAACCAGGATACTGAGTGCCGCCATACGCTCCAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGACCAGGCCGCCGTA‑‑‑AACCAGCCATTCGCCGCCGAAATGCGCCGAGGCCCCGGCCTGCCCGGAGAACAGCAGCGTCCACAGGCGCAGCACGGCGTACAGGCCGACTTTGGTCATGATGGCG  >  NC_002947/2534745‑2535027
                                                                                                                                                    |                                                                                                                                         
aaGGGTCGAGTTGACCAGGTAGAACAGCGCAGCGCCGGTCAGGGCCGCCTGGGCAAAACCCACCGCACCCAGCAGCGTGCCGGCCGAAACCAGGATACTGAGTGCCGCCATACGCTCCAGGCGCTGAGCTGCCAGGATCGATACCGCTg                                                                                                                                           >  5:45698/1‑149 (MQ=255)
       gAGTTGACCAGGTAGAACAGCGCAGCGCCGGTCAGGGCCGCCTGGGCAAAACCCACCGCACCCAGCAGCGTGCCGGCCGAAACCAGGATACTGAGTGCCGCCATACGCTCCAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCACCGc                                                                                                                                    <  3:138122/149‑1 (MQ=255)
                  gTAGAACAGCGCAGCGCCGGTCAGGGCCGCCTGGGCAAAACCCACCGCACCCAGCAGCGTGCCGGCCGAAACCAGGATACTGAGTGCCGCCATACGCTCCAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGACCa                                                                                                                         <  3:291747/149‑1 (MQ=255)
                   tAGAACAGCGCAGCGCCGGTCAGGGCCGCCTGGGCAAAACCCACCGCACCCAGCAGCGTGCCGGCCGAAACCAGGAGACTGAGTGCCGCCATACGCTCCAGGAGCTGAGCAGCCAGGATCGATACCGCTGCCCCCGCCAGGGTGACCAg                                                                                                                        >  8:262155/1‑149 (MQ=255)
                      aaCAGCGCAGCGCCGGTCAGGGCCGCCTGGGCAAAACCCACCGCACCCAGCAGCGTGCCGGCCGAAACCAGGATACTGAGTGCCGCCATACGCTCCAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGACCAGgcc                                                                                                                     >  2:88989/1‑149 (MQ=255)
                                          ggCCGCCTGGGCAAAACCCACCGCACCCAGCAGCGTGCCGGCCGAAACCAGGATACTGAGTGCCGCCATACGCTCCAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGACCAGGCCGCCGTA‑‑‑AACCAGCCATTcgc                                                                                              >  3:208228/1‑149 (MQ=255)
                                                     cAAAACCCACCGCACCCAGCAGCGTGCCGGCCGAAACCAGGATACTGAGTGCCGCCATACGCTCCAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGGCCAGGCCGCCGTA‑‑‑AACCAGCCATTCGCCGCCGAAATgc                                                                                   >  6:191433/1‑149 (MQ=255)
                                                       aaaCCCACCGCACCCAGCAGCGTGCCGGCCGAAACCAGGATACTGAGTGCCGCCATACGCTCCAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGACCAGGccgccg                                                                                                                 >  4:94164/1‑120 (MQ=255)
                                                       aaaCCCACCGCACCCAGCAGCGTGCCGGCCGAAACCAGGATACTGAGTGCCGCCATACGCTCCAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGACCAGGccgccg                                                                                                                 <  3:94164/120‑1 (MQ=255)
                                                          cccACCGCACCCAGCAGCGTGCCGGCCGAAACCAGGATACTGAGTGCCGCCATACGCTCCAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGACCAGGCCGCCGTA‑‑‑AACCAGCCATTCGCCGCCGAAATGCGcc                                                                                >  7:299598/1‑147 (MQ=255)
                                                            cACCGCACCCAGCAGCGTGCCGGCCGAAACCAGGATACTGAGTGCCGCCATACGCTCCAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGACCAGGCCGCCGTA‑‑‑AACCAGCCATTCGCCGCCGAAATGCGCCGAgg                                                                            >  6:74033/1‑149 (MQ=255)
                                                                  aCCCAGCAGCGTGCCGGCCGAAACCAGGATACTGAGTGCCGCCATACGCTCCAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGACCAGGCCGCCGTA‑‑‑AACCAGCCATTCGCCGCCGAAATGCGCCGAGGCCCCgg                                                                      >  5:239514/1‑149 (MQ=255)
                                                                   cccAGCAGCGTGCCGGCCGAAACCAGGATACTGAGTGCCGCCATACGCTCCAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGGCCAGGCCGCCGTA‑‑‑AACCAGCCATTCGCCGCCGAAATGCGCCGAGGCCCCGGc                                                                     <  5:191433/149‑1 (MQ=255)
                                                                     cagcagCGTGCCGGCCGAAACCAGGATACTGAGTGCCGCCATACGCTCCAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGACCAGGCCGCCGTA‑‑‑AACCAGCCATTCGCCGCCGAAATGCGCCGAGGCCCCGGCCt                                                                   >  3:201381/1‑149 (MQ=255)
                                                                          gcGTGCCGGCCGAAACCAGGATACTGAGTGCCGCCATACGCTCCAGGCGCTGAGCTGCCAGGATCGATACCGCTACCACCGCCAGGGTGACCAGGCCGCCGTA‑‑‑AACCAGCCATTCGCCGCCGAAATGCGCCGAGGCCCCGGCCTGCCCg                                                              >  6:103182/1‑149 (MQ=255)
                                                                          gcGTGCCGGCCGAAACCAGGATACTGAGTGCCGCCATACGCTCCAGGCGCTGAGCTGCCAGGATCGATACCGCTACCACCGCCAGGGTGACCAGGCCGCCGTA‑‑‑AACCAGCCATTCGCCGCCGAAATGCGCCGAGGCCCCGGCCTGCCCg                                                              >  4:260523/1‑149 (MQ=255)
                                                                          gcGTGCCGGCCGAAACCAGGATACTGAGTGCCGCCATACGCTCCAGGCGCTGAGCTGCCAGGATCGATACCGCTACCACCGCCAGGGTGACCAGGCCGCCGTA‑‑‑AACCAGCCATTCGCCGCCGAAATGCGCCGAGGCCCCGGCCTGCCCg                                                              <  5:103182/149‑1 (MQ=255)
                                                                          gcGTGCCGGCCGAAACCAGGATACTGAGTGCCGCCATACGCTCCAGGCGCTGAGCTGCCAGGATCGATACCGCTACCACCGCCAGGGTGACCAGGCCGCCGTA‑‑‑AACCAGCCATTCGCCGCCGAAATGCGCCGAGGCCCCGGCCTGCCCg                                                              <  3:260523/149‑1 (MQ=255)
                                                                            gTGCCGGCCGAAACCAGGATACTGAGTGCCGCCATACGCTCCAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGACCAGGCCGCCGTA‑‑‑AACCAGCCATTCGCCGCCGAAATGCGCCGAGGCCCCGGCCTGCCCGga                                                            <  4:208228/149‑1 (MQ=255)
                                                                             tGCCGGCCGAAACCAGGATACTGAGTGCCGCCATACGCTCCAGGCGCTTAGCTGCCAGGATCGATACCGCTGCGACCGCCAGGGTGACCAGGCCGCCGTA‑‑‑AACCAGCCATTCGCCGCCGAAATGCGCCGAGGCCCCGGCCTGCGCGgag                                                           <  7:262155/149‑1 (MQ=255)
                                                                             tGCCGGCCGAAACCAGGATACTGAGTGCCGCCATACGCTCCAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGACCAGGCCGCCGTA‑‑‑AACCAGCCATTCGCCGCCGAAATGCGCCGAGGCCCCGGCCTGCCCGgag                                                           <  2:105538/149‑1 (MQ=255)
                                                                             tGCCGGCCGAAACCAGGATACTGAGTGCCGCCATACGCTCCAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGACCAGGCCGCCGTA‑‑‑AACCAGCCATTCGCCGCCGAAATGCGCCGAGGCCCCGGCCTGCCCGgag                                                           <  6:45698/149‑1 (MQ=255)
                                                                                  gccgAAACCAGGATACTGAGTGCCGCCATACGCTCCAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGACCAGGCCGCCGTA‑‑‑AACCAGCCATTCGCCGCCGAAATGCGCCGAGGCCCCgg                                                                      >  3:222293/1‑133 (MQ=255)
                                                                                  gccgAAACCAGGATACTGAGTGCCGCCATACGCTCCAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGACCAGGCCGCCGTA‑‑‑AACCAGCCATTCGCCGCCGAAATGCGCCGAGGCCCCgg                                                                      <  4:222293/133‑1 (MQ=255)
                                                                                           aGGATACTGAGTGCCGCCATACGCTCCAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGACCAGGCCGCCGTA‑‑‑AAc                                                                                                         <  7:54401/89‑1 (MQ=255)
                                                                                           aGGATACTGAGTGCCGCCATACGCTCCAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGACCAGGCCGCCGTA‑‑‑AAc                                                                                                         >  8:54401/1‑89 (MQ=255)
                                                                                                  tGAGTGCCGCCATACGCTCCAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGACCAGGCCGCCGTA‑‑‑AACCAGCCATTCGCCGCCGAAATGCGCCGAGGCCCCAGCCTGCCCGga                                                            >  1:191765/1‑127 (MQ=255)
                                                                                                  tGAGTGCCGCCATACGCTCCAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGACCAGGCCGCCGTA‑‑‑AACCAGCCATTCGCCGCCGAAATGCGCCGAGGCCCCAGCCTGCCCGga                                                            <  2:191765/127‑1 (MQ=255)
                                                                                                   gAGTGCCGCCATACGCTCCAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGACCAGGCCGCCGTA‑‑‑AACCAGCCATTCGCCGCCGAAATGCGCCGAGGCCCCGGCCTGCCCGGAGAACAGCAGCGTCCACAGgcg                                       >  2:233387/1‑147 (MQ=255)
                                                                                                   gAGTGCCGCCATACGCTCCAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGACCAGGCCGCCGTA‑‑‑AACCAGCCATTCGCCGCCGAAATGCGCCGAGGCCCCGGCCTGCCCGGAGAACAGCAGCGTCCACAGgcg                                       >  6:155406/1‑147 (MQ=255)
                                                                                                     gTGCCGCCATACGCTCCAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGACCAGGCCGCCGTA‑‑‑AACCAGCCATTCGCCGCCGAAATGCGCCGAGGCCCCGGCCTGCCCGGAGAACAGCAGCGTCCACAGGCgcag                                    >  4:213812/1‑148 (MQ=255)
                                                                                                         cgccATACGCTCCAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGACCAGGCCGCCGTA‑‑‑AACCAGCCATTCGCCGCCGAAATGCGCCGAGGCCCCGGCCTGCCCGGAGAACAGCAGCGTCCACAGGCGCAGCACgg                               >  1:188716/1‑149 (MQ=255)
                                                                                                           ccATACGCTCCAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGACCAGGCCGCCGTATGGAGCCAGCCATTCGCCGGCCAAATGCCCCGACGCCCCCGCCGGCCCCGACAACCGCAACCTACACAGGCCCAacacc                                >  1:98394/1‑148 (MQ=255)
                                                                                                           ccATACGCTCCAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGACCAGGCCGCCGTA‑‑‑AACCAGCCATTCGCCGCCGAAATGCGCCGAGGCCCCGGCCTGCCCGGAGAACAGCAGCGTCCACAGGCGCAGCACGGCg                             >  2:78909/1‑149 (MQ=255)
                                                                                                              tACGCTCCAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGACCAGGCCGCCGTA‑‑‑AACCAGCCATTCGCCGCCGAAATGCGCCGAGGCCCCGGCCTGCCCGGAGAACAGCAGCGTCCACAGGCGCAGCACGGCGTAc                          <  1:88989/149‑1 (MQ=255)
                                                                                                                cGCTCCAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGACCAGGCCGCCGTA‑‑‑AACCAGCCATTCGCCGCCGAAATGCGCCGAGGCCCCGGCCTGCCCGGAGAACAGCAGCGTCCACAGGCGCAGCACGGCGTACa                         <  5:74033/148‑1 (MQ=255)
                                                                                                                     cAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGACCAGGCCGCCGTA‑‑‑AACCAGCCATTCGCCGCCGAAATgcg                                                                                  >  2:89493/1‑86 (MQ=255)
                                                                                                                     cAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGACCAGGCCGCCGTA‑‑‑AACCAGCCATTCGCCGCCGAAATgcg                                                                                  <  1:89493/86‑1 (MQ=255)
                                                                                                                     cAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGACCAGGCCGCCGTA‑‑‑AACCAGCCATTCGCCGCCGAAATgc                                                                                   >  2:296167/1‑85 (MQ=255)
                                                                                                                     cAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGACCAGGCCGCCGTA‑‑‑AACCAGCCATTCGCCGCCGAAATgc                                                                                   <  1:296167/85‑1 (MQ=255)
                                                                                                                     cAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGACCAGGCCGCCGTA‑‑‑AACCAGCCATTCGCCGCCGAAATGCGCCGAGGCCCCGGCCTGCCCGGAGAACAGCAGCGTCCACAGGCGCAGCACGGCGTACAGGCCGa                   <  6:115350/149‑1 (MQ=255)
                                                                                                                                 tGCCAGGATCGATACCGCTGCCACCGCCAGGGTGACCAGGCCGCCGTA‑‑‑AACCAGCCATTCGCCGCCGAAATGCGCCGAGGCCCCGGCCTGCCCGGAGAACAGCAGCGTCCACAGGCGCAGCACGGCGTACAGGCCGACTTTGGTCatga       <  4:201381/149‑1 (MQ=255)
                                                                                                                                   ccAGGATCGATACCGCTGCCACCGCCAGGGTGACCAGGCCGCCGTA‑‑‑AACCAGCCATTCGCCGCCGAAATGCGCCGAGGCCCCGGCCTGCCCGGAGAACAGCAGCGTCCACAGGCGCAGCACGGCGTACAGGCCGACTTTGGTCatgatg     <  8:299598/149‑1 (MQ=255)
                                                                                                                                       gATCGATACCGCTGCCACCGCCAGGGTGACCAGGCCGCCGTA‑‑‑AACCAGCCATTCGCCGCCGAAATGCGCCGAGGCCCCGGCCTGCCCGGGGAACAGCAGCGTCCACAGGCGCAGCACGGCGTACAGGCCGACTTTGGTCATGATGGCg  >  3:203884/1‑148 (MQ=255)
                                                                                                                                            aTACCGCTGCCACCGCCAGGGTGACCAGGCCGCCGTA‑‑‑AACCAGCCATTCGCCGCCGAAATGCGCCGAGGCCCCGGCCTGCCCGGAGAACAGCAGCGTCCACAGGCGCAGCACGGCGTACAGGCCGACTTTg              >  7:208515/1‑131 (MQ=255)
                                                                                                                                            aTACCGCTGCCACCGCCAGGGTGACCAGGCCGCCGTA‑‑‑AACCAGCCATTCGCCGCCGAAATGCGCCGAGGCCCCGGCCTGCCCGGAGAACAGCAGCGTCCACAGGCGCAGCACGGCGTACAGGCCGACTTTg              <  8:208515/131‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                         
AAGGGTCGAGTTGACCAGGTAGAACAGCGCAGCGCCGGTCAGGGCCGCCTGGGCAAAACCCACCGCACCCAGCAGCGTGCCGGCCGAAACCAGGATACTGAGTGCCGCCATACGCTCCAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGACCAGGCCGCCGTA‑‑‑AACCAGCCATTCGCCGCCGAAATGCGCCGAGGCCCCGGCCTGCCCGGAGAACAGCAGCGTCCACAGGCGCAGCACGGCGTACAGGCCGACTTTGGTCATGATGGCG  >  NC_002947/2534745‑2535027

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: ‑