Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,488,314	A→G	N28S (AAC→AGC)	aldA →	aldehyde dehydrogenase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,488,314	0	A	G	100.0%	193.0 / NA	54	N28S (AAC→AGC)	aldA	aldehyde dehydrogenase A
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (31/23);  total (31/23)

TGTTAATTAACAATGTATTCACCGAAAACAAACATATAAATCACAGGAGTCGCCCATGTCAGTACCCGTTCAACATCCTATGTATATCGATGGACAGTTTGTTACCTGGCGTGGAGACGCATGGATTGATGTGGTAAACCCTGCTACAGAGGCTGTCATTTCCCGCATACCCGATGGTCAGGCCGAGGATGCCCGTAAGGCAATCGATGCAGCAGAACGTGCACAACCAGAATGGGAAGCGTTGCCTGCTATTGAACGCGCCAGTTGGT  >  NC_000913/1488177‑1488445                                                                                                                                          |                                                                                                                                    tGTTAATTAACAATGTATTCACCGAAAACAAACATATAAATCACAGGAGTCGCCCATGTCAGTACCCGTTCAACATCCTATGTATATCGATGGACAGTTTATTACCTGTCGTGGAGATGCATTGATTGATGTGGTAAGcc                                                                                                                                   >  2:283890/1‑140 (MQ=255)         aaCAATGTATTCACCTAAAACATACGTATTAATCACACGAGTCGCCCATGTCAGTACCCGTTCAACATCCTATGTATATCGATGGACAGTTTGTTACCTGGCGTGGAGACGCATGGATTGATGTGGTAAGCCCTGCTACa                                                                                                                           <  1:259220/140‑1 (MQ=255)         aaCAATGTATTCACCGAAAACAAACGTATAAATCACAGGAGTCGCCCATGTCAGTACCCGTTCAACATCCTATGTATATCGATGGACAGTTTGTTACCTGGCGTGGAGACGCATGGATTGATGTGGTAAGCCCTGCTACa                                                                                                                           <  1:283890/140‑1 (MQ=255)           cAATGTATTCACCGAAAACAAACATATAAATCACAGGAGTCGCCCATGTCAGTACCCGTTCAACATCCTATGTATATCGATGGACAGTTTGTTACCTGGCGTGGAGACGCATGGATTGATGTGGTAAGCCCTGCTACaga                                                                                                                         >  2:487711/1‑140 (MQ=255)               gTATTCACCGAAAACAAACATATAAATCACAGGAGTCGCCCATGTCAGTACCCGTTCAACATCCTATGTATATCGATGGACAGTTTGTTACCTGGCGTGGAGACGCATGGATTGATGTGGTAAGCCCTGCTACAGAGGCt                                                                                                                     >  1:125832/1‑140 (MQ=255)                   tCACCGAAAACAAACATATAAATCACAGGAGTCGCCCATGTCAGTACCCGTTCAACATCCTATGTATATCGATGGACAGTTTGTTACCTGGCGTGGAGACGCGTGGATTGATGTAGTGAGCCCTGCTACAGAGGCTGTCa                                                                                                                 >  1:249977/1‑140 (MQ=255)                        gAAAACAAACATATAAATCACAGGAGTCGCCCATGTCAGTACCCGTTCAACATCCTATGTATATCGATGGACAGTTTGTTACCTGGCGTGGAGACGCATGGATTGATGTGGTAAGCCCTGCTACAGAGGCTGTCATTTcc                                                                                                            <  2:125832/140‑1 (MQ=255)                         aaaaCAAACATATAAATCACAGGAGTCGCCCATGTCAGTACCCGTTCAACATCCTATGTATATCGATGGACAGTTTGTTACCTGGCGTGGAGACGCATGGATTGATGTGGTAAGCCCTGCTACAGAGGCTGTCATTTccc                                                                                                           >  1:657067/1‑140 (MQ=255)                          aaacaaacATATAAATCACAGGAGTCGCCCATGTCAGTACCCGTTCAACATCCTATGTATATCGATGGACAGTTTGTTACCTGGCGTGGAGACGCATGGATTGATGTGGTAAGCCCTGCTACAGAGGCTGTCATTTCCCg                                                                                                          >  1:720974/1‑140 (MQ=255)                                       aaTCACAGGAGTCGCCCATGTCAGTACCCGTTCAACATCCTATGTATATCGATGGACAGTTTGTTACCTGGCGTGGAGACGCATGGATTGATGTGGTAAGCCCTGCTACAGAGGCTGTCATTTCCCGCATACCCAATGGt                                                                                             <  2:114989/140‑1 (MQ=255)                                        aTCACAGGAGTCGCCCATGTCAGTACCCGTTCAACATCCTATGTATATCGATGGACAGTTTGTTACCTGGCGTGGAGACGCATGGATTGATGTGGTAAGCCCTGCTACAGAGGCTGTCATTTCCCGCATACCCGATGGTc                                                                                            >  1:613927/1‑140 (MQ=255)                                            caGGAGTCGCCCATGTCAGTACCCGTTCAACATCCTATGTATATCGATGGACAGTTTGTTACCTGGCGTGGAGACGCATGGATTGATGTGGTAAGCCCTGCTGCAGAGGCAGTCATTTCCCGCAGACCCGATGGTCAGGc                                                                                        >  1:373903/1‑140 (MQ=255)                                              ggAGTCGCCCATGTCAGTACCCGTTCAACATCCTATGTATATCGATGGACAGTTTGTTACCTGGCGTGGAGACGCATGGATTGATGTGGTAAGCCCTGCTACAGAGGCTGTCATTTCCCGCATACCCGATGGt                                                                                             >  2:81396/1‑133 (MQ=255)                                              ggAGTCGCCCATGTCAGTACCCGTTCAACATCCTATGTATATCGATGGACAGTTTGTTACCTGGCGTGGAGACGCATGGATTGATGTGGTAAGCCCTGCTACAGAGGCTGTCATTTCCCGCATACCCGATGGt                                                                                             <  1:81396/133‑1 (MQ=255)                                              ggAGTCGCCCATGTCAGTACCCGTTCAACATCCTATGTATATCGATGGACAGTTTGTTACCTGGCGTGGAGACGCATGGATTGATGTGGTAAGCCCTGCTACAGAGGCTGTCATTTCCCGCATACCCGATGGTCAGGCCg                                                                                      <  1:487711/140‑1 (MQ=255)                                                  tCGCCCATGTCAGTACCCGTTCAACATCCTATGTATATCGATGGACAGTTTGTTACCTGGCGTGGAGACGCATGGATTGATGTGGTAAGCCCTGCTACAGAGGCTGTCATTTCCCGCATACCCGATGGTCAGGCCGAGGa                                                                                  <  1:734662/140‑1 (MQ=255)                                                   cGCCCATGTCAGTACCCGTTCAACATCCTATGTATATCGATGGACAGTTTGTTACCTGGCGTGGAGACGCATGGATTGATGTGGTAAGCCCTGCTACAGAGGCTGTCATTTCCCGCATACCCGATGGTCAGGCCGAGGAt                                                                                 >  2:25348/1‑140 (MQ=255)                                                   cGCCCATGTCAGTACCCGTTCAACATCCTATGTATATCGATGGACAGTTTGTTACCTGGCGTGGAGACGCATGGATTGATGTGGTAAGCCCTGCTACAGAGGCTGTCATTTCCCGCATACCCGATGGTCAGGCCGAGGAt                                                                                 <  1:426443/140‑1 (MQ=255)                                                   cGCCCATGTCAGTACCCGTTCAACATCCTATGTATATCGATGGACAGTTTGTTACCTGGCGTGGAGACGCATGGATTGATGTGGTAAGCCCTGCTACAGAGGCTGTCATTTCCCGCATACCCGATGGTCAGGCCGAGGAt                                                                                 <  2:71906/140‑1 (MQ=255)                                                    gCCCATGTCAGTACCCGTTCAACATCCTATGTATATCGATGGACAGTTTGTTACCTGGCGTGGAGACGCATGGATTGATGTGGTAAGCCCTGCTACAGAGGCTGTCATTTCCCGCATACCCGATGGTCAGGCCGAGGATg                                                                                >  2:726558/1‑140 (MQ=255)                                                          gTCAGTACCCGTTCAACATCCTATGTATATCGATGGACAGTTTGTTACCTGGCGTGGAGACGCATGGATTGATGTGGTAAGCCCTGCTACAGAGGCTGTCATTTCCCGCATACCCGATGGTCAGGCCGAGGATGCCCGTa                                                                          <  1:853229/140‑1 (MQ=255)                                                                atccGTTCAACATCCTATGGATATCGATGGACAGTTTGTTACCTGGCGTGGAGACGCATGGATTGTTGTGGTAAGCCGTGATACCGAGGCTGTCATTTCCAGAAGCCCCGATGGTCGAGCCGGGGATGCCCGTAGGGCaa                                                                    >  2:301513/3‑140 (MQ=255)                                                                aCCCGTTCAACATCCTATGTATATCGATGGACAGTTTGTTACCTGGCGTGGAGACGCATGGATTGATGTGGTAAGCCCTGCTACAGAGGCTGTCATTTCCCGCATACCCGATGGTCAGGCCGAGGATGCCCGTAAGGCaa                                                                    >  2:5133/1‑140 (MQ=255)                                                                aCCCGTTCAACATCCTATGTATATCGATGGACAGTTTGTTACCTGGCGTGGAGACGCATGGATTGATGTGGTAAGCCCTGCTACAGAGGCTATCATTTCCGGCAAACCTGATGGTCAAGCCGCGGATGCCCGTAAGGCaa                                                                    >  2:315152/1‑140 (MQ=255)                                                                 cccGTTCAACATCCTATGTATATCGATGGACAGTTTGTTACCTGGCGTGGAGACGCATGGATTGATGTGGTAAGCCCTGCTACAGAGGCTGTCATTTCCCGCATACCCGATGGTCAGGCCGAGGATGCCCGTAAGGCAAt                                                                   <  2:720974/140‑1 (MQ=255)                                                                   cGTTCAACATCCTATGTATATCGATGGACAGTTTGTTACCTGGCGTGGAGACGCATGGATTGATGTGGTAAGCCCTGCTACAGAGGCTGTCATTTCCCGCATACCCGATGGTCAGGCCGAGGa                                                                                  >  1:774280/1‑123 (MQ=255)                                                                   cGTTCAACATCCTATGTATATCGATGGACAGTTTGTTACCTGGCGTGGAGACGCATGGATTGATGTGGTAAGCCCTGCTACAGAGGCTGTCATTTCCCGCATACCCGATGGTCAGGCCGAGGa                                                                                  <  2:774280/123‑1 (MQ=255)                                                                     ttCAACATCCTATGTTTATCGAGGGCCAGATTGTTACCTGGCGTGGAGACGCATGGATTGATGTGGTAAGCCCTGCTACAGAGGCTGTCATTTCCCGCATACCCGATGGTCAGGCCGAGGATGCCCGTAAGGCAATCGAt                                                               <  1:301513/140‑1 (MQ=255)                                                                     ttCAACATCCTATGTGTATCGACGGCCAGTTTGTTACCTGGCGTGGAGACGCATGGATTGATGTGGTAAGCCCTGCTACAGAGGCTGTCATTTCCCGCATACCCGATGGTCAGGCCGAGGATGCCCGTAAGGCAATCGAt                                                               <  1:315152/140‑1 (MQ=255)                                                                     ttCAACATCCTATGTATATCGATGGACAGTTTGTTACCTGGCGTGGAGACGCATGGATTGATGTGGTAAGCCCTGCTACAGAGGCTGTCATTTCCCGCATACCCGATGGTCAGGCCGAGGATGCCCGTAAGGCAATCGAt                                                               <  1:5133/140‑1 (MQ=255)                                                                     ttCAACATCCTATGTATATCGATGGACAGTTTGTTACCTGGCGTGGAGACGCATGGATTGATGTGGTAAGCCCTGCTACAGAGGCTGTCATTTCCCGCATACCCGATGGTCAGGCCGAGGATGCCCGTAAGGCAATCGAt                                                               >  1:696542/1‑140 (MQ=255)                                                                             ccTATGTATATCGATGGACAGTTTGTTACCTGGCGTGGAGACGCATGGATTGATGTGGTAAGCCCTGCTACAGAGGCTGTCATTTCCCGCATACCCGATGGTCAGACCGAGGATGCCCGTAAGGCAATCGATGCAGCAGa                                                       <  2:613927/140‑1 (MQ=255)                                                                                      atCGATGGACAGTTTGTTACCTGGCGTGGAGACGCATGGATTGATGTGGTAAGCCCTGCTACAGAGGCTGTCATTTCCCGCATACCCGATGGTCAGGCCGAGGATGCCCGTAAGGCAATCGATGCAGCAGAACGTGcaca                                              >  2:176371/1‑140 (MQ=255)                                                                                      atCGATGGACAGTTTGTTACCTGGCGTGGAGACGCATGGATTGATGTGGTAAGCCCTGCTACAGAGGCTGTCATTTCCCGCATACCCGATGATCAGGCCGAGGATGCCCGTAAGGCAATCGATGCAGCAGAACGTg                                                  <  2:610385/136‑1 (MQ=255)                                                                                      atCGATGGACAGTTTGTTACCTGGCGTGGAGACGCATGGATTGATGTGGTAAGCCCTGCTACAGAGGCTGTCATTTCCCGCATACCCGATGATCAGGCCGAGGATGCCCGTAAGGCAATCGATGCAGCAGAACGTg                                                  >  1:610385/1‑136 (MQ=255)                                                                                        cGATGGACAGTTTGTTACCTGGCGTGGAGACGCATGGATTGATGTGGTAAGCCCTGCTACAGAGGCTGTCATTTCCCGCATACCCGATGGTCAGGCCGAGTATGCCCGGAAGGCAATTGATGGAGCAGAACGTGCACAAc                                            >  2:218710/1‑140 (MQ=255)                                                                                        cGATGGACAGTTTGTTACCTGGCGTGGAGACGCATGGATTGATGTGGTAAGCCCTGCTACAGAGGCTGTCATTTCCCGCATACCCGATGGTCAGGCCGAGGATGCCCGTAAGGCAATCGATGCAGCAGAACGTGCACAAc                                            >  2:201406/1‑140 (MQ=255)                                                                                           tGGACAGTTTGTTACCTGGCGTGGAGACGCATGGATTGATGTGGTAAGCCCTGCTACAGAGGCTGTCATTTCCCGCATACCCGATGGTCAGGCCGAGGATGCCCGTAAGGCAATCGATGCAGCAGAACGTGCACAACCAg                                         >  2:114674/1‑140 (MQ=255)                                                                                           tGGACAGTTTGTTACCTGGCGTGGAGACGCATGGATTGATGTGGTAAGCCCTGCTACAGAGGCTGTCATTTCCCGCATACCCGATGGTCAGGCCGAGGATGCCCGTAAGGCAATCGATGCAGCAGAACGTGCACAACCAg                                         >  1:544510/1‑140 (MQ=255)                                                                                           tGGACAGTTTGTTACCTGGCGTGGAGACGCATGGATTGATGTGGTAAGCCCTGCTACAGAGGCTGTCATTTCCCGCATACCCGATGGTCAGGCCGAGGATGCACGTAAGGGAATCGACGCAGCATAACGTGCACAACCAg                                         >  1:334020/1‑140 (MQ=255)                                                                                           tGGACAGTTTGTTACCTGGCGTGGAGACGCATGGATTGATGTGGTAAGCCCTGCTACAGAGGCTGTCATTTCCCGCATACCCGATGGTCAGGCCGAGGATGCACGTAAGGCAATCGATGCAGCAGAACGTGCACAACCAg                                         >  1:302613/1‑140 (MQ=255)                                                                                           tGGACAGTTTGTTACCTGGCGTGGAGACGCATGGATTGATGTGGTAAGCCCTGCTACAGAGGCTGTCATTTCCCGCATACCCGATGGTCAGGCCGAGGATACCCCTAAGGCAATCGACGCAGCAGAACGTGCACAACCAg                                         >  2:267001/1‑140 (MQ=255)                                                                                             gACAGTTTGTTACCTGGCGTGGAGACGCATGGATTGATGTGGTAAGCCCTGCTACAGAGGCTGTCATTTCCCGCATACCCGATGGTCAGGCCGAGGATGCCCGTAAGGCACTCGGTGCAGCAGAACGTGCACAACCAGaa                                       >  1:333922/1‑140 (MQ=255)                                                                                               cAGTTTGTTACCTGGCGTGGAGACGCATGGATTGATGTGGTAAGCCCTGCTACAGAGGCTGTCATTTCCCGCATACCCGATGGTCAGGCCGCGGATGCCCGTAAGGCAATCGATGCAGCAGAACGCGCACAACCAGAATg                                     >  2:682562/1‑140 (MQ=255)                                                                                                 gTTTGTTACCTGGCGTGGAGACGCATGGATTGATGTGGTAAGCCCTGCTACAGAGGCTGTCATTTCCCGCATACCCGATGGTCAGGCCGAGGATGCCCGTAAGGCAATCGATGCAGCAGAACGTGCACAACCAGAATggg                                   <  1:726558/140‑1 (MQ=255)                                                                                                 gTTTGTTACCTGGCGAGGAGACTCACGGACTGATGTGGTAAGCCCTGCTACAGAGGCTGTCATTTCCCGCATACCCGATGGTCAGGCCGAGGATGCCCGTAAGGCAATCGATGCAGCAGAACGTGCACAACCAGAATggg                                   <  1:267001/140‑1 (MQ=255)                                                                                                  tttGTTACCTGGCGTGGAGACGCATGGATTGATGTGGTAAGCCCTGCTACAGAGGCTGTCATTTCCCGCATACCCGATGGTCAGGCCGAGGATGCCCGTAAGGCAATCGATGCAGCAGAACGTGCACAACCAGAATGGGa                                  >  1:37096/1‑140 (MQ=255)                                                                                                          cTGGCGTGGAGACGCATGGATTGATGTGGTAAGCCCTGCTACAGAGGCTGTCATTTCCCGCATACCCGATGGTCAGGCCGAGGATGCCCGTAAGGCAATCGATGCAGCAGAACGTGCACAACCAGAATGGGAAGCGTTGc                          <  2:37096/140‑1 (MQ=255)                                                                                                            ggCGTGGAGACGCATGGATTGATGTGGTAAGCCCTGCTACAGAGGCTGTCATTTCCCGCATACCCGATGGTCAGGCCGAGGATGCCCGTAAGGCAATCGATGCAGCAGAACGTGCACAACCAGAATGGGAAGCGTTGCCt                        >  1:593468/1‑140 (MQ=255)                                                                                                                tGGAGACGCATGGATTGATGTGGTAAGCCCTGCTACAGAGGCTGTCATTTCCCGCATACCCGATGGTCAGGCCGAGGATGCCCGTAAGGCAATCGATGCAGCAGAACGTGCACAACCAGAATGGGAAGCGTTGCCTGCTa                    >  2:56307/1‑140 (MQ=255)                                                                                                                    gaCGCATGGATTGATGTGGTAAGCCCTGCTACAGAGGCTGTCATTTCCCGCATACCCGATGGTCAGGCCGAGGATGCCCGTAAGGCAATCGATGCAGCAGAACGTGCACATCCAGAAAGGGAAGCGTTGCCTGCTATTGa                >  1:316662/1‑140 (MQ=255)                                                                                                                     aCGCATGGATTGATGTGCTAAGGCCTGCTACCGAGCCTGCCATTTCCCACATACCCGATGGTCAGGCCGAGGATGCCCGTAAGGCAATCGATGCAGCAGAACGTGCACAACCAGAATGGGAAGCGTTGCCTGCTATTGaa               <  2:249977/140‑1 (MQ=255)                                                                                                                       gCATGGATTGATGTGGTAAGCCCTGCTACAGAGGCTGTCATTTCCCGCATACCCGATGGTCAGGCCGAGGATGCCCGTAAGGCAATCGATGCAGCAGAACGTGCACAACCAGAATGGGAAGCGTTGCCTGCTATTGAAcg             <  2:696542/140‑1 (MQ=255)                                                                                                                                 aTGTGGTAAGCCCTGCTACAGAGGCTGTCATTTCCCGCATACCCGATGGTCAGGCCGAGGATGCGCGTAAGGCAATCGATGCAGCAGAACGTGCACAACCAGAATGGGAAGCGTTGCCTGCTATTGAACGCGCCAgttgg   >  2:862283/1‑140 (MQ=255)                                                                                                                                  tgtgGTAAGCCCTGCGACAGAGCCTCTCATTTCCCGCCTACCCGATGGTCAGGCCGAGGATGCCCGTAAGGCAATCGATGCAGCAGAACGTGCACAACCAGAATGGGAAGCGTTGCCTGCTATTGAACGCGCCAgttggt  <  1:224980/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                    TGTTAATTAACAATGTATTCACCGAAAACAAACATATAAATCACAGGAGTCGCCCATGTCAGTACCCGTTCAACATCCTATGTATATCGATGGACAGTTTGTTACCTGGCGTGGAGACGCATGGATTGATGTGGTAAACCCTGCTACAGAGGCTGTCATTTCCCGCATACCCGATGGTCAGGCCGAGGATGCCCGTAAGGCAATCGATGCAGCAGAACGTGCACAACCAGAATGGGAAGCGTTGCCTGCTATTGAACGCGCCAGTTGGT  >  NC_000913/1488177‑1488445

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