Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,646,078 G→A 19.0% Q199Q (CAG→CAA fusB → elongation factor G 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,646,0780GA19.0% 36.2 / 4.3 21Q199Q (CAG→CAAfusBelongation factor G 2
Reads supporting (aligned to +/- strand):  ref base G (10/7);  new base A (2/2);  total (12/9)
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.99e-01

GCCGACTTCCTGCGCGTGGTCAAGCAGATCGATCAACGTCTGGGGCACCACCCGGTGCCCATCCAGCTGGCCATCGGCAGCGAAGAAAACTTCATGGGCCAGATCGACCTGGTGAAGATGAAGGCCATCTACTGGAACGATGCCGACCAGGGCACCAGCTACCGCGAAGAAGAAATCCCCGCCGAACTGAAGGCTCTGGCCGACGAATGGCGCGCGCACATGATCGAAGCTGCGGCCGAGGCCAATGACGAGCTGACCATGAAGTTTCTCGATGGCGAGGAGCTGAGCATCGAAGAG  >  NC_002947/4645929‑4646225
                                                                                                                                                     |                                                                                                                                                   
gCCGACTTCCTGCGCGTGGTCAAGCAGATCGATCAACGTCTGGGGCACCACCCGGTGCCCATCCAGCTGGCCATCGGCAGCGAAGAAAACTTCATGGGCCAGATCGACCTGGTGAAGATGAAGGCCATCTACTGGAACGATGCCGACCAg                                                                                                                                                     <  1:577011/150‑1 (MQ=255)
         cTGCGCGTGGTCAAGCAGATCGATCAACGTCTGGGGCACCACCCGGTGCCCATCCAGCCGGCCATCGGCAGCGAAGAAAACTTCATGGGCCAGATCGACCTGGTGAAGATGAAGGCCATCTACTGGAACGATGCCGACCAGGGCACCAGc                                                                                                                                            >  1:172347/1‑150 (MQ=255)
                     aaGCAGATCGATCAACGTCTGGGGCACCACCGGGTGCCCAAGCAGCTGGCCATCGGCAGAGAAAAAAACTACATCCGCTGGAACGGCAGCGACCAGTTCAAGTCTATCTACTGGAACGATGCCGACCAGGGCACCAGCAACAGCgaagaa                                                                                                                                >  1:408939/1‑150 (MQ=255)
                                                   ccGGTGCCCATCCAGCTGGCCATCGGCAGCGAAGAAAACTTCATGGGCCAGATCGACCTGGTGAAGATGAAGGCCATCTACTGGAACGATGCCGACCAGGGCACCAGCTACCGCGAAGAAGAAATCCCCGCCGAACTGAAGGCTCTGGcc                                                                                                  >  2:722837/1‑150 (MQ=255)
                                                      gTGCCCATCCAGCTGGCCATCGGCAGCGAAGAAAACTTCATGGGCCAGATCGACCTGGTGAAGATGAAGGCCATCTACTGGAACGATGCCGACCAGGGCACCCGCTACCGCGAAGAAGAAATCCCCGCCGAACTGAAGGCTCTGGCcgac                                                                                               >  1:199619/1‑150 (MQ=255)
                                                       tGCCCATCCAGCTGGCCATCGGCAGCGAAGAAAACTTCATGGGCCAGATCGACCTGGTGAAGATGAAGGCCATCTACTGGAACGATGCCGACCAGGGCACCAGCTACCGCGAAGAAGAAATCCCCGCCGAACTGAAGGCTCTGGCcgacg                                                                                              <  2:794061/150‑1 (MQ=255)
                                                                              aGCGAAGAAAACTTCATGGGCCAGATCGACCTGGTGAAGATGAAGGCCATCTACTGGAACGATGCCGACCAGGGCACCAGCTACCGCGAAGAAGAAATCCCCGCCGAACTGAAGGCTCTGGCCGACGAATGGCGCGCGCACATGATCGaa                                                                       >  2:200410/1‑150 (MQ=255)
                                                                                 gaagaaAACTTCATGGGCCAGATCGACCTGGTGAAGATGAAGGCCATCTACTGGAACGATGCCGACCAGGGCACCAGCTACCGCGAAGAAGAAATCCCCGCCGAACTGAAGGCTCTGGCCGACGAATGGCGCGCGCACATGATCGAAGCt                                                                    <  1:387751/150‑1 (MQ=255)
                                                                                   agaaAACTTCATGGGCCAGATCGACCTGGTGAAGATGAAGGCCATCTACTGGAACGATGCCGACCAGGGCACCAGCTACCGCGAAGAAGAAATCCCCGCCGCACTGAAGGCTCTGGCCGACGAATGGCGCGCGCACATGATCGAAGCTGc                                                                  >  1:845063/1‑150 (MQ=255)
                                                                                          ttCATGGGCCAGATCGACCTGGTGAAGATGAGGCCATCTACTGGAACGATGCCGACCAGGGCACCAGCTACCGCGAAGAAGAAATCCCCGCCGAACTGAAGGCTCTGGCCGACGAATGGCGCGCGCCCATGACCGAAGCTGCGGCCGAgg                                                          >  2:773769/1‑150 (MQ=255)
                                                                                            cATGGGCCAGATCGACCTGGTGAAGATGAAGGCCATCTACTGGAACGATGCCGACCAGGGCACCAGCTACCGCGAAGAAGAAATCCCCGCCGAACTGAAGGCTCTGGCCGACGAATGGCGCGCGCACATGATCGAAGCTGCGGCCGAGGc                                                         <  1:722837/150‑1 (MQ=255)
                                                                                                      aTCGACCTGGTGAAGATGAAGGCCATCTACTGGAACGATGCCGACCAAGGCACCAGCTACCGCGAAGAAGAAATCCCCGCCGAACTGAAGGCTCTGGCCGACGAATGgc                                                                                        >  2:802411/1‑109 (MQ=255)
                                                                                                      aTCGACCTGGTGAAGATGAAGGCCATCTACTGGAACGATGCCGACCAAGGCACCAGCTACCGCGAAGAAGAAATCCCCGCCGAACTGAAGGCTCTGGCCGACGAATGgc                                                                                        >  2:596855/1‑109 (MQ=255)
                                                                                                      aTCGACCTGGTGAAGATGAAGGCCATCTACTGGAACGATGCCGACCAAGGCACCAGCTACCGCGAAGAAGAAATCCCCGCCGAACTGAAGGCTCTGGCCGACGAATGgc                                                                                        <  1:802411/109‑1 (MQ=255)
                                                                                                      aTCGACCTGGTGAAGATGAAGGCCATCTACTGGAACGATGCCGACCAAGGCACCAGCTACCGCGAAGAAGAAATCCCCGCCGAACTGAAGGCTCTGGCCGACGAATGgc                                                                                        <  1:596855/109‑1 (MQ=255)
                                                                                                                tGAAGATGAAGGCCATCTACTGGAACGATGCCGACCAGGGCACCAGCTACCGCGAAGAAGAAATCCCCGCCGAACTGAAGGCTCTGGCCGACGAATGGCGCGCGCACATGATCGAAGCTGCGGCCGAGGCCAATGACGAGCTGACCATGa                                     <  2:199619/150‑1 (MQ=255)
                                                                                                                 gAAGATGAAGGCCATCTACTGGAACGATGCCGACCAGGGCACCAGCTACCGCGAAGAAGAAATCCCCGCCGAACTGAAGGCTCTGGCCGACGAATGGCGCGCGCACATGATCGAAGCTGCGGCCGAGGCCAATGACGAGCTGACCATGaa                                    >  2:378469/1‑150 (MQ=255)
                                                                                                                        aaGGCCATCTACTGGAACGATGCCGACCAGGGCACCAGCTACCGCGAAGAAGAAAGCCCCGCCGAACTGAAGGCTCTGGCCGACGAATGGCGCGCGCACATGATCGAAGCTGCGGCCGAGGCCAATGACGAGCTGACCATGAAGTTtctc                             <  2:408939/150‑1 (MQ=255)
                                                                                                                                      gAACGATGCCGACCAGGGCACCAGCTACCGCGAAGAAGAAATCCCCGCCGAACTGAAGGCTCTGGCCGACGAATGGCGCGCGCACATGATCGAAGCTGCGGCCGAGGCCAATGACGAGCTGACCATGAAGttt                                >  2:474417/1‑133 (MQ=255)
                                                                                                                                      gAACGATGCCGACCAGGGCACCAGCTACCGCGAAGAAGAAATCCCCGCCGAACTGAAGGCTCTGGCCGACGAATGGCGCGCGCACATGATCGAAGCTGCGGCCGAGGCCAATGACGAGCTGACCATGAAGttt                                <  1:474417/133‑1 (MQ=255)
                                                                                                                                      gAACGATGCCGACCAGGGCACCAGCTACCGCGAAGAAGAAATCCCCGCCGAACTGAAGGCTCTGGCCGACGAATGGCGCGCGCACATGATCGAAGCTGCGGCCGAGGCCAATGACGAGCTGACCATGAAGTTTCTCGATGGCGAGGAGCt               <  2:172347/150‑1 (MQ=255)
                                                                                                                                                   cAGGGCACCAGCTACCGCGAAGAAGAAATCCCCGCCGAACTGAAGGCTCTGGCCGACGAATGGCGCGCGCACATGATCGAAGCTGCGGCCGAGGCCAATGACGAGCTGACCATGAAGTTTCTCGATGGCGAGGAGCTGAGCATCGAagag  >  2:457729/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GCCGACTTCCTGCGCGTGGTCAAGCAGATCGATCAACGTCTGGGGCACCACCCGGTGCCCATCCAGCTGGCCATCGGCAGCGAAGAAAACTTCATGGGCCAGATCGACCTGGTGAAGATGAAGGCCATCTACTGGAACGATGCCGACCAGGGCACCAGCTACCGCGAAGAAGAAATCCCCGCCGAACTGAAGGCTCTGGCCGACGAATGGCGCGCGCACATGATCGAAGCTGCGGCCGAGGCCAATGACGAGCTGACCATGAAGTTTCTCGATGGCGAGGAGCTGAGCATCGAAGAG  >  NC_002947/4645929‑4646225

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: