Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_bme_pgi	3,440,086	C→A	G315V (GGC→GTC)	rpoA ←	RNA polymerase, alpha subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_bme_pgi	3,440,086	0	C	A	100.0%	122.7 / NA	35	G315V (GGC→GTC)	rpoA	RNA polymerase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (17/18);  total (17/18)

TTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTT  >  NC_000913_3_bme_pgi/3439962‑3440228                                                                                                                             |                                                                                                                                               ttGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAg                                                                                                                       <  1:92750/150‑1 (MQ=255)      gTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGt                                                                                                                  <  4:116224/150‑1 (MQ=255)       ttCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCagag                                                                                                                                          >  2:277357/1‑125 (MQ=255)       ttCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCagag                                                                                                                                          <  1:277357/125‑1 (MQ=255)                         ttACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCa                                                                                                       <  2:282805/142‑1 (MQ=255)                         ttACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCa                                                                                                       >  1:282805/1‑142 (MQ=255)                                      tGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAgg                                                                                  >  2:208629/1‑150 (MQ=255)                                             aTCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAAtc                                                                                                          >  4:261592/1‑119 (MQ=255)                                             aTCCTTCTCAGTAAAACCTTAACATGTGAGCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAAtc                                                                                                          <  3:261592/119‑1 (MQ=255)                                                 ttCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAg                                                                                                                       <  2:333684/102‑1 (MQ=255)                                                 ttCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAg                                                                                                                       >  1:333684/1‑102 (MQ=255)                                                    tcAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAg                                                                    <  1:158238/150‑1 (MQ=255)                                                      aGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGa                                                                  >  2:9845/1‑150 (MQ=255)                                                      aGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGa                                                                  <  2:162722/150‑1 (MQ=255)                                                              cTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAAtc                                                                                                          <  2:256045/102‑1 (MQ=255)                                                              cTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAAtc                                                                                                          >  1:256045/1‑102 (MQ=255)                                                                               ttACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAg                                                                 <  3:69934/126‑1 (MQ=255)                                                                               ttACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAg                                                                 >  4:69934/1‑126 (MQ=255)                                                                               ttACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAg                                         <  1:208629/150‑1 (MQ=255)                                                                               ttACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAg                                         <  1:9845/150‑1 (MQ=255)                                                                                 aCTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGAt                                       <  4:121766/150‑1 (MQ=255)                                                                                        aGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCa                                                                                                       <  4:28347/79‑1 (MQ=255)                                                                                        aGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCa                                                                                                       >  3:28347/1‑79 (MQ=255)                                                                                                  gCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAc                                    >  1:456682/1‑136 (MQ=255)                                                                                                  gCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAc                                    <  2:456682/136‑1 (MQ=255)                                                                                                      gTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCtt                  >  1:509648/1‑150 (MQ=255)                                                                                                          ccAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGc              <  2:509648/150‑1 (MQ=255)                                                                                                                    aGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAg    >  4:177636/1‑150 (MQ=255)                                                                                                                     ggCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTc                                                              <  3:141875/91‑1 (MQ=255)                                                                                                                     ggCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTc                                                              >  4:141875/1‑91 (MQ=255)                                                                                                                      gcgcATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGtt  >  4:17806/1‑150 (MQ=255)                                                                                                                      gcgcATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGtt  >  1:226002/1‑150 (MQ=255)                                                                                                                      gcgcATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGtt  >  4:4837/1‑150 (MQ=255)                                                                                                                        gcATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATc                                      >  1:467354/1‑112 (MQ=255)                                                                                                                        gcATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATc                                      <  2:467354/112‑1 (MQ=255)                                                                                                                             |                                                                                                                                               TTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTT  >  NC_000913_3_bme_pgi/3439962‑3440228

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