Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,632,887 A→C 8.7% A533A (GCA→GCC lepA → elongation factor 4

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,632,8870AC 8.7% 101.8 / 2.3 42A533A (GCA→GCClepAelongation factor 4
Reads supporting (aligned to +/- strand):  ref base A (19/19);  new base C (2/2);  total (21/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 = 6.96e-01

TCAAGCTGGACGTGCTGATCAACGGCGACAAGGTCGATGCCCTGGCATTGATCGTGCACCGCGACAACGCGGCCTATAAAGGCCGTGCGTTGACCGAGAAGATGAAGGAACTGATCCCTCGGCAGATGTTCGACGTAGCGATCCAGGCAGCCATTGGCGGCCAGATCATTGCGCGGACAACCGTCAAGGCGCTCAGAAAGAACGTACTGGCCAAGTGCTACGGCGGTGACGTCAGCCGTAAGAAGAAGCTGCTGGAGAAGCAGAAGGCCGGTAAGAAACGCATGAAACAGGTGG  >  NC_002947/1632739‑1633032
                                                                                                                                                    |                                                                                                                                                 
tCAAGCTGGACGTGCTGATCAACGGCGACAAGGTCGATGCCCTGGCATTGATCGTGCACCGCGACAACGCGGCCTATAAAGGCCGTGCGTTGACCGAGAAGATGAAGGAACTGATCCCTCGGCAGATGTTCGACGTAGCGATCCAGGCAg                                                                                                                                                  >  1:168384/1‑150 (MQ=255)
              cTGATCAACGGCGACAAGGTCGATGCCCTGGCATTGATCGCGCACCGCGACAACGCGGCCTATAAAGGCCGTGCGTTGACCGAGAAGATGAAGGAACTGATCCCTCGGCAGATGTTCGACGTAGCGATCCAGGCAGCCATTGGCGGCCAg                                                                                                                                    <  1:385169/150‑1 (MQ=255)
                gATCAACGGCGACAAGGTCGATGCCCTGGCATTGATCGTGCACCGCGACAACGCGGCCTATAAAGGCCGTGCGTTGACCGAGAAGATGAAGGAACTGATCCCTCGGCAGATGTTCGACGTAGCGATCCAGGCCGCCAATGGCTGcctgga                                                                                                                                  <  1:390013/150‑5 (MQ=255)
                                 tCGATGCCCTGGCATTGATCGTGCACCGCGACAACGCGGCCTATAAAGGCCGTGCGTTGACCGAGAAGATGAAGGAACTGATCCCTCGGCAGATGTTCGACGTAGCGATCCAGGCAGCCATTGGCGGCCAGATCATTGCGCGGACAACCg                                                                                                                 >  2:453012/1‑150 (MQ=255)
                                        ccTGGCATTGATCGTGCACCGCGACAACGCGGCCTATAAAGGCCGTGCGTTGACCGAGAAGATGAAGGAACTGATCCCTCGGCAGATGTTCGACGTAGCGATCCAGGCAGCCATTGGCGGCCAGATCATTGCGCGGACAACCGTCAAGgc                                                                                                          <  2:168384/150‑1 (MQ=255)
                                              aTTGATCGTGCACCGCGACAACGCGGCCTATAAAGGCCGTGCGTTGACCGAGAAGATGAAGGAACTGATCCCTCGGCAGATGTTCGACGTAGCGATCCAGGCAGCCATTGGCGGCCAGATCATTGCGCGGACAACCGTCAAGGCGCTCag                                                                                                    >  1:186569/1‑150 (MQ=255)
                                              aTTGATCGTGCACCGCGACAACGCGGCCTATAAAGGCCGTGCGTTGACCGAGAAGATGAAGGAACTGATCCCTCGGCAGATGTTCGACGTAGCGATCCAGGCAGCCATTGGCGGCCAGATCATTGCGCGGACAACCGTCAAGGCGCTCag                                                                                                    >  2:749283/1‑150 (MQ=255)
                                                 gATCGTGCACCGCGACAACGCGGCCTATAAAGGCCGTGCGTTGACCGAGAAGATGAAGGAACTGATCCCTCGGCAGATGTTCGACGTAGCGATCCAGGCCGCCAATGGCTGcct                                                                                                                                     <  2:281116/114‑2 (MQ=255)
                                                 gATCGTGCACCGCGACAACGCGGCCTATAAAGGCCGTGCGTTGACCGAGAAGATGAAGGAACTGATCCCTCGGCAGATGTTCGACGTAGCGATCCAGGCCGCCAATGGCTGcct                                                                                                                                     >  1:281116/1‑113 (MQ=255)
                                                     gTGCACCGCGACAACGCGGCCTATAAAGGCCGTGCGTTGACCGAGAAGATGAAGGAACTGATCCCTCGGCAGATGTTCGACGTAGCGATCCAGGCCCCAACAGGCGGCCAGATCATTGCGCGGACAACCGTCAAGGCGCTCAGAAAGAAc                                                                                             >  1:480368/1‑150 (MQ=255)
                                                         aCCGCGACAACGCGGCCTATAAAGGCCGTGCGTTGACCGAGAAGATGAAGGAACTGATCCCTCGGCAGATGTTCGACGTAGCGATCCAGGCAGCCATTGGCGGCCAGATCATTGCGCGGACAACCGTCAAGGCGATCAGAAAGAACGTAc                                                                                         >  1:731803/1‑150 (MQ=255)
                                                          ccGCGACAACGCGGCCTATAAAGGCCGTGCGTTGACCGAGAAGATGAAGGAACTGATCCCTCGGCAGATGTTCGACGTAGCGATCCAGGCAGCCATTGGCGGCCAGATCATTGCGCGGACAACCGTCAAGGCGCTCAGAAAGAACGTACt                                                                                        <  2:186569/150‑1 (MQ=255)
                                                          ccGCGACAACGCGGCCTATAAAGGCCGTGCGTTGACCGAGAAGATGAAGGAACTGATCCCTCGGCAGATGTTCGACGTAGCGATCCAGGCAGCCATTGGCGGCCAGATCATTGCGCGGACAACCGTCAAGGCGCTCAGAAAGAACGTACt                                                                                        >  1:34714/1‑150 (MQ=255)
                                                           cgcgACAACGCGGCCTATAAAGGCCGTGCGTTGACCGAGAAGATGAAGGAACTGATCCCTCGGCAGATGTTCGACGTAGCGATCCAGGCAGCCATTGGCGGCCAGATCATTGCGCGGACAACCGTCAAGGCGCTCAGAAAGAACGTACTg                                                                                       <  2:34714/150‑1 (MQ=255)
                                                               aacgcacgGCCTTTATAGGCC‑‑GCGTCGACCGAGAAGATGAAGGAACTGATCCCTCGGCAGATGTTCGGCGTAGCGATCCAGGCAGCCAATGGCGGCCAGATCATTGGGGGGACAACCGTCAAGGCGCCCAGAAAGAACGTCCTGGCCAAg                                                                                 >  2:411075/7‑150 (MQ=255)
                                                                 aaCGCGGCCTATAAAGGCCGTGCGTTGACCGAGAAGATGAAGGAACTGATCCCTCGGCAGATGTTCGACGTAGCGATCCAGGCAGCCATTGGCGGCCAGATCATTGCGCGGACAACCGTCAAGGCGCTCAGAAAGAACGTACTGGCCAAg                                                                                 >  2:181022/1‑150 (MQ=255)
                                                                  aCGCGGCCTATAAAGGCCGTGCGTTGACCGAGAAGATGACGAAACTGATCCCTCGGCAGATGTTCGACGTAGCGATCCAGGCAGCCATTGGCGGCCAGATCATTGCGCGGACAACCGTCAAGGCGCTCAGAAAGAACGTACTGGCCAAGt                                                                                <  1:649498/150‑1 (MQ=255)
                                                                                       cGTTGACCGAGAAGATGAAGGAACTGATCCCTCGGCAGATGTTCGACGTAGCGCTCCAGGCAGCCATTGGCGGCCAGATCATTGCGCGGACAACCGTCAAGGCGCTCAGAAAGAACGTACTGGCCAAGTGCTACGGCGGTGACGTCAGcc                                                           <  1:295466/150‑1 (MQ=255)
                                                                                              cGAGAAGATGAAGGAACTGATCCCTCGGCAGATGTTCGACGTAGCGATCCAGGCAGCCATTGGCGGCCAGATCATTGCGCGGACAACCGTCAAGGCGCTCAGAAAGAACGTACTGGCCAAGTGCTACGGCGGTGACGt                                                                <  1:238331/138‑1 (MQ=255)
                                                                                              cGAGAAGATGAAGGAACTGATCCCTCGGCAGATGTTCGACGTAGCGATCCAGGCAGCCATTGGCGGCCAGATCATTGCGCGGACAACCGTCAAGGCGCTCAGAAAGAACGTACTGGCCAAGTGCTACGGCGGTGACGt                                                                >  2:238331/1‑138 (MQ=255)
                                                                                                    gaTGAAGGAACTGATCCCTCGGCAGATGTTCGACGTAGCGATCCAGGCAGCCATTGGCGGCCAGATCATTGCGCGGACAACCGTCAAGGCGCTCAGAAAGAACGTACTGGCCAAGTGCTACGGCGGTGACGTCAGCCGTAAGAAGAAgct                                              >  2:481359/1‑150 (MQ=255)
                                                                                                      tGAAGGAACTGATCCCTCGGCAGATGTTCGACGTAGCTATCCAGGCAGCCATTggc                                                                                                                                          <  2:737450/56‑1 (MQ=255)
                                                                                                      tGAAGGAACTGATCCCTCGGCAGATGTTCGACGTAGCGATCCAGGCAGCCATTggc                                                                                                                                          >  1:737450/1‑56 (MQ=255)
                                                                                                      tGAAGGAACTGATCCCTCGGCAGATGTTCGACGTAGCGATCCAGGCAGCCATTGGCGGCCAGATCATTGCGCGGACAACCGTCAAGGCGCTCAGAAAGAACGTACTGGCCAAGTGCTACGGCGGTGACGTCAGCCGTAAGAAGAAcctgc                                            <  2:260453/150‑1 (MQ=255)
                                                                                                        aaGGAACTGATCCCTCGGCAGATGTTCGACGTAGCGATCCAGGCAGCCAGTGGCGGCCAGATCATTGCGCGGACAACCGTCAAGGCGCTCAGAAAGAACGTACTGGCCAAGTGCTACGGCGGTGACGTCAGCCGTAAGAAGAAGctgctg                                          >  2:32837/1‑150 (MQ=255)
                                                                                                              cttaccccTCGGCAGATGTTCGACGTAGCGATCCAGGCAGCAATTGGCGGCCAGATCATTGCGCGGACAACCGTCAAGGCGCTAAGAAAGAACGTACTGGCCAAGTGCTACGGCGGTGACGTCAGCCGTAAGAAGAAGCTGCTGGAGAAg                                    <  2:731803/145‑1 (MQ=255)
                                                                                                              cTGATCCCTCGGCAGATGTTCGACGTAGCGATCCAGGCAGCCATTGGCGGCCAGATCATTGCGCGGACAACCGTCAAGGCGCTCAGAAAGAACGTACTGGCCAAGTGCTACGGCGGTGACGTCAGCCGTAAGAAGAAGCTGCTGGAGAAg                                    <  1:749283/150‑1 (MQ=255)
                                                                                                                         gCAGATGTTCGACGTAGCGATCCAGGCAGCCATTGGCGGCCAGATCATTGCGCGGACAACCGTCAAGGCGCTCAGAAAGAACGTa                                                                                          <  2:798350/85‑1 (MQ=255)
                                                                                                                         gCAGATGTTCGACGTAGCGATCCAGGCAGCCATTGGCGGCCAGATCATTGCGCGGACAACCGTCAAGGCGCTCAGAAAGAACGTa                                                                                          >  1:798350/1‑85 (MQ=255)
                                                                                                                         gCAGATGTTCGACGTAGCGATCCAGGCAGCCATTGGCGGCCAGATCATTGCGCGGACAACCGTCAAGGCGCTCAGAAAGAACGTACTGGCCAAGTGCTACGGCGGTGACGTCAGCCGTAAGAAGAAGCTGCTGGAGAAGCAGAAGGCCgg                         <  1:481359/150‑1 (MQ=255)
                                                                                                                                ttCGACGTAGCGATCCAGGCAGCCATTGGCGGCCAGATCATTGCGCGGACACCCGTCAAGGCGCACAGCACGCACGTACAGGCCAACTGCTACGCCGGTGACGACAGCCCCAAGAAGAAGCTGCTGGAGAAGCAGAAGGCCGGTAAGaaa                  >  2:312787/1‑150 (MQ=255)
                                                                                                                                   gACGTAGCGATCCAGGCAGCCATTGGCGGCCAGATCATTGCGCGGACAACCGTCAAGGCGCTCAGAAAGAACGTACTGGCCAAGTGCTACGGCGGTGACGTCAGCCGTAAGAAGAAGCTGCTGGAGAAGCAGAAGGCCGGTAAGAAACGc               <  1:453012/150‑1 (MQ=255)
                                                                                                                                    aCGTAGCGATCCAGGCAGCCATTGGCGGCCAGATCATTGCGCGGACAACCGTCAAGGCGCTCAGAAAGAACGTACTGGCCAAGTGCTACGGCGGTGACGTCAGCCGTAAGAAGAAGCTGCTGGAGAAGCAGAAGGCCGGTAAGAAACGCa              >  2:768805/1‑150 (MQ=255)
                                                                                                                                      gTAGCGATCCAGGCAGCCATTGGCGGCCAGATCATTGCGCGGACAACCGTCAAGGCGCTCAGAAAGAACGTACTGGCCAAGTGCTACGGCGGTGACGTCAGCCGTaagaa                                                    >  2:27860/1‑110 (MQ=255)
                                                                                                                                      gTAGCGATCCAGGCAGCCATTGGCGGCCAGATCATTGCGCGGACAACCGTCAAGGCGCTCAGAAAGAACGTACTGGCCAAGTGCTACGGCGGTGACGTCAGCCGTaagaa                                                    <  1:27860/110‑1 (MQ=255)
                                                                                                                                            aTCCAGGCAGCCATTGGCGGCCAGATCATTGCGCGGACAACCGTCAAgg                                                                                                           <  2:243870/49‑1 (MQ=255)
                                                                                                                                            aTCCAGGCAGCCATTGGCGGCCAGATCATTGCGCGGACAACCGTCAAgg                                                                                                           >  1:243870/1‑49 (MQ=255)
                                                                                                                                              ccAGGCAGCCATTGGCGGCCAGATCATTGCGCGGACAACCGTCAAGGCGCTCAGAAAGAACGTACTGGCCAAGTGCTACGGCGGTGACGTCAGCCGTAAGAAGAAGCTGCTGGAGAAGCAGAAGGCCGGTAAGAAACGCATGAAACAGGt    >  1:523656/1‑150 (MQ=255)
                                                                                                                                               cAGGCAGCCATTGGCGGCCAGATCATTGCGCGGACAACCGTCAAGGCGCTCAGAAAGAACGTACTGGCCAAGTGCTACGGCGGTGACGTCAGCCGTAAGAAGAAGCTGCTGGAGAAGCAGAAGGCCGGTAAGAAACGCATGAAACAGGTg   <  2:549304/150‑1 (MQ=255)
                                                                                                                                                aGGCAGCCATTGGCGGCCAGATCATTGCGCGGACAACCGTCAAGGCGCTCAGAAAGAACGTACTGGCCAAGTGCTATGGAGGTGACGTCATACGGAAGGATAAGCTGCTGGAGAAGCAGAAGGCCGGTAAGAAACGCATGAAACAGGTgg  <  1:531711/150‑1 (MQ=255)
                                                                                                                                                aGGCAGCCATTGGCGGCCAGATCATTGCGCGGACAACCGTCAAGGCGCTCAGAAAGAACGTACTGGCCAAGTGCTACGGCGGTGACGTCAGCCGTAAGAAGAAGCTTCTGGAGAAGCAGAAGGCCGGTAAGAAACGCATGAAACAGGTgg  <  1:181022/150‑1 (MQ=255)
                                                                                                                                                aGGCAGCCATTGGCGGCCAGATCATTGCGCGGACAACCGTCAAGGCGCTCAGAAAGAACGTACTGGCCAAGTGCTACGGCGGTGACGTCAGCCGTAAGAAGAAGCTGCTGGAGAAGCAGAAGGCCGGTAAGAAACGCATGAAACAGGTgg  >  2:374667/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
TCAAGCTGGACGTGCTGATCAACGGCGACAAGGTCGATGCCCTGGCATTGATCGTGCACCGCGACAACGCGGCCTATAAAGGCCGTGCGTTGACCGAGAAGATGAAGGAACTGATCCCTCGGCAGATGTTCGACGTAGCGATCCAGGCAGCCATTGGCGGCCAGATCATTGCGCGGACAACCGTCAAGGCGCTCAGAAAGAACGTACTGGCCAAGTGCTACGGCGGTGACGTCAGCCGTAAGAAGAAGCTGCTGGAGAAGCAGAAGGCCGGTAAGAAACGCATGAAACAGGTGG  >  NC_002947/1632739‑1633032

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: