Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,188,003	T→C	V885A (GTT→GCT)	rpoC →	RNA polymerase subunit beta'

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,188,003	0	T	C	100.0%	146.4 / NA	42	V885A (GTT→GCT)	rpoC	RNA polymerase subunit beta'
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (17/25);  total (17/25)

CGCGCAACACGCTGCTGCACGAACAGTGGTGTGACCTGCTGGAAGAGAACTCTGTCGACGCGGTTAAAGTACGTTCTGTTGTATCTTGTGACACCGACTTTGGTGTATGTGCGCACTGCTACGGTCGTGACCTGGCGCGTGGCCACATCATCAAC  >  NC_000913/4187925‑4188079                                                                               |                                                                             cgcgCAACACGCTGCTGCACGAACAGTGGTGTGACCTGCTGGAAGAGAACTCTGTCGACGCGGTTAAAGTACGTTctgctg                                                                            <  2:852807/81‑1 (MQ=255)    gcAACACGCTGCTGCACGAACAGTGGTGTGACCTGCTGGAAGAGAACTCTGTCGACGCGGTTAAAGTACGTTCTGCTGTAt                                                                         >  2:852863/1‑81 (MQ=255)     cAACACGCTGCTGCACGAACAGTGGTGTGACCTGCTGGAAGAGAACTCTGTCGACGCGGTTAAAGTACGTTCTGCTGTATc                                                                        <  2:852808/81‑1 (MQ=255)      aaCACGCTGCTGCACGAACAGTGGTGTGACCTGCTGGAAGAGAACTCTGTCGACGCGGTTAAAGTACGTTCTGCTGTATCt                                                                       <  1:852809/81‑1 (MQ=255)       acacGCTGCTGCACGAACAGTGGTGTGACCTGCTGGAAGAGAACTCTGTCGACGCGGTTAAAGTACGTTCTGCTGTATCtt                                                                      >  2:852870/1‑81 (MQ=255)       acacGCTGCTGCACGAACAGTGGTGTGACCTGCTGGAAGAGAACTCTGTCGACGCGGTTAAAGTACGTTCTGCTGTATCtt                                                                      <  1:852810/81‑1 (MQ=255)       acacGCTGCTGCACGAACAGTGGTGTGACCTGCTGGAAGAGAACTCTGTCGACGCGGTTAAAGTACGTTCTGCTGTATCtt                                                                      <  1:852811/81‑1 (MQ=255)        cacGCTGCTGCACGAACAGTGGTGTGACCTGCTGGAAGAGAACTCTGTCGACGCGGTTAAAGTACGTTCTGCTGTATCTtg                                                                     >  1:852871/1‑81 (MQ=255)          cGCTGCTGCACGAACAGTGGTGTGACCTGCTGGAAGAGAACTCTGTCGACGCGGTTAAAGTACGTTCTGCTGTATCTtgtg                                                                   <  2:852812/81‑1 (MQ=255)            ctgctgCACGAACAGTGGTGTGACCTGCTGGAAGAGAACTCTGTCGACGCGGTTAAAGTACGTTCTGCTGTATCTTGTGa                                                                  >  2:852866/1‑80 (MQ=255)            ctgatgCACGAACAGTGGTGTGAACTGATGGAAGAGAACGCTGTAGACGCGGTTAAAGTACGTTCTGCTGTATCTTGTGaa                                                                 <  1:1038017/81‑2 (MQ=255)               ctgcACGAACAGTGGTGTGACCTGCTGGAAGAGAACTCTGTCGACGCGGTTAAAGTACGTTCTGCTGTATCTTGTGaca                                                                <  2:852813/79‑1 (MQ=255)               ctgcACGAACAGTGGTGTGACCTGCTGGAAGAGAACTCTGTCGACGCGGTTAAAGTACGTTCTGCTGTATCTTGTGACAcc                                                              >  1:852874/1‑81 (MQ=255)                    cGAACAGTGGTGTGACCTGCTGGAAGAGAACTCTGTCGACGCGGTTAAAGTACGTTCTGCTGTATCTTGTGACACCGACtt                                                         <  2:852815/81‑1 (MQ=255)                    cGAACAGTGGTGTGACCTGCTGGAAGAGAACTCTGTCGACGCGGTTAAAGTACGTTCTGCTGTATCTTGTGACACCGACtt                                                         <  2:852814/81‑1 (MQ=255)                     gAACAGTGGTGTGACCTGCTGGAAGAGAACTCTGTCGACGCGGTTAAAGTACGTTCTGCTGTATCTTGTGACACCGACttt                                                        >  2:852873/1‑81 (MQ=255)                                tgACCTGCTGGAAGAGAACTCTGTCGACGCGGTTAAAGTACGTTCTGCTGTATCTTGTGACACCGACTTTGGTGTATGTgc                                             <  1:852816/81‑1 (MQ=255)                                  aCCTGCTGGAAGAGAACTCTGTCGACGCGGTTAAAGTACGTTCTGCTGTATCTTGTGACACCGACTTTGGTGTATGTgcg                                            <  1:852817/80‑1 (MQ=255)                                      gctgGAAGAGAACTCTGTCGACGCGGTTAAAGTACGTTCTGCTGTATCTTGTGACACCGACTTTGGTGTATGTGCGCACTg                                       >  2:852884/1‑81 (MQ=255)                                       ctgGAAGAGAACTCTGTCGACGCGGTTAAAGTACGTTCTGCTGTATCTTGTGACACCGACTTTGGTGTATGTGCGCACt                                        <  1:852818/79‑1 (MQ=255)                                       ctgGAAGAGAACTCTGTCGACGCGGTTAAAGTACGTTCTGCTGTATCTTGTGACACCGACTTTGGTGTATGTGCGCACt                                        <  2:852819/79‑1 (MQ=255)                                       ctgGAAGAGAACTCTGTCGACGCGGTTAAAGTACGTTCTGCTGTATCTTGTGACACCGACTTTGGTGTATGTGCGCACTGc                                      >  1:852885/1‑81 (MQ=255)                                        tgGAAGAGAACTCTGTCGACGCGGTTAAAGTACGTTCTGCTGTATCTTGTGACACCGACTTTGGTGTATGTGCGCACTGc                                      >  1:852888/1‑80 (MQ=255)                                          gAAGAGAACTCTGTCGACGCGGTTAAAGTGCGTTCTGCTGTATCTTGTGACACCGACTTTGGTGTATGTGCGCACTGCt                                     <  2:852820/79‑1 (MQ=255)                                             gagaACTCTGTCGACGCGGTTAAAGTACGTTCTGCTGTATCTTGTGACACCGACTTTGGTGTATGTGCGCACTGCTACGGt                                <  1:852821/81‑1 (MQ=255)                                             gagaACTCTGTCGACGCGGTTAAAGTACGTTCTGCTGTATCTTGTGACACCGACTTTGGTGTATGTGCGCACTGCTACGGt                                <  2:852822/81‑1 (MQ=255)                                              agaACTCTGTCGACGCGGTTAAAGTACGTTCTGCTGTATCTTGTGACACCGACTTTGGTGTATGTGCGCACTGCTACGGTc                               >  1:852895/1‑81 (MQ=255)                                              agaACTCTGTCGACGCGGTTAAAGTACGTTCTGCTGTATCTTGTGACACCGACTTTGGTGTATGTGCGCACTGCTACGGTc                               <  2:852823/81‑1 (MQ=255)                                               gaACTCTGTCGACGCGGTTAAAGTACGTTCTGCTGTATCTTGTGACACCGACTTTGGTGTATGTGCGCACTGCTACGGTCg                              >  1:852853/1‑81 (MQ=255)                                                aaCTCTGTCGACGCGGTTAAAGTACGTTCTGCTGTATCTTGTGACACCGACTTTGGTGTATGTGCGCACTGCTACGGTCGt                             <  2:852824/81‑1 (MQ=255)                                                 aCTCTGTCGACGCGGTTAAAGTACGTTCTGCTGTATCTTGTGACACCGACTTTGGTGTATGTGCGCACTGCTACGGTCGTg                            <  2:852827/81‑1 (MQ=255)                                                 aCTCTGTCGACGCGGTTAAAGTACGTTCTGCTGTATCTTGTGACACCGACTTTGGTGTATGTGCGCACTGCTACGGTCGTg                            <  1:852826/81‑1 (MQ=255)                                                 aCTCTGTCGACGCGGTTAAAGTACGTTCTGCTGTATCTTGTGACACCGACTTTGGTGTATGTGCGCACTGCTACGGTCGTg                            <  1:852825/81‑1 (MQ=255)                                                     tGTCGACGCGGTTAAAGTACGTTCTGCTGTATCTTGTGACACCGACTTTGGTGTATGTGCGCACTGCTACGGTCGTGAc                          <  2:852828/79‑1 (MQ=255)                                                      gTCGACGCGGTTAAAGTACGTTCTGCTGTATCTTGTGACACCGACTTTGGTGTATGTGCGCACTGCTACGGTCGTGACCTg                       >  2:852881/1‑81 (MQ=255)                                                       tCGACGCGGTTAAAGTACGTTCTGCTGTATCTTGTGACACCGACTTTGGTGTATGTGCGCACTGCTACGGTCGTGAcc                         <  1:852829/78‑1 (MQ=255)                                                       tCGACGCGGTTAAAGTACGTTCTGCTGTATCTTGTGACACCGACTTTGGTGTATGTGCGCACTGCTACGGTCGTGACCTgg                      <  1:852830/81‑1 (MQ=255)                                                         gACGCGGTTAAAGTACGTTCTGCTGTATCTTGTGACACCGACTTTGGTGTATGTGCGCACTGCTACGGTCGTGACCTGgcg                    >  2:852883/1‑81 (MQ=255)                                                         gACGCGGTTAAAGTACGTTCTGCTGTATCTTGTGACACCGACTTTGGTGTATGTGCGCACTGCTACGGTCGTGACCTGgc                     <  2:852831/80‑1 (MQ=255)                                                                   aaGTACGTTCTGCTGTATCTTGTGACACCGACTTTGGTGTATGTGCGCACTGCTACGGTCGTGACCTGGCGCGTGGCcaca          >  1:852875/1‑81 (MQ=255)                                                                    aGTACGTTCTGCTGTATCTTGTGACACCGACTTTGGTGTATGTGCGCACTGCTACGGTCGTGACCTGGCGCGTGGCCAcat         >  1:852889/1‑81 (MQ=255)                                                                       aCGTTCTGCTGTATCTTGTGACACCGACTTTGGTGTATGTGCGCACTGCTACGGTCGTGACCTGGCGCGTGGCCAcatcat      >  1:852901/1‑81 (MQ=255)                                                                           tCTGCTGTATCTTGTGACACCGACTTTGGTGTATGTGCGCACTGCTACGGTCGTGACCTGGCGCGTGGCCACATCATcaa   >  2:852899/1‑80 (MQ=255)                                                                            ctgctgTATCTTGTGACACCGACTTTGGTGTATGTGCGCACTGCTACGGTCGTGACCTGGCGCGTGGCCACATCATcaac  <  2:852832/80‑1 (MQ=255)                                                                            ctgctgTATCTTGTGACACCGACTTTGGTGTATGTGCGCACTGCTACGGTCGTGACCTGGCGCGTGGCCACATCATcaac  >  2:852900/1‑80 (MQ=255)                                                                               |                                                                             CGCGCAACACGCTGCTGCACGAACAGTGGTGTGACCTGCTGGAAGAGAACTCTGTCGACGCGGTTAAAGTACGTTCTGTTGTATCTTGTGACACCGACTTTGGTGTATGTGCGCACTGCTACGGTCGTGACCTGGCGCGTGGCCACATCATCAAC  >  NC_000913/4187925‑4188079

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