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% 102.5 / ‑6.7 39P1100Q (CCG→CAG) rpoBRNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (15/23);  total (16/23)
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.

GGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCA  >  NC_000913/4184406‑4184681
                                                                                                                                         |                                                                                                                                          
tccagTTATACGCCGTAGCCAGCGTGTTGTCAAGATGGAAGGTGGTCACGGTGACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTgg                                                                                                                                       <  2:254858/139‑1 (MQ=255)
   ggTTAAACGCGGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGCACGCCGGTAGACATCGCACTGAACCAGCCGGGCg                                                                                                                                    >  2:784822/1‑143 (MQ=255)
           gCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGtct                                                                                                                            <  1:262265/143‑1 (MQ=255)
            ccGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGtctc                                                                                                                           >  2:229824/1‑143 (MQ=255)
               tATCCGGCCTGGTGACAGGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTa                                                                                                                        <  2:320046/143‑1 (MQ=255)
                 tCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATg                                                                                                                      <  2:743294/143‑1 (MQ=255)
                   cAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGaa                                                                                                                    >  1:370088/1‑143 (MQ=255)
                    aGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGagc                                                                                                                   >  2:193365/1‑141 (MQ=255)
                      ccTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGGACCGTCTGGTATGAACAt                                                                                                                 >  2:208491/1‑143 (MQ=255)
                             aCAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAg                                                                                                          >  2:616254/1‑143 (MQ=255)
                             aCAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAg                                                                                                          >  2:570191/1‑143 (MQ=255)
                              cAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGa                                                                                                         <  1:124973/143‑1 (MQ=255)
                              cAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGa                                                                                                         <  2:290546/143‑1 (MQ=255)
                                     gccggTCGTCACGGTAACAAGGGTGTACTTTCTAAGATCAACCCGATCGAAGCTATGCCTTACGCTGAAACCGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGa                                                                                                  <  2:641006/140‑1 (MQ=255)
                                      cAGGTCGTCCTGGTAACAAGGGTGTCATTGCTAACACCAACCCAATCGAAGAGATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGaa                                                                                                 <  1:156680/143‑1 (MQ=255)
                                      cAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGaa                                                                                                 <  1:406089/143‑1 (MQ=255)
                                       agttcgtcACGGTAACAAGGGTGTACTCTCTAAGATCACCCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGaaa                                                                                                <  1:608403/140‑1 (MQ=255)
                                         gtcgtcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATCAACCTCGGTCAGATACTCGaaagg                                                                                              >  1:174925/1‑141 (MQ=255)
                                             tcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACAACCGTCAGGTCCTCGAAACCCAcc                                                                                          >  1:606515/1‑143 (MQ=255)
                                                 ggTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGAACCTCGCAACCCACCCggg                                                                                      >  1:377965/1‑143 (MQ=255)
                                                     aCAAGGGTGTAATTTCTACGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATg                                                                                  <  1:298994/143‑1 (MQ=255)
                                                             gtAATTTCCAACATCAACCCGATCGAAGAGATGCCTTACGATGCAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAGCATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGaa                                                                          <  2:399344/143‑1 (MQ=255)
                                                              tAACTTCTAAGCTCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGaaa                                                                         <  1:440545/143‑1 (MQ=255)
                                                                aTTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCGGCGACAgg                                                                       >  2:574946/1‑143 (MQ=255)
                                                                 gttgtAAGATCAACCCGATCGAGGATATGCCTTACGATGGAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGt                                                                      <  2:603036/139‑1 (MQ=255)
                                                                      aaGATCAACCCAATCGAGGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCgg                                                                 <  2:305682/143‑1 (MQ=255)
                                                                         ctCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGCAGACATCGTACTGAACCAGCCCGGCGTACCGTCTCGTCTGAACATCGGCCAGCGCCCCGAACCCCCCCCCGGCACGGCCGCGCAAGGTACCGGccc                                                              >  2:241323/2‑141 (MQ=255)
                                                                            aaCCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACaa                                                           <  1:229824/143‑1 (MQ=255)
                                                                                   tCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAAc                                                    <  2:105002/143‑1 (MQ=255)
                                                                                    cGAAGATACGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACg                                                   <  1:418426/143‑1 (MQ=255)
                                                                                                   cgctGAAAGCGGTACGCCGGTAGACCTCGTACTGAACCAGCTGGGCGTTCCGTCTCGTATGAACATCGGTCAGATCCCCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAAcagc                                    <  2:447550/140‑1 (MQ=255)
                                                                                                           aCGGTCCCCCGTTACACATCGTACTGACCCAGCTGGGCGTACCGACTCGTATGAACATCGGTCAGATCTTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCaagaa                            <  1:28750/143‑1 (MQ=255)
                                                                                                                  gCCGGTAGACATCGTACTGCAGCAGCTGGGCGTGCTGTATCACCTGAACAGCTGTCAGATCCTCGAAACCCTCCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGa                     <  1:570191/143‑1 (MQ=255)
                                                                                                                  gCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCGCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGa                     <  1:733903/143‑1 (MQ=255)
                                                                                                                   ccGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGGCATGCTGAAACAGCAGCAGGAAGGCGCGaa                    >  1:31959/1‑143 (MQ=255)
                                                                                                                             tCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGaa          >  1:735836/1‑143 (MQ=255)
                                                                                                                               gTACTGACCCACCTGGGCGGACCGTCTCGTACGAGCATCGGTCAAATCCTCGAAACCCACCTGGTTATCGCTGCGAAAGCTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAAtt        <  2:373045/143‑1 (MQ=255)
                                                                                                                                     aaCCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCa  >  2:757595/1‑143 (MQ=255)
                                                                                                                                     aaCCAGCTGGGCGTACCGTCTCGTATGAACATCAGTCAGATCCTCGAAACCCACGTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAGGAAGTCGCGAAGCTGCGCGAATTcctgca  >  1:159464/1‑138 (MQ=255)
                                                                                                                                         |                                                                                                                                          
GGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCA  >  NC_000913/4184406‑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: