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%	137.8 / NA	39	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 (17/22);  total (17/22)

AATTTTGTCCTAAACTTGATCTCGACGAAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGAGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCTGGTGGAAT  >  NC_000913/1862503‑1862768                                                                                                                                    |                                                                                                                                       aaTTTTGTCCTAAACTTGATCTCGACGAAATGGCTGCACCTAAATCGTGATGAAAACCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGagt                                                                                                                                 <  2:1287985/139‑1 (MQ=255)         ccTAAACTTGATCTCGACGAAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTcgcg                                                                                                                         <  2:1271027/139‑1 (MQ=255)                      tcGACGAAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGca                                                                                                            <  2:1174560/139‑1 (MQ=255)                      tcGACGAAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGca                                                                                                            <  1:874595/139‑1 (MQ=255)                       cgacgaAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGcac                                                                                                           >  2:81639/1‑139 (MQ=255)                                    gCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGt                                                                                              >  2:31692/1‑139 (MQ=255)                                            aTCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGa                                                                                      <  2:158761/139‑1 (MQ=255)                                                  atgaAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGaca                                                                                >  2:1068194/1‑139 (MQ=255)                                                               ttttATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGa                                                                   <  1:1068194/139‑1 (MQ=255)                                                                               cTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTcg                                                                                                                           <  2:32368/67‑1 (MQ=255)                                                                               cTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTcg                                                                                                                           >  1:32368/1‑67 (MQ=255)                                                                                   aaaaTTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTa                                               <  1:1252031/139‑1 (MQ=255)                                                                                             gcgcCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAgg                                     <  2:685981/139‑1 (MQ=255)                                                                                             gcgcCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAgg                                     <  1:81639/139‑1 (MQ=255)                                                                                               gcCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCa                                   <  1:1050133/139‑1 (MQ=255)                                                                                                 cGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAAc                                 <  1:632614/139‑1 (MQ=255)                                                                                                               ctcAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGca                                                                                                            >  2:668478/1‑50 (MQ=255)                                                                                                               ctcAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGca                                                                                                            <  1:668478/50‑1 (MQ=255)                                                                                                               ctcAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTaa                   >  2:751463/1‑139 (MQ=255)                                                                                                                 cAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACa                 >  2:747059/1‑139 (MQ=255)                                                                                                                  aaGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGca                                                                                                            <  2:693313/47‑1 (MQ=255)                                                                                                                  aaGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGca                                                                                                            >  1:693313/1‑47 (MQ=255)                                                                                                                  aaGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTAtt                         >  1:177769/1‑130 (MQ=255)                                                                                                                  aaGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTAtt                         <  2:177769/130‑1 (MQ=255)                                                                                                                   aGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTAtt                         >  2:1269380/1‑129 (MQ=255)                                                                                                                   aGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTAtt                         <  1:1269380/129‑1 (MQ=255)                                                                                                                    gCCGAAAGGAAGGGTGTGGATAGTCAGTCGGGTAATTCTTAGGCACAGGATTGATTTGTCGCAATGATTGGCACGATTCCGCTTGACGCTGGGTAAGGTTTTTGGAATTTTACAGGCAACCTTTTATTCACTAACAAAt              <  1:31692/139‑1 (MQ=255)                                                                                                                       cAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGc                                                                                                             >  2:1338341/1‑41 (MQ=255)                                                                                                                       cAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGc                                                                                                             <  1:1338341/41‑1 (MQ=255)                                                                                                                       cAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACg                                                                 <  2:149709/85‑1 (MQ=255)                                                                                                                       cAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACg                                                                 >  1:149709/1‑85 (MQ=255)                                                                                                                        aaaGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGAc                                                                  <  2:1095433/83‑1 (MQ=255)                                                                                                                        aaaGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGAc                                                                  >  1:1095433/1‑83 (MQ=255)                                                                                                                        aaaGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCt          >  2:476061/1‑139 (MQ=255)                                                                                                                         aaGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAAc                                 >  2:1370734/1‑115 (MQ=255)                                                                                                                         aaGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAAc                                 <  1:1370734/115‑1 (MQ=255)                                                                                                                            gAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCtggtg      <  1:747059/139‑1 (MQ=255)                                                                                                                                agTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCa                                                                                          >  2:1322435/1‑51 (MQ=255)                                                                                                                                agTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCa                                                                                          <  1:1322435/51‑1 (MQ=255)                                                                                                                                agTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCTGGTGGAat  >  1:111391/1‑139 (MQ=255)                                                                                                                                    |                                                                                                                                       AATTTTGTCCTAAACTTGATCTCGACGAAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGAGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCTGGTGGAAT  >  NC_000913/1862503‑1862768

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