Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,187,214	A→G	D622G (GAC→GGC)	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,187,214	0	A	G	100.0%	150.1 / NA	42	D622G (GAC→GGC)	rpoC	RNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (21/21);  total (21/21)

GATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTC  >  NC_000913/4187089‑4187358                                                                                                                              |                                                                                                                                                 gatgatTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTacac                                                                                                                                   <  5:105877/141‑1 (MQ=255) gatgatTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTacac                                                                                                                                   >  6:105877/1‑141 (MQ=255) gatgatTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCg                                                                                                                           >  6:44325/1‑149 (MQ=255)     atTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTa                                                                                                                       >  2:129547/1‑149 (MQ=255)     atTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTa                                                                                                                       <  2:15937/149‑1 (MQ=255)           ccGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAgc                                                                                                                 >  4:68464/1‑149 (MQ=255)           ccGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAgc                                                                                                                 <  3:68464/149‑1 (MQ=255)              aaaGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGcgcg                                                                                                              >  4:56220/1‑149 (MQ=255)                       ccTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTgg                                                                                                        <  5:44325/146‑1 (MQ=255)                             tCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGt                                                                                               >  3:69304/1‑149 (MQ=255)                                   gTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTAt                                                                                         >  4:65895/1‑149 (MQ=255)                                    tCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACAcc                                                                                                                                  >  4:33347/1‑107 (MQ=255)                                    tCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACAcc                                                                                                                                  <  3:33347/107‑1 (MQ=255)                                       aCCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGAt                                                                                     <  3:56220/149‑1 (MQ=255)                                             cgcTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATg                                                                               <  3:65895/149‑1 (MQ=255)                                                   gTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATc                                                                         >  3:121617/1‑149 (MQ=255)                                                   gTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATc                                                                         >  5:6956/1‑149 (MQ=255)                                                    tAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATcc                                                                        >  7:16789/1‑149 (MQ=255)                                                        aaaGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGga                                                                    <  6:6956/149‑1 (MQ=255)                                                             aaTCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGagaaga                                                               <  1:129547/149‑1 (MQ=255)                                                              aTCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGATCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCTGagaaga                                                               >  4:93697/1‑148 (MQ=255)                                                                   cAAACTGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGGTGGTATCGATGACATGGTCATCCCGGAGAAGaaa                                                             <  3:93697/145‑1 (MQ=255)                                                                   cAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAaca                                                           <  4:121617/147‑1 (MQ=255)                                                                           tGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTg                                                                                                         >  2:73782/1‑93 (MQ=255)                                                                           tGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTg                                                                                                         <  1:73782/93‑1 (MQ=255)                                                                             aaCACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCCTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCAtctc                                               <  8:16789/149‑1 (MQ=255)                                                                              acacCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGaaa                                                      <  2:5826/141‑1 (MQ=255)                                                                              acacCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGaaa                                                      >  1:5826/1‑141 (MQ=255)                                                                                  cTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTAt                                                                                                                      <  1:156340/73‑1 (MQ=255)                                                                                  cTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTAt                                                                                                                      >  2:156340/1‑73 (MQ=255)                                                                                  cTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTgg                                                                                                        <  2:27884/87‑1 (MQ=255)                                                                                  cTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTgg                                                                                                        >  1:27884/1‑87 (MQ=255)                                                                                  cTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGAt                                                                                     >  5:62975/1‑106 (MQ=255)                                                                                  cTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGAt                                                                                     <  6:62975/106‑1 (MQ=255)                                                                                        ccGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGgcagaag                                    >  4:108057/1‑149 (MQ=255)                                                                                             ttCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGgcagaagcag                                 >  8:65314/1‑147 (MQ=255)                                                                                                          accgaccgTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCa                                                                           >  5:9255/1‑92 (MQ=255)                                                                                                          accgaccgTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCa                                                                           <  6:9255/92‑1 (MQ=255)                                                                                                            cgaccgTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAg                <  4:63448/149‑1 (MQ=255)                                                                                                              accgTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAg                >  6:71676/1‑147 (MQ=255)                                                                                                                       tttGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTcc      <  7:65314/148‑1 (MQ=255)                                                                                                                           cGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAgtc  <  4:69304/148‑1 (MQ=255)                                                                                                                              |                                                                                                                                                 GATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTC  >  NC_000913/4187089‑4187358

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