Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	344,072	G→A	intergenic (+139/‑104)	yahK → / → yahL	aldehyde reductase, NADPH‑dependent/uncharacterized protein YahL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	344,072	0	G	A	100.0%	72.0 / NA	21	intergenic (+139/‑104)	yahK/yahL	aldehyde reductase, NADPH‑dependent/uncharacterized protein YahL
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (12/9);  total (12/9)

AAATATCGTTTTGTTATCGATAATCGCACACTAACAGACTGAAAAAATTAATAAATACCCTGTGGTTTAACATATTAACTTCGCTCTCCACTTAACTTTTTAGTTAAGGAGAGCGAAATAATATCAAAGTAGCAGTAAAACCTATAACGTAAATTTAAATTGTTAAATTAACGCCCTCCAGTACACAATACTTCACACGTTAGTTATGAGCGATTTCTGATAGTGCCTGGTTTAATCAGAGCTTTATTATCTGCGACGTTTATTTTTATTTAAGAGAGATCGCGATGATATCATTAAAAGCTCCGCACAATAATTTGATGCCATATACGCAACAAAGCATACTTAATACGGTTAAAAACAATCAGTTACCAGAGGATATCAAAAGCTCCCTGGTTTCCTGTGTCGATATATTCAAGGTTTTGAT  >  NC_000913/343892‑344315                                                                                                                                                                                     |                                                                                                                                                                                                                                                    aaaTATCGTTTTGTTATCGATAATCGCACACTAACAGACTGAAAAAATTAATAAATACCCTGTGGTTTAACATATTAACTTCGCTCTCCACTTAACTTTTTAGTTAAGGAGAGCGAAATAATATCAAAGTAGCAGTAAAACCTATAACGTAAATTTAAATTGTTAAATTAACGCCCTCCAATACACAATACTTCACACGTTAGTTATGAGCGATTTCTGATAGTGCCTGGTTTAATCAGAGCTTTATTATc                                                                                                                                                                               >  2:33703/1‑251 (MQ=255)                             cacTAACAGACTGAAAAAATTAATAAATACCCTGTGGTTTAACATATTAACTTCGCTCTCCACTTAACTTTTTAGTTAAGGAGAGCGAAATAATATCAAAGTAGCAGTAAAACCTATAACGTAAATTTAAATTGTTAAATTAACGCCCTCCAATACACAATACTTCACACGTTAGTTATGAGCGATTTCTGATAGTGCCTGGTTTAATCAGAGCTTTATTATCTGCGACGTTTATTTTTATTTAagag                                                                                                                                                      >  1:172493/1‑248 (MQ=255)                              acTAACAGACTGAAAAAATTAATAAATACCCTGTGGTGTAACATATTAACTTCGCTCTCCACTTAACTTTTTAGTTAAGGAGAGCGAAATAATATCAAAGTAGCAGGAAAACCTATAACGTAAATTTAAATTGTTAAATTAACGCCCTCCAATACACAATACTTCACACTTTAGTTATGAGCGATTTCTGATAGTGCCTGGTTTAATCAGAGCTTTATTATCGGCGACGTTTATTTTTATTTAAGAGAGAt                                                                                                                                                  <  2:131306/251‑1 (MQ=255)                                      aCTGACAATATTAACAAATACAATATGGTTTAACATATTAACTTCACTATCCACTTAACTTTTTAGTTAAGAAGAGCGAACTAATATCAAAGACACAGTATCAACTATAACGTAAATTTAAATTGTTAAATTAACGCCATCCAATACACAATACTTCACACGATCGTTATGACCGATATCTGATAGTGCCTGGTTTACTCACAGCTTTATTATATGCGACGTTTCTTTTTATTTAATAGAGATCGCgatgt                                                                                                                                          >  2:126077/1‑250 (MQ=255)                                                 taataaATACCCTGTGGTTTAACATATTAACTTCGCTCTCCACTTAACTTTTTAGTTAAGGAGAGCGAAATAATATCAAAGTAGCAGTAAAACCTATAACGTAAATTTAAATTGTTAAATTAACGCCCTCCAATACACAATACTTCACACGTTAGTTATGAGCGATTTCTGATAGTGCCTGGTTTAATCAGAGCTTTATTATCTGCGACGTTTATTTTTATTTAAGAGAGATCGCGATGATATCATTaaa                                                                                                                                >  1:186789/1‑250 (MQ=255)                                                     aaataCCCTGTGGTTTAACATATTAACTTCGCTCTCCACTTAACTTTTTAGTTAAGGAGAGCGAAATAATATCAAAGTAGCAGTAAAACCTATAACGTAAATTTAAATTGTTAAATTAACGCCCTCCAATACACAATACTTCACACGTTAGTTATGAGCGATTTCTGATAGTGCCTGGTTTAATCAGAGCTTTATTATCTGCGACGTTTATTTTTATTTAAGAGAGATCGCGATGATATCATTAAAAGCTc                                                                                                                           >  1:81188/1‑251 (MQ=255)                                                         aCCCTGTGGTTTAACATATTAACTTCGCTCTCCACTTAACTTTTTAGTTAAGGAGAGCGAAATAATATCAAAGTAGCAGTAAAACCTATAACGTAAATTTAAATTGTTAAATTAACGCCCTCCAATACACAATACTTCACACGTTAGTTATGAGCGATTTCTGATAGTGCCTGGTTTAATCAGAGCTtta                                                                                                                                                                                    >  2:283161/1‑190 (MQ=255)                                                         aCCCTGTGGTTTAACATATTAACTTCGCTCTCCACTTAACTTTTTAGTTAAGGAGAGCGAAATAATATCAAAGTAGCAGTAAAACCTATAACGTAAATTTAAATTGTTAAATTAACGCCCTCCAATACACAATACTTCACACGTTAGTTATGAGCGATTTCTGATAGTGCCTGGTTTAATCAGAGCTtta                                                                                                                                                                                    <  1:283161/190‑1 (MQ=255)                                                                              aCTTCGCTCTCCACTTAACTTTTTAGTTAAGGAGAGCGAAATAATATCAAATTAGCAGTAAAACCTATAACGTAAATTGAAATTGTTAAAGTAACTCCCTCCAATACACAATTCTTCACACGTTAGTTATGATCGATTTCTGATAGTGCTTGGTTTAATCAGAGCTTTATTATTTGCGACGTCTATTTTTATTTAATAGAGATCGCGATGATATCATTAAAAGCTCCGCACAATAATCTGATGCGATATAc                                                                                                  <  2:120713/251‑1 (MQ=255)                                                                                  cGCTCTCCACTTAACTTTTTAGTTAAGGAAAGCGAAATAATATCAAAGTAGCAGTAAAACCTATATCGTAAATTTAAATTGTTAAATTAACGCCCTCCAATACACAATGCTTCACACGTTAGTTATGAGCGATTTCTGATAGTGCCGTGTTTAATCAGAGCTTTAGTATCTGCTACGTTTATTTTTATTTAAGAGAGATCGCGATGATATCATTAAAAGCTCCGCACAATAATTTGATGCTTTATACGcat                                                                                              <  2:108584/251‑2 (MQ=255)                                                                                                             gagaGCGAAATAATCTCAAAGTAGCAGTAAAACCTATAACGTACATTTAAATTGTTAAATTAACGCCCTCCAATACACAATACTTCACACGTTAGTTATGAGCGATTTCTGATAGTGCCTGGTTTAATCAGAGCTTTATTATCTGCGACGTTTATTTTTATTTAAGAGAGATCGCGCTGATATCATTAAAAGCTCCGCACAATAATTTGATGCCATATACGCAACAAAGCATACTTAATACGGTTAAAAAc                                                                   >  2:29429/1‑251 (MQ=255)                                                                                                             gagaGCGAAATAATATCAAAGTAGCAGTAAAACCTATAACGTAAATTTAAATTGTTAAATTAACGCCCTCCAATACACAATACTTCACACGTTAGTTATGAGCGATTTCTGATAGTGCCTGGTTTAATCAGAGCTTTATTATCTGCGACGTTTATTTTTATTTAAGAGAGATCGCGATGATATCATTAAAAGCTCCGCACAATAATTTGATGCCATATACGCAACAAAGCATACTTAATACGGTTAAAAAc                                                                   <  1:33703/251‑1 (MQ=255)                                                                                                                agCGAAATAATATCAAAGTAGCAGTAAAACCTATAACGTAAATTTAAATTGTTAAATTAACGCCCTCCAATACACAATACTTCACACGTTAGTTATGAGCGATTTCTGATAGTGCCTGGTTTAATCAGAGCTTTATTATCTGCGACGTTTATTTTTATTTAAGAGAGATCGCGATGATATCATTAAAAGCTcc                                                                                                                          >  1:263655/1‑193 (MQ=255)                                                                                                                agCGAAATAATATCAAAGTAGCAGTAAAACCTATAACGTAAATTTAAATTGTTAAATTAACGCCCTCCAATACACAATACTTCACACGTTAGTTATGAGCGATTTCTGATAGTGCCTGGTTTAATCAGAGCTTTATTATCTGCGACGTTTATTTTTATTTAAGAGAGATCGCGATGATATCATTAAAAGCTcc                                                                                                                          <  2:263655/193‑1 (MQ=255)                                                                                                                                   aGCAGTAAAACCTATAACGTAAATTTAAATTGTTAAATTAACGCCCTCCAATACACAATACTTCACACGTTAGTTATGAGCGATTTCTGATAGTGCCTGGTTTAATCAGAGCTTTATTATCTGCGACGTTTATTTTTATTTAAGAGAGATCGCGATGATATCATTAAAAGCTCCGCACAATAATTTGATGCCATATACGCAACAAAGCATACTTAATACGGTTAAAAACAATCAGTTACCAGAGGATATc                                              >  1:94511/1‑250 (MQ=255)                                                                                                                                   aGCAGTAAAACCTATAACGTAAATTTAAATTGTTAAATTAACGCCCTCCAATACACAATACTTCACACGTTAGTTATGAGCGATTTCTGATAGTGCCTGGTTTAATCAGAGCTTTATTATCTGCGACGGTTATTTTTATTTAAGAGAGATCGCGATGATATCATTAAAAGCTCCGCACAATAATTTGATGCCATATACGCAACAAAGCATACTTAATACGGTTAAAAACAATCAGTTACCAGAGGATATCa                                             <  1:264926/251‑1 (MQ=255)                                                                                                                                       gTAAAACCTATAACGTAAATTTAAATTGTTAAATTAACGCCCTCCAATACACAATACTTCACACGTTAGTTATGAGCGATTTCTGATAGTGCCTGGTTTAATCAGAGCTTTATTATCTGCGACGTTTATTTTTATTTAAGAGAGATCGCGATGATATCATTAAAAGCTCCGCACAATAATTTGATGCCATATACGCAACAAAGCATACTTAATACGGTTAAAAACAATCAGTTACCAGAGGATATCaaa                                           >  1:206690/1‑249 (MQ=255)                                                                                                                                                ataACGTAAATTTAAATTGTTAAATTAACGCCCTCCAATACACAATACTTCACACGTTAGTTATGAGCGATTTCTGATAGTGCCTGGTTTAATCAGAGCTTTATTATCTGCGACGTTTATTTTTATTTAAGAGAGATCGCGATGATATCATTAAAAGCTCCGCACAATAATTTGATGCCATATACGCAACAAAGCATACTTAATACGGTTAAAAACAATCAGTTACCAGAGGATATCAAAAGCTCCCTGGt                                >  1:118878/1‑251 (MQ=255)                                                                                                                                                    cGTAAATTTAAATTGTTAAATTAACGCCCTCCAATACACAATACTTCACACGTTAGTTATGAGCGATTTCTGATAGTGCCTGGTTTAATCAGAGCTTTATTATCTGCGACGtttattttt                                                                                                                                                               >  1:126747/1‑120 (MQ=255)                                                                                                                                                    cGTAAATTTAAATTGTTAAATTAACGCCCTCCAATACACAATACTTCACACGTTAGTTATGAGCGATTTCTGATAGTGCCTGGTTTAATCAGAGCTTTATTATCTGCGACGtttattttt                                                                                                                                                               <  2:126747/120‑1 (MQ=255)                                                                                                                                                                              cccTCCAATACACAATACTTCACACGTTAGTTATGAGCGATTTCTGATAGTGCCTGGTTTAATCAGAGCTTTATTATCTGCGACGTTTATTTTTATTTAAGAGAGATCGCGATGATATCATTAAAAGCTCCGCACAATAATTTGATGCCATATACGCAACAAAGCATACTTAATACGGTTAAAAACAATCAGTTACCAGAGGATATCAAAAGCTCCCTGGTTTCCTGTGTCGATATATTCAAGGTTTTGAt  <  1:88174/251‑1 (MQ=255)                                                                                                                                                                                     |                                                                                                                                                                                                                                                    AAATATCGTTTTGTTATCGATAATCGCACACTAACAGACTGAAAAAATTAATAAATACCCTGTGGTTTAACATATTAACTTCGCTCTCCACTTAACTTTTTAGTTAAGGAGAGCGAAATAATATCAAAGTAGCAGTAAAACCTATAACGTAAATTTAAATTGTTAAATTAACGCCCTCCAGTACACAATACTTCACACGTTAGTTATGAGCGATTTCTGATAGTGCCTGGTTTAATCAGAGCTTTATTATCTGCGACGTTTATTTTTATTTAAGAGAGATCGCGATGATATCATTAAAAGCTCCGCACAATAATTTGATGCCATATACGCAACAAAGCATACTTAATACGGTTAAAAACAATCAGTTACCAGAGGATATCAAAAGCTCCCTGGTTTCCTGTGTCGATATATTCAAGGTTTTGAT  >  NC_000913/343892‑344315

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