Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_vch_pgi 4,183,379 C→T S712F (TCC→TTC)  rpoB → RNA polymerase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_vch_pgi4,183,3790CT97.4% 114.1 / NA 39S712F (TCC→TTC) rpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (19/19);  minor base A (1/0);  total (20/19)
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.88e-01

GTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGTTCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGAAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTT  >  NC_000913_3_vch_pgi/4183239‑4183519
                                                                                                                                            |                                                                                                                                            
gTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTTCGGTGTaa                                                                                                                                      >  2:160092/1‑149 (MQ=255)
   cccgATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGACCGTGCTGTTGCCGTTGACTTCGGTGTAACTg                                                                                                                                   <  1:16721/146‑1 (MQ=255)
    cTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTTCGGTGTAACTGc                                                                                                                                  <  1:298030/149‑1 (MQ=255)
        tCCCGTTCCGGGAACACGATGACGCCAACCGGGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTTCGGTGTAACTGCGGTa                                                                                                                              <  2:134531/149‑1 (MQ=255)
         cccGTTCCTGGAACACGATGACGCCAACCGTTAATATATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTTCGGTGTAACTGCGGTAg                                                                                                                             >  1:189975/1‑149 (MQ=255)
          ccGTTCCTGGAACACGATGACGCCAACCGTGCATTGCTGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTTCGGTGTACCGGCGGTAGc                                                                                                                            >  1:131971/1‑149 (MQ=255)
                cTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTTCGGTGTAACTGCGGTAGCTAACCg                                                                                                                      >  1:1107/1‑149 (MQ=255)
                 tGGAACACGATGACGCCAACCGCGCATTCAGCAGCGCGACCCTGCCACGACAGGCCGTTCCGACTCTTCGCGCTGATAAGCCCCTGGTTTGTACCGGTATTTAACGCGCTGTTGCCGTTCGCTACAGTGAAACTGCCGTAGCCAAAcat                                                                                                                     >  2:259173/1‑147 (MQ=255)
                     acacGATAACGCCAACCGTGCATTGATGGGTCGGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTTCGGTGTAACTGCGGTAGCTAAACgtggtg                                                                                                                 <  2:117087/149‑1 (MQ=255)
                       acGGTGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGGCTCTGCGGGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTTCGGTGTAACTGCGGTAGCTAAACggggtggt                                                                                                               <  2:367788/149‑1 (MQ=255)
                             aCGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTTCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGtt                                                                                                         <  2:346310/149‑1 (MQ=255)
                               gCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTTCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGTTCa                                                                                                       <  1:240803/149‑1 (MQ=255)
                                         cATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTTCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGTTCAGTAcatacga                                                                                             <  2:339760/149‑7 (MQ=255)
                                         cATTGATGGGTGCGAACAGGGAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAATGTGCTGTTGCCGTTGACTTCGGTGTAACTGCGGTAGCTAACAGTAGTAGTCTCCTTCCGTAACTGGAt                                                                                             <  2:85290/149‑1 (MQ=255)
                                            tGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTTCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGTTCAGTACGTGGATGCt                                                                                          >  1:27944/1‑149 (MQ=255)
                                             gATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTTCGGTGTAACTGCGGTAGCTAAACGTGTTGGTTTCGTTCAGTACGTGTATGCtt                                                                                         >  2:187622/1‑149 (MQ=255)
                                                gggTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTTCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGTTCAGTACGTGGATGCTTccc                                                                                      >  2:190685/1‑149 (MQ=255)
                                                      gAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGAATGGAACGTGCTGTTGCCGTTGACTTCAGTGTAACTGCGGTACCTAAAAGTGGCGGTGTCGCGCAGCACGTGGATGCTGCCCGTATCg                                                                                >  2:298140/1‑149 (MQ=255)
                                                             cAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTTCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGTTCAGTACGTGGATGCTTCCCGTATCGTAATCaa                                                                         >  2:192618/1‑149 (MQ=255)
                                                              aaCGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTTCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGTTCAGTACGTGGATGCTTCCCGTATCGTTATCaaa                                                                        >  1:218386/1‑149 (MQ=255)
                                                              aaCGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTTCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGTTCAGTACGTGGATGCTTCCCGTATCGTTATCaaa                                                                        >  1:70930/1‑149 (MQ=255)
                                                                                 tctGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTTCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGTCCAGTACGTGGATGCTTACCGTATCGTTATCAACGCTAACCAAGCCGAGATGt                                                     >  2:168839/1‑149 (MQ=255)
                                                                                          tGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTTCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGTTCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTg                                            >  1:43467/1‑149 (MQ=255)
                                                                                            aTAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTTCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGTTCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGaa                                          <  1:190685/149‑1 (MQ=255)
                                                                                             tAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTTCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGTTCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGAAg                                         >  1:344160/1‑149 (MQ=255)
                                                                                              aaGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTTCGGTGTAACTGCGGTAGCGAAACGTGGTGGTGTCGTTCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGAAGc                                        >  2:195064/1‑149 (MQ=255)
                                                                                                gCCGCTGGTTGGTCCTGGTATGGAACGTGCTGTTGCCGTTGACTTCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGTTCAGTACTTGTTTCTTTGTGACAACCATAGCAAAGTTAACGAAGACGAGTTGTATCCGGGTGAAGCAg                                      <  1:286810/149‑1 (MQ=255)
                                                                                                       gttggtACTGGTATGGAACGGGCTGTTGCCGTTGACTTCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGTTCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGAAGCAGGTATCGa                               <  2:197133/149‑1 (MQ=255)
                                                                                                          ggtACTGGTATGGAACGTGCTGTTGCCGTTGACTTCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGTTCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGAAGCAGGTATCGACAt                            <  1:160092/149‑1 (MQ=255)
                                                                                                          gggaCTGGTATGGAACGTGCTGTTGCCGTTGACTTCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGTTCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGAAGCAGGTATCGACAt                            <  1:4172/146‑1 (MQ=255)
                                                                                                            tACTGGTATGGAACGTGCTGTTGCCGTTGACTTCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGTTCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGAAGCAGGTATCGACATCt                          <  2:86548/149‑1 (MQ=255)
                                                                                                            tACTGGTATGGAACGTGCTGTTGCCGTTGACTTCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGTTCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGAAGCAGGTATCGACATCt                          >  1:86548/1‑149 (MQ=255)
                                                                                                              cgggTATGGAACGTGCTGTTGCCGTGTACTTCGGTGTAACTGCGGCAGCTAAACGTGGTGGTGTCGTTCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGGAGCAGGTATCGACATCTAc                        <  2:327757/147‑1 (MQ=255)
                                                                                                                 gTATGGAACGTGCTGTTGTCGTTGACTTCGGGGTAACTGCGGTAGCTAAACGTGGTGGTGTCGTTCAGTACGTGGATGCTTCCCTTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGAAGCAGGTATCGACATCTACAAc                     <  2:76816/149‑1 (MQ=255)
                                                                                                                 gTATGGAACGTGCTGTTGCCGTTGACTTCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGTTCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGAAGCAGGTATCGACATATAAAAc                     >  2:226767/1‑149 (MQ=255)
                                                                                                                     gtaaCGTGCTGTTGCCGTTGACTTCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGTTCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGAAGCAGGTATCGACATCTACAACCTGa                 <  2:242501/147‑1 (MQ=255)
                                                                                                                      gAACGTGCTGTTGCCGTTGACTTCGGCATAACAGCAGTGGCTAAAAGTGGTGGTGTCGTTAAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACTAAGACGAGATGTATCCGGATGCAGCAGGTATCGACATCTACAACCTGAc                >  1:132887/1‑149 (MQ=255)
                                                                                                                                 tggcGTTGACTTCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGTTCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGAAGCAGGTATCGACATCTACAACCTGACCAAATACAccc     <  2:364504/146‑1 (MQ=255)
                                                                                                                                    cGTTGCCTTCGGTGTAACTACGGTAGCTAAACGTGGTGGTGTCGTTCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGTTGTATCCGGGTGAAGCAGGTATCGACATCTACAACCTGACTAAATACACCCGtt  >  2:226384/1‑149 (MQ=255)
                                                                                                                                            |                                                                                                                                            
GTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGTTCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGAAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTT  >  NC_000913_3_vch_pgi/4183239‑4183519

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: ‑