Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,184,543 C→A P1100Q (CCG→CAG)  rpoB → RNA polymerase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,184,5430CA89.3% 70.0 / NA 29P1100Q (CCG→CAG) rpoBRNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base A (12/13);  minor base G (0/3);  total (13/16)
Fisher's exact test for biased strand distribution p-value = 2.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCA  >  NC_000913/4184403‑4184681
                                                                                                                                             |                                                                                                                                          
tCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGc                                                                                                                                          >  1:745456/1‑143 (MQ=255)
   ggCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTgg                                                                                                                                       <  1:99119/143‑1 (MQ=255)
           aaCGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCATACTGAACCAGATGGGCGTACCg                                                                                                                               >  2:708226/1‑143 (MQ=255)
            aCGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTCGCCACAGTAGTGACACAGCGGGGCGGACCGt                                                                                                                              >  2:364333/1‑143 (MQ=255)
            aCGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGt                                                                                                                              >  1:307590/1‑143 (MQ=255)
             cGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTAAGGCGGTCGACATCGTACTGAACCAGCTGTGCGTAcagga                                                                                                                             >  2:656224/1‑139 (MQ=255)
               ccGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGactc                                                                                                                           <  1:442106/143‑1 (MQ=255)
                cGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAGAAAGGCACGCCGGTAGACCTCGTACTGAACCAGCTGGGCGTACCGTCTCg                                                                                                                          >  1:476472/1‑143 (MQ=255)
                  gagccAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTa                                                                                                                        <  1:197376/140‑1 (MQ=255)
                      cAGCCTGGGGGCAAGTTGGCAGGTCGTAACGGTAACAAGGGTGGAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGaa                                                                                                                    <  1:362088/143‑1 (MQ=255)
                             gggagaaGGTGGCAGGTCGTGAGGGGAACAAGGGTGGAATGTCTAAGATCAACCCGATGGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACAGCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGt                                                                                                             <  2:660483/138‑1 (MQ=255)
                                            gtcgtcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAAcc                                                                                              >  2:688292/1‑143 (MQ=255)
                                            gtcgtcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGAGCCTGGAAAcc                                                                                              >  1:208148/1‑143 (MQ=255)
                                                      tAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTa                                                                                    <  1:346009/143‑1 (MQ=255)
                                                                            aTCCACCTGATCGAAGATATGCCTTCCGGTGGAAGCCGTACGCCGGTAGACATCGTACTGACCCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGa                                                              <  1:616161/143‑1 (MQ=255)
                                                                            aTCCACCCGATCGAAGATATGCCTTACGATGAAAACGGTCCGCCGGTTGCTTTATGACCGACACAGAATTTTTGTTCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGa                                                              <  2:339481/143‑1 (MQ=255)
                                                                               ggcccGATCGAAGTTGTGTGTTACGAGGGATACGGGAGGCTGGTAGACATCGTACTTAGCCAGCAGGGCGAACCGTCTCGTATGAACATCGGTCAGATCCTAGAAACCCACCTGGGTATGGGTGCGAAAGGTATCGGCGACaa                                                           <  2:328973/141‑1 (MQ=255)
                                                                                 cccGATCGAAGATATGCCTGACGATGAAAACGGTACGCCATTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGa                                                         <  1:733153/143‑1 (MQ=255)
                                                                                       cGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGCCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGACCGGTATCGGCGACCAGATCAACg                                                   >  2:461158/1‑143 (MQ=255)
                                                                                            atatGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTAGCGTCTCGTATGAGCATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATg                                              <  2:265974/143‑1 (MQ=255)
                                                                                             tatGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGACCATCGGTCAGATCCTCGAAACCCACCTGGGTCTGGCTGCGCAAGGTATCGGCGACAAGATCAGCGccgcgc                                             >  1:53514/1‑139 (MQ=255)
                                                                                             tacgCCTTGCGATGAAAACGGTACGCCGGTAGACATCGTACTGAGCCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGc                                             <  1:213877/140‑1 (MQ=255)
                                                                                               tGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTg                                           <  1:119066/143‑1 (MQ=255)
                                                                                                tgcTTTCGATGAAACCGGTACGCCGGTAGACATCGTCCTGACCCGGCTGGGCCTACCGTCTGGTATGAGCATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGa                                          <  2:47308/141‑1 (MQ=255)
                                                                                                            aaaCGGTACGCCCGTAGACATCGTACTGAACCAGCTGGGCTTACCGTCTCGTATGAACATCGGTCAGATCCTCGAATCCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCaag                              <  2:255355/143‑1 (MQ=255)
                                                                                                              cccgcaCGCCGGTAGACGCCGTACTGAACCGGCCGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCaagaa                            <  1:242301/138‑1 (MQ=255)
                                                                                                                                tCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATGGGTGAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCCCGaa          >  1:110981/1‑143 (MQ=255)
                                                                                                                                        aaCCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCa  >  1:572032/1‑143 (MQ=255)
                                                                                                                                        aaCCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATGCATccc  >  1:654522/1‑142 (MQ=255)
                                                                                                                                             |                                                                                                                                          
TCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCA  >  NC_000913/4184403‑4184681

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