Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,183,238	C→T	A665V (GCG→GTG)	rpoB →	RNA polymerase, beta subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,183,238	0	C	T	100.0%	119.6 / NA	35	A665V (GCG→GTG)	rpoB	RNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (16/19);  total (16/19)

GAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGT  >  NC_000913/4183117‑4183362                                                                                                                          |                                                                                                                             gaagaaGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGTGTCCCTGATCCCGTTcc                                                                                                             <  1:365773/139‑1 (MQ=255)   agaagGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGgtgtgt                                                                                                                            >  1:623416/1‑122 (MQ=255)   agaagGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGgtgtgt                                                                                                                            <  2:623416/122‑1 (MQ=255)                 tAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGTGTCCCTGATCCCGTTCCTGGAACACGATGACGc                                                                                             <  2:468733/139‑1 (MQ=255)                    aagaCCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGTGTCCCTGATCCCGTTCCTGGAACACGATGACGCCaa                                                                                          <  1:269023/139‑1 (MQ=255)                       aCCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGTGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCg                                                                                       >  1:585319/1‑139 (MQ=255)                        ccTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGTGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGt                                                                                      <  2:652127/139‑1 (MQ=255)                         cTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGTGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTg                                                                                     >  1:50414/1‑139 (MQ=255)                                   gCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGTGTCCCTGATCCCGTTCCTGGAac                                                                                                       >  1:699790/1‑111 (MQ=255)                                   gCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGTGTCCCTGATCCCGTTCCTGGAac                                                                                                       <  2:699790/111‑1 (MQ=255)                                         gCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGTGTCCCTGATCCCGTTCCTGGAACACGATGACGCCaa                                                                                          >  1:368763/1‑118 (MQ=255)                                         gCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGTGTCCCTGATCCCGTTCCTGGAACACGATGACGCCaa                                                                                          <  2:368763/118‑1 (MQ=255)                                              ggCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGTGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATgg                                                                >  1:570939/1‑138 (MQ=255)                                                      cAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGTGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAg                                                        <  2:79467/139‑1 (MQ=255)                                                               cAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGTGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCg                                               <  1:365111/139‑1 (MQ=255)                                                                       aCCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGTGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGCTGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTgcg                                       >  2:109792/1‑139 (MQ=255)                                                                                aCTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGTGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATg                                                                 >  1:281417/1‑104 (MQ=255)                                                                                aCTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGTGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATg                                                                 <  2:281417/104‑1 (MQ=255)                                                                                    cATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGTGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCg                          <  1:109792/139‑1 (MQ=255)                                                                                    cATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGTGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCg                          <  1:223018/139‑1 (MQ=255)                                                                                      tGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGTGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCt                        >  1:300767/1‑139 (MQ=255)                                                                                       ggACGTATCCACCCAGCAGGTGGTATCCGTCGGTGTGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTCCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCtg                       >  1:573088/1‑139 (MQ=255)                                                                                         aCGTATCCACCTAGCAGGTGGTATCCGTCGGTGTGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCtggt                     <  2:520478/139‑1 (MQ=255)                                                                                         aCGTATCCACCCAGCAGGTGGTATCCGTCGGTGTGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAAcc                                                                                        >  2:716970/1‑72 (MQ=255)                                                                                         aCGTATCCACCCAGCAGGTGGTATCCGTCGGTGTGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAAcc                                                                                        <  1:716970/72‑1 (MQ=255)                                                                                         aCGTATCCACCCAGCAGGTGGTATCCGTCGGTGTGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCtggt                     <  2:99985/139‑1 (MQ=255)                                                                                         aCGTATCCACCCAGCAGGTGGTATCCGTCGGTGTGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCtggt                     <  2:156014/139‑1 (MQ=255)                                                                                          cGTATCCACCCAGCAGGTGGTATCCGTCGGTGTGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCtggtt                    <  1:310033/139‑1 (MQ=255)                                                                                          cGTATCCACCCAGCAGGTGGTATCCGTCGGTGTGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCtggtt                    >  1:2482/1‑139 (MQ=255)                                                                                            tATCCACCCAGCAGGTGGTATCCGTCGGTGTGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCtggttgg                  >  2:40752/1‑139 (MQ=255)                                                                                              tCCACCCAGCAGGTGGTATCCGTCGGTGTGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTa                <  1:669253/139‑1 (MQ=255)                                                                                              tCCACCCAGCAGGTGGTATCCGTCGGTGTGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTa                <  1:741003/139‑1 (MQ=255)                                                                                              tCCACCCAGCAGGTGGTATCCGTCGGTGTGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTa                >  2:741003/1‑139 (MQ=255)                                                                                               ccACCCAGCAGGTGGTATCCGTCGGTGTGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACg                                                            >  2:105604/1‑94 (MQ=255)                                                                                               ccACCCAGCAGGTGGTATCCGTCGGTGTGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACg                                                            <  1:105604/94‑1 (MQ=255)                                                                                                     agcagGTGGTATCCGTCGGTGTGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTAt         >  1:590857/1‑139 (MQ=255)                                                                                                            ggtATCCGTCGGTGTGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGt  >  1:727869/1‑139 (MQ=255)                                                                                                                          |                                                                                                                             GAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGT  >  NC_000913/4183117‑4183362

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