Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_vch_pgi 4,187,214 A→C D622A (GAC→GCC)  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_pgi4,187,2140AC97.2% 103.1 / NA 36D622A (GAC→GCC) rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base C (20/15);  minor base G (0/1);  total (20/16)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGGCCGTGCCATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAG  >  NC_000913_3_vch_pgi/4187070‑4187347
                                                                                                                                                |                                                                                                                                     
ttGGCCGTGCCATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAg                                                                                                                                   >  1:106637/1‑149 (MQ=255)
       tGCCATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTGTTTTTGCGGCCCAGATCATGt                                                                                                                            >  1:162639/1‑149 (MQ=255)
         ccATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGGAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTAACGCATTCTCGGTCTGAAACCGTCGGTTATTTTTGCTGCCCACATCagggac                                                                                                                          >  1:212845/1‑144 (MQ=255)
           agtCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCTAGGGCATGTacac                                                                                                                        <  2:344566/147‑1 (MQ=255)
           aTTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTacac                                                                                                                        >  1:459786/1‑149 (MQ=255)
            ttCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACAcc                                                                                                                       >  2:95949/1‑149 (MQ=255)
             tCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTATACCg                                                                                                                      <  2:263588/149‑1 (MQ=255)
                 tgGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCtt                                                                                                                  >  1:384362/1‑149 (MQ=255)
                 tgGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCtt                                                                                                                  >  1:534318/1‑149 (MQ=255)
                      gatTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGCTCATGTACACCGGCTTCGCCt                                                                                                             >  2:157580/1‑149 (MQ=255)
                                                 cATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTg                                                                                  >  1:203329/1‑149 (MQ=255)
                                                 cATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTg                                                                                  <  1:499636/149‑1 (MQ=255)
                                                  aTCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTgg                                                                                 >  2:408313/1‑149 (MQ=255)
                                                   tCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGt                                                                                <  1:370159/149‑1 (MQ=255)
                                                         aCCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGAt                                                                          <  1:228267/149‑1 (MQ=255)
                                                               cgcTGGGTAAACAAGCACTCTCCAAAATGCTGAACACCTGCCACCGCATTCTCTGTCTGAAACCGACCGTTATTCTTGCGGCCCAGTTCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATg                                                                    <  1:187840/149‑1 (MQ=255)
                                                               cgcTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATg                                                                    <  2:384362/149‑1 (MQ=255)
                                                                    ggTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCAt                                                               >  1:330117/1‑149 (MQ=255)
                                                                         aaaaGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGACCATGTAAACCGGCTTCGCCTCTGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACCTGGTTATCCCgg                                                          >  2:80816/1‑149 (MQ=255)
                                                                                    ccAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAacac                                               >  2:71519/1‑149 (MQ=255)
                                                                                      aaaaTGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGa                                             >  1:387408/1‑149 (MQ=255)
                                                                                       aaaTGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTTCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACAAGaa                                            >  2:62355/1‑149 (MQ=255)
                                                                                             tGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCTAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTGATCTCGGTGCGGAAACACGAGatcatt                                      >  2:245101/1‑148 (MQ=255)
                                                                                              gAACACCGGCTCCCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCAtct                                     <  2:74947/149‑1 (MQ=255)
                                                                                                  acCTGCTGCCTCATTCTCCGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGa                                 <  2:366005/149‑1 (MQ=255)
                                                                                                         aCCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGgcagaa                          <  2:203329/149‑1 (MQ=255)
                                                                                                             cATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCCTCCCGGAGAAGAAACACGAAATCATCTCCGCGgcagaagcag                      >  2:497186/1‑149 (MQ=255)
                                                                                                                         gAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGaaa          <  2:459786/149‑1 (MQ=255)
                                                                                                                             ccgaccgTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAAATCa      >  1:387459/1‑149 (MQ=255)
                                                                                                                              ctaacgTTATTTTTGCGGGCCAGATCATGTACAACGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAg     <  2:25970/145‑1 (MQ=255)
                                                                                                                              cgaccgTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAg     <  2:373649/149‑1 (MQ=255)
                                                                                                                               gaccgTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAgg    >  2:35897/1‑149 (MQ=255)
                                                                                                                                accgGTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGa   <  1:508987/149‑1 (MQ=255)
                                                                                                                                 ccgTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGGAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAg  <  1:71519/149‑1 (MQ=255)
                                                                                                                                         tttGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCAtctc                                    <  2:174469/107‑1 (MQ=255)
                                                                                                                                         tttGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCAtctc                                    >  1:174469/1‑107 (MQ=255)
                                                                                                                                                |                                                                                                                                     
TTGGCCGTGCCATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAG  >  NC_000913_3_vch_pgi/4187070‑4187347

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