Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,187,647	T→C	G766G (GGT→GGC)	rpoC →	RNA polymerase subunit beta'

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,187,647	0	T	C	100.0%	141.4 / NA	41	G766G (GGT→GGC)	rpoC	RNA polymerase subunit beta'
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (17/24);  total (17/24)

CATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCA  >  NC_000913/4187512‑4187763                                                                                                                                        |                                                                                                                     cATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGCCTGa                                                                                                                  <  1:498895/140‑1 (MQ=255)  aTCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGCCTGaa                                                                                                                 >  1:45912/1‑140 (MQ=255)  aTCTACATGATGGCCGACTCCGGTGAGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGCCTGaa                                                                                                                 >  2:555568/1‑140 (MQ=255)    cTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGCCTGAACg                                                                                                               <  2:92668/140‑1 (MQ=255)      aCATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGCCTGAACGTa                                                                                                             <  2:45912/140‑1 (MQ=255)      aCATGATGGCCGACTCCGGTGCCCGCGGTCCTGCGGCTCAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGCCTGAACGTa                                                                                                             <  1:371844/140‑1 (MQ=255)           atgGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGCCTGAACg                                                                                                               >  1:193735/1‑133 (MQ=255)           atgGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGCCTGAACg                                                                                                               <  2:193735/133‑1 (MQ=255)           atgGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGCCTGAACGTACTCCa                                                                                                        >  1:317326/1‑140 (MQ=255)           atgGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCACCGAAACGCTAATGACCGCGAACTTCGGTGAGGGCCTGAACGTACTCCa                                                                                                        >  2:232429/1‑140 (MQ=255)            tgGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGCCTGAACg                                                                                                               >  2:779430/1‑132 (MQ=255)            tgGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGCCTGAACg                                                                                                               <  1:779430/132‑1 (MQ=255)            tgGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGCCTGAACGTACTCCAg                                                                                                       >  2:151953/1‑140 (MQ=255)                 gACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGCCTGAACGTACTCCAGTACtt                                                                                                  <  1:55233/140‑1 (MQ=255)                     ccGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGCCTGAACGTACTCCAGTACTTCAtc                                                                                              <  1:151953/140‑1 (MQ=255)                     ccGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGCCTGAACGTACTCCAGTACTTCAtc                                                                                              <  2:317326/140‑1 (MQ=255)                       ggTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGCCTGAACGTACTCCAGTACTTCAtctc                                                                                            <  2:764072/140‑1 (MQ=255)                             cgTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGGAACGCCAATCACCGCGAACTTCCGTGAAGGCCTGAACGTACTCCAGTACTTCATCTccaccca                                                                                      >  1:696350/1‑140 (MQ=255)                               tGGTTCTGCGGCACAGATTCGTCAGATTGCTGGTATGTGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGCCTGAACGTACTCCAGTACTTCATCTCCACCCACg                                                                                    <  1:395944/140‑1 (MQ=255)                                  ttCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGCCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTg                                                                                 <  1:555568/140‑1 (MQ=255)                                  ttCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGCCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTg                                                                                 <  2:726374/140‑1 (MQ=255)                                    cTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGCCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCt                                                                               >  1:654462/1‑140 (MQ=255)                                        ggCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGCCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTa                                                                           <  1:550793/140‑1 (MQ=255)                                         gCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGCCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTaa                                                                          <  2:469787/140‑1 (MQ=255)                                             aGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGCCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGt                                                                      <  2:621112/140‑1 (MQ=255)                                             aGATTCGTCAGCTTGCTGGTATGCGCGGTCTGATGGCCAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGCCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGt                                                                      <  1:220580/140‑1 (MQ=255)                                                         ttGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGCCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATAcc                                                          <  2:517331/140‑1 (MQ=255)                                                           gCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGCCTGAACGTACTCCAGTACTTCATCGCCACCCACGGTACTGGTAAAGGTCTGGCGGATACCGc                                                        >  1:367048/1‑140 (MQ=255)                                                           gCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGCCTGAACGTACTCCAGTACTTCATCCCCACCCAGGGTACTCGTAAAGGTCTGGCGGATACCGc                                                        >  1:299712/1‑140 (MQ=255)                                                                  tGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGCCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGaaa                                                 >  1:15171/1‑140 (MQ=255)                                                                    cGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGCCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAAc                                               >  2:136121/1‑140 (MQ=255)                                                                         tCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGCCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCGCGTAAAGGTCTAGCTGATACCGCACTGAAAACTGCGa                                          >  1:235793/1‑140 (MQ=255)                                                                          cTGATGGCGAAGCCGGACGGCTGCATCATCGCAACCCCAATCACCGCGTACTTCCGTGAAGGCCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGaa                                         <  2:235793/140‑1 (MQ=255)                                                                              tGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGCCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTcc                                     <  2:15171/140‑1 (MQ=255)                                                                                gCGAAGCCGGATGGCACCATCAACGCAACCCCAATCACCGCGAACTTCCGTGAAGGCCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCgg                                   <  1:232429/140‑1 (MQ=255)                                                                                     gCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGCCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTAcc                              >  1:805614/1‑140 (MQ=255)                                                                                     gCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGCCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTAGTAAAGGTTTGGGGGAAACCGCACTGAAAACTGCGAACTCCGGTTAcc                              >  2:289359/1‑140 (MQ=255)                                                                                                          aaCGCCAATCACCGCGAACTTCCGTGAAGGCCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACg         <  2:207923/140‑1 (MQ=255)                                                                                                          aaCGCCAATCACCGCGAACTTCCGTGAAGGCCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACg         <  1:327097/140‑1 (MQ=255)                                                                                                                 aTCACCGCGAACTTCCGTGAAGGCCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCa  >  1:696402/1‑140 (MQ=255)                                                                                                                 aTCACCGCGAACTTCCGTGAAGACCCGAACGTACTCCCGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCa  <  1:267247/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                     CATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCA  >  NC_000913/4187512‑4187763

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