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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,184,5420CA97.4% 112.4 / NA 38P1100Q (CCG→CAG) rpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (18/19);  total (18/20)
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

TCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCA  >  NC_000913/4184402‑4184680
                                                                                                                                            |                                                                                                                                          
tCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGACCCAGc                                                                                                                                          <  1:666900/143‑1 (MQ=255)
   ggCGGTTAAACGCCGTATACAGCCTGGTGACAAGATGGCAGGTCGCCACGGGAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTgg                                                                                                                                       <  2:62151/143‑1 (MQ=255)
      ggTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAAAACAATCTAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCg                                                                                                                                    <  1:506236/143‑1 (MQ=255)
       gTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGt                                                                                                                                   >  1:483336/1‑143 (MQ=255)
            aCGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCgg                                                                                                                              >  1:191821/1‑142 (MQ=255)
              gCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAAATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGAACGCCGGGAGACATCGTACTGAACAAGCGGGGCGTagaagct                                                                                                                            >  2:781156/1‑137 (MQ=255)
               ccGTATCCAGCCTGGTGACACGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGtctc                                                                                                                           <  2:552622/143‑1 (MQ=255)
                   aTCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTAt                                                                                                                       >  2:553186/1‑143 (MQ=255)
                           tGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCg                                                                                                               <  1:563241/143‑1 (MQ=255)
                                   aGATGGCAGGTCGTCACGGTAACAAGGGGGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGAACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGAATGAACATCGGTCAGATc                                                                                                       >  2:386911/1‑143 (MQ=255)
                                          gcgtcgtcGCGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGCTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGaaa                                                                                                <  2:413106/141‑1 (MQ=255)
                                                 cACGGTAACGAGGGTGGAACTGCTAAGAGCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCt                                                                                         <  2:407398/143‑1 (MQ=255)
                                                 cACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCt                                                                                         <  1:262627/143‑1 (MQ=255)
                                                 cACGGTAACAAGGGTGTAATTTCTAAGATCAACCAGACAGAGAATCGGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGACGTACAGTCTCGTATGAACATCGCGCAGATCCTCGAAACCCACCt                                                                                         >  2:439509/1‑143 (MQ=255)
                                                  aCGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGGCGAAAACGACACGCCGATAGACATCGTCCTGAACCAGCTGGGCGCACCGTGTCGTATGAAAATCGGTCGGAACCCCGAAcccagcccc                                                                                        >  2:700252/1‑137 (MQ=255)
                                                               tgtAATTTCTAAGATGAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCGCGATACCCACCTGGGTATGGccggga                                                                           >  2:677478/1‑138 (MQ=255)
                                                                        tAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCGTCGAAACCCACCTGGGTGTGGCTGCGAAGGGGATCg                                                                  >  1:4880/1‑143 (MQ=255)
                                                                         cagctCAACCCGATCGAGGATATGTCTTAGCATGTAAGCGGTACGTCGGTAGAGATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCgg                                                                 <  2:512338/139‑1 (MQ=255)
                                                                         agggtCAACCCGATCGAGGATATGCCTTACGAGGGATACGGTACGCCGGTAGACATTGTGCTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAGCCCACCTGGGTATGGCTGCGAAAGGTATCgg                                                                 <  2:725607/139‑1 (MQ=255)
                                                                         aacaTCACCCCGATCGAAGCTATGCCCTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCgg                                                                 <  1:527395/140‑1 (MQ=255)
                                                                            aTCAACCCGATCGAAGATATGCCTTACGATGAAAACTGGACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTCTGGCTGCGAAAGGTATCGGCgc                                                              >  1:340586/1‑142 (MQ=255)
                                                                              cAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTCTGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGAcc                                                            >  2:513514/1‑142 (MQ=255)
                                                                              cAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACa                                                            <  1:553186/143‑1 (MQ=255)
                                                                                 cccGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGa                                                         <  1:774913/143‑1 (MQ=255)
                                                                                  ccGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCGGGGTATGGCTGCGAAAGGTATCGGGGACAAggt                                                        >  1:154757/1‑141 (MQ=255)
                                                                                      tCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCAGCGGCAAGATCAAc                                                    >  1:383637/1‑143 (MQ=255)
                                                                                          agagATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCa                                                <  1:111365/143‑1 (MQ=255)
                                                                                          aGATATGCCTTAAGATGAAAACGGTGCGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGGCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCa                                                <  1:755015/143‑1 (MQ=255)
                                                                                                  cTTGCGACGAAAACGGTACGCGGGTGGTCGACGTACTAAACCCGTAGTGCGTTCCTGCTCGTATGAACAGCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGaaa                                        <  2:739465/143‑1 (MQ=255)
                                                                                                      cgacgaAAACGGTACGCCGGTAGACAGCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCCGCGAAAGATATCGGCGACAAGATCAACGCCATGCTGAAAcagc                                    <  2:366201/143‑1 (MQ=255)
                                                                                                         tGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACagcagc                                 <  1:32927/143‑1 (MQ=255)
                                                                                                          gAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAgcagca                                >  2:371212/1‑143 (MQ=255)
                                                                                                             aaCGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATAAACGCCATGCTGCAACAGCAGCaagc                             >  2:522544/1‑142 (MQ=255)
                                                                                                              aCGTTACGCCATTAGACATCGTACTTTACCATCTTGGCGTACCGTCACGGATGAACGACGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCaagaa                            <  1:307766/143‑1 (MQ=255)
                                                                                                                      ccGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCCCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGaa                    <  2:683261/143‑1 (MQ=255)
                                                                                                                                tCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCCGATCCTCGAAACCCACCTGGGTACGGACGGGAAATGTAGCGGCCCCAACATCCACGCCATGCTGAAACAGCAGCAAGAAGTCGGGAAACTGCGCGaa          >  2:240627/1‑143 (MQ=255)
                                                                                                                                   tACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTc       >  1:658353/1‑143 (MQ=255)
                                                                                                                                        aaCCAGCTGGGCGTACAGTCTCGTATGAACATCGGTCAGATCCTCGAAAACCACCTGGGTATGGCTGCGAAAGATATCGGCGACAAGATCAACGCCATGCTGCAACAGCAGCAAGAAGTCGCGACACTGCGCGAAGTCCTCCa  >  1:170131/1‑143 (MQ=255)
                                                                                                                                            |                                                                                                                                          
TCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCA  >  NC_000913/4184402‑4184680

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