Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,184,543	C→A	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_000913	4,184,543	0	C	A	100.0%	122.9 / NA	39	P1100Q (CCG→CAG)	rpoB	RNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (16/23);  total (16/23)

GGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAG  >  NC_000913/4184406‑4184682                                                                                                                                          |                                                                                                                                            ggtggTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTgg                                                                                                                                        <  1:569640/140‑1 (MQ=255) ggCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGCTCGTCACGGCAGCAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTgg                                                                                                                                        <  1:243218/143‑1 (MQ=255)   cGGTTAATCGCCGTATTCAGTCTTGTGACAAGATGGCAGGACGTCACGGTAACAAGGGTGTAATGTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGc                                                                                                                                      <  1:328526/143‑1 (MQ=255)      ttAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAGCCAGCTGGGCGTa                                                                                                                                   >  1:924583/1‑143 (MQ=255)        acatgCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCACCCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTAcc                                                                                                                                 <  1:863906/139‑1 (MQ=255)          aCGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGCCATCGTACTGAACCAGCTGGGCGTACCGt                                                                                                                               <  1:936605/143‑1 (MQ=255)          aCGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGt                                                                                                                               >  2:933643/1‑143 (MQ=255)                  tacAGCCGGGTGACAACATGGCAGGTCGTCACGGTAACAAGCGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATg                                                                                                                       <  1:683512/141‑1 (MQ=255)                         tGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCGACCCGATCGACGATATGCCTTACGCTGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTCTGAACATCg                                                                                                                >  1:321772/1‑143 (MQ=255)                                aaGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCATCTGGGCGTACCGTCTCGTATGAACATCGGTCAGAt                                                                                                         >  1:52528/1‑143 (MQ=255)                                        aGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGaaa                                                                                                 <  1:287289/143‑1 (MQ=255)                                           tcgtcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAAccc                                                                                              <  2:913801/143‑1 (MQ=255)                                            cgtcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTCTGAACGTCGGTCAGATCCTCGGACCCCa                                                                                             >  1:481969/1‑143 (MQ=255)                                              tcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCGcc                                                                                           >  2:819146/1‑143 (MQ=255)                                                       cacgggTGTAATTTCTACGACCAACCGGATCGAAGGTATGCCTGGCGATGAAGACGGTGCTCCTGTAGGCATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATgg                                                                                  <  2:930732/140‑1 (MQ=255)                                                               tAATTTCTAAGCTCACCCCGATCGAAGATATGCCGTACGACGACAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGaaa                                                                          <  2:55213/143‑1 (MQ=255)                                                                          gtCACCCCGATCGAAGATATGCCTTACGCTGAGAACGGTACGCGGGTAGACATCGTACTGAACCAGCTGGGCGCACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGa                                                               <  2:764213/142‑1 (MQ=255)                                                                          gccAACCCGAGCGAGGACAGGCCTGGCGATGAAGACGGTACGCCGGTAGACAGCTGACTGAACCAGCTGGGCGTACCGTCGCGTATGAGCATCGGTCAGAGCCTGGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGa                                                               <  2:731572/141‑1 (MQ=255)                                                                          accAACCCGATCGAAGATAGGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGa                                                               <  2:209961/141‑1 (MQ=255)                                                                          aTCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGGACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCCACTGGGTAGGGCTGCGAAAGGTGTCGGCGa                                                               >  1:413164/1‑143 (MQ=255)                                                                             aaCCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACaa                                                            <  1:56105/143‑1 (MQ=255)                                                                                           gatGCCTTACGATGAAAACGGTACGCCGGTAGACAGCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGc                                              <  1:534833/142‑1 (MQ=255)                                                                                              gCCTTACGATGAAAACGGTACGCCGGTAGGCATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCAACTGGGTATGGCTGCGAAAGGTATCGGCGACAGGATCAACGACATGCTGa                                           >  1:628/1‑143 (MQ=255)                                                                                              gCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCGTGCTGa                                           >  2:789252/1‑143 (MQ=255)                                                                                                cTTACGTTGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGTGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGaaa                                         <  1:616961/143‑1 (MQ=255)                                                                                                cTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGaaa                                         <  1:131249/143‑1 (MQ=255)                                                                                                   aCGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATACTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAAcag                                      >  1:493535/1‑143 (MQ=255)                                                                                                    cGATGACATCGGTACGCCGGTAGACATCGTACTGATCCAGCTGGGCGTACCGGCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAAcagc                                     <  2:417825/143‑1 (MQ=255)                                                                                                       tGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACagcagc                                  <  1:49159/143‑1 (MQ=255)                                                                                                            ccGGTTCGCCGGTGGACATCGCTCTGAACCAGCTGGGCGTACCGCCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCaagaa                             <  2:903310/142‑1 (MQ=255)                                                                                                                 aCGCAGCTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTcgc                        <  1:657375/143‑1 (MQ=255)                                                                                                                      ggTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAAc                   <  1:933643/143‑1 (MQ=255)                                                                                                                          gACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGGATCGGCGACAAGATCAACACCATGCTGAAACAGCAGAAAGAAGGCGGGAAACGgcg               >  2:316146/1‑143 (MQ=255)                                                                                                                             aTCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCGACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGGTGAAACAGCAGCAAGAAGTCGCGAAACTGGGCGa            >  1:409274/1‑143 (MQ=255)                                                                                                                             aTCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGa            >  1:639160/1‑143 (MQ=255)                                                                                                                              tCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGaa           <  1:113849/143‑1 (MQ=255)                                                                                                                                  aCTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCCAGACGTCGCGCAACTGCGCGACTTCa       >  1:329725/1‑143 (MQ=255)                                                                                                                                    tGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCACGAAACCCACCTGGGTATGGCGGCGAAAGGGATCGGCGACAAGATCAACGCCATGCGGAAACAGCAGCAAGAAGGCGCGAAACTGCGCGAGTTGATc     >  1:327855/1‑143 (MQ=255)                                                                                                                                       aCCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAg  >  1:717739/1‑143 (MQ=255)                                                                                                                                          |                                                                                                                                            GGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAG  >  NC_000913/4184406‑4184682

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