Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,862,634	A→G	intergenic (‑205/‑137)	msrB ← / → gapA	methionine sulfoxide reductase B/glyceraldehyde‑3‑phosphate dehydrogenase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,862,634	0	A	G	100.0%	109.1 / NA	31	intergenic (‑205/‑137)	msrB/gapA	methionine sulfoxide reductase B/glyceraldehyde‑3‑phosphate dehydrogenase A
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (13/18);  total (13/18)

GTCCTAAACTTGATCTCGACGAAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGAGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCT  >  NC_000913/1862509‑1862760                                                                                                                              |                                                                                                                               gTCCTAAACTTGATCTCGACGAAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGagtcag                                                                                                                      <  2:103818/136‑1 (MQ=255) gTCCTAAACTTGATCTCGACGAAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGagtcag                                                                                                                      >  1:103818/1‑136 (MQ=255)          ttGATCTCGACGAAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCt                                                                                                          <  1:723689/139‑1 (MQ=255)           tGATCTCGACGAAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCtt                                                                                                         <  1:334648/139‑1 (MQ=255)                 cgacgaAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGcac                                                                                                   >  2:602120/1‑139 (MQ=255)                   acgaAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAg                                                                                                 >  2:691160/1‑139 (MQ=255)                   acgaAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGACGAGTCAGTCGCGTAATGCTTAGGCACAg                                                                                                 >  2:212679/1‑139 (MQ=255)                    cgaAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAgg                                                                                                >  2:349678/1‑139 (MQ=255)                          ggCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGgattgat                                                                                          <  1:400605/139‑1 (MQ=255)                               cACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTc                                                                                     >  1:523069/1‑139 (MQ=255)                                 ccTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGc                                                                                   >  1:757175/1‑139 (MQ=255)                                               aaaaTCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGa                                                                     <  2:429126/139‑1 (MQ=255)                                                      caTTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGATACGATTCCGCt                                                              <  1:349678/139‑1 (MQ=255)                                                          tttATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGAc                                                          <  2:523069/139‑1 (MQ=255)                                                           ttATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTgggg                                                                                                                              <  1:65373/70‑1 (MQ=255)                                                           ttATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTgggg                                                                                                                              >  2:65373/1‑70 (MQ=255)                                                                         cTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCt                                                                                                          >  1:602354/1‑76 (MQ=255)                                                                         cTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCt                                                                                                          <  2:602354/76‑1 (MQ=255)                                                                                       gcgcCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGACGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAgg                             <  1:212679/139‑1 (MQ=255)                                                                                         gcCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCa                           <  1:691160/139‑1 (MQ=255)                                                                                         gcCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCa                           <  2:757175/139‑1 (MQ=255)                                                                                         gcCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCa                           <  1:227870/139‑1 (MQ=255)                                                                                         gcCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCa                           <  1:224006/139‑1 (MQ=255)                                                                                                      ggTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCAc              <  1:602120/139‑1 (MQ=255)                                                                                                       gTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACt             >  2:600783/1‑139 (MQ=255)                                                                                                        tctcAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTa            >  1:787568/1‑139 (MQ=255)                                                                                                        tctcAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTa            >  1:799396/1‑139 (MQ=255)                                                                                                         ctcAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTaa           >  2:550421/1‑139 (MQ=255)                                                                                                            aaGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTAtt                 <  1:783282/130‑1 (MQ=255)                                                                                                            aaGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTAtt                 >  2:783282/1‑130 (MQ=255)                                                                                                              gCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAAt      <  2:201612/139‑1 (MQ=255)                                                                                                              gCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAAt      <  1:600783/139‑1 (MQ=255)                                                                                                                  aaaGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCt  <  2:390173/139‑1 (MQ=255)                                                                                                                              |                                                                                                                               GTCCTAAACTTGATCTCGACGAAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGAGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCT  >  NC_000913/1862509‑1862760

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