Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,175,524 G→T 61.5% V83F (GTC→TTC)  rpoC → RNA polymerase, beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,175,5240GT61.5% 11.7 / 12.9 23V83F (GTC→TTC) rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base G (4/5);  new base T (6/8);  total (10/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 = 9.46e-01

TTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCTG  >  CP009273/4175380‑4175655
                                                                                                                                                |                                                                                                                                   
ttCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAGGCGCCTGAAACACCGTGGCTTCAt                                                                                                                                 <  2:202202/149‑1 (MQ=255)
              aaaGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGTGGTAAGTACAAGCGCCTGAAACACCGTGGCTTCATCTGTGAGAAGTgcg                                                                                                                   <  2:7640/149‑1 (MQ=255)
                    ggAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCTTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCTTCATCTGTGAGAAGGgcgg                                                                                                                  >  1:650/1‑144 (MQ=255)
                         ccATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTg                                                                                                             >  1:233195/1‑144 (MQ=255)
                          cATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTACGCCCGAAACTTGAAGCCGGGAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCTTCATCTGTGAGAAGTGCGGCGTTGAAGTGa                                                                                                       >  2:136774/1‑149 (MQ=255)
                                cTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTATAAGCGCCTGAAACATCGTGGCTTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGa                                                                                                 <  2:85008/149‑1 (MQ=255)
                                   ccGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCTTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCa                                                                                                   >  1:103871/1‑144 (MQ=255)
                                           tCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCTTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAg                                                                                           <  1:220207/144‑1 (MQ=255)
                                              aaCCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGt                                                                                   <  2:268200/149‑1 (MQ=255)
                                                         gACGGCCTTTTCGGCGCCCGTGTCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCTTCATCTGTGAGAATTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATggg                                                                        <  2:27007/149‑1 (MQ=255)
                                                                    ccgcgcCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCTTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACAt                                                                  <  1:123589/142‑1 (MQ=255)
                                                                           cGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGTCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTCAGCGTATGGGCCACATCGAACTGGCTTc                                                      <  2:202677/149‑1 (MQ=255)
                                                                              aTCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCTTCATCTGTGAGAAGTGCGGCGTTGAAGTGTCCCAGTCTAAAGTACGCCGTGAGCGTATGGGCCACAGCGAACTGGCTTcccc                                                   <  2:801/149‑1 (MQ=255)
                                                                                    gggCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCTTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCg                                                  <  1:162534/144‑1 (MQ=255)
                                                                                           tAAAAGATTACGAGTGCCTGGGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTgcgc                                           <  1:110783/144‑1 (MQ=255)
                                                                                                 aTTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGACGTCATCTGTGAGAAGTGCGGCGTTGAAGTAACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGCTTc                                >  2:8998/1‑149 (MQ=255)
                                                                                                       aGTGCCTGTGCGGGTAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGaaa                          <  2:127903/149‑1 (MQ=255)
                                                                                                        gTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGACTTCATCTGTGAGAAGTGCGGCATTGAAGTGACACAGACTAAAGTACGCCATGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTACTGAAAt                         >  2:119054/1‑149 (MQ=255)
                                                                                                            cTGTGCGGTAAGTACAAGCGCCTGAAACACAGTGGCTTCATCTGTGAGAAGTGCGGCGTTGAAGTGACGCAGACGAAAGTACGCCGTGAGCGTATGGGCCAAATCGACCTGGCTTCCACGACTGCGACAATCAGGTTGCTaaaa                          >  1:50633/1‑144 (MQ=255)
                                                                                                             tgtgCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATAAGTGAGAAGTGAAACGTTAAAGTGAACCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTg                    >  2:99689/1‑149 (MQ=255)
                                                                                                                             gcgcCTGAAACACCGTGGCTTCCTCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTc    >  2:76243/1‑149 (MQ=255)
                                                                                                                               gcCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGTTTCCCCGACTGCGTACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTctg  <  2:255232/149‑1 (MQ=255)
                                                                                                                                   gAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTct   >  1:113546/1‑144 (MQ=255)
                                                                                                                                                |                                                                                                                                   
TTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCTG  >  CP009273/4175380‑4175655

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: