Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,719,426 T→C 42.3% 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,4260TC42.3% 7.0 / 25.1 26intergenic (‑321/+2)kgtP/rrfGalpha‑ketoglutarate transporter/5S ribosomal RNA of rrnG operon
Reads supporting (aligned to +/- strand):  ref base T (8/7);  new base C (4/7);  total (12/14)
Fisher's exact test for biased strand distribution p-value = 4.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

ACGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTAT  >  CP009273/2719280‑2719560
                                                                                                                                                  |                                                                                                                                      
acGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATAAGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGa                                                                                                                                      >  2:2349/1‑149 (MQ=255)
     tGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTCGa                                                                                                                                      >  1:83080/1‑144 (MQ=255)
               cAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTccc                                                                                                                       >  2:245078/1‑149 (MQ=255)
                     tgggAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCATCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGTCATCTGTCAGGATGGCCTTTCGCATAATTCGATGCCTGGCAGTTCCCTActct                                                                                                                 <  2:59802/148‑1 (MQ=255)
                       ggAAGCCGTATCCGTTGCTGAATCTGGCCTATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTCGATGCCTGGCAGTTCCCTAc                                                                                                                    <  1:24855/144‑1 (MQ=255)
                            acGTATCCGTTGCTGAATCTGGCAAATGTGGGCGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTCGCTGCCTGGCAGTTCCCTACTCTCg                                                                                                               >  1:144067/2‑144 (MQ=255)
                               tATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGTATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGga                                                                                                       <  2:5899/149‑1 (MQ=255)
                                       gCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTCGTTGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCcac                                                                                               <  2:89577/149‑1 (MQ=255)
                                         tGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGTTGGTCTTTCGCATAATTCGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCacac                                                                                             <  2:91800/149‑1 (MQ=255)
                                          gAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTTAGGTTGTCCTTTCGCATTATTTGTTGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACt                                                                                            <  2:90204/149‑1 (MQ=37)
                                          gAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCAACAGGCAATATAAACAACACCAATACAAAAAGGCCAACCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACt                                                                                            >  2:122272/1‑149 (MQ=25)
                                           aaTCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTGCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTCGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTa                                                                                           >  2:147463/1‑149 (MQ=255)
                                            aTCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCcac                                                                                               <  1:223807/144‑1 (MQ=255)
                                                          ggAGTATAAGACGCGCAGCTTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTTTCGCATGGGGAGACCCCACACTACCATCGGCGCTACgg                                                                            <  2:166405/149‑1 (MQ=255)
                                                              tataAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCATTGCAAAAAGGCCATCCGTCAGGATGGACTTTCGCATAATTTGATGCCTGGCAGTTGCCTACTCTCGCATCGGGAGACCGCACACTACTATCGGCGCTCCg                                                                             >  1:6896/1‑144 (MQ=21)
                                                                     cgcgCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTCGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGtgta                                                                      >  1:153101/1‑141 (MQ=255)
                                                                             gTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGAATAGAAAATCTTTACTCTCGCTTGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGa                                                              >  1:79989/1‑144 (MQ=25)
                                                                                          tttttttCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATgggg                                                 <  1:2057/144‑1 (MQ=255)
                                                                                                    cgcCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCACATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGaccccc                                  >  2:8699/1‑149 (MQ=255)
                                                                                                    cgcCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAAATTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGaccacc                                  >  2:140462/1‑149 (MQ=255)
                                                                                                     gcCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACcg                                 >  2:186409/1‑149 (MQ=255)
                                                                                                           gCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTCGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTAc                           <  2:134032/149‑1 (MQ=255)
                                                                                                            cataaaGTCCACCCGTCAGGATGGCCTTTCGCATAATTCGATGCCTGGCTGTTCCCTCCTCTCGCATGGGGAGACCCCACAGTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACcgcg                               <  1:177101/141‑1 (MQ=17)
                                                                                                                               gATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCCTGGGGAGACCCCACACTACCATCGGCGCTCCGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAAt            <  1:217418/144‑1 (MQ=255)
                                                                                                                                    ccTTTCGCATAATTCGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTTAGGTGGGACCATCGTGCTACGGCCGCCAGGCAAATTCTGTTTTAt  <  2:61687/149‑1 (MQ=255)
                                                                                                                                        tCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTa   <  1:109652/144‑1 (MQ=33)
                                                                                                                                                  |                                                                                                                                      
ACGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTAT  >  CP009273/2719280‑2719560

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