Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,130,161	G→A	L19F (CTT→TTT)	yghO ←	putative DNA‑binding transcriptional regulator YghO

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,130,161	0	G	A	100.0%	73.1 / NA	24	L19F (CTT→TTT)	yghO	putative DNA‑binding transcriptional regulator YghO
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (15/9);  total (15/9)

CCGTAACGCTCGCGGTGCAAGGTATCTATTTGCGCGGTAATACGCCCAACTATCTGCCCCGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACATGGTGGAAGTCATCGAGCCGTTTTACCCCAAGGCTGGTAATGGCCGGCGACCTTATCCGCTG  >  NC_000913/3129923‑3130384                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                ccGTAACGCTCGCGGTGCAAGGTATCTATTTGCGCGGTAATACGCCCAACTATCTGCCCCGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAATCCCTGATATg                                                                                                                                                                                                                      >  1:204937/1‑250 (MQ=255)                                gcgcGGTAATACGCCCAACTATCTGCCCCGTTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAATCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACa                                                                                                                                                                                      <  1:339708/251‑1 (MQ=255)                                         tACGCCCAACTATCTGCCCCGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGAGCGAACAAATCCCTGATATGAGATCCTGTTTGTCATCTGGAGCCATAGAACAGGGTtcatc                                                                                                                                                                             >  2:195553/1‑251 (MQ=255)                                               cAACTATCTGCCCCGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAATCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGt                                                                                                                                                                       <  1:357850/251‑1 (MQ=255)                                                       tGCCCCGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGAGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAATCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATACGtcatcacat                                                                                                                                                               >  1:329278/1‑248 (MQ=255)                                                               tttttttGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAATCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCg                                                                                                                                                       >  1:219599/1‑251 (MQ=255)                                                                                                  ccgtccCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAATCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGAc                                                                                                                     <  2:390005/250‑1 (MQ=255)                                                                                                                  cGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAATCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCt                                                                                                    <  1:195553/251‑1 (MQ=255)                                                                                                                      gACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAATCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTc                                                                                                <  2:158124/251‑1 (MQ=255)                                                                                                                       taaaaTGCTCGTTGCGTGCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGACCAAATCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTAATCAACTGACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGAACAAAAACGAGATTTTCTTGTcc                                                                                               >  2:42195/3‑251 (MQ=21)                                                                                                                             gCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAATCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCAt                                                                                          <  1:81501/250‑1 (MQ=255)                                                                                                                              cTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAATCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATgg                                                                                        >  1:332276/1‑251 (MQ=255)                                                                                                                                    gCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGAGGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAATCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTGTCTTGTCCCGCATGGAGCAg                                                                                   <  2:288686/250‑1 (MQ=255)                                                                                                                                        tCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAATCCATGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCt                                                                              >  1:59426/1‑251 (MQ=255)                                                                                                                                                     ggggTATCCAGTTTGGCTAGTCGGGATAAAGTGATGACGGAAAAGCGATAACTGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAATCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTAAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAAc                                                                 >  2:253150/1‑251 (MQ=255)                                                                                                                                                       ggTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAATCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACAt                                                               <  1:196685/251‑1 (MQ=255)                                                                                                                                                                                 aaGTGATGACGGAAAAGCGATTAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAATCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGGCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACATGGTGGAAGTAATCGGGCCGtatgacc                                     >  2:313598/1‑245 (MQ=14)                                                                                                                                                                                  aGTGATGACGGAAAAGCGATAACTGCCTTGCGGTCATTTTTATTAAGGAAGGTGAGAACAAATACCTGATATGAGATACTGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACCGAGAATTCAGCAGTAAGCGCCGTCCGACCCGAAAAGAGATTTTCTTGTCACGCATGGAGCAGATTCTGCAATGGCAAAACATGGTGGAAGTCATCGAGCCGTTTTcccc                                    >  2:175322/1‑251 (MQ=14)                                                                                                                                                                                    tgatgaCGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAATCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACATGGTGGAAGTCATCGAGCCGTTTTAcccc                                   >  1:220502/1‑250 (MQ=255)                                                                                                                                                                                          cGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAATCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACATGGTGGAAGTCATCGAGCCGTTTTACCCCAAGGCTg                            >  1:94268/1‑251 (MQ=255)                                                                                                                                                                                                   gATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAATCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGaaaa                                                                                                              <  1:364469/160‑1 (MQ=255)                                                                                                                                                                                                   gATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAATCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGaaaa                                                                                                              >  2:364469/1‑160 (MQ=255)                                                                                                                                                                                                    aTAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAATCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACATGGTGGAAGTCATCGAGCCGTTTTACCCCAAGGCTGGTAATGGCCg                  >  1:264333/1‑251 (MQ=255)                                                                                                                                                                                                                    cATTTTTATTAAGGAAGGTGCGAACAAATCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACATGGTGGAAGTCATCGAGCCGTTTTACCCCAAGGCTGGTAATGGCCGGCGACCTTATCCGCTg  >  1:62866/1‑251 (MQ=17)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                CCGTAACGCTCGCGGTGCAAGGTATCTATTTGCGCGGTAATACGCCCAACTATCTGCCCCGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACATGGTGGAAGTCATCGAGCCGTTTTACCCCAAGGCTGGTAATGGCCGGCGACCTTATCCGCTG  >  NC_000913/3129923‑3130384

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