Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,244,402 G→C 14.6% S221R (AGC→AGG aruH ← arginine‑‑pyruvate transaminase AruH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,244,4020GC14.6% 67.0 / 3.2 26S221R (AGC→AGGaruHarginine‑‑pyruvate transaminase AruH
Reads supporting (aligned to +/- strand):  ref base G (11/11);  new base C (2/2);  total (13/13)
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

CGTAGAGCATGCACAGGGCGAGGTTTTCCAGATGCGTGGCCAGTGCTGTCGAGCCGACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTGGCGCGCGGCAGGCTCGCCCCTGAAGGGTTGTGC  >  NC_002947/4244260‑4244543
                                                                                                                                              |                                                                                                                                              
cGTACAGCATGCACCGGGCGAAGGTTCTTGGGTTTGTTGTCATTGCTGTGGAGGTGGCTACGCACCCCATGCGCCAGCCGGTAATGTCATGCGATTTGGACAGGCTGTTGAACGTTACAGTGCGCTCAGCTATGCCCGGCAGCCTGCCgg                                                                                                                                         <  2:222525/150‑1 (MQ=255)
  tAGAGCATGCACAGGGCGAGGTTTTCCAGATGCGTGGCCAGTGCTGTCGAGCCGACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCgggg                                                                                                                                       >  2:166969/1‑150 (MQ=255)
   agagCATGCACAGGGCGAGGTTTTTTTTTTTCGTGGCCAGTGCTGTCGAGCCGACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGCCTGCCgggca                                                                                                                                      <  1:21/150‑3 (MQ=255)
        aTGCACAGGGCGAGGTTTTCCAGATGCGTGGCCAGTGCTGTCGAGCCGACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACa                                                                                                                                 >  1:302004/1‑150 (MQ=255)
            acaGGGCGAGGTTTTCCAGATGCGTGGCCAGTGCTGTCGAGCCGACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGCGACTGACATGCt                                                                                                                             >  1:401024/1‑150 (MQ=255)
             caGGGCGAGGTTTTCCAGATGCGTGGCCAGTGCTGTCGAGCCGACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTc                                                                                                                            <  2:84800/150‑1 (MQ=255)
                                          gTGCTGTCGAGCCGACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCt                                                                                              <  2:409867/150‑1 (MQ=255)
                                             cTGTCGAGCCGACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGt                                                                                           >  1:454241/1‑150 (MQ=255)
                                             cTGTCGAGCCGACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGt                                                                                           <  1:609172/150‑1 (MQ=255)
                                                  gAGCCGACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCagag                                                                                      <  2:490731/150‑1 (MQ=255)
                                                    gCCGACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGAt                                                                                    >  1:286105/1‑150 (MQ=255)
                                                       gACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCGAGCTGCTGTatgaag                                                                                              >  1:306475/1‑132 (MQ=255)
                                                       gACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCGAGCTGCTGTatgaag                                                                                              <  2:306475/137‑6 (MQ=255)
                                                        accaccCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGc                                                                                <  1:531049/150‑1 (MQ=255)
                                                             ccacccCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAg                                                                           >  2:522994/1‑150 (MQ=255)
                                                                    aTGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATgg                                                                    <  2:587609/150‑1 (MQ=255)
                                                                       cgcCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCa                                                                 >  1:426583/1‑150 (MQ=255)
                                                                                                       gCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGa                                                                                                                                   >  1:157136/1‑53 (MQ=255)
                                                                                                       gCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGa                                                                                                                                   <  2:157136/53‑1 (MQ=255)
                                                                                                           ttCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCa                                       >  1:401683/1‑140 (MQ=255)
                                                                                                           ttCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCa                                       <  2:401683/140‑1 (MQ=255)
                                                                                                              aaCGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATgg                                                                    <  2:600426/108‑1 (MQ=255)
                                                                                                              aaCGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATgg                                                                    >  1:600426/1‑108 (MQ=255)
                                                                                                                                  cATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTGGCGCGCGGCAGGCTCGCCCCTGAAGGGtt      <  1:70559/150‑1 (MQ=255)
                                                                                                                                    gccccGGCAGCCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTGGCGCGCGGCAGGCTCGCCCCTGAAGGGTtgt    >  1:238576/3‑150 (MQ=255)
                                                                                                                                      ccaggCAGCCTTCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTGGCGCGCGGCAGGCTCGCCCCTGAAGGGTTGTGc  >  2:609389/4‑150 (MQ=255)
                                                                                                                                              |                                                                                                                                              
CGTAGAGCATGCACAGGGCGAGGTTTTCCAGATGCGTGGCCAGTGCTGTCGAGCCGACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTGGCGCGCGGCAGGCTCGCCCCTGAAGGGTTGTGC  >  NC_002947/4244260‑4244543

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