Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_vch_tpiA 4,188,312 G→A R978H (CGT→CAT)  rpoC → RNA polymerase, beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_vch_tpiA4,188,3120GA96.8% 79.2 / NA 31R978H (CGT→CAT) rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (20/10);  minor base T (0/1);  total (20/11)
Fisher's exact test for biased strand distribution p-value = 3.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CTGACCATGCGTACGTTCCACATCGGTGGTGCGGCATCTCGTGCGGCTGCTGAATCCAGCATCCAAGTGAAAAACAAAGGTAGCATCAAGCTCAGCAACGTGAAGTCGGTTGTGAACTCCAGCGGTAAACTGGTTATCACTTCCCGTAATACTGAACTGAAACTGATCGACGAATTCGGTCGTACTAAAGAAAGCTACAAAGTACCTTACGGTGCGGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCG  >  NC_000913_3_vch_tpiA/4188167‑4188445
                                                                                                                                                 |                                                                                                                                     
cTGACCCTGCGTACGTTCCAAAGAGGTGGTGCGGCATCTCGTGCGGCTGCTGGATACAGCATACAAGTGAGAATTAAAGGTATCATCAAGCTTATCAACGTGTATTCAGTGATAAACTCCACCGCTACACTGGCTATCACTTCCCATaa                                                                                                                                    >  2:365454/1‑149 (MQ=255)
   aCCATGCGTACGTTCCACATCGGTGGTGCGGCATCTCGTGCGGCTGCTGAATCCAGCATCCAAGTGAAAAACAAAGGTAGCATCAAGCTCAGCAACGTGAAGTCGGTTGTGAACTCCAGCGGTAAAATGGTTATCACTTCCCATAATac                                                                                                                                 >  1:251555/1‑149 (MQ=255)
     cATGCGTACGTTCCACATCGGTGGTGCGGCATCTCGTGCGGCTGCTAAATACAACATCCAAGAGAAAAACAAAGGTAGCATCAAGATCAGCAACGTGAAGTCGGTTGTGAACTCCAGCGGTAAACTGGTTACCACTTCCCATAATactg                                                                                                                               >  1:58273/1‑149 (MQ=255)
     cATGCGGACGTTCCACCTCGGTGGTGCGGCATCTCGTGCAGCTGCTGAATCCAGCATCCAAGGGAAAAACAAAGGTAGCATCAAGCTCAGCAACGTGAAGTCGGTTGTGAGCTCCAGCGGTAAACTGGTTATCACTTCCCATAATactg                                                                                                                               <  2:477254/149‑1 (MQ=255)
             cGTTCCACCTCGTTGGGGCGGCAGCTCGTGCGGCTGCTGAATCCAGCATCCAAGTGAAAAACAAAGGTAGCATCAAGCTCAGCAATGTGAAGTCGGTTGTGAACTCCAGCGGTAAACTGGTTATCACTTCCCATAATACTGaactgaaa                                                                                                                       <  1:240295/149‑1 (MQ=255)
             cGTTCCACATCGGTGGTGCGGCATCTCGTGCGGCTGCTGAATCCAGCATCCAAGTGAAAAACAAAGGTAGCATCAAGCTCAGCAACGTGAAGTCGGTTGTGAACTCCAGCGGTATATTGGTTATCACTTCCCAGAATACTGaactgaaa                                                                                                                       <  1:387725/149‑1 (MQ=255)
               ttCCACATCGGTGGTGCGGCATCTCGTGCGGCTGCTGAATCCAGCAACCAAGTGAAAAACAAAGGTAACATCAAGCTCAGCAACGTGAAGTCGGTTGTGAACTCCAGCGGTAAACTGGTTATCACTTCCCATAATACTGaactgaaact                                                                                                                     >  1:198666/1‑149 (MQ=255)
                     aTCGGTGGTACAGCTGCTCGTGCGGCTGGTGGATACAGCATCAAAGAGAAAAACAAAGGTAGCATCACGCTCAGCAACGTGAAATCGGTGGTGAACTCCAGCGGTAAACTGGTTATCACTACCCATAATGCTGaactgcaaatgcacgc                                                                                                               >  2:144069/1‑144 (MQ=255)
                                ggCATCTCGTTCTGCTGCTGAATCCAGCATCCAAGTGAAAAACAAAGGTAGCATCAAGCTCAGCAACGTGAAGTCGGTTGTGAACTCCAGCGGTAAACTGGTTATCACTTCCCATAATACTGAACTGAAACTGATCGACGAATTCGGTc                                                                                                    <  2:361400/149‑1 (MQ=255)
                                      tcGTGCGGCTGCTGAATCCAGCATCCAAGTGAAAAACAAAGGTAACATAAAGCTCAACAACGTGAAGTCGGTTGTGAACTCCAGCGGTAACCTAGATATCCCGTCCCATAAGACTGAACTGAAACTGATCGACGAATTCGGTCGTACTa                                                                                              >  1:13915/1‑149 (MQ=255)
                                             gctgctGAATCCAGCATCCAAGTGAAAAACAAAGGTAGCATCAAGCTCAGCAACGTAAAGTAGGTTGTGAACTCCAGCGGTAAACTGGTTATCACTTCCCATAATACTGAACTGAAACTAAACGACGAATTCGGCCGTACTaaagaaag                                                                                       >  2:502400/1‑149 (MQ=255)
                                                    aaTCCAGCATCCAAGTGAAAAACAAAGGTAGCATCAAGCTCAGCAACGTGAAGTCGGTTGCGAACTCCAGCGGTACACTGGTTATCACTTCCCATAATACTGAACTGAAACTGATCGACAAATTCGGTCGGACTAAAGAAAGCTaacaa                                                                                >  2:25406/1‑145 (MQ=255)
                                                      tCCAGCATCCAAGTGAAAAACAAAGGTAGCATCAAGCTCAGCAACGTGAAGTCGGTTGTGAACTCCAGCGGTAAACTGGTTATCACTTCCCATAATACTGAACTGAAACTGATCGACGAATTCGGTc                                                                                                    >  1:484242/1‑127 (MQ=255)
                                                      tCCAGCATCCAAGTGAAAAACAAAGGTAGCATCAAGCTCAGCAACGTGAAGTCGGTTGTGAACTCCAGCGGTAAACTGGTTATCACTTCCCATAATACTGAACTGAAACTGATCGACGAATTCGGTc                                                                                                    <  2:484242/127‑1 (MQ=255)
                                                        cAGCATCCAAGTGAAAAACAAAGGTAGCATCAATCTCAGCAACGTGAAGTCGGTTGTGAACTTCAGCGGTAAACTGGTTATCACTTTCCTTTATTCTTATCTTTTTTTGATTGATGAATTCGGTTGTACTAAAGAAAGCTACAAAGTAc                                                                            <  1:247729/149‑1 (MQ=255)
                                                        cAGCATCCAAGTGAAAAACAAAGGTAGCATCAAGCTCAGCAACGTGAAGTCGGTTGTGAACTCCAGCGGTAAACTGGTTATCACTTCCCATAATACTGAACTGAAACTGATCGACGAATTCGGTCGTACTAAAGAAAGCTACAAAGTAc                                                                            >  1:583070/1‑149 (MQ=255)
                                                                    gataaaCAAAGGTAGCATCAAGCTCAGCAACGTGAAGTCGGTTGTGAACTCCAGCGGTAAACTGGTTATCACTTCCCATAATACTGAACTGAAACTGATCGACGAATTCGGTCGTACTAAAGAAAGCTAAAAAGTACCTTACGGTGCgg                                                                >  2:380610/4‑149 (MQ=255)
                                                                    gAAAAACAAAGGTAGCATCAAGCTCAGCAACGTGAAGTCGGTTGTGAACTCCAGCGGTAAACTGGTTATCACTTCCCATAATACTGAACTGAAACTGATCGACGAATTCGGTCGTACTACAAACAGCTACAAAGTACATTACGGTGCgg                                                                >  1:253572/1‑149 (MQ=255)
                                                                           aaaGGTAGCATCAAGCTCAGCAAAGAGAAGTCGGTTGTGAACTACAGCGGTAAACTGGTTATAACTACCCATAATACTGAACTGAAACTGATAGACGAATTCGGACGTACTAAAGAAAGCTACACAATACATTACGTTGCGGTACTGGc                                                         >  2:124002/1‑149 (MQ=255)
                                                                               gTAGCATCAAGCTCAGCAACGTGACGTCGGTTGTGAACTCCAGCGGTAAACTGGTTATCACTTCCCATAATACTGAACTGAAACTGATCGACGAATTCGGTCGTACTAAAGAAAGCTACAAAGTACCTTACGGTGCGGTACTGGCGaaa                                                     >  1:488726/1‑149 (MQ=255)
                                                                               gTAGCATCAAGCTCAGCAACGTGAAGTCGGTTGTGAACTCCAGCGGTAAACTAGTTATCACTTCCCATAATACTGAACTGAAACTGATCGACGAATTCGGTCGTACTAAAGAAAGCTACAAAGTACCTTACGGTGCGGTACTGGCGaaa                                                     >  1:486360/1‑149 (MQ=255)
                                                                                       aaGCTCAGCAACGTGAAGTCGGTTGTGAACTCCAGCGGTAAACTGGTTATCCCTTCCTATAATACTGAACTGAAACTGATCGACGAATTCGTTCGTACTAAAGAAAGCTTTAAAGTACCTTACGGTGCGGTACTGGCGAAAGGCGATgg                                             <  1:170010/149‑1 (MQ=255)
                                                                                           tCAGCAACGTGAAGTCGGTTGTGAACTCCAGCGGTAAACTGGTTATCACTTCCCATAATACTGAACTGAAACTGATCGACGAATTCGGTCGTACTAAAGAAAGCTACAAAGTACCTTACGGTGCGGTACTGGCGAAAGGCGATGGCGaa                                         >  1:304271/1‑149 (MQ=255)
                                                                                                      aaGTGGGTTGTGAACATCAGTGGTAAACGGGTTATCAGTTCCCAGAATACTGAACTGAAACTGATCTACGAATTCTGTCTTACTAAAGAAAGCTACAAAGTACCTTACGGTGCGGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTgg                              <  2:13915/149‑1 (MQ=255)
                                                                                                               gtgAACTCCAGCGGTAAACTGGTTATCACTTCCCATAATACTGAACTGAAACTGATCGACGAATTCGGTCGTACTAAAGAAAGCTACAAAGTACCTTACGGTGCGGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAAc                     >  1:61810/1‑149 (MQ=255)
                                                                                                                  aaCTCCAGCGGTAAACTGGTTATCACTTCCCATAATACTGAACTGAAACTGATCGACGAATTCGGTCGTACTAAAGAAAGCTACAAAGTACCTTACGGTGCGGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGt                  <  1:581534/149‑1 (MQ=255)
                                                                                                                     tCCAGCGGTAAACTGGTTATCACTTCCCATAATACTGAACTGAAACTGATCGACGAATTCGGTCGTACTAAAGAAAGCTACAAAGTACCTTACGGTGCGGTACTGGCGAAAGGCGATGGCGAACAGGTTGGTGGCGAGGAAACCGTTGc               >  2:395653/1‑149 (MQ=255)
                                                                                                                           tgtcaccTGTTTTTCACTTCCCATAATACTGAACTGAAACTGATCGACGAATTCGGTCGTACTAAAGAAAGCTACAAAGTACCTTACGGTGCGGTACTTGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTg         <  1:170961/143‑1 (MQ=255)
                                                                                                                                 cTGGTTATCACTTCCCATAATACTGAACTGAAACTGATCGACGAATTCGGTCGTACTAAAGAAAGCTACAAAGTACCTTACg                                                                      <  2:336738/82‑1 (MQ=255)
                                                                                                                                 cTGGTTATCACTTCCCATAATACTGAACTGAAACTGATCGACGAATTCGGTCGTACTAAAGAAAGCTACAAAGTACCTTACg                                                                      >  1:336738/1‑82 (MQ=255)
                                                                                                                                  tGGTTATCACTTCCCATAATACTGAACTGAAACTGATCGACGAATTCGGTGGTACTCAACAAAGCATTAAAGTATCTAACAGTCCGGGACTGCCCAAAGACGATGCCGAAAAGGGTGCTGACGGCGAACCCCTTCCAACCggggcgcac  >  2:430740/1‑144 (MQ=255)
                                                                                                                                                 |                                                                                                                                     
CTGACCATGCGTACGTTCCACATCGGTGGTGCGGCATCTCGTGCGGCTGCTGAATCCAGCATCCAAGTGAAAAACAAAGGTAGCATCAAGCTCAGCAACGTGAAGTCGGTTGTGAACTCCAGCGGTAAACTGGTTATCACTTCCCGTAATACTGAACTGAAACTGATCGACGAATTCGGTCGTACTAAAGAAAGCTACAAAGTACCTTACGGTGCGGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCG  >  NC_000913_3_vch_tpiA/4188167‑4188445

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