Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,278,488	C→T	P4S (CCA→TCA)	ghxP →	guanine/hypoxanthine transporter GhxP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,278,488	0	C	T	100.0%	91.9 / NA	33	P4S (CCA→TCA)	ghxP	guanine/hypoxanthine transporter GhxP
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (12/21);  total (12/21)

CAAGACAATAAAAGAGAAAAAAGCAGCAAACTTCGGTTGAAAAAGCCGCTATGATCGCCGGATAATCGTTTGCTTTTTTTACCACCCGTTTTGTATGCGCGGAGCTAAACGTTTGCTTTTTTGCGACGCAGCAAATTGTCGCAAACCTGGAGCAGGAAGATAACGTTTCGCTGGCAGGGGATTGTCCGCCACGCATCTTGACGAAAATTAAACTCTCAGGGGATGTTTTCTTATGTCTACGCCATCAGCGCGTACCGGCGGTTCACTCGACGCCTGG‑TTTAAAATTTCACAACGTG‑GAAGCACTGTCCGTCAGG‑AAGTGG‑TTGCCGGGTTAACAACGTTTCTGGCGATGGTCTACTCGGTCATCGTCGTTCCAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGAT  >  NC_000913/4278247‑4278732                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                         cAAGACAATAAAAGAGAAAAAAGCAGCAAACTTCGGTTGAAAAAGCCGCTATGATCGCCGGATAATCGTTTGCTTTTTTTACCACCCGTTTTGTATGCGCGGAGCTAAACGTTTGCTTTTTTGCGACGCAGCAAATTGTCGCAAACCTGGAGCAGGAAGATAACGTTTCGCTGGCAGGGGATTGTCCGCCACGCAGCTTGACGAAAATTAAACTCTCCGGGGATGTTTTCTTATGTCTACGTCATCAgcgc                                                                                                                                                                                                                                                 >  2:192150/1‑251 (MQ=255)   aGACAATAAAAGAGAAAAAAGCAGCAAACTTCGGTTGACAAAGCCGCTATGATCGCCGGATAATCGTTTGCTTTTTTTACCACCCGTTTTGTATGCGCGGAGCTAAACGTTTGCTTTTTTGCGACGCAGCAAATTGTCGCAAACCTGGAGCAGGAAGATAACGTTTCGCTGGCAGGGGATTGTCCGCCACGCCTCTTGACGAAAATTAAACACTCAGGGGATGTTTTCTTCTGTATACGTCATCAAcgcgg                                                                                                                                                                                                                                               >  2:167784/1‑250 (MQ=255)                     aaGCAGCAAACTTCGGTTGAAAAAGCCGCTATGATCGCCGGATAATCGTTTGCTTTTTTTACCTCCCGTTTTGTAGGCGCGGAGCTAAACGTTTGCTTTTTTGCGACGCAGCAAATTGTCGCAAACCTGGAGCAGGAAGATAACGTTTCGCTGGCAGGGGATTGTCCGCCACGCATCTTGACGAAAATTAAACTCTCAGGGGATGTTTTCTTATGTCTACGTCATCAGCGCGTACCGGCGGTTCACTCGAc                                                                                                                                                                                                                             <  2:411580/251‑1 (MQ=255)                       gcagcaAACTTCGGTTGAAAAAGCCGCTATGATCGCCGGATAATCGTTTGCTTTTTTTACCACCCGTTTTGTATGCGCGGAGCTAAACGTTTGCTTTTTTGCGACGCAGCAAATTGTCGCAAACCTGGAGCAGGAAGATAACGTTTCGCTGGCAGGGGATTGTCCGCCACGCATCTTGACGAAAATTAAACTCTCAGGGGATGTTTTCTTATGTCTACGTCATCAGCGCGTACCGGCGGTTCAccccactc                                                                                                                                                                                                                           >  1:162442/1‑244 (MQ=255)                          gcaAACTTCGGTTGAAAAAGCCGCTATGATCGCCGGATAATCGTTTGCTTTTTTTACCACCCGTTTTGTATGCGCGGAGCTAAACGTTTGCTTTTTTGCGACGCAGCAAATTGTCGCAAACCTGGAGCAGGAAGATAACGTTTCGCTGGCAGGGGATTGTCCGCCACGCATCTTGACGAAAATTAAACTCTCAGGGGATGTTTTCTTATGTCTACGTCATCAGCGCGTACCGGCGGTTCACTCGACGCCTg                                                                                                                                                                                                                        <  1:383346/251‑1 (MQ=255)                                   ggTTGAAAAAGCCGCTATGATCGCCGGATAATCGTTTGCTTTTTTTACCACCCGTTTTGTATGCGCGGAGCTAAACGTTTGCTTTTTTGCGACGCAGCAAATTGTCGCAAACCTGGAGCAGGAAGATAACGTTTCGCTGGCAGGGGATTGTCCGCCACGCATCTTGACGAAAATTAAACTCTCAGGGGATGTTTTCTTATGTCTACGTCATCAGCGCGTACCGGCGGTTCACTTGACGCCTGGGTTTaaat                                                                                                                                                                                                               >  1:80851/1‑250 (MQ=255)                                             gCCGCTATGATCGCCGGATAATCGTTTGCTTTTTTTACCACCCGTTTTGTATGCGCGGAGCTAAACGTTTGCTTTTTTGCGACGCAGCAAATTGTCGCAAACCTGGAGCAGGAAGATAACGTTTCGCTGGCAGGGGATTGTCCGCCACGCATCTTGACGAAAATTAAACTCTCAGGGGATGTTTTCTTATGTCTACGTCATCAGCGCGTACCGGCGGTTCACTCGACGCCTGG‑TTTTAAAtttatcaaacg                                                                                                                                                                                                    >  1:15213/1‑243 (MQ=255)                                                     gATCGCCGGATAAGCGTTTTCTTTTTTTACCACCCGTTTTGTATGCGCGGAGCTAAACGTTTGCTTTTTTGCGACGCAGCAAATTGTCGCAAACCTGGAGCAGGAAGATAACGTTTCGCTGGCAGGGGATTGTCCGCCACGCATCTTGACGAAAATTAAACTCTCAGGGGATGTTTTCTTATGTCTACGTCATCAGCGCGTACCGGCGGTTCACTCGACGCCTGG‑TTTAAAATTTCACAACGTG‑GAAGCa                                                                                                                                                                                            <  1:531685/250‑1 (MQ=255)                                                                   cGTTTGCTTTTTTTACCACCCGTTTTGTATGCGCGGAGCTAAACGTTTGCTTTTTTGCGACGCAGCAAATTGTCGCAAACCTGGAGCAGGAAGATAACGTTTCGCTGGCAGGGGATTGTCCGCCACGCATCTTGACGAAAAGTAAACTCTCAGGGGATGTTTTCTTATGTCTACGTCATCAGCGCGTACCGGCGGTTCACTCGACGCCTGG‑TTTAAAATTTCACAACGTG‑GAAGCACTGTCCTTCAGG‑aa                                                                                                                                                                             <  2:132800/250‑1 (MQ=255)                                                                            tttttACCACCCGTTTTGTATGCGCGGAGCTAAACGTTTGCTTTTTTGCGACGCAGCAAATTGTCGCAAACCTGGAGCAGGAAGATAACGTTTCGCTGGCAGGGGATTGTCCGCCACGCATCTTGACGAAAATTAAACTCTCAGGGGATGTTTTCTTATGTCTACGTCATCAGCGCGTACCGGCGGTTCACTCGACGCCTGG‑TTTAAAATTTCACAACGTG‑GAAGCACTGTCCGTCAGG‑AAGTGG‑TTGcc                                                                                                                                                                   <  1:167784/250‑1 (MQ=255)                                                                                                              cgtgtGCTTTTTTGCGACGCAGCAAATTGTCGCAAACCTGGAGCAGGAAGATAACGTTTCGCTGGCAGGGGATTGTCCGCCACGCATCTTGACGAAAATTAAACTCTCAGGGGATGTTTTCTTATGTCTACGTCATCAGCGCGTACCGGCGGTTCACTCGACGCCTGG‑TTTAAAATTTCACAACGTG‑GAAGCACTGTCCGTCAGG‑AAGTGG‑TTGCCGGGTTAACAACGTTTCTGGCGATTGTCTACTCGGt                                                                                                                                <  2:146675/247‑1 (MQ=255)                                                                                                              cGTTTGCTTTTTTGCGACGCAGCAAATTGTCGCAAACCTGGAGCAGGAAGATAACGTTTCGCTGGCAGGGGATTGTCCGCCACGCATCTTGACGAAAATTAAACTCTCAGGGGATGTTTTCTTATGTCTACGTCATCAGCGCGTACCGGCGGTTCACTCGACGCCTGG‑TTTTAAATTTCACAAAGTG‑GAAGAACTGTCCGTTAGGGAAG‑GG‑GTGCCGGGGTAACAAAGTTTTTGGCGATGGTGTACACGGt                                                                                                                                >  1:503710/1‑251 (MQ=255)                                                                                                                 ttGCTTTTTGGCGACGCAGCAAATTGTCGCAAACCTGGAGCAGGAAGATAACGTTTCGCTGGCAGGGGATTGTCCGCCACGCATCTTGACGAAAATTAAACTCTCAGGGGATGTTTTCTTATGTCTACGTCATCAGCGCGTACCGGCGGTTCACTCGACGCCTGG‑TTTAAAATTTCACAACGTG‑GAAGCACTGTCCGTCAGG‑AAGTGG‑TTGCCGGGTTAACAACGTTTCTGGCGATGGTCTACTCGGTCa                                                                                                                              <  1:253973/250‑1 (MQ=255)                                                                                                                  gtcTTTTTTGTGAGGCAGTAAATTGTCGCAAACCTGGAGCAGGAAGAATACGTTTGTCTGGCAGGGTATGGTCCTCCATGCATCCTTAATAAAATTAAACTCTCAGTGGATGTTTTCTTATGTCTTCGTCATCAGCGCGTACCGGCGGTTCACTCGACGCCTGG‑TTTAAAATTTCCCAACGTT‑GAAGCACGGTCCGTCAGG‑AAGTGT‑TTGTCGGGTTAACAATGTTTCTGGCGTTGGTCTTCGCGGTCAt                                                                                                                             <  2:479002/248‑1 (MQ=255)                                                                                                                                                aaCCTGGAGCAGGAAGATAACGTTTCGCTGGCAGGGGATTGTCCGCCACGCATCTTGACGAAAATTAAACTCTCAGGGGATGTTTTCTTATGTCTACGTCATCAGCGCGTACCGGCGGTTTACTCGGCGCCTTG‑TTTTAAATTTTTCACCGTG‑GGAGCACTGTCCGTCAGG‑GAATGG‑TTGCCGGGTTAACAACGTTTCTGGGGGTGGTCTTCTCCGTAATCTTCGTTCCCGGTTTGTTTGGGAAAACAggg                                                                                              >  1:333778/1‑250 (MQ=255)                                                                                                                                                        gCAGGAAGATAACGTTTCGCTGGCAGGGGATTGTCCGCCACGCATCTTGACGAAAATTAAACTCTCAGGGGATGTTTTCTTATGTCTACGTCATCAGCGCGTACCGGCGGTTCACTCGACGCCTGG‑TTTAAAATTTCACAACGTG‑GAAGCACTGTCCGTCAGG‑AAGTGG‑TTGCCGGGTTAACAACGTTTCTGGCGATGGTCTACTCGGTCATCGTCGTTCCAGGTATGTTGGGTAAAGCAGGCTTCCCGcc                                                                                      <  1:14160/251‑1 (MQ=255)                                                                                                                                                                    cGTTCCGCTGGCAGGGGATTGTCCGCCATGCATCTTGACGAAAATTAAACTCTCAGGGGATGTTTTCTTAGGTCTACGTCATCAGCGCGTACCGGCGGTTCACCCGACGCCTGG‑TTAAAAATTTCACAACGGG‑GAAGCACTGTCCGTCAGG‑AAGGGG‑TTGCCGGGTTAACAACGTTTCTGGGGATGGTCTACTGGGTCATCGTCGTTCCAGGTATGTGGGTTAAAGCAGGCTCCCCGCCTGCGGCAGtttt                                                                          <  1:411948/251‑1 (MQ=255)                                                                                                                                                                         cGCTGGCAGGGGATTGTCCGCCACGCATCTTGACGAAAATTAAACTCTCAGGGGATGTTTTCTTATGTCTACGTCATCAGCGCGTACCGGCGGTTCACTCGACGCCTGG‑TTTAAAATTTCACAACGTG‑GAAGCACTGTCCGTCAGG‑AAGTGG‑TTGCCGGGTTAACAACGTTTCTGGCGATGGTCTACTCGGTCATCGTCGTTCCAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTg                                                                     <  1:497572/251‑1 (MQ=255)                                                                                                                                                                                          ccgtcACGCATCTTGACGAAAATTAAACCCTCAGGGTATGTTTTCTTATGTCTGCGTCCTCAGCGCGTACCGGCGGTTCACTCGACGCCTGG‑TTTAAAATTTCACAACGTG‑GAAGCGCTGTCCGCCAGG‑AAGTGT‑TTGCCGGGTTAACATCGTTTCTGGCGATGTTCTTCTCGCTCATCGTCGTTCCAGGTATGTTGGGTAAGGCAGGCTGCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGcc                                                    <  2:166270/247‑1 (MQ=255)                                                                                                                                                                                          ccccccataattttGACTAAAATTAAACTCTCAGGGGATGTTTTCGTCTGTCTACGTCATCAGCGCGTCCCGGGGTTTCACTCGACGCCTGG‑TTTAAAATTGCACAACGTG‑GAATCACTGTCCGTCAGG‑AATTGG‑TTGCCGGGTTAACAACGTTTGTGGCGATGTTCTCCTCGGTCATCGTCGTTCCAGGTATGTGGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGcc                                                    <  2:162442/239‑1 (MQ=255)                                                                                                                                                                                              cACGCATCTTGACGAAAATTAAACTCCCAGGGGATGTTTTCTTATGTCTACGTCACCAGCGCGTACCGGCGGTTCACTCGACGCCTGG‑TTTAAAATTTCACAACGTG‑GAAGCACTGTCCGTCAGG‑AAGTGG‑TTGCCGGGTTAACAACGTTTCTGGCGATGGTCTACTCGGTCATCGTCGTTCCAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGAc                                                <  1:19271/251‑1 (MQ=255)                                                                                                                                                                                                   aTCTTGACGAAAATTAAACTCTCAGGGGATGTTTTCTTATGTCTACGTCATCAGCGCGTACCGGCGGTTCCCCCGCCCCCCGG‑GTTAAAAATTCCCAACGAG‑GAAGCACTGTCCGTCAGG‑GAGTGG‑TTGCCGGGTTTACAACGTTTTTTGCGGTGGTCTACTCGGTCATCGTCGTTCCAGGTATGTTGGGGAAAACAGGCTTCCCGCCCGCGGCCGTTTTCGTTGCAACCTGTCTGGTTGCCGGGCTCggg                                           >  1:544608/1‑250 (MQ=255)                                                                                                                                                                                                    tCTTGACGAAAATTAAACTCTCAGGGGATGTTTTCTTATGTCTACGTCATCAGCGCGTACCGGCGGTTCACTCGGCGCCTGG‑TTTAAAAATTCACAACGTG‑GAAGCACTGTCCGTCAGG‑GAGTGG‑TTGCCGGGGTAACAACGGTTCTGGCGGTGGTCTACTCGGTCCTCCTCGTTCCAGGGGTGTTGGGGAAAGCAGGGCTCCCCCCCGCGGGAGGTTTTGTTTCAAACTGTCTGGTTTCCCGCCTCGGtt                                          >  1:411976/1‑251 (MQ=255)                                                                                                                                                                                                                 tAAACTCTCAGGGGATGTTTTCTTATGTCTACGTCATCAGCGCGTACCGGCGGTTCAATCCACGCCTGG‑GTTAAAATTTCACAACCTG‑GGAGCACAGTCCGGCCGG‑GAGTGGGTTGC‑GGGTTTACAAAGTTTTTGGCGATGGTCTAATCGGTCATCGTCGTTCCAGGGATGTTGGGTTAAGCAGGCCTCCCGCCTGCGGCAGATTTTGTTGCTAACTATCTTGTTTTCGGGCTCGGTTCtttttgtatggg                             >  1:386481/1‑240 (MQ=255)                                                                                                                                                                                                                   aaCTCTCAGGGAAGGTTTTCTGATGCCTACGTCAACAGCCCTTCCCGGCGGTTCACTCGCCGCCTGG‑TTTAAAATTTCACAACGTG‑GAAGCACTGTCCGGCAGG‑AAGTGT‑TTGCCGGTTTAACAACGTTTCTGGCGATGGCCTACTCGGTCATCGTCGTTCCAGGTATGTTGTGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTc                           <  1:299512/251‑1 (MQ=255)                                                                                                                                                                                                                       ctcAGGGGATGTTTTCTTATGTCTACGTCATCAGCGCGTCCCGTCGGTTCACTCGCCGCCGTG‑TTTAAAAGTTCATACAGTG‑GAAGCACGGTCCGTCAGG‑AAGTGG‑TTGCCGGGTTAAAAACGTTGTTGGCGATGGTCTACTCGGTCATCGTCGTTCCAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCtg                         <  2:464173/249‑1 (MQ=255)                                                                                                                                                                                                                         cAGGGGATGTTTTCTTATGTCTACGTCATCAGCGCGTACCGGCGGTTCACTCGACGCCTGG‑TTTAAAATTTCACAACGTG‑GAAGCACTGTCCGTCAGG‑AAGTGG‑TTGCCGGGTTAACAACGTTTCTGGCGATGGTCTACTCGGTCATCGTCGTTCCAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTggg                     <  1:197698/251‑1 (MQ=255)                                                                                                                                                                                                                            gggATTTTTTCTTATGTCTTCGTCATCAGCGCTAACCGGCGGTTCACTCGACGCCTGT‑TTTAAAATTTCACTAC‑TGTGAAGCACTGTCCGTCAGG‑AAGTGG‑TTGCCGGGTTAACAACGGTTCTGTCGATGGTCTCCTCGGTCAGCGTCTTTCCAGCTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTa                  <  2:80851/251‑1 (MQ=255)                                                                                                                                                                                                                            gggATGTTTTCTTATGTCTACGTCATCAGCGCGTACCGGCGGTTCACTCGACGCCTGG‑TTTAAAATTTCACAACGTG‑GAAGCACTGTCCGTCAGG‑AAGTGG‑TTGCCGGGTTAACAACGTTTCTGGCGATGGTCTACTCGGTCATCGTCGTTCCa                                                                                                                   <  2:80100/154‑1 (MQ=255)                                                                                                                                                                                                                                    ttCTTATGTCTACGTCATCAGCGCGTACCGGCGGTTCACTCGACGCCTGG‑TTTAAAATTTCACAACGTG‑GAAGCACTGTCCGTCAGG‑AAGTGG‑TTGCCGGGTTAACAACGTTTCTGGCGATGGTCTACTCGGTCATCGTCGTTCCAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGcc           <  1:233763/250‑1 (MQ=255)                                                                                                                                                                                                                                      cTTATGTCTACGTCATCAGCGCGTACCGGCGGTTCACTCGACGCCTGG‑TTTAAAATTTCCCAACGTG‑GGAGCACTGTCCGTCCGG‑GAGTGG‑TTTCCGGGTTTACAACGTTTCTTGCGATGGTCTACTCGGTCATCGTCGTTCCCGGTATGTTTGGTAAAGCCAGCTTCCCCCCCGCGGCCGTTTTTGTTTCAACCTGTCTGGTTGCCGGACTCCGTTTTTTCGTGGTGGGGCTGTTGGCTTATCTGccctt        >  1:396166/1‑248 (MQ=255)                                                                                                                                                                                                                                         aTGTCTACGTCATCAGCGCGTACCGGCGGTTCACTCGACGCCTGG‑TTTAAAATTTCACAACGTG‑GAAGCACTGCCCGTAAGG‑AAGTGG‑TTGCCGGGTTAACAACGTTTCTGGGGATGGTCTACTCGGTCATCGTCGTTCCAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGc     <  1:90441/251‑1 (MQ=255)                                                                                                                                                                                                                                            tCTACGTCATCAGCGCGTACCGGCGGTTCACTCGGCGCCTGG‑TTTAAAATTTCACAACGTG‑GAAGCACTGTCCCTCAGG‑GAGGGG‑GTGCCGGGTTAACAACGTTTCTGGCGGTGGGCTACTCGGTCATCGTCGTTCCCGGGATGTTGGGTAAAACAGGCTTCCCGCCTGCGGCAGGTTTCGTTTCAACCTGTCTTGTTGCCGGGCTCGGGTCTATCGTGATGGGTCTGTGGGCTGATCTTGCGTTgggggt  >  1:497559/1‑247 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                         CAAGACAATAAAAGAGAAAAAAGCAGCAAACTTCGGTTGAAAAAGCCGCTATGATCGCCGGATAATCGTTTGCTTTTTTTACCACCCGTTTTGTATGCGCGGAGCTAAACGTTTGCTTTTTTGCGACGCAGCAAATTGTCGCAAACCTGGAGCAGGAAGATAACGTTTCGCTGGCAGGGGATTGTCCGCCACGCATCTTGACGAAAATTAAACTCTCAGGGGATGTTTTCTTATGTCTACGCCATCAGCGCGTACCGGCGGTTCACTCGACGCCTGG‑TTTAAAATTTCACAACGTG‑GAAGCACTGTCCGTCAGG‑AAGTGG‑TTGCCGGGTTAACAACGTTTCTGGCGATGGTCTACTCGGTCATCGTCGTTCCAGGTATGTTGGGTAAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGGTTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGAT  >  NC_000913/4278247‑4278732

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