Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,214,757 T→C Y571H (TAC→CAC)  rpoD → RNA polymerase, sigma 70 (sigma D) factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,214,7570TC100.0% 85.9 / NA 25Y571H (TAC→CAC) rpoDRNA polymerase, sigma 70 (sigma D) factor
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (12/13);  total (12/13)

ACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCG  >  NC_000913/3214622‑3214891
                                                                                                                                       |                                                                                                                                      
accaccCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACcaca                                                                                                                                     >  2:668280/1‑139 (MQ=255)
accaccCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACcaca                                                                                                                                     >  2:101484/1‑139 (MQ=255)
  caccCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACCACACg                                                                                                                                   <  1:930996/139‑1 (MQ=255)
          aGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACCACACGCTGgaaga                                                                                                                           >  1:1066190/1‑139 (MQ=255)
                          tGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACCACACGCTGGaagaag                                                                                                                         >  2:546160/1‑125 (MQ=255)
                          tGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACCACACGCTGGaagaag                                                                                                                         <  1:546160/125‑1 (MQ=255)
                                caccGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACCACACGCTGGAAGAAGTGGGTAAACAGTTCGACGtt                                                                                                     <  2:1066190/139‑1 (MQ=255)
                                   cGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGTCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACCACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTAcc                                                                                                  >  1:1207686/1‑139 (MQ=255)
                                        gCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACCACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGa                                                                                             <  1:101484/139‑1 (MQ=255)
                                        gCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACCACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGa                                                                                             <  2:1107107/139‑1 (MQ=255)
                                             cgtgcgGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACCacac                                                                                                                                    >  2:624499/1‑95 (MQ=255)
                                             cgtgcgGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACCacac                                                                                                                                    <  1:624499/95‑1 (MQ=255)
                                                    cAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACCACACGCTGGAAGAAGTGGGTaa                                                                                                                  >  2:902552/1‑106 (MQ=255)
                                                    cAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACCACACGCTGGAAGAAGTGGGTaa                                                                                                                  <  1:902552/106‑1 (MQ=255)
                                                                  ctggctggCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACCACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAgg                                                                   >  2:562926/1‑139 (MQ=255)
                                                                                                   cTGCGTATGCGTTTCGGTATCGATATGAACACCGACCACACGCTGGTAGAAGTGGGTAAACAGTTCGACTTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCta                                  <  1:875807/139‑2 (MQ=255)
                                                                                                      cGTATGCGTTTCGGTATCGATATGAACACCGACCACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGa                                                                                             >  2:404646/1‑77 (MQ=255)
                                                                                                      cGTATGCGTTTCGGTATCGATATGAACACCGACCACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGa                                                                                             <  1:404646/77‑1 (MQ=255)
                                                                                                              tttCGGTATCGATATGAACTCCGACCACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGt                       <  2:807492/139‑1 (MQ=255)
                                                                                                                tCGGTATCGATATGAACACCGACCACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAg                     >  2:107251/1‑139 (MQ=255)
                                                                                                                       cGATATGAACACCGACCACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAAc                                                                                           >  2:1006312/1‑62 (MQ=255)
                                                                                                                       cGATATGAACACCGACCACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAAc                                                                                           <  1:1006312/62‑1 (MQ=255)
                                                                                                                       cGATATGAACACCGACCACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTg              <  1:284826/139‑1 (MQ=255)
                                                                                                                                cacCGACCACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCg                                            <  2:103412/100‑1 (MQ=255)
                                                                                                                                cacCGACCACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCg                                            >  1:103412/1‑100 (MQ=255)
                                                                                                                                  ccGACCACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATc   >  1:862778/1‑139 (MQ=255)
                                                                                                                                   cGACCACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCg  >  2:1076155/1‑139 (MQ=255)
                                                                                                                                       |                                                                                                                                      
ACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCG  >  NC_000913/3214622‑3214891

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: