Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,184,543 C→A 100% 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,5430CA97.9% 163.1 / ‑4.8 54P1100Q (CCG→CAG) rpoBRNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base A (30/22);  total (31/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCAT  >  NC_000913/4184403‑4184678
                                                                                                                                            |                                                                                                                                       
tCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGc                                                                                                                                       >  2:648452/1‑143 (MQ=255)
  tGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTg                                                                                                                                     >  2:787896/1‑143 (MQ=255)
     cGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGc                                                                                                                                  <  1:383666/143‑1 (MQ=255)
      ggTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCg                                                                                                                                 >  1:301677/1‑143 (MQ=255)
          aaaCGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGATCCAGCTGGGCGTAcc                                                                                                                             <  2:934521/143‑1 (MQ=255)
               ccGTATCCAGAGTGGTGACAAGAAGGCAGGGGGGCACGGAAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGAAATGCCTTTCGTTGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGtctc                                                                                                                        <  1:43190/143‑1 (MQ=255)
                 gTATCTAGCCTGGTGACAAGATGGCAGGGCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGt                                                                                                                      <  2:978723/143‑1 (MQ=255)
                       aGCCTGGTGACAAGATGGCAGGTCGTCACGGTAAGAAGGGTGTAATTTCTGAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCCGGGCGTACCGTGTCAGATGAAc                                                                                                                >  2:566424/1‑143 (MQ=255)
                             gTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGt                                                                                                          <  1:786150/143‑1 (MQ=255)
                                   aGATGGCAGGTCGTCACGGTAACAAGGGGGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATc                                                                                                    >  1:868213/1‑143 (MQ=255)
                                   aGATGGCAGGTCGCCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATc                                                                                                    <  1:382573/143‑1 (MQ=255)
                                       ggCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCg                                                                                                >  1:978282/1‑143 (MQ=255)
                                       ggCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCg                                                                                                >  1:693154/1‑143 (MQ=255)
                                          cggTCGTCACGGTAACAACGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGaaa                                                                                             <  1:175824/142‑1 (MQ=255)
                                          aGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGaaa                                                                                             >  1:324574/1‑143 (MQ=255)
                                               gtcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGTGCGTACCGTCTCGTATGAACATCGGTCAGCTCCTCGAAACCCAc                                                                                        >  2:96739/1‑143 (MQ=255)
                                                    ggTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGCCCCCGCGTTGTTTTTCGGCTCGTATGAACATCGGTCAGATCCTCGAAACCCAACTggg                                                                                   >  1:509485/1‑143 (MQ=255)
                                                    ggTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTggg                                                                                   >  2:387944/1‑143 (MQ=255)
                                                    ggTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCGCGAAACCCACCTggg                                                                                   >  2:267779/1‑143 (MQ=255)
                                                        aCAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATg                                                                               >  1:230568/1‑143 (MQ=255)
                                                          aaGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATCAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGc                                                                             >  2:462829/1‑143 (MQ=255)
                                                           aGGGTTTAAGGTCTAAGACCAGCCCGATCGAAGCTATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCt                                                                            <  2:1082111/143‑1 (MQ=255)
                                                            gggTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTg                                                                           <  1:648452/143‑1 (MQ=255)
                                                             gggttAATTTCTGAGTTCACCCCGATCGAAGGGATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCACCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGc                                                                          <  2:315805/139‑1 (MQ=255)
                                                              gtgtAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCCCGAAACCCACCCGGGGATGgcggcg                                                                         >  2:2941/1‑143 (MQ=255)
                                                                  aactacgaaGATCAACCCGATCGAAGATATGCCTTACGATGAAAGCGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGAGCCTCGAAACCCACCTGGGTATGGCTGCGAAAg                                                                     <  2:944679/136‑1 (MQ=255)
                                                                  aaattCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTNCGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAg                                                                     <  1:558151/140‑1 (MQ=255)
                                                                  aaTTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTCTGAACATCGGTCAGATCCTCGCTACCCACCTGGGTTTGGCTGCGAAAg                                                                     >  1:694334/1‑143 (MQ=255)
                                                                   aTTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGAGCGTACCGTCTCGTATGAACATCGGTCAGATCCGCGAAACCCACCCGGGGATGGCTGCGAAAgg                                                                    >  2:891126/1‑143 (MQ=255)
                                                                     ttCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTa                                                                  >  1:520694/1‑143 (MQ=255)
                                                                      tCTAAGATCACCCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTAt                                                                 <  1:785724/143‑1 (MQ=255)
                                                                      tCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTCCCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTAGGGCTGCGAAAGGTAt                                                                 >  1:556586/1‑143 (MQ=255)
                                                                         aaGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGGATCgg                                                              >  2:214865/1‑143 (MQ=255)
                                                                            aTCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGa                                                           <  1:620615/143‑1 (MQ=255)
                                                                               aaCCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTCTCGGCGACag                                                        >  1:241806/1‑142 (MQ=255)
                                                                               aaCCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACaa                                                        <  1:536434/143‑1 (MQ=255)
                                                                                  ccGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAgct                                                     >  2:196133/1‑141 (MQ=255)
                                                                                         aaGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACACACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGcc                                              >  2:722011/1‑143 (MQ=255)
                                                                                                ggcTTATGATGAAAGCGCTACGCCGGTAGACATCGTACTGAACCCGGTGGGCGTACCGTCTCGTATGAACATTGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAATGGTATCGGCGACAAGATCAACGCCATGCTGa                                       <  2:296359/141‑1 (MQ=255)
                                                                                                    tACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAAca                                   <  1:624486/143‑1 (MQ=255)
                                                                                                       gATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTCTGGCTGCGAAAGGTATCGGCGACAAGATCACCGCCATGCTGAAAcagca                                >  2:259295/1‑143 (MQ=255)
                                                                                                        aTGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAAcagccg                               >  2:618278/1‑141 (MQ=255)
                                                                                                         tGAAAACGGTACGCCGGTAGACATCGTACCGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGACAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGACCAACCCCATGCTGACAaagcaga                              >  2:335860/1‑142 (MQ=255)
                                                                                                         cgAATACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCTTACCGGCTCGTATGAACCTCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACagcagc                              <  2:644615/142‑1 (MQ=255)
                                                                                                                ggTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGt                       <  1:787896/143‑1 (MQ=255)
                                                                                                                  tACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCGATGCTGAAACAGCAGCAAGAAGTcg                     >  1:720869/1‑143 (MQ=255)
                                                                                                                   tctccGCTATACACCGTACTCGCCCAGCTGGGCGTACCGTCTCTTCTGACCATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTcgc                    <  2:755732/140‑1 (MQ=255)
                                                                                                                   acgacgCTAGACGTCGTACTGACCCAGCGGGGCGTACCGTCTCGTATGAACATCGGTCAAATCCTCGAAACCCACCTGGGTATGCCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTTAAACAGCAGCAAGAAGTcgc                    <  2:548389/143‑1 (MQ=255)
                                                                                                                       cGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGGAAGCAGTCGCGaaa                >  1:647705/1‑143 (MQ=255)
                                                                                                                       cGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGaaa                <  1:267779/143‑1 (MQ=255)
                                                                                                                       cGGTAGACATCGTACTGAACCAGCTGGGCGTAAAGTATCACATGAACATCGGTGAGATCCTCGAAAACCACCTGGGTATGTGGGCGGAAGGGAGCCGCTAGCAGATCAACGCCATGCTCAAAGCGCAGCAAGAAGAcgcgcga                >  2:19899/1‑140 (MQ=255)
                                                                                                                                 cGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAAt      >  1:1010676/1‑143 (MQ=255)
                                                                                                                                    aCTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCa   >  2:1040061/1‑143 (MQ=255)
                                                                                                                                     cTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCAt  <  2:1010676/143‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                       
TCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCAT  >  NC_000913/4184403‑4184678

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