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%	177.2 / NA	50	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 (23/27);  total (23/27)

GGCAAAATTTTGTCCTAAACTTGATCTCGACGAAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGAGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCTGGTGGAAT  >  NC_000913/1862498‑1862768                                                                                                                                         |                                                                                                                                       ggCAAAATTTTGTCCTAAACTTGATCTCGACGAAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTgggg                                                                                                                                      <  2:455435/139‑1 (MQ=255)          ttGTCCTAAACTTGATCTCGACGAAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGagtcagt                                                                                                                             <  2:889663/139‑1 (MQ=255)               cTAAACTTGATCTCGACGAAATGGCTGCACCTAAAACGTGATTAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGACGAGTCAGTCGCGt                                                                                                                        <  2:1468504/139‑1 (MQ=255)                        aTCTCGACGAAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCtt                                                                                                                 >  2:571446/1‑137 (MQ=255)                        aTCTCGACGAAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCtt                                                                                                                 <  1:571446/137‑1 (MQ=255)                            cgacgaAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGcac                                                                                                           >  2:823635/1‑139 (MQ=255)                                            ccTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGc                                                                                           >  1:769869/1‑139 (MQ=255)                                             cTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCa                                                                                          <  2:1344732/139‑1 (MQ=255)                                              tAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCaa                                                                                         >  2:1306673/1‑139 (MQ=255)                                                     tgatgaAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGgattgat                                                                                                  <  2:710646/123‑1 (MQ=255)                                                     tgatgaAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGgattgat                                                                                                  >  1:710646/1‑123 (MQ=255)                                                          aaaaTCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGa                                                                             <  1:163126/139‑1 (MQ=255)                                                             aTCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGAt                                                                                     <  2:741111/128‑1 (MQ=255)                                                             aTCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGAt                                                                                     >  1:741111/1‑128 (MQ=255)                                                                      ttATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACg                                                                 <  1:823635/139‑1 (MQ=255)                                                                         tCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTcgcg                                                                                                                         >  2:89125/1‑80 (MQ=255)                                                                         tCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTcgcg                                                                                                                         <  1:89125/80‑1 (MQ=255)                                                                            tAATTGCCCTTTAAAATTCGGGTCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACg                                                                 >  2:1418102/1‑133 (MQ=255)                                                                            tAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACg                                                                 <  1:1418102/133‑1 (MQ=255)                                                                              aTTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGTCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACg                                                                 <  1:999219/131‑1 (MQ=255)                                                                              aTTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACg                                                                 >  2:999219/1‑131 (MQ=255)                                                                                      ttAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGca                                                                                                            >  1:1287186/1‑80 (MQ=255)                                                                                      ttAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGca                                                                                                            <  2:1287186/80‑1 (MQ=255)                                                                                               ggggCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTAc                                        >  2:403403/1‑139 (MQ=255)                                                                                               ggggCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTAc                                        >  1:15530/1‑139 (MQ=255)                                                                                                 ggCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGt                                                           <  1:1275697/118‑1 (MQ=255)                                                                                                 ggCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGt                                                           >  2:1275697/1‑118 (MQ=255)                                                                                                 ggCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTaa                                                         <  2:909921/120‑1 (MQ=255)                                                                                                 ggCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCAGCTTGACGCTGCGTaa                                                         >  1:909921/1‑120 (MQ=255)                                                                                                             aTGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGca                                                                                                            <  1:959186/57‑1 (MQ=255)                                                                                                             aTGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGca                                                                                                            >  2:959186/1‑57 (MQ=255)                                                                                                              tgtgGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTAtt                         <  2:15530/139‑1 (MQ=255)                                                                                                                 ggTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCAc                      <  1:1306673/139‑1 (MQ=255)                                                                                                                  gTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACt                     <  2:769869/139‑1 (MQ=255)                                                                                                                   tctcAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGca                                                                                                            <  1:262233/51‑1 (MQ=255)                                                                                                                   tctcAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGca                                                                                                            >  2:262233/1‑51 (MQ=255)                                                                                                                   tctcAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGacac                                                                               >  2:934151/1‑80 (MQ=255)                                                                                                                   tctcAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGacac                                                                               <  1:934151/80‑1 (MQ=255)                                                                                                                    ctcAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCt                                                                                                                  <  1:286613/44‑1 (MQ=255)                                                                                                                    ctcAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCt                                                                                                                  >  2:286613/1‑44 (MQ=255)                                                                                                                      cAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGa                                                                             <  2:1076284/79‑1 (MQ=255)                                                                                                                      cAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGa                                                                             >  1:1076284/1‑79 (MQ=255)                                                                                                                        aGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACg                                                                 >  2:1246276/1‑89 (MQ=255)                                                                                                                        aGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACg                                                                 <  1:1246276/89‑1 (MQ=255)                                                                                                                            taagagacaggTGGGGCGCGTCAGTCGCGTAATGCTGAAGCACAGGTTTGATTTGTCGCAATGATTGACACGATTCAGCTTGACGCTGCGTAAGGTTTTTGCATTTTTACAGGCAACCTTTTATTCACGAACAAATAGc           <  2:1007259/129‑1 (MQ=255)                                                                                                                                 gAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCtggtg      <  1:403403/139‑1 (MQ=255)                                                                                                                                   agagTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGa                                                                             >  2:430403/1‑66 (MQ=255)                                                                                                                                   agagTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGa                                                                             <  1:430403/66‑1 (MQ=255)                                                                                                                                     agTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTg                                                                                               >  2:748064/1‑46 (MQ=255)                                                                                                                                     agTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTg                                                                                               <  1:748064/46‑1 (MQ=255)                                                                                                                                     agTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCTGGTGGAat  >  1:1474197/1‑139 (MQ=255)                                                                                                                                         |                                                                                                                                       GGCAAAATTTTGTCCTAAACTTGATCTCGACGAAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGAGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCTGGTGGAAT  >  NC_000913/1862498‑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: ‑