Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,552,855 C→A 21.0% A113A (GCG→GCT rhtA ← threonine/homoserine exporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,552,8550CA21.0% 49.8 / 2.2 19A113A (GCG→GCTrhtAthreonine/homoserine exporter
Reads supporting (aligned to +/- strand):  ref base C (9/6);  new base A (2/2);  total (11/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCATAGGCGGCGCCCACTGGGTCGAGTGCCTGGCCGCTGTGACCTGCCGGCAGCAACAGCAGCAGGCCTGCGATGGCCAGGGCAATCCACAGGAAATCGATCGCTCGGCGCGAGGAGAAGATCGCCACGGCCAACGGGCCGGTGAATTCCAGGGCGACGGCGATACCGATTGGTACCGTTTGCAGGGCCATATAGAAGAGGAAGTTCATGCCGCCCAGTGCCATGCCGTAGATGACCACATTGC  >  NC_002947/2552733‑2552976
                                                                                                                          |                                                                                                                         
gcATAGGCGGCGCCCACTGGGTCGAGTGCCTGGCCGCTGTGACCTGCCGGCAGCAACAGCAGCAGGCCTGCGATGGCCAGGGCAATCCACAGGAAATCGATCGCTCGGCGCGAGGAGAAGATCGCCACGGCCAACg                                                                                                              >  2:265715/1‑136 (MQ=255)
      gcggcgCCCACTGGGTCGAGTGCCTGGCCGCTGTGACCTGCCGGCAGCAACAGCAGCAGGCCTGCGATGGCCAGGGCAATCCACAGGAAATCGATCGCTCGGCGCGAGGAGAAGATCGCCACGGCCAACGGGCCgg                                                                                                        >  1:477418/1‑136 (MQ=255)
          cgcCCACTGGGTCGAGTGCCTGGCCGCTGTGACCTGCCGGCAGCAACAGCAGCAGGCCTGCGATGGCCAGGGCAATCCACAGGAAATCGATCGCTCGGCGCGAGGAGAAGATCGCCACGGCCAACGGGCCGGTGaa                                                                                                    <  2:213678/136‑1 (MQ=255)
            cccACTGGGTCGAGTGCCTGGCCGCTGTGACCTGCCGGCAGCAACAGCAGCAGGCCTGCGATGGCCAGGGCAATCCACAGGAAATCGATCGCTCGGCGCGAGGAGAAGATCGCCACGGCCAACGGGCCGGTGAAtt                                                                                                  >  2:306435/1‑136 (MQ=255)
                                 ccGCTGTGACCTGCCGGCAGCAACAGCAGCAGGCCTGCGATGGCCAGGGCAATCCACAGGAAATCGATCGCTCGGCGCGAGGAGAAGATCGCCACGGCCAACGGGCCGGTGAATTCCAGGGCGACGGCGATACCGa                                                                             >  1:265109/1‑136 (MQ=255)
                                   gCTGTGACCTGCCGGCAGCAACAGCAGCAGGCCTGCGATGGCCAGGGCAATCCACAGGAAATCGATCGCTCGGCGCGAGGAGAAGATCGCCACGGCCAACGGGCCGGTGAATTCCAGGGCGACGGCGATACCGAtt                                                                           <  1:265715/136‑1 (MQ=255)
                                             gCCGGCAGCAACAGCAGCAGGCCTGCGATGGCCAGGGCAATCCACAGGAAATCGATCGCTCGGCGCGAGGAGAAGATCGCCACGGCCAACGGGCCGGTGAATTCCAGGGCGACGGCGATACCGATTGGTACCGttt                                                                 <  1:256556/136‑1 (MQ=255)
                                             gCCGGCAGCAACAGCAGCAGGCCTGCGATGGCCAGGGCAATCCACAGGAAATCGATCGCTCGGCGCGAGGAGAAGATCGCCACGGCCAACGGGCCGGTGAATTCCAGGGCGACGGCGATACCGATTGGTACCGttt                                                                 <  2:59183/136‑1 (MQ=255)
                                                 gcagcaACAGCAGCAGGCCTGCGATGGCCAGGGCAATCCACAGGAAATCGATCGCTCGGCGCGAGGAGAAGATAGCCACGGCCAACGGGCCGGTGAATTCCAGGGCGACGGCGATACCGATTGGTACCGTTTGCAg                                                             >  2:348096/1‑136 (MQ=255)
                                                 gcagcaACAGCAGCAGGCCTGCGATGGCCAGGGCAATCCACAGGAAATCGATCGCTCGGCGCGAGGAGAAGATAGCCACGGCCAACGGGCCGGTGAATTCCAGGGCGACGGCGATACCGATTGGTACCGTTTGCAg                                                             >  2:350938/1‑136 (MQ=255)
                                                      aaCAGCAGCAGGCCTGCGATGGCCAGGGCAATCCACAGGAAATCGATCGCTCGGCGCGAGGAGAAGATCGCCACGGCCAACGGGCCGGTGAATTCCAGGGCGACGGCGATACCGATTGGTACCGTTTGCAGGGCCa                                                        >  1:520494/1‑136 (MQ=255)
                                                                   cTGCGATGGCCAGGGCAATCCACAGGAAATCGATCGCTCGGCGCGAGGAGAAGATCGCCACGGCCAACGGGCCGGTGAATTCCAGGGCGACGGCGATACCGATTGGTACCGTTTGCAGGGCCATATAGAAGAGGaa                                           >  2:56637/1‑136 (MQ=255)
                                                                      cGATGGCCAGGGCAATCCACAGGAAATCGATCGCTCGGCGCGAGGAGAAGATCGCCACGGCCAACGGGCCGGTGAATTCCAGGGCGACGGCGATACCGATTGGTACCGTTTGCAGGGCCATATAGAAGAGGAAGtt                                        >  1:513265/1‑136 (MQ=255)
                                                                         tGGCCAGGGCAATCCACAGGAAATCGATCGCTCGGCGCGAGGAGAAGATCGCCACGGCCAACGGGCCGGTGAATTCCAGGGCGACGGCGATACCGATTGGTACCGTTTGCAGGGCCATATAGAAGAGGAAGTTCAt                                     >  1:429767/1‑136 (MQ=255)
                                                                                   aaTCCACAGGAAATCGATCGCTCGGCGCGAGGAGAAGATCGCCACGGCCAACGGGCCGGTGAATTCCAGGGCGACGGCGATACCGATTGGTACCGTTTGCAGGGCCATATAGAAGAGGAAGTTCATGCCGCCCAGt                           >  1:149377/1‑136 (MQ=255)
                                                                                           ggAAATCGATCGCTCGGCGCGAGGAGAAGATCGCCACGGCCAACGGGCCGGTGAATTCCAGGGCGACGGCGATACCGATTGGTACCGTTTGCAGGGCCATATAGAAGAGGAAGTTCATGCCGCCCAGTGCCATGcc                   <  2:265109/136‑1 (MQ=255)
                                                                                                            cgcgAGGAGAAGATCGCCACGGCCAACGGGCCGGTGAATTCCAGGGCGACGGCGATACCGATTGGTACCGTTTGCAGGGCCATATAGAAGAGGAAGTTCATGCCGCCCAGTGCCATGCCGTAGATGACCACATTgc  <  2:403973/136‑1 (MQ=255)
                                                                                                            cgcgAGGAGAAGATAGCCACGGCCAACGGGCCGGTGAATTCCAGGGCGACGGCGATACCGATTGGTACCGTTTGCAGGGCCATATAGAAGAGGAAGTTCATGCCGCGCAGTGCCATGCCGTAGATGACCACATTgc  <  1:350938/136‑1 (MQ=255)
                                                                                                            cgcgAGGAGAAGATAGCCACGGCCAACGGGCCGGTGAATTCCAGGGCGACGGCGATACCGATTGGTACCGTTTGCAGGGCCATATAGAAGAGGAAGTTCATGCCGCGCAGTGCCATGCCGTAGATGACCACATTgc  <  1:348096/136‑1 (MQ=255)
                                                                                                                          |                                                                                                                         
GCATAGGCGGCGCCCACTGGGTCGAGTGCCTGGCCGCTGTGACCTGCCGGCAGCAACAGCAGCAGGCCTGCGATGGCCAGGGCAATCCACAGGAAATCGATCGCTCGGCGCGAGGAGAAGATCGCCACGGCCAACGGGCCGGTGAATTCCAGGGCGACGGCGATACCGATTGGTACCGTTTGCAGGGCCATATAGAAGAGGAAGTTCATGCCGCCCAGTGCCATGCCGTAGATGACCACATTGC  >  NC_002947/2552733‑2552976

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: