Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,186,086	C→A	P246Q (CCG→CAG)	rpoC →	RNA polymerase, beta prime subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,186,086	0	C	A	100.0%	131.4 / NA	45	P246Q (CCG→CAG)	rpoC	RNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (21/24);  total (21/24)

CGTGAAGAGCTGAACGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATC  >  NC_000913/4185953‑4186222                                                                                                                                      |                                                                                                                                         cattaaGAGCTGAACGAAAACAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCAGCCAGATCt                                                                                                                                 <  1:841362/139‑1 (MQ=255)  gTGAAGAGCTGAACGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGCATGATCCTGACCGTTCTGCCGGTACTGCAGCAAGATCTg                                                                                                                                >  2:90898/1‑143 (MQ=255)   tGAAGAGCTGAACGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCAGCCAGATCTGc                                                                                                                               <  2:850607/143‑1 (MQ=255)          cTGAACGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCAGCCAGATCTGCGTCCGCt                                                                                                                        <  1:940575/143‑1 (MQ=255)             aaCGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTATGGTAACAAACCAGAGCGGATGATCCTGACCGTTCTGCCGGTACCGCAGCCAGATCTGCCTACGCTGGt                                                                                                                     >  2:133632/1‑143 (MQ=255)                   aCCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCAGCCAGATCTGCGTCCGCTGGTTCCGCt                                                                                                               <  1:151276/143‑1 (MQ=255)                     cgaCTCCGAAGCCAAGCGTAAAAAGCTGACCAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGGAACAAACCCGAGTACTTTATCCTGTCCGTTCTGCCGGTACTGCAGCCAGATCTGCGTCCGCTGGTTCCGCTgg                                                                                                             <  1:548651/141‑1 (MQ=255)                          ccGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCAGCCAGATCTGCGTCCGCTGGTTCCGCTGGAtggt                                                                                                        >  1:925564/1‑143 (MQ=255)                          ccGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCAGCCAGATCTGCGTCCGCTGGGTCCGCTGGAtggt                                                                                                        >  1:173437/1‑143 (MQ=255)                            gAAACCAAGCGTAAAAAGCTGACCAATCGTATCAAACTGCTAGAAGAATTCGTTAAGTCTGGTAACAAACAAGAGTGGATGATCCTGACCGTCCTGCAGGAACTGCAGCCAGATCTCCGTCCCCGCGTTCCTCTGGAtgcggg                                                                                                      >  2:67028/1‑139 (MQ=255)                             aaaCCAAGCGTAAAAAGCTGACCAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCAGCCAGATCTGCGTCCGCTGGTTCCGCTGGATggtggt                                                                                                     >  1:674206/1‑143 (MQ=255)                                     cGTAAAAAGCTGACCAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCAGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTcgc                                                                                             <  1:430231/143‑1 (MQ=255)                                       tAAAAAGCTGACCAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCAGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGa                                                                                           >  1:240028/1‑143 (MQ=255)                                       tAAAAAGCTGACCAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCAGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGa                                                                                           <  1:799834/143‑1 (MQ=255)                                       gaaaaaGCTGACCAAGCGTATCACACTGCTGGAAGCGTTCTTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCAGCCACAGCTTCTTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGa                                                                                           <  2:291918/142‑1 (MQ=255)                                              cTGACCAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCAGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGa                                                                                    <  1:50889/143‑1 (MQ=255)                                              cTGACCAAGCGTATCAAACTGCTGGAAGCGTTCGTGCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCAGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGa                                                                                    <  1:513946/143‑1 (MQ=255)                                                gACCAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGCCCGTTCTGCCGGTACTGCAGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGAcc                                                                                  >  1:209482/1‑143 (MQ=255)                                                 aCCAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGGTCTGCCGGTACTGCAGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGAGACTTGTGGCCt                                                                                 >  2:301872/1‑143 (MQ=255)                                                   cAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCAGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCGGa                                                                               >  1:921867/1‑143 (MQ=255)                                                    aaGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCCGACCGTTCTGCCGGTACTGCAGCCAGATCTGCGTCCGCTGGTTCCGCTGGAGGGTGGTCGTTCCGCGACTTCTGGccatga                                                                              >  2:88599/1‑139 (MQ=255)                                                           tCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCAGCCAGATCTGCGTCCGCTGGTTCCGCTGGCTGGTGGTCGTTTCGCGAATTCTGACCTGAACGATCTg                                                                       >  2:750568/1‑143 (MQ=255)                                                            cAAACTGCCGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCAGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGt                                                                      <  1:92228/143‑1 (MQ=255)                                                                    tgGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCAGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTcgc                                                              <  1:545799/143‑1 (MQ=255)                                                                     ggAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCAGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTcgcg                                                             <  2:407192/143‑1 (MQ=255)                                                                        aGCGTTCTTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCAGCCGGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCa                                                          <  1:940115/143‑1 (MQ=255)                                                                              cgttcAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCAGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAAcc                                                    <  1:335798/143‑1 (MQ=255)                                                                               gttcAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCAGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGCCCTGCACGATCTGTATCGTCGCGTCATGAACCg                                                   >  2:235127/1‑143 (MQ=255)                                                                                ttcAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCAGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGt                                                  >  1:621534/1‑143 (MQ=255)                                                                                  cattCTGGTAACAAGTTAGAGTGGTTGATGCTGATCGTTCTGCCGGTACTGCAGCCAGATCTGCGTTCGTTTGTTCTGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTaa                                                <  2:541076/140‑1 (MQ=255)                                                                                    gTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGAAACAAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCATTTCGCGACTTCTGACCTGAACGATCTGTACCGGCACGTCATGAAcctgagaa                                              >  2:99017/1‑137 (MQ=255)                                                                                     tCTGGTAACGAGCCTGAGTGGATGAGCCTGACCGTTCTGCGGGTACTGCAGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTaacaa                                             <  2:300557/143‑1 (MQ=255)                                                                                         gTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCAGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGt                                         >  1:826416/1‑143 (MQ=255)                                                                                                aCCGGAGTGGATGATCCGGACCGTTCTGCCGGTACTGCAGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAAc                                  <  1:492157/143‑1 (MQ=255)                                                                                                aCCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCAGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAAc                                  <  1:722455/143‑1 (MQ=255)                                                                                                  cAGAGTGGATGTGCTTGACCGGTCTGCCGTTACTGCAGCCAGAGCTGCGTCCGCGGGTTCCGCTGGATGGTGGTCGGTTCGCGAGTTCTGACCGGAGCGATCTGTATCGTCGCGTCATTAACCGTGACAACCGTCTGAAACGt                                <  2:6187/143‑1 (MQ=255)                                                                                                     agTGGATGATCCTGTCCGTTCTGCCGGTACTGCAGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTctg                             <  2:393375/143‑1 (MQ=255)                                                                                                           tgatCCTGACCGTTCTGCCGGTACTGCAGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCAACTTCTGACCTGCAGGATCTGTATCGTCGCCTCCTTAAACGTAACAACCGTCTGAAACGTCTGCAGGaa                       >  2:61804/1‑142 (MQ=255)                                                                                                            gatCCTGACCGTTCTGCCGGTACTGCAGCCAGAACTGCGTCCGCTGGTTCCGCTGGATGGTGATCGTTTCGCGACTTCTGACCTTAACGATCTGTCTCGTCGCGTCATTAACCGTCACAAGCGTTTGCAACGTCTGCTGGATc                      >  1:631030/1‑143 (MQ=255)                                                                                                                cTGACCGTTCTGCCGGTACTGCAGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAATGATCTGTATCGTCGCGTCATTAACCGAAACAACCGTCTGAAACGTCTGCGGGAttttgc                  >  1:671143/1‑138 (MQ=255)                                                                                                                       ttCTGCCGGTACTGCAGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCg           <  2:842999/143‑1 (MQ=255)                                                                                                                         cggCCGGTACTGCAGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGa         <  1:525575/141‑1 (MQ=255)                                                                                                                           gCCGGTACTGCAGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCACGACTTCTGACCTGAACGATCTGTATCGTCGAGTCATTAACCGTAACAACAGTCTGAAACGTATGCTGTGTCTTAATGAGCCGGAca       >  2:727825/1‑143 (MQ=255)                                                                                                                              ggTACTGCAGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGCACGCCCTGTATCGTCGCGTCATTCACTGAAAAAATCGTCCGAACCGACTGCTGGCTCTGGCTGCGCCGGAcatca    >  2:567280/1‑143 (MQ=255)                                                                                                                                tACTGCAGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTCTCGCTACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACatcatc  <  2:674206/143‑1 (MQ=255)                                                                                                                                      |                                                                                                                                         CGTGAAGAGCTGAACGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATC  >  NC_000913/4185953‑4186222

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