Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,630,547	T→C	T123A (ACA→GCA)	ydfI ←	putative oxidoreductase YdfI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,630,547	0	T	C	100.0%	90.9 / NA	26	T123A (ACA→GCA)	ydfI	putative oxidoreductase YdfI
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (14/12);  total (14/12)

GACAACGTCACGCATAACATGACCGTTTTCTGGCATGTTGTCACATGACATGACGGTAAATGCGGGAAGTCCTGCCGCTTTACGGCGAGCCAGCGCCTCAACAATCACCCCTGTTGCTGTTTTCGGCTGGTGGGGATTTTGCACGTCGGCAGCTACCATCGGGTGATCGAGCATTAACTGTCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGTCAGAGAGACAATCGCGATTTGCGGTTCACACATCGCTGCCAACACGGTTTCTAAGCCATCTATCTGTACGTGCAAGGCTTTTTTAACGACGCCAACGACGCGAGCCGTCCACACATCGGCCGACATTTCCGCAACGGTATAAAGATTATCTTGCTGTTGTAAATCGGCAATTTGCTGTTCGCCGCCGATTAAGTTGACCTCATAATATCCCCAGTCACTGAAATGTTCCGTA  >  NC_000913/1630326‑1630779                                                                                                                                                                                                                              |                                                                                                                                                                                                                                         gACAACGTCACGCATAACATGACCGTTTTCTGGCATGTTGTCACATGACATGACGGTAAATGCGGGAAGTCCTGCCGCTTTACGGCGAGCCAGCGCCTCAACAATCACCCCTGTTGCTGTTTTCGGCTGGTGGGGATTTTGCACGTCGGCAGCTACCATCGGGTGATCGAGCATTAACTGTCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGCCAGAGAGACAATCGCGAt                                                                                                                                                                                                                        >  1:399866/1‑240 (MQ=255)    aaCGTCACGCATAACATGACCGTTTGCTGGCATGTTGTCACATGACATGACGGTAAATGCGGGAAGTCCTGCCGCTTTACGGCGAGCCAGCGCCTCAACAATCACCCCTGTTGCTGTTTTCGGCTGGTGGGGATTTTGCACGTCGGCAGCTACCATCGGGTGATCGAGCATTAACTGTCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGCCAGAGAGACAATCGCGAttt                                                                                                                                                                                                                      <  2:597446/239‑1 (MQ=255)      cGTCACGCATAACATGACCGTTTTCTGGCATGTTGTCACATGACATGACGGTAAATGCGGGAAGTCCTGCCGCTTTACGGCGAGCCAGCGCCTCAACAATCACCCCTGTTGCTGTTTTCGGCTGGTGGGGATTTTGCACGTCGGCAGCTACCATCGGGTGATCGAGCATTAACTGTCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGCCAGAGAGACAATCGCGATTTGCg                                                                                                                                                                                                                   >  2:205262/1‑240 (MQ=255)       gTCACGCATAACATGACCGTTTTCTGGCATGTTGTCACATGACATGACGGTAAATGCGGGAAGTCCTGCCGCTTTACGGCGAGCCAGCGCCTCAACAATCACCCCTGTTGCTGTTTTCGGCTGGTGGGGATTTTGCACGTCGGCAGCTACCATCGGGTGATCGAGCATTAACTGTCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGCCAGAGAGACAATCGCGATTTg                                                                                                                                                                                                                     >  1:475388/1‑237 (MQ=255)            gCATAACATGACCGTTTTCTGGCATGTTGTCACATGACATGACGGTAAATGCGGGAAGTCCTGCCGCTTTACGGCGAGCCAGCGCCTCAACAATCACCCCTGTTGCTGTTTTCGGCTGGTGGGGATTTTGCACGTCGGCAGCTACCATCGGGTGATCGAGCATTAACTGTCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGCCAGAGAGACAATCGCGATTTGCGGt                                                                                                                                                                                                                 <  2:364571/236‑1 (MQ=255)                                                                       ccTGCCGCTTTACGGCGAGCCAGCGCCTCAACAATCACCCCTGTTGCTGTTTTCGGCTGGTGGGGATTTTGCACGTCGGCAGCTACCATCGGGTGATCGAGCATTAACTGTCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGCCAGAGAGACAATCGCGATTTGCGGTTCACACATCGCTGCCAACACGGTTTCTAAGCCATCTATCTGTACGTGCAAGGCTTTTTTAac                                                                                                                                                   <  1:166447/239‑1 (MQ=255)                                                                            cGCTTTACGGCGAGCCAGCGCCTCAACAATCACCCCTGTTGCTGTTTTCGGCTGGTGGGGATTTTGCACGTCGGCAGCTACCATCGGGTGATCGAGCATTAACTGTCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGCCAGAGAGACAATCGCGATTTGCGGTTCACACATCGCTGCCAaca                                                                                                                                                                                              <  2:431305/191‑1 (MQ=255)                                                                                               gcCTCAACAATCACCCCTGTTGCTGTTTTCGGCTGGTGGGGATTTTGCACGTCGGCAGCTACCATCGGGTGATCGAGCATTAACTGTCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGCCAGAGAGACAATCGCGATTTGCGGTTCACACATCGCTGCCAacac                                                                                                                                                                                             <  2:450196/173‑1 (MQ=255)                                                                                                         tCACCCCTGTTGCTGTTTTCGGCTGGTGGGGATTTTGCACGTCGGCAGCTACCATCGGGTGATCGAGCATTAACTGTCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGCCAGAGAGACAATCGCGATTTGCGGTTCACACATCGCTGCCAACACGGTTTCTAAGCCATCTATCTGTACGTGCAAGGCTTTTTTAACGACGCCAACGACGCGAGCCGTCCACACATCGGcc                                                                                                                 >  1:504028/1‑239 (MQ=255)                                                                                                         tCACCCCTGTTGCTGTTTTCGGCTGGTGGGGATTTTGCACGTCGGCAGCTACCATCGGGTGATCGAGCATTAACTGTCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGCCAGAGAGACAATCGCGATTTGCGGTTCACACATCGCTGCCAACACGGTTTCTAAGCCATCTATCTGTACGTGCAAGGCTTTTTTAACGACGCCAACGACGCGAGCCGTCCACACATCGGCCg                                                                                                                >  1:267574/1‑240 (MQ=255)                                                                                                                       gTTTTCGGCTGGTGGGGATTTTGCACGTCGGCAGCTACCATCGGGTGATCGAGCATTAACTGTCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGCCAGAGAGACAATCGCGATTTGCGGTTCACACATCGCTGCCAACACGGTTTCTAAGCCATCTATCTGTACGTGCAAGGCTTTTTTAACGACGCCAACGACGCGAGCCGTCCACACATCGGCCGACATTTCCGCAAc                                                                                                   >  1:362322/1‑239 (MQ=255)                                                                                                                         tttCGGCTGGTGGGGATTTTGCACGTCGGCAGCTACCATCGGGTGATCGAGCATTAACTGTCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGCCAGAGAGACAATCGCGATTTGCGGTTCACACATCGCTGCCAACACGGTTTCTa                                                                                                                                                                                     <  2:37782/155‑1 (MQ=255)                                                                                                                         tttCGGCTGGTGGGGATTTTGCACGTCGGCAGCTACCATCGGGTGATCGAGCATTAACTGTCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGCCAGAGAGACAATCGCGATTTGCGGTTCACACATCGCTGCCAACACGGTTTCTAAGCCATCTATCTGTACGTGCAAGGCTTTTTTAACGACGCCAACGACGCGAGCCGTCCACACATCGGCCGa                                                                                                               <  2:365049/225‑1 (MQ=255)                                                                                                                                            tGCACGTCGGCAGCTACCATCGGGTGATCGAGCATTAACTGTCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGCCAGAGAGACAATCGCGATTTGCGGTTCACACATCGCTGCCAacac                                                                                                                                                                                             <  2:186543/128‑1 (MQ=255)                                                                                                                                                      cAGCTACCATCGGGTGATCGAGCATTAACTGTCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGCCAGAGAGACAATCGCGATTTGCGGTTCACACATCGCTGCCAACACGGTTTCTAAGCCATCTATCTGTACGTGCAAGGCTTTTTTAACGACGCCAACGACGCGAGCCGTCCACACATCGGCCGACATTTCCGCaa                                                                                                    <  2:312213/207‑1 (MQ=255)                                                                                                                                                           aCCATCGGGTGATCGAGCATTAACTGTCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGCCAGAGAGACAATCGCGATTTGCGGTTCACACATCGCTGCCAACACGGTTTCTAAGCCATCTATCTGTACGTGCAAGGCTTGTTTAACGACGCCAACGACGCGAGCCGTCCACACATCGGCCGACATTTCCGCAACGGTATAAAGAtt                                                                                       >  1:93354/1‑215 (MQ=255)                                                                                                                                                            ccATCGGGTGATCGAGCATTAACTGTCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGCCAGAGAGACAATCGCGATTTGCGGTTCACACATCGCTGCCAACACGGTTTCTAAGCCATCTATCTGTACGTGCAAGGCTTTTTTAACGACGCCAACGACGCGAGCCGTCCACACATCGGCCGACATTTCCGCAACGGTATAAAGATTATCTTGCTGTTGTAAATCGGCAAttt                                                             >  2:549052/1‑240 (MQ=255)                                                                                                                                                                  ggTGATCGAGCATTAACTGTCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGCCAGAGAGACAATCGCGATTTGCGGTTCACACATCGCTGCCAACACGGTTTCTAAGCCATCTATCTGTACGTGCAAGGCTTTTTTAACGACGCCAACGACGCGAGCCGTCCACACATCGGCCGACATTTCCGCAACGGTATAAAGATTATCTTGCTGTTGTAAATCGGCAATTTGCTGtt                                                       >  2:514809/1‑240 (MQ=255)                                                                                                                                                                   gTGATCGAGCATTAACTGTCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGCCAGAGAGACGATcgc                                                                                                                                                                                                                           >  2:72419/1‑75 (MQ=255)                                                                                                                                                                   gTGATCGAGCATTAACTGTCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGCCAGAGAGACAATCGCGATTTGCGGTTCACACATCGCTGCCAACACGGTTTCTAAGCCATCTATCTGTACGTGCAAGGCTTTTTTAACGACGCCAACGACGCGAGCCGTCCACACATCGGCCGACATTTCCGCAACGGTATAAAGATTATCTTGCTGTTGTAAATCGGCAATTTGCTGtt                                                       >  1:471910/1‑239 (MQ=255)                                                                                                                                                                                  cTGTCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGCCAGAGAGACAATCGCGATTTGCGGTTCACACATCGCTGCCAACACGGTTTCTAAGCCATCTATCTGTACGTGCAAGGCTTTTTTAACGACGCCAACGACGCGAGCCGTCCACACATCGGCCGACATTTCCGCAACGGTATAAAGATTATCTTGCTGTTGTAAATCGGCAATTTGCTGTTCGCCGCCGATTAAGtt                                       <  2:428198/240‑1 (MQ=255)                                                                                                                                                                                     tCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGCCAGAGAGACAATCGCGATTTGCGGTTCACACATCGCTGCCAACACGGTTTCTAAGCCATCTATCTGTACGTGCAAGGCTTTTTTAACGACGCCAACGACGCGAGCCGTCCACACATCGGCCGACATTTCCGCAACGGTATAAAGATTATCt                                                                                   >  2:379792/1‑193 (MQ=255)                                                                                                                                                                                      ccGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGCCAGAGAGACAATCGCGATTTGCGGTTCACACATCGCTGCCAACACGGTTTCTAAGCCATCTATCTGTACGTGCAAGGCTTTTTTAACGACGCCAACGACGCGAGCCGTCCACACATCGGCCGACATTTCCGCAACGGTATAAAGATTATCTTGCTGTTGTAAATCGGCAATTTGCTGTTCGCCGCCGATTAAGTTGAc                                    <  2:405324/239‑1 (MQ=255)                                                                                                                                                                                               ggAGAGTGGAAATACCCTTTTTCGGTGATTGCCAGAGAGACAATCGCGATTTGCGGTTCACACATCGCTGCCAACACGGTTTCTAAGCCATCTATCTGTACGTGCAAGGCTTTTTTAACGACGCCAACGACGCGAGCCGTCCACACATCGGCCGACATTTCCGCAACGGTATAAAGATTATCTTGCTGTTGTAAATCGGCAATTTGCTGTTCGCCGCCGATTAAGTTGACCTCataata                           >  2:382136/1‑239 (MQ=255)                                                                                                                                                                                                                 ttttCGGTGATTGCCAGAGAGACAATCGCGATTTGCGGTTCACACATCGCTGCCAACACGGTTTCTAAGCCATCTATCTGTACGTGCAAGGCTTTTTTAACGACGCCAACGACGCGAGCCGTCCACACATCGGCCGACATTTCCGCAACGGTATAAAGATTATCTTGCTGTTGTAAATCGGCAATTTGCTGTTCGCCGCCGATTAAGTTGACCTCATAATATCCCCAGTCa                 <  2:399858/231‑1 (MQ=255)                                                                                                                                                                                                                        tGATTGCCAGAGAGACAATCGCGATTTGCGGTTCACACATCGCTGCCAACACGGTTTCTAAGCCATCTATCTGTACGTGCAAGGCTTTTTTAACGACGCCAACGACGCGAGCCGTCCACACATCGGCCGACATTTCCGCAACGGTATAAAGATTATCTTGCTGTTGTAAATCGGCAATTTGCTGTTCGCCGCCGATTAAGTTGACCTCATAATATCCCCAGTCACTGAAATGGTCCGTa  >  2:366138/1‑239 (MQ=255)                                                                                                                                                                                                                              |                                                                                                                                                                                                                                         GACAACGTCACGCATAACATGACCGTTTTCTGGCATGTTGTCACATGACATGACGGTAAATGCGGGAAGTCCTGCCGCTTTACGGCGAGCCAGCGCCTCAACAATCACCCCTGTTGCTGTTTTCGGCTGGTGGGGATTTTGCACGTCGGCAGCTACCATCGGGTGATCGAGCATTAACTGTCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGTCAGAGAGACAATCGCGATTTGCGGTTCACACATCGCTGCCAACACGGTTTCTAAGCCATCTATCTGTACGTGCAAGGCTTTTTTAACGACGCCAACGACGCGAGCCGTCCACACATCGGCCGACATTTCCGCAACGGTATAAAGATTATCTTGCTGTTGTAAATCGGCAATTTGCTGTTCGCCGCCGATTAAGTTGACCTCATAATATCCCCAGTCACTGAAATGTTCCGTA  >  NC_000913/1630326‑1630779

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