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%	148.9 / NA	41	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 (22/19);  total (22/19)

AACTTGATCTCGACGAAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGAGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCTGGTGGAATA  >  NC_000913/1862515‑1862769                                                                                                                        |                                                                                                                                        aaCTTGATCTCGACGAAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAAt                                                                                                                      >  1:221157/1‑139 (MQ=255)        tctcGACGAAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAgg                                                                                                               <  1:120850/139‑1 (MQ=255)                  aTGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGgattg                                                                                                     >  1:265184/1‑139 (MQ=255)                   tGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGgattga                                                                                                    >  2:822499/1‑139 (MQ=255)                          aCCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCt                                                                                                                   <  2:978701/117‑1 (MQ=255)                          aCCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCt                                                                                                                   >  1:978701/1‑117 (MQ=255)                                   gTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCt                                                                                                                   <  1:508526/108‑1 (MQ=255)                                   gTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCt                                                                                                                   >  2:508526/1‑108 (MQ=255)                                    tgatgaAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGa                                                                                   <  1:557516/139‑1 (MQ=255)                                      atgaAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGt                                                                                                                         <  2:940725/99‑1 (MQ=255)                                      atgaAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGt                                                                                                                         >  1:940725/1‑99 (MQ=255)                                           aaTCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGAtt                                                                            <  2:265184/139‑1 (MQ=255)                                                     ttATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACg                                                                  <  2:221157/139‑1 (MQ=255)                                                             aTTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGca                                                                                                             >  1:69638/1‑88 (MQ=255)                                                             aTTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGca                                                                                                             <  2:69638/88‑1 (MQ=255)                                                                  ccTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGttt                                                     >  1:744311/1‑139 (MQ=255)                                                                   cTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGtttt                                                    >  1:326832/1‑139 (MQ=255)                                                                      tAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGAtt                                                                            <  2:975622/112‑1 (MQ=255)                                                                      tAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGAtt                                                                            >  1:975622/1‑112 (MQ=255)                                                                           ttCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGc                                                                 <  2:528587/118‑1 (MQ=255)                                                                           ttCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGc                                                                 >  1:528587/1‑118 (MQ=255)                                                                             cGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTa                                          >  2:291539/1‑139 (MQ=255)                                                                                   gcCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCa                                    <  2:326832/139‑1 (MQ=255)                                                                                                 gTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGt                                                                                                                         <  1:329239/40‑1 (MQ=255)                                                                                                 gTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGt                                                                                                                         >  2:329239/1‑40 (MQ=255)                                                                                                 gTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACg                                                                  <  2:906202/95‑1 (MQ=255)                                                                                                 gTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACg                                                                  >  1:906202/1‑95 (MQ=255)                                                                                                  tctcAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAgg                                      >  1:243988/1‑122 (MQ=255)                                                                                                  tctcAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAgg                                      <  2:243988/122‑1 (MQ=255)                                                                                                  tctcAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTa                     >  1:938115/1‑139 (MQ=255)                                                                                                   ctcAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTaa                                                          <  2:510240/101‑1 (MQ=255)                                                                                                   ctcAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTaa                                                          >  1:510240/1‑101 (MQ=255)                                                                                                       aGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAgg                                                        >  2:728253/1‑99 (MQ=255)                                                                                                       aGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAgg                                                        <  1:728253/99‑1 (MQ=255)                                                                                                            aaaGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCt           <  2:295772/139‑1 (MQ=255)                                                                                                              aGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAgg                                      <  2:895290/110‑1 (MQ=255)                                                                                                              aGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAgg                                      >  1:895290/1‑110 (MQ=255)                                                                                                              aGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGT‑AGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCtggt        >  2:304982/1‑139 (MQ=255)                                                                                                                gAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCtggtg       <  2:744311/139‑1 (MQ=255)                                                                                                                   gagTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCTGGTGGaa    >  1:420117/1‑139 (MQ=255)                                                                                                                    agTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACC‑TTTATTCACTAACAAATAGCTGGTGGAata  >  1:771071/1‑139 (MQ=255)                                                                                                                        |                                                                                                                                        AACTTGATCTCGACGAAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGAGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCTGGTGGAATA  >  NC_000913/1862515‑1862769

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