Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,187,597 C→T 100% P750L (CCG→CTG)  rpoC → RNA polymerase, beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,187,5970CT100.0% 161.8 / NA 50P750L (CCG→CTG) rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (28/22);  total (28/22)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTC  >  NC_000913/4187451‑4187740
                                                                                                                                                  |                                                                                                                                               
gCAAACTGAAACCGTGATGAACCGTAACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGTTGGTATGCGTGGTCTGATGGCGAAGCTgg                                                                                                                                               <  1:447347/149‑1 (MQ=255)
  aaaCTGAAACCCTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGAt                                                                                                                                             <  1:340276/149‑1 (MQ=255)
        aaaCCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTcc                                                                                                                                       >  2:376616/1‑149 (MQ=255)
                        tGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCACGTTGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCaa                                                                                                                       >  1:110742/1‑149 (MQ=255)
                            ggTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCAc                                                                                                                   <  1:275619/149‑1 (MQ=255)
                                  gAAGAGGAGCCGGTTCCCTTCAACAGCAGCTACATGGGGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGaa                                                                                                             <  2:269285/149‑1 (MQ=255)
                                  gAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGaa                                                                                                             >  1:283339/1‑149 (MQ=255)
                                  gAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACAGCGaa                                                                                                             <  1:215294/149‑1 (MQ=255)
                                        aaGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATTCGTGGTCTGATGGCGAAGCTGAGCGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCg                                                                                                       >  1:497489/1‑149 (MQ=255)
                                        aaGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCCCCATCATCGAAACGCCAATCACCGCGAACTTCCg                                                                                                       >  2:117493/1‑149 (MQ=255)
                                               tttCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGt                                                                                                <  1:429593/149‑1 (MQ=255)
                                               tttCCGTCAACAGCATCTACATGATGGCCGACTCCGGTGCCCGGGGTGCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGt                                                                                                <  2:153912/149‑1 (MQ=255)
                                                   cTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGa                                                                                            >  2:291280/1‑149 (MQ=255)
                                                    ttCAGCAGCATCTACTTGATGTCCGATTCCGGTGCGCGTGGGTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGaa                                                                                           >  2:68468/1‑149 (MQ=255)
                                                    ttCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGaa                                                                                           >  1:343240/1‑149 (MQ=255)
                                                    ttCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGaa                                                                                           >  1:75183/1‑149 (MQ=255)
                                                     tCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGCAACGCCAATCACCGCGAACTGCCGTGAAGGTCTGAAc                                                                                          <  1:465521/149‑1 (MQ=255)
                                                     tCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAAc                                                                                          <  1:422404/149‑1 (MQ=255)
                                                        aCAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGCTGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGGGAAGGTCTGAACGTa                                                                                       >  2:506589/1‑149 (MQ=255)
                                                           gCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTc                                                                                    <  1:76189/149‑1 (MQ=255)
                                                             aTCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATCCGTCAGCTGGCTGGTATGCGTGGTCTGATGGCGAAGCTGGAGGGCTCGGTCATCGCAACGCCAATAACGGCGAACTTCCCTGACGGACTGAACGTCCGCCa                                                                                  >  2:258576/1‑149 (MQ=255)
                                                                tACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAGACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTa                                                                               >  1:540959/1‑149 (MQ=255)
                                                                 accggatgGCCGACTCCGGTGCGCGTGGTGCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCACCAGCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTAc                                                                              <  1:688567/145‑1 (MQ=255)
                                                                 aCATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCACTCACCGCGAGCTTCCGTGAAGGTCTGAACGTACTCCAGTAc                                                                              >  2:584873/1‑149 (MQ=255)
                                                                 aCATGATGCCCGACTCCGGTGCGCCGGGTTCCGCGGCACAGACTCCTCAGCTTGCCGGTATGCGTTGTCTGAGGGCGAAGCTGGATGGCTCCCTAAACGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTAc                                                                              <  1:688492/149‑1 (MQ=255)
                                                                        ggCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCAtct                                                                       >  1:691952/1‑149 (MQ=255)
                                                                          ccGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCAGAACGTACTCCAGTACTTCATCTcc                                                                     <  1:365974/149‑1 (MQ=255)
                                                                           cGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGCAGGTCTGAACGTACTCCAGTACTTCATCTcca                                                                    >  1:652648/1‑149 (MQ=255)
                                                                           cGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTcca                                                                    <  1:74668/149‑1 (MQ=255)
                                                                                 cGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGGTGGCTCCATCATCGAAAGGCCAATCACCGCGGACTTCCGTGAGGGTCTGACCGTACTCCAGTAGTTAGTCGCCAACCACg                                                              >  2:169098/1‑149 (MQ=255)
                                                                                 cGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCGGTGAGGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACg                                                              >  1:158906/1‑149 (MQ=255)
                                                                                 cGGTGCGCCTTGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTCTGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACg                                                              <  2:615083/149‑1 (MQ=255)
                                                                                    ggcgcGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTg                                                           <  2:562923/148‑1 (MQ=255)
                                                                                           ggTTCGGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTaa                                                    <  1:7562/149‑1 (MQ=255)
                                                                                            gTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTaca                                                   >  2:576755/1‑147 (MQ=255)
                                                                                                   ggCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTTTgg                                            >  2:29097/1‑149 (MQ=255)
                                                                                                   ggCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTgg                                            <  1:376616/149‑1 (MQ=255)
                                                                                                    gCACAGGTGCGTCAGCTTGCTGGTAGGCGTGGCCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGc                                           <  2:661230/149‑1 (MQ=255)
                                                                                                              gTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCAGGGTGCTCGTAAAGGTATGGCGGTTACCGCa                                 >  1:408922/1‑149 (MQ=255)
                                                                                                              gTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCa                                 <  1:131183/149‑1 (MQ=255)
                                                                                                                        tGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTg                       >  1:364734/1‑149 (MQ=255)
                                                                                                                              gCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCTCCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAGCt                 >  1:527673/1‑149 (MQ=255)
                                                                                                                                 tGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCg              >  1:401536/1‑149 (MQ=255)
                                                                                                                                   gTCTGTTGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGt            <  1:416357/149‑1 (MQ=255)
                                                                                                                                    tCTGATGGCGAAGCTGGATGGCTCCATCACCGAAACGCCAATCACCGCGAACTGCCGTGAAAGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGCCGGATACTGCACTGAAAAGTGCGAACTCCGGtt           >  1:447862/1‑149 (MQ=255)
                                                                                                                                      tGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTAc         >  1:545446/1‑149 (MQ=255)
                                                                                                                                       gATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGACCTCCGGTTAcc        >  2:492482/1‑149 (MQ=255)
                                                                                                                                        ctGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCACTACTTCATCCCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTCAccc       <  2:549543/148‑2 (MQ=255)
                                                                                                                                          ggCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCGGATTACCTGa     >  1:37564/1‑149 (MQ=255)
                                                                                                                                             gAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAGGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACtc  >  1:690007/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
GCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTC  >  NC_000913/4187451‑4187740

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑