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.6% 109.9 / ‑6.5 42P1100Q (CCG→CAG) rpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (16/25);  minor base T (1/0);  total (17/25)
Fisher's exact test for biased strand distribution p-value = 4.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATC  >  NC_000913/4184407‑4184678
                                                                                                                                       |                                                                                                                                        
cGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGc                                                                                                                                   <  1:648751/143‑1 (MQ=255)
        cGCCGTAACCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCTAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGtc                                                                                                                           <  1:358468/143‑1 (MQ=255)
          ccGTATCCAGTCTTGGTACAGAATGCCAGGTTGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCTATCGAAGATAGGCCTTACGATGAAAACGGTACGCCGTTAGACATCGTACTGAACCAGCTGGGCGTACCGtctc                                                                                                                         <  1:612452/143‑1 (MQ=255)
           cGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATAGGCCTTACGATGAAAACGATACGCCGGTAGACAACGTAGTGAACAAGCTGGGCGGACCGTCACg                                                                                                                        >  1:386352/1‑143 (MQ=255)
             tATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTGTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTa                                                                                                                      <  2:328682/143‑1 (MQ=255)
                 cAGCCTGGTGACAAGATGGCGCGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGaa                                                                                                                  <  1:480254/143‑1 (MQ=255)
                 cAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGaa                                                                                                                  <  2:538247/143‑1 (MQ=255)
                  aGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGACGATATGCCTTACGATGAAAACGGTACGCCGGTAGACAGCGTGCTGAACCAGCTGGGCGTACCGTCTCGTCTGAAc                                                                                                                 >  2:235432/1‑143 (MQ=255)
                   gCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGAGCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGCACTGAACCAGCTGGGCGTACCGTCTCGTATGAGCa                                                                                                                >  1:634771/1‑143 (MQ=255)
                   gCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGGGTAATTTATAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGGACCGTCCCGTATGGACa                                                                                                                >  1:306303/1‑143 (MQ=255)
                        gTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGACTCCTATGAACATCgct                                                                                                           >  1:236310/1‑141 (MQ=255)
                        gTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGGACCGTCTCGTGTGAACACCGGt                                                                                                           >  1:360969/1‑143 (MQ=255)
                              aGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATG‑AAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATcc                                                                                                    <  2:3529/143‑1 (MQ=255)
                               gATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCGGGTAGACATCGTCCTTAACCAGATGGGCGTAACGCCTCGTATGGACATCGGTCCGATcc                                                                                                    >  2:645380/1‑143 (MQ=255)
                               gATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGACCCAGCTGGGTGTATCGCCTCGTATGAACATCGGTCAGTTcc                                                                                                    >  2:264729/1‑143 (MQ=255)
                               gATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTAGCGGCTCGTATGAACATAGGTCAGATcc                                                                                                    >  2:43341/1‑143 (MQ=255)
                               gATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGGCCCGTCGGACCATCGGGCAGATcc                                                                                                    >  1:559162/1‑143 (MQ=255)
                                        tcgtcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACCTCGGTCAGATCCTCGAAAccc                                                                                           >  2:443462/1‑143 (MQ=255)
                                                    cAACGGTCTAATTTCTAAGATCAACCCGATCGAATATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATgg                                                                               <  2:475460/143‑1 (MQ=255)
                                                      aGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAGAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGTCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCt                                                                             <  1:332042/143‑1 (MQ=255)
                                                         gtgtAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCg                                                                          >  2:99568/1‑143 (MQ=255)
                                                           ctaacacatAAAATCATCTCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGaa                                                                        <  2:94177/135‑1 (MQ=255)
                                                            tAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACAAACTGGTCGTATTGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGaaa                                                                       <  2:269177/143‑1 (MQ=255)
                                                             actttCCACGATCAGCGCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAg                                                                      <  2:274397/141‑1 (MQ=255)
                                                                     tgATCAACCCGATGTACGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGc                                                              <  1:598703/142‑1 (MQ=255)
                                                                                 tCGAAGATATGCCTTACGAGGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAAc                                                  <  2:702198/143‑1 (MQ=255)
                                                                                 gcGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAAc                                                  <  1:97138/142‑1 (MQ=255)
                                                                                      ggtatGCCTGACGAGGAAAACGGTACGCCGGCAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCAt                                             <  2:236310/141‑1 (MQ=255)
                                                                                      gATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCAt                                             <  2:98949/143‑1 (MQ=255)
                                                                                         atGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCt                                          >  2:503344/1‑143 (MQ=255)
                                                                                                 cGAGGAAAACGGCACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCGCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAAcagc                                  <  1:607900/143‑1 (MQ=255)
                                                                                                  gacgAAAACGGTACGCCGGTCGACATCGTACTGAACCAGCTGGCCGTACCGTCTCGTATGAACACCGGTCAGTTCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAAcagca                                 <  1:53411/140‑1 (MQ=255)
                                                                                                  gATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAAcagca                                 <  2:315507/143‑1 (MQ=255)
                                                                                                    tGAAAACGGTACGCCGGTAGACATCGTACGGCACATACTTTGCGTGCCGGCTCGTATGAACATCGGTCAGGTCCTCTATACTCAGCTGGGTTTGGCTGCGGAGTGTGTCGGCGACAAGAACAAAGCCATGCTGAAAcaggatg                               >  1:474011/1‑139 (MQ=255)
                                                                                                             tACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGAACAACGCCATGCTGAAACAGCAGCAAGAAGTcg                      >  2:86053/1‑143 (MQ=255)
                                                                                                                gCCGGTAGACATCGTACTGGACCAGCTGGGCGTATGCTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGa                   <  1:503344/143‑1 (MQ=255)
                                                                                                                       gacggcGTACTGAACCAGCTGGGCGTACCGGCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTgcg            <  1:408149/138‑1 (MQ=255)
                                                                                                                           tCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGACACTGCGCGaa        >  1:24545/1‑143 (MQ=255)
                                                                                                                           tCGTACCGAACCAGCTGGGCGTACCGTCTCGTATGAAGATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGaa        <  2:141480/143‑1 (MQ=255)
                                                                                                                            cGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTTGGTATGTCTGCGAAATGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAAt       <  2:704222/143‑1 (MQ=255)
                                                                                                                               aCTGAACCAGCTGTGCTTCCCGGCTCGTATGAACATCGGTCAGACCCTCGAAAGCCACCTGGGTATGGCGGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCa    <  1:643730/143‑1 (MQ=255)
                                                                                                                                 ttaaCCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCCCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGAAAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCaac  >  2:684500/3‑141 (MQ=255)
                                                                                                                                       |                                                                                                                                        
CGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATC  >  NC_000913/4184407‑4184678

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