Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,719,426 T→C 38.5% intergenic (‑321/+2) kgtP ← / ← rrfG alpha‑ketoglutarate transporter/5S ribosomal RNA of rrnG operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,719,4260TC38.5% 5.9 / 15.9 24intergenic (‑321/+2)kgtP/rrfGalpha‑ketoglutarate transporter/5S ribosomal RNA of rrnG operon
Reads supporting (aligned to +/- strand):  ref base T (9/4);  new base C (4/5);  total (13/11)
Fisher's exact test for biased strand distribution p-value = 3.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01

CGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATCAGACCGC  >  CP009273/2719281‑2719568
                                                                                                                                                 |                                                                                                                                              
cGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCTTTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGAt                                                                                                                                             <  1:153012/149‑1 (MQ=255)
                aGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGAGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGGTCCCTa                                                                                                                             >  1:203332/1‑149 (MQ=255)
                      gtaaGCCCTATCCGTTGGTGAATTTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTTAGGATGGCCTTTCGCATAATTGGATGCCTGGCAGTTCCCTACTCTCg                                                                                                                       <  1:80336/147‑1 (MQ=25)
                                tCCGTTGCTGAATCTGGCAAATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGACAAATGCCTGTATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGgaga                                                                                                             >  1:241490/1‑149 (MQ=25)
                                     tGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGGGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCca                                                                                                        <  2:21605/149‑1 (MQ=255)
                                                 cATATGTGGGAGTATAGGTCGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCTATCCGTCAGGATGGCCTTTCGCATAATTCGATGTCTGGTAGTTACCTACTCTCGCATGGGGAGACCCCACACTACCATCgg                                                                                            <  1:242333/149‑1 (MQ=25)
                                                      gtgGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATCATTTGATGCCTGGCAGTTCCCTACTCTCGCCTGGGGAGCCCCCACACTACCCTCggacgta                                                                                       >  2:202706/1‑144 (MQ=21)
                                                       tgGGAGTATAAGACGCGCAGCGTCGCATCCGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTAc                                                                                      >  1:199072/1‑149 (MQ=255)
                                                       tgGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTTTTTTTTTTTGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTAc                                                                                      >  2:50674/1‑149 (MQ=25)
                                                          gattctcacaCGCGCAGCGTCGCATCAGGCCTTTTTTTCTCCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTCGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGc                                                                                   <  2:254169/140‑1 (MQ=25)
                                                           aGTATAAGACGCGCAGCGTCGCATCAGGCTTGTCTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTCGATGCCTTGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCg                                                                                  <  1:128485/149‑1 (MQ=37)
                                                                   aCGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTCGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACtt                                                                          <  2:138600/149‑1 (MQ=255)
                                                                            gTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAAGATGGCCTTTAGCATAATACGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTca                                                                 >  1:72507/1‑148 (MQ=255)
                                                                                  tcgggCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGGCATGATGGCCTTTCGCATAATTCGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATgg                                                           <  2:91173/146‑1 (MQ=255)
                                                                                         tttttttCTGCGCCAATGCAACAAGGCCATCCGTCAGGATGACATTTAGCATAATTTGATGCCTGACAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTAAgg                                                    >  1:171735/1‑149 (MQ=18)
                                                                                         tttttttCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTCGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCGGAGTTCGGCATGGGGTCAgg                                                    >  2:220144/1‑149 (MQ=255)
                                                                                                      cAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACcgcg                                       <  1:262603/149‑1 (MQ=255)
                                                                                                          gCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGCCAGGTGGGACCACCGCGCTAc                                   >  1:268809/1‑149 (MQ=255)
                                                                                                                 ggCCATCCGTCAGGATGGCCTTTCGCATAATTCGATGCCTGGCAGTTCCCTACTCTCGCATAGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGgccgcc                            >  1:112399/1‑149 (MQ=255)
                                                                                                                            ggtaTGTCCTTTCGCATAAGTTGATGCCTGGCAGTTCCATTTTCGAGCATGGGGATACCCCACACTACCACCGGCGCTACGGCGTTTCTCTTTTGAGTTCGGCATGCGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCt                 <  2:90671/146‑1 (MQ=11)
                                                                                                                            agtggggCCTTTCGCATAAGTGGATGCCGGGCAGTTCGCTGCACGCGCATGGGGAGACCCCACACTACTAGCGGGGCTACGGCGTGTCACTTCTGAGTGCGGCATGGGGTCAGGTGGGACCCCCGCGCTACGGCCGGCAGGCTAATTCt                 <  2:90661/144‑1 (MQ=11)
                                                                                                                                    cTTTCGCATAATTCGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGGCAGGCAAATTCTGTTTTATc         >  1:316256/1‑149 (MQ=255)
                                                                                                                                      ttCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATCAg       >  1:35765/1‑149 (MQ=34)
                                                                                                                                           aTAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATCAGACCGc  >  1:11817/1‑149 (MQ=11)
                                                                                                                                                 |                                                                                                                                              
CGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATCAGACCGC  >  CP009273/2719281‑2719568

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