Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,186,578 A→T 100% D410V (GAT→GTT)  rpoC → RNA polymerase subunit beta'

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,186,5780AT90.6% 50.5 / ‑2.7 21D410V (GAT→GTT) rpoCRNA polymerase subunit beta'
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base T (9/9);  total (11/10)
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.89e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TACCTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCT  >  NC_000913/4186427‑4186764
                                                                                                                                                       |                                                                                                                                                                                          
tACCTGCGTCTGCATCAGAGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGTGTGGTCTTGCTACCACCATTAAAACTGCGAAGAAAATGGTTGAGCCCGACGAAGCTGTCGTTTGGGTTATGCTGGACGAAGTTATCCGCGAacac                                                                                                                                                                >  1:135037/1‑180 (MQ=255)
      cGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGTTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCAcc                                                                                                                                         >  2:90310/1‑197 (MQ=255)
      cGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGTTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAgg                                                                                                               >  1:366468/1‑223 (MQ=255)
      cGTCTGCAGCAGTGCGGTCTGCCGAAGAAAATGGCACTGGGGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGTTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAgg                                                                                                               <  2:366467/223‑1 (MQ=255)
         cTGCATCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGGCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGTTATCCTGGACGAAGTTATCCGCGAACAcc                                                                                                                                                               <  1:245324/172‑1 (MQ=255)
          tGCATCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGTTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGAct                                                                                                                                     <  1:90310/197‑1 (MQ=255)
                    cGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACAGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGACGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCAGGTACTGCTGAACCGTGCACCGACTCTGCACAGTCCGGGTACCCAGGCATTTGAACCGGTACTGATCGATGGTAAAGCGATCGAGCTGCACCCGCGGGGTTTTGTgggggaaa                                             >  2:141904/1‑269 (MQ=255)
                                                   ttCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGTTATCCTGGACGAAGTTATCCGCGAAAACCCGGTACTGCTGAACCGt                                                                                                                                              >  1:179587/1‑147 (MQ=255)
                                                       aaCCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGTTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCAc                                                                                                                                          <  2:179586/147‑1 (MQ=255)
                                                             tCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGTTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATaa                                            >  2:117969/1‑235 (MQ=255)
                                                               acctccGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCt  <  1:141905/270‑1 (MQ=255)
                                                                 cTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGTGATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGcc                                        <  1:117969/235‑1 (MQ=255)
                                                                                         tCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGTTATCCTGGAAGAAGTTATCCGCGAACACCCGGTACTGCCGAACCGTGCCCCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGCTCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATaa                                            >  1:357072/1‑207 (MQ=255)
                                                                                         tCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTaaa                                                                                  >  1:245485/1‑169 (MQ=255)
                                                                                             gCTACCACCATTAAAGCAGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGTTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTa                                                                              <  2:245484/169‑1 (MQ=255)
                                                                                                  caccaTTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGTTATCCTGGACGAAGTTATCCGCGAACAcc                                                                                                                                                               <  1:179071/83‑1 (MQ=255)
                                                                                                   accaTTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGTTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGt                                                                                                                                              >  2:317116/1‑99 (MQ=255)
                                                                                                       ttAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGTTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCAc                                                                                                                                          <  1:317117/99‑1 (MQ=255)
                                                                                                               gCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGTTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCg                                                                                                    >  2:123982/1‑129 (MQ=255)
                                                                                                               gCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGTTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTtg                                                        >  1:244222/1‑173 (MQ=255)
                                                                                                                   agaaAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGTTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTAc                                                                                                <  1:123982/129‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                                                          
TACCTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCT  >  NC_000913/4186427‑4186764

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