Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,828,362	+A	coding (1373/1479 nt)	astD ←	aldehyde dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,828,362	1	.	A	100.0%	107.0 / NA	33	A458V (GCC→GTC)	astD	aldehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base A (18/15);  total (18/15)

GTCGAAATTGACTTCCCAGGCGTTCATCGCACCACCTCATCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGTCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGG‑CATACCAGGCGCTGGGGCGATGGTTACCGGAAGCACCAATGCCGCCGAATGGCGCGGTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGTTAACAATCCCCGCCCGCGCCTCC  >  NC_000913/1828233‑1828475                                                                                                                                   |                                                                                                                  gTCGAAATTGACTTCCCAGGCGTTCATCGCACCACCTCATCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGTCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGACATACCAgg                                                                                                          >  2:1224186/1‑140 (MQ=255)  tCGAAATTGACTTCCCAGGCGTTCATCGCACCACCTCATCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGTCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGACATACCAGgc                                                                                                         >  1:267452/1‑140 (MQ=255)  tCGAAATTGACTTCCCAGGCGTTAATCGCACCACCTCATCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGTCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGACATACCAGgc                                                                                                         >  2:532595/1‑140 (MQ=255)     aaaTTGACTTCCCAGGCGTTCATCGCACCACCTCATCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGTCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGACATACCAGGCGCt                                                                                                      <  1:302588/140‑1 (MQ=255)     aaaTTGACTTCCCAGGCGTTCATCGCACCACCTCATCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGTCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGACATACCAGGCGCt                                                                                                      <  1:522132/140‑1 (MQ=255)         tGACTTCCCAGGCGTTCATCGCACCACCTCATCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGTCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGACATACCAGGCGCTgggg                                                                                                  <  2:868407/140‑1 (MQ=255)            cTTCCCAGGCGTTCATCGCACCACCTCATCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGTCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGACATACCAGGCGCTGGGGCGa                                                                                               <  1:905103/140‑1 (MQ=255)                         cATCGCACCACCTCATCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGTCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGACATACCAGGCGCTGGGGCGATGGTTACCGGAAg                                                                                  >  1:1295923/1‑140 (MQ=255)                         cATCGCACCACCTCATCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGTCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGACATACCAGGCGCTGGGGCGATGGTTACCGGAAg                                                                                  >  1:300834/1‑140 (MQ=255)                                ccaccTCATCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGTCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGACATACCAGGCGCTGGGGCGATGGTTACCGGAAGCACCAAt                                                                           >  2:1214672/1‑140 (MQ=255)                                      cATCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGTCCGACTCCAGAATCGCCATCGGCCATGCGCAGTAATCTGCGGACATACCAGGCGCTGGGGCGATGGTTACCGGAAGCACCAATgccgcc                                                                     >  2:1268609/1‑140 (MQ=255)                                             aaaaTCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGTCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGACATACCAGGCGCTGGGGCGATGGTTACCGGAAGCACCAATGCCGCGGAATGgc                                                              <  1:763275/140‑1 (MQ=255)                                                aTCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGTCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGACATACCAGGCGCTGGGGCGATGGTTACCGGAAGCACCAATGCCGCCGAATGGcgcg                                                           <  1:532595/140‑1 (MQ=255)                                                                    tGGCGGGCAATGTTAACGAGTCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGACATACCAGGCGCTGGGGCGATGGTTACCGGAAGCACCAATGCCGCCGAATGGCGCGGTACTGGCAGCACCGGTAAg                                       <  1:1268609/140‑1 (MQ=255)                                                                     ggCGGGCAATGTTAACGAGTCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGACATACCAGGCGCTGGGGCGATGGTTACCGGAAGCACCAATGCCGCCGAATGGCGCGGTACTgg                                                    >  2:523593/1‑126 (MQ=255)                                                                     ggCGGGCAATGTTAACGAGTCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGACATACCAGGCGCTGGGGCGATGGTTACCGGAAGCACCAATGCCGCCGAATGGCGCGGTACTgg                                                    <  1:523593/126‑1 (MQ=255)                                                                     ggCGGGCAATGTTAACGAGTCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGACATACCAGGCGCTGGGGCGATGGTTACCGGAAGCACCAATGCCGCCGAATGGCGCGGTACTGGCAGCACCGGTAAGc                                      >  1:301348/1‑140 (MQ=255)                                                                     ggCGGGCAATGTTAACGAGTCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGACATACCAGGCGCTGGGGCGATGGTTACCGGAAGCACCAATGCCGCCGAATGGCGCGGTACTGGCAGCACCGGTAAGc                                      >  1:267162/1‑140 (MQ=255)                                                                     ggCGGGCAATGTTAACGAGTCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGACATACCAGGCGCTGGGGCGATGGTTACCGGAAGCACCAATGCCGCCGAATGGCGCGGTACTGGCAGCACCGGTAAGc                                      >  2:1416374/1‑140 (MQ=255)                                                                     ggCGGGCAATGTTAACGAGTCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGACATACCAGGCGCTGGGGCGATGGTTACCGGAAGCACCAATGCCGCCGAATGGCGCGGTACTGGCAGCACCGGTAAGc                                      >  2:16942/1‑140 (MQ=255)                                                                     ggCGGGCAATGTTAACGAGTCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGACATACCAGGCGCT‑GGGCGATGGTTACCGGAAGCACCAATGCCGCCGAATGGCGCGGTACTGGCAGCACCGGTAAGCg                                     >  2:806704/1‑140 (MQ=255)                                                                      gCGGGCAATGTTAACGAGTCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGACATACCAGGCGCTGGGGCGATGGTTACCGGAAGCACCAATGCCGCCGAATGGCGCGGTACTGGCAGCACCGGTAAGCg                                     <  2:1124042/140‑1 (MQ=255)                                                                       cGGGCAATGTTAACGAGTCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGACATACCAGGCGCTGGGGCGATGGTTACCGGAAGCACCAATGCCGCCGAATGGCGCGGTACTGGCAGCACCGGTAAGCgg                                    <  1:656059/140‑1 (MQ=255)                                                                               gTTAACGAGTCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGACATACCAGGCGCT‑GGGCGATGGTTACCGGAAGCACCAATGCCGCCGAATGGCGCGGTACTGGCAGCACCGGTAAGCGGTTTGTTCCa                           <  1:806704/140‑1 (MQ=255)                                                                                   aCGAGTCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGACATACCAGGCGCTGGGGCGATGGTTACCGGAAGCACCAATGCCGCCGAATGGCGCGGTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGtt                        >  2:473514/1‑140 (MQ=255)                                                                                        tCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGACATACCAGGCGCTGGGGCGATGGTTACCGGAAGCACCAATGCCGCCGAATGGCGCGGTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGTTAACaa                   >  1:569521/1‑140 (MQ=255)                                                                                            aCTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGACATACCAGGCGCTGGGGCGATGGTTACCGGAAGCACCAATGCCGCCGAATGGCGCGGTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGTTAACAATccc               >  2:578511/1‑140 (MQ=255)                                                                                                 aGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGACATACCAGGCGCTGGGGCGATGGTTACCGGAAGCACCAATGCCGCCGAATGGCGCGGTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGTTAACAATCcccgccc          >  2:922935/1‑140 (MQ=255)                                                                                                    cTCGCCATCGGCCATGCGCAGTAATCTGCGGACATACCAGGCGCTGGGGCGATGGTTACCGGAAGCACCAATGCCGCCGAATGGCGCGGTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGTTAACAATCCCCGCCcgcg       <  1:1416374/140‑1 (MQ=255)                                                                                                       gCCATCGGCCATGCGCAGTAATCTGCGGACATACCAGGCGCTGGGGCGATGGTTACCGGAAGCACCAATGCCGCCGAATGGCGCGGTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGTTAACAATCCCCGCCCGCGCCt    <  1:16942/140‑1 (MQ=255)                                                                                                       gCCATCGGCCATGCGCAGTAATCTGCGGACATACCAGGCGCTGGGGCGATGGTTACCGGAAGCACCAATGCCGCCGAATGGCGCGGTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGTTAACAATCCCCGCCCGCGCCt    <  1:1224186/140‑1 (MQ=255)                                                                                                       gCCATCGGCCATGCGCAGTAATCTGCGGACATACCAGGCGCTGGGGCGATGGTTACCGGAAGCACCAATGCCGCCGAATGGCGCGGTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGTTAACAATCCCCGCCCGCGCCt    >  1:1127760/1‑140 (MQ=255)                                                                                                            cGGCCATGCGCAGTAATCTGCGGACATACCAGGCGCTGGGGCGATGGTTACCGGAAGCACCAATGCCGCCGAATGGCGCGGTACTGGCAGCCCCGGTAGGCGGTTTGTTCCAGTTAACAATCCCCGCCCGCGCCTcc  <  2:661661/137‑1 (MQ=255)                                                                                                            cGGCCATGCGCAGTAATCTGCGGACATACCAGGCGCTGGGGCGATGGTTACCGGAAGCACCAATGCCGCCGAATGGCGCGGTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGTTAACAATCCCCGCCCGCGCCTcc  >  1:661661/1‑137 (MQ=255)                                                                                                                                   |                                                                                                                  GTCGAAATTGACTTCCCAGGCGTTCATCGCACCACCTCATCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGTCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGG‑CATACCAGGCGCTGGGGCGATGGTTACCGGAAGCACCAATGCCGCCGAATGGCGCGGTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGTTAACAATCCCCGCCCGCGCCTCC  >  NC_000913/1828233‑1828475

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