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,5430CA100.0% 104.5 / NA 34P1100Q (CCG→CAG) rpoBRNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (16/18);  total (16/18)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTC  >  NC_000913/4184409‑4184676
                                                                                                                                      |                                                                                                                                     
gcttAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGCCGTCACGGTAACAAGGGTGTAATTTCTAAGATCATCCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCg                                                                                                                               <  2:493141/141‑1 (MQ=255)
         cggtggccAGACTGGTGCTGAGATGGGGGGTCGTAAGCGTAACAAGGGTGTAATGTCGAAAATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGCTAGACATCGTACTGAACCAGCTGGGCGTACCGtctc                                                                                                                      <  2:204821/137‑1 (MQ=255)
              tCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGGTATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGACCCAGCTGGCCGTACCGTGTCGTATg                                                                                                                 >  1:150504/1‑143 (MQ=255)
                cAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGaa                                                                                                               <  2:350167/143‑1 (MQ=255)
                 aGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGCAc                                                                                                              >  2:317356/1‑143 (MQ=255)
                      ggTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCgg                                                                                                         >  2:44059/1‑143 (MQ=255)
                            aaGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGAt                                                                                                   <  1:219473/143‑1 (MQ=255)
                            aaGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGACGTACCGTCTCGTATGAACATCGCGCAGAt                                                                                                   >  1:82709/1‑143 (MQ=255)
                                  gCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCCATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGa                                                                                             <  1:500623/143‑1 (MQ=255)
                                           cACGGTAACAAGGGTGTAATTGCTAAGGTCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCt                                                                                    <  2:48143/143‑1 (MQ=255)
                                             cGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTgg                                                                                  >  1:402874/1‑143 (MQ=255)
                                                   cAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATgg                                                                            >  1:221242/1‑143 (MQ=255)
                                                         tggactttCTGAGATCGAGCCGATGGAAGATATGCCGTACGATTAAAACGGTACGCCGGTAGACATTGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGTGa                                                                      <  2:495321/138‑1 (MQ=255)
                                                           tAATTTCTAAGACCAACCCGATCGAAGATATGCCTTACGTTGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGaaa                                                                    <  2:214979/143‑1 (MQ=255)
                                                              tttCTAAGATCAACCCGATCGAAGATATGCCTTACGATCAAAACGGTACGCCGGTAAACAACGTACGGAACCAGGTGGGCGTACCGTCTCGGATGCAGATCGGTCAGACCGTCCACACCCCCCTGGGTAGTGCTGCGAGAGGt                                                                 >  2:190327/1‑143 (MQ=255)
                                                                       tCAACCCGATCGAAGATATGCTTTACTATGAAAACGGTACGTTGGTATACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGAc                                                        <  1:419251/143‑1 (MQ=255)
                                                                        cAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACa                                                       <  1:44059/143‑1 (MQ=255)
                                                                            ccGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGGCTCGTCTGAACATCGGTCAGATCCTCGAAACCCACCGGAGCATGGCTGCGAAGGGTCTCGGCGAccaaat                                                   >  2:416464/1‑138 (MQ=255)
                                                                             cGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGCTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAACGGTATCGGCGGCAAGATc                                                  >  1:339737/1‑143 (MQ=255)
                                                                               aTCGAAGATATGCCTTACGATGAGGACGGTACGCCGGTAGACAGCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCaa                                                <  2:239935/143‑1 (MQ=255)
                                                                                      atcggCCGTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATg                                         <  2:328219/139‑1 (MQ=255)
                                                                                        atGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAAACAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGACTGGCCACGAAAGGTATCGACGCGCAGATACACGCCATgat                                       >  1:37699/1‑141 (MQ=255)
                                                                                         tGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGACCATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGCAAGGTATCGGCGACAAGATCAACGCCCCGCTg                                      >  2:98881/1‑143 (MQ=255)
                                                                                            cTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGaaa                                   <  1:28822/143‑1 (MQ=255)
                                                                                             ttACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAGGATCAACGCCATGCTGGAAc                                  >  1:316550/1‑143 (MQ=255)
                                                                                                 gATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAAcagca                              <  1:417351/143‑1 (MQ=255)
                                                                                                     caaaCGGCAGCCCGGGAGACATCGGACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAATCCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCaa                          <  1:373176/142‑1 (MQ=255)
                                                                                                          tgTAAGCCGGTAGACATCGGACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGt                     <  1:317228/142‑1 (MQ=255)
                                                                                                               gCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGAATGGCTGCGAAAGGTATCGGCGGCAAGATCAACGCCATGCTGAAACAGCAGCACGAAGCcgggg                >  1:514907/1‑140 (MQ=255)
                                                                                                                   gTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACt            >  1:187101/1‑143 (MQ=255)
                                                                                                                   gTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCATAGAAACCCACCTGGGTATGGCTGCGAAAGGTAACGGCGAAAAGATCAACGCCCAGGGGAAACAGCAGAAAGAAGAcgacaacca            >  2:143717/1‑136 (MQ=255)
                                                                                                                         aTCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCAGCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAGGTCGCGAAACTGCGCGa      >  1:400187/1‑143 (MQ=255)
                                                                                                                          tCGCACGGAACCAGCCGGGCGTACCGTCGCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGaa     <  1:209091/143‑1 (MQ=255)
                                                                                                                             tACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTc  <  1:455994/143‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                     
GGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTC  >  NC_000913/4184409‑4184676

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