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%	138.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 (16/23);  total (16/23)

TTTTGTCCTAAACTTGATCTCGACGAAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGAGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCTGGTG  >  NC_000913/1862505‑1862764                                                                                                                                  |                                                                                                                                   ttttGTCCTAAACTTGATCTCGACGAAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGagtca                                                                                                                           <  2:1356342/139‑1 (MQ=255)        cTAAACTTGATCTCGACGAAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGt                                                                                                                    <  1:480734/139‑1 (MQ=255)                     cgacgaAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTcgcg                                                                                                                     >  2:1451607/1‑125 (MQ=255)                     cgacgaAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTcgcg                                                                                                                     <  1:1451607/125‑1 (MQ=255)                       acgaAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGcac                                                                                                       <  2:1506155/137‑1 (MQ=255)                       acgaAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGcac                                                                                                       >  1:1506155/1‑137 (MQ=255)                                         aaTCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATg                                                                                   >  2:126041/1‑139 (MQ=255)                                           tCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGAt                                                                                 <  2:389740/139‑1 (MQ=255)                                              tgatgaAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCtt                                                                                                             >  2:322379/1‑108 (MQ=255)                                              tgatgaAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCtt                                                                                                             <  1:322379/108‑1 (MQ=255)                                                atgaAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTaa                                                                                                                  >  1:1427558/1‑101 (MQ=255)                                                atgaAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTaa                                                                                                                  <  2:1427558/101‑1 (MQ=255)                                                atgaAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGaca                                                                            <  1:455802/139‑1 (MQ=255)                                                   aaaaTCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGa                                                                         <  1:126041/139‑1 (MQ=255)                                                           aTTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCtt                                                                 >  2:1427946/1‑139 (MQ=255)                                                                     tAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGt                                                       >  1:324897/1‑139 (MQ=255)                                                                       aTTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTaa                                                     >  1:804153/1‑139 (MQ=255)                                                                               ttAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTg                                             <  2:804153/139‑1 (MQ=255)                                                                                tAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACg                                                             <  1:2684/122‑1 (MQ=255)                                                                                tAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACg                                                             >  2:2684/1‑122 (MQ=255)                                                                                tAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGt                                            >  1:1339985/1‑139 (MQ=255)                                                                                        ggggCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTAc                                    >  1:296980/1‑139 (MQ=255)                                                                                         gggCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACa                                   <  2:296980/139‑1 (MQ=255)                                                                                             gcCGACCCCATGTGGTCTCAAGCCCAAAGGAAGGGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCa                               <  1:679758/139‑1 (MQ=255)                                                                                                       tgtgGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTAtt                     <  1:30219/139‑1 (MQ=255)                                                                                                           gTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACt                 <  2:324897/139‑1 (MQ=255)                                                                                                             ctcAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCa                               <  1:867336/123‑1 (MQ=255)                                                                                                             ctcAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCa                               >  2:867336/1‑123 (MQ=255)                                                                                                               cAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACa             <  2:204094/139‑1 (MQ=255)                                                                                                                aaGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGt                                                                                                                    >  2:1427612/1‑35 (MQ=255)                                                                                                                aaGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGt                                                                                                                    <  1:1427612/35‑1 (MQ=255)                                                                                                                       aaGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCtg     >  1:48409/1‑139 (MQ=255)                                                                                                                          gAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCtggtg  <  2:427094/139‑1 (MQ=255)                                                                                                                          gAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCtggtg  <  1:1427946/139‑1 (MQ=255)                                                                                                                          gAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCtggtg  <  2:95935/139‑1 (MQ=255)                                                                                                                              agTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGgattgat                                                                                              >  2:277244/1‑43 (MQ=255)                                                                                                                              agTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGgattgat                                                                                              <  1:277244/43‑1 (MQ=255)                                                                                                                              agTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAgg                                 <  2:101719/104‑1 (MQ=255)                                                                                                                              agTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAgg                                 >  1:101719/1‑104 (MQ=255)                                                                                                                                  |                                                                                                                                   TTTTGTCCTAAACTTGATCTCGACGAAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGAGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCTGGTG  >  NC_000913/1862505‑1862764

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