Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,244,739	A→G	S113P (TCC→CCC)	malF ←	maltose transporter subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,244,739	0	A	G	100.0%	136.6 / NA	38	S113P (TCC→CCC)	malF	maltose transporter subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (17/21);  total (17/21)

CAGTTTTTGCTCGCCGCCAAATTTAAAAGCGTCGGAGAGGTAATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCTGCTTGCCAGGAGCGATCTAACAACACTTCCTGCGCACGTTCAAAAGTCAGCTGGTTAGTGCTGCTGTAGTTGGTGAAGGCAATGGCGATGGTGCAGACCAGAGGGAAGAGGACGAATAATCCCATTCCAGCCAT  >  NC_000913/4244605‑4244862                                                                                                                                       |                                                                                                                            cAGTTTTTGCTCGCCGCCAAATTTAAAAGCGTCGGAGAGGTAATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCTGCTTGCCAGGGGCGa                                                                                                                         <  2:785569/139‑1 (MQ=255) cAGTTTTTGCTCGCCGCCAAATTTAAAAGCGTCGGAGAGGTAATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCTGCTTGCCAGGGGCGa                                                                                                                         <  1:1218374/139‑1 (MQ=255) cAGTTTTTGCTCGCCGCCAAATTTAAAAGCGTCGGAGAGGTAATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCTGCTTGCCAGGGGCGa                                                                                                                         >  1:635978/1‑139 (MQ=255)         gCTCGCCGCCAAATTTAAAAGCGTCGGAGAGGTAATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCTGCTTGCCAGGGGCGATCTaacaa                                                                                                                 <  2:591482/139‑1 (MQ=255)                                tCGGAGAGGTAATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCTGCTTGCCAGGGGCGATCTaacaa                                                                                                                 <  2:21509/116‑1 (MQ=255)                                tCGGAGAGGTAATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCTGCTTGCCAGGGGCGATCTaacaa                                                                                                                 >  1:21509/1‑116 (MQ=255)                                tCGGAGAGGTAATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCTGCTTGCCAGGGGCGATCTAACAACACTTCCTGCGCACGTTCAAAAg                                                                                          >  2:1120005/1‑139 (MQ=255)                                    agagGTAATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCTGCTTGCCAGGGGCGATCTAACAACACTTCCTGCGCACGTTCAAAAGTCAg                                                                                      <  1:1163562/139‑1 (MQ=255)                                    agagGTAATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCTGCTTGCCAGGGGCGATCTAACAACACTTCCTGCGCACGTTCAAAAGTCAg                                                                                      <  2:1172563/139‑1 (MQ=255)                                            tttttGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCTGCTTGCCAGGGGCGATCTAACAACACTTCCTGCGCACGTTCAAAAGTCAGCTGGTTAg                                                                              <  1:985260/139‑1 (MQ=255)                                             ttttGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCTGCTTGCCAGGGGCGATCTAACAACACTTCCTGCGCACGTTCAAAAGTCAGCTGGTTAGt                                                                             <  2:635978/139‑1 (MQ=255)                                                  ccGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCTGCTTGCCAGGGGCGATCTAACAACACTTCCTGCGCACGTTCAAAAGTCAGCTGGTTAGtgctgc                                                                        >  1:428822/1‑139 (MQ=255)                                                   cGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCTGCTTGCCAGGGGCGATCTAACAACACTTCCTGCGCACGTTCAAAAGTCAGCTGGTTAGTgctgct                                                                       >  2:814961/1‑139 (MQ=255)                                                       ttCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCTGCTTGCCAGGGGCGATCTAACAACACTTCCTGCGCACGTTTAAAAGTCAGCTGGTTAGTGCTGCTGTAg                                                                   <  1:814961/139‑1 (MQ=255)                                                                gCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCTGCTTGCCAGGGGCGATCTAACAACACTTCCTGCGCACGTTCAAAAGTCAGCTGGTTAGTgctgct                                                                       >  1:1144519/1‑126 (MQ=255)                                                                gCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCTGCTTGCCAGGGGCGATCTAACAACACTTCCTGCGCACGTTCAAAAGTCAGCTGGTTAGTgctgct                                                                       <  2:1144519/126‑1 (MQ=255)                                                                    aGCGCCAGTAGCCACTCATCGCCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCTGCTTGCCAGGGGCGATCTAACAACACTTCCTGCGCACGTTCAAAAGTCAGCTGGTTAGTGCTGCTGTAGTTGGTGAAGGCaa                                                      >  2:629596/1‑139 (MQ=255)                                                                          aGTTGCCACTCATCGCCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCTGCTTGCCAGGGGCGATCTAACAACACTTCCTGCGCACGTTCAAAAGTCAGCTGGTTAGTGCTGCTGTAGTTGGTGAAGGCAATGGCGa                                                <  1:468170/139‑1 (MQ=255)                                                                              gCCACTCATCGCCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCC‑‑‑‑TGCCAGGGGCGATCTAACAACACTTCCTGCGCACGTTCAAAAGTCAg                                                                                      >  2:1210765/1‑93 (MQ=255)                                                                              gCCACTCATCGCCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCC‑‑‑‑TGCCAGGGGCGATCTAACAACACTTCCTGCGCACGTTCAAAAGTCAg                                                                                      <  1:1210765/93‑1 (MQ=255)                                                                               ccACTCATCGCCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCTGCTTGCCAGGGGCGATCTAACAACACTTCCTGCGCACGTTCAAAAGTCAGCTGGTTAGTGCTGCTGTAGTTGGTGAAGGCAATGGCGATGGTg                                           <  1:629596/139‑1 (MQ=255)                                                                                      tCGCCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCTGCTTGCCAGGGGCGATCTAACAACACTTCCTGCGCACGTTCAAAAGTCAGCTGGTTAGTGCTGCTGTAGTTGGTGAAGGCAATGGCGATGGTGCAGACCa                                    <  2:428822/139‑1 (MQ=255)                                                                                       cgcccgccGGGTAAAGACCAAAGTTATAGGTTTTGCCTGCTTGCCAGGGGCGATCTAACAACACTTCCTGCGCACGTTCAAAAGTCAGCTGGTTAGt                                                                             >  2:1152714/1‑97 (MQ=255)                                                                                       cgcccgccGGGTAAAGACCAAAGTTATAGGTTTTGCCTGCTTGCCAGGGGCGATCTAACAACACTTCCTGCGCACGTTCAAAAGTCAGCTGGTTAGt                                                                             <  1:1152714/97‑1 (MQ=255)                                                                                       cgcccgccGGGTAAAGACCAAAGTTATAGGTTTTGCCTGCTTGCCAGGGGCGATCTAACAACACTTCCTGCGCACGTTCAAAAGTCAGCTGGTTAGTGCTGCTGTAGTTGGTGAAGGCAATGGCGATGGTGCAGACCag                                   >  1:101496/1‑139 (MQ=255)                                                                                                 gTAAAGACCAAAGTTATAGGTTTTGCCTGCTTGCCAGGGGCGATCTAACAACACTTCCTGCGCACGTTCAAAAGTCAGCTGGTTAg                                                                              <  2:131888/86‑1 (MQ=255)                                                                                                 gTAAAGACCAAAGTTATAGGTTTTGCCTGCTTGCCAGGGGCGATCTAACAACACTTCCTGCGCACGTTCAAAAGTCAGCTGGTTAg                                                                              >  1:131888/1‑86 (MQ=255)                                                                                                                      tttGCCTGCTTGCCAGGGGCGATCTAACAACACTTCCTGCGCACGTTCAAAAGTCAGCTGGTTAg                                                                              <  1:673238/65‑1 (MQ=255)                                                                                                                      tttGCCTGCTTGCCAGGGGCGATCTAACAACACTTCCTGCGCACGTTCAAAAGTCAGCTGGTTAg                                                                              <  2:714735/65‑1 (MQ=255)                                                                                                                      tttGCCTGCTTGCCAGGGGCGATCTAACAACACTTCCTGCGCACGTTCAAAAGTCAGCTGGTTAg                                                                              >  2:673238/1‑65 (MQ=255)                                                                                                                      tttGCCTGCTTGCCAGGGGCGATCTAACAACACTTCCTGCGCACGTTCAAAAGTCAGCTGGTTAg                                                                              >  1:714735/1‑65 (MQ=255)                                                                                                                        tGCCTGCTTGCCAGGGGCGATCTAACAATACTTCCTGCGCACGTTCAAAAGTCAGATGGTTAGTGCTGCTGTAGTTGGTGAAGGCAATGGCGATGGTGCa                                         >  2:1190133/1‑100 (MQ=255)                                                                                                                        tGCCTGCTTGCCAGGGGCGATCTAACAATACTTCCTGCGCACGTTCAAAAGTCAGATGGTTAGTGCTGCTGTAGTTGGTGAAGGCAATGGCGATGGTGCa                                         <  1:1190133/100‑1 (MQ=255)                                                                                                                        tGCCTGCTTGCCAGGGGCGATCTAACAACACTTCCTGCGCACGTTCAAAAGTCAGCTGGTTAGt                                                                             <  2:324842/64‑1 (MQ=255)                                                                                                                        tGCCTGCTTGCCAGGGGCGATCTAACAACACTTCCTGCGCACGTTCAAAAGTCAGCTGGTTAGt                                                                             >  1:324842/1‑64 (MQ=255)                                                                                                                        tGCCTGCTTGCCAGGGGCGATCTAACAACACTTCCTGCGCACGTTCAAAAGTCAGCTGGTTAGTGCTGCTGTAg                                                                   <  2:1286349/74‑1 (MQ=255)                                                                                                                        tGCCTGCTTGCCAGGGGCGATCTAACAACACTTCCTGCGCACGTTCAAAAGTCAGCTGGTTAGTGCTGCTGTAg                                                                   >  1:1286349/1‑74 (MQ=255)                                                                                                                        tGCCTGCTTGCCAGGGGCGATCTAACAACACTTCCTGCGCACGTTCAAAAGTCAGCTGGTTAGTGCTGCTGTAGTTGGTGAAGGCAATGGCGATGGTGCAGACCAGAGGGAAGAGGACGAATAATCCCATTCCAGCCAt  >  1:770994/1‑139 (MQ=255)                                                                                                                                       |                                                                                                                            CAGTTTTTGCTCGCCGCCAAATTTAAAAGCGTCGGAGAGGTAATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCTGCTTGCCAGGAGCGATCTAACAACACTTCCTGCGCACGTTCAAAAGTCAGCTGGTTAGTGCTGCTGTAGTTGGTGAAGGCAATGGCGATGGTGCAGACCAGAGGGAAGAGGACGAATAATCCCATTCCAGCCAT  >  NC_000913/4244605‑4244862

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