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% 146.3 / ‑6.7 52P1100Q (CCG→CAG) rpoBRNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (29/22);  minor base T (0/1);  total (29/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.

GCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATT  >  NC_000913/4184407‑4184675
                                                                                                                                         |                                                                                                                                    
gCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTggg                                                                                                                                >  1:16010/1‑143 (MQ=255)
gCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCAAGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACCGAACCAGCTggg                                                                                                                                >  2:687294/1‑143 (MQ=255)
     tAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTAc                                                                                                                           <  1:205676/143‑1 (MQ=255)
       aaggCCGTATCCAGCCTGGTTACAAGATGGCAGGTCGTACCGATCAACTACTGTTTATTTCTAAGATCAACCCGATCGAAGATATGCCTTACTATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCg                                                                                                                         <  1:735571/140‑1 (MQ=255)
          gCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGAAGACATCGTACTGAACCAGCGGGGCGTACCGtac                                                                                                                      >  2:521415/1‑141 (MQ=255)
           ccGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGGACCGtctc                                                                                                                     >  2:830847/1‑143 (MQ=255)
                tCCAGCCTGGTGCCAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAATCGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATg                                                                                                                <  1:281459/143‑1 (MQ=255)
                     ccTGTTGACACGATGTCACGTCGTCACGGTAACAAGGGTGTAATTTCTAAGACCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACAt                                                                                                           <  2:745658/143‑1 (MQ=255)
                          tGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTc                                                                                                      >  2:644921/1‑143 (MQ=255)
                             caggATGGCAGGTCGTCACGGTAACAAGGGTGTAATGTCTCAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACAGCGTACTGAACCAGCTGGGCGTACCGTCTCGTATAAACATCGGTCAGa                                                                                                   <  2:395213/140‑1 (MQ=255)
                                gATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATcc                                                                                                >  2:79533/1‑143 (MQ=255)
                                 aTGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAGACGGTACGCCGGTAGACAGCGTACTGAACCAGCTGGGCGTACCGTCGCGTATGAACATCGGTCAGATCCt                                                                                               <  2:511833/143‑1 (MQ=255)
                                  ttgCATGTCGTGAGGTTAATAAGGGTGTAATTTCTAATATCAAGCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTc                                                                                              <  1:777487/141‑1 (MQ=255)
                                    gCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGa                                                                                            <  1:560194/143‑1 (MQ=255)
                                      aggaccccACGGTAACAAGGGTGGAATTTCTAAGAGCAACCCAATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACGTCGTACTGAACCAGCTGGGAGCACCGCATCGAAAGGACATCGGTCAGCTAGTCGaaa                                                                                          >  2:338732/8‑143 (MQ=255)
                                      aGGTCGTCACGGTAACATGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGaaa                                                                                          <  1:653284/143‑1 (MQ=255)
                                      aGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGATCAGATCCTCGaaa                                                                                          >  2:383408/1‑143 (MQ=255)
                                      aGGGCGTCACGGTAACAAGGGTGTAATTTCTAAGATCGACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGaaa                                                                                          <  2:398942/143‑1 (MQ=255)
                                            tcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCAcc                                                                                    >  1:463875/1‑143 (MQ=255)
                                              aCGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTg                                                                                  >  1:859303/1‑143 (MQ=255)
                                                ggTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTggg                                                                                <  1:493906/143‑1 (MQ=255)
                                                cgctgagaaGGGTGTCATTTCTACGAGCATCCCCATCGATGTTACTCCTTACGGCGACCACGGTACGCCGCTAGACATCGTACTCACCCTGCGTGCCCTACCCTCTCGTATGATCATCGGTCAGATCCTCGAAACCCACCTgg                                                                                 <  2:847034/136‑1 (MQ=255)
                                                       aaGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCGCGAAACCCGCCTGGGGATGGc                                                                          >  2:86795/1‑143 (MQ=255)
                                                          ggTGTATTTTCTAGCATCAACCCGATCGTACATATGTCTTAGGATGGAAACGCTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAATCCCACCTGGGTATGGCTGc                                                                       <  2:722626/143‑1 (MQ=255)
                                                          cgtgtAATCTCGAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGc                                                                       <  1:864131/142‑1 (MQ=255)
                                                            tgtAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGATCAGATCCTCGCAACCCACCTGGCTCCGGCTGCgc                                                                     >  2:207934/1‑142 (MQ=255)
                                                              tAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGCACATCGGTCAGATCCTCGAAACCCACCCGGGTATGGCTGCGaaa                                                                   >  2:425953/1‑143 (MQ=255)
                                                              tAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGCGTCTGGCTGCGaca                                                                   >  1:100551/1‑141 (MQ=255)
                                                              caaTTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACGGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGaaa                                                                   <  1:540758/142‑1 (MQ=255)
                                                               aaTTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTCTGGCTGCGACAg                                                                  >  1:279393/1‑143 (MQ=255)
                                                                 ttgcTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTATTCATCGTATTGAACCAGCTGGGCGTACCGGCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGt                                                                <  1:572371/140‑1 (MQ=255)
                                                                    cTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATc                                                             >  1:707567/1‑143 (MQ=255)
                                                                     tAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTCGACATCGTACTGGACCAGCGGGGCGTACCGTCGCGTAGGAACGTCGGTCAGACCCGCGAAACCCACCTGGGCATGGCTGCGAAAGGTAACg                                                            >  2:43844/1‑143 (MQ=255)
                                                                       aGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCggg                                                          >  1:883422/1‑142 (MQ=255)
                                                                       aGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGc                                                          >  2:838508/1‑143 (MQ=255)
                                                                        gATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGCATCGGCg                                                         >  1:753450/1‑143 (MQ=255)
                                                                            agcccactCGAAGATATGCCTTCCGGTGCAAGCGGTCCGCCGGTAGACCTCGTACTGAACCACCTGGGCGTACCGCCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACaa                                                     <  2:308502/136‑1 (MQ=255)
                                                                                   tCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGGGTAGCGTGTGGTATGAACATCGGTCAGATCCTCCAACCCCACCTGGGTATGGCTGCGAAACGTATCGGCGACCAGACCaag                                              >  2:871748/1‑142 (MQ=255)
                                                                                   tCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGAACCTCGAAACCCACCTGGGTATGGCTGCGAAAGCGCTCGGCAACAAGAGCacc                                              >  2:868402/1‑141 (MQ=255)
                                                                                      aaGAGATGCCTTAGGAGGAAAACGGTACGCCGGTAGACAGCGGACTGAACGAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGcc                                           <  1:57621/143‑1 (MQ=255)
                                                                                               cTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGaaa                                  <  2:707567/143‑1 (MQ=255)
                                                                                                    gatgatAGCGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAAcagca                             <  1:424868/143‑1 (MQ=255)
                                                                                                     aTGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAAcagcag                            >  1:665829/1‑143 (MQ=255)
                                                                                                      tGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACagcagc                           <  1:79533/143‑1 (MQ=255)
                                                                                                       gAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCGGAAACAgcagca                          >  1:484999/1‑143 (MQ=255)
                                                                                                        aaaaCGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATAGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGGGAGGTATCGGCGACAAGATCAACGCCATGCGGCAACAGGGGCaa                         >  1:101853/1‑143 (MQ=255)
                                                                                                         aaaCGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCaag                        >  1:918021/1‑143 (MQ=255)
                                                                                                                   ccGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGaa              >  1:881460/1‑143 (MQ=255)
                                                                                                                     ggTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAAc            <  1:796299/143‑1 (MQ=255)
                                                                                                                            aTCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGa     >  1:804435/1‑143 (MQ=255)
                                                                                                                             tCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGGATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCCGCAACAGCAGCAAGAAGTCCCGCACCCgcgcaga    >  2:553599/1‑140 (MQ=255)
                                                                                                                               gTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAAtt  <  2:209587/143‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                    
GCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATT  >  NC_000913/4184407‑4184675

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: