Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,266,937	C→A	intergenic (+76/‑177)	alr → / → tyrB	alanine racemase 1/tyrosine aminotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,266,937	0	C	A	100.0%	160.6 / NA	47	intergenic (+76/‑177)	alr/tyrB	alanine racemase 1/tyrosine aminotransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (25/22);  total (25/22)

CGCTTACGAACTTATTACGCGCCTGACTTCAAGGGTCGCGATGAAATACGTGGATTAATCGTTCTGTAATATTTGATTGTCTGTGCCGGATGCGGCGTGAATGCCTTATCCGGCCAATAAAATCCTAAAAATTCAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCGTGCCTCTGGCGGATGTACGTTTGTCATGAGTCTCACTCTGTTGCTAATTGCCGTTCGCTCCTGAACATCCACTCGATCTTCGCCTTCTTCCGGTTTATT  >  NC_000913/4266804‑4267073                                                                                                                                      |                                                                                                                                         cgcTTACGAACTTATTACGCGCCTGACTTCAAGGGTCGCGATGAAATACGTGGATTAATCGTTCTGTAATATTTGATTGTCTGTGCCGGATGCGGCGTGAATGCCTTATCCGGCCAATAAAATCCTAAAAATTAAATAAg                                                                                                                                    <  2:583326/140‑1 (MQ=255)      aCGAACTTATTACGCGCCTGACTTCAAGGGTCGCGATGAAATACGTGGATTAATCGTTCTGTAATATTTGATTGTCTGTGCCGGATGCGGCGTGAATGCCTTATCCGGCCAATAAAATCCTAAAAATTAAATAAGTTGAt                                                                                                                               >  1:836793/1‑140 (MQ=255)       cGAACTTATTACGCGCCTGACTTCAAGGGTCGCGATGAAATACGTGGATTAATCGTTCTGTAATATTTGATTGTCTGTGCCGGATGCGGCGTGAATGCCTTATCCGGCCAATAAAATCCTAAAAATTAAATAAGTTGATg                                                                                                                              <  2:836793/140‑1 (MQ=255)              attaCGCGCCTGACTTCAAGGGTCGCGATGAAATACGTGGATTAATCGTTCTGTAATATTTGATTGTCTGTGCCGGATGCGGCGTGAATGCCTTATCCGGCCAATAAAATCCTAAAAATTAAATAAGTTGATGTTCTTTc                                                                                                                       >  1:746678/1‑140 (MQ=255)              attaCGCGCCTGACTTCAAGGGTCGCGATGAAATACGTGGATTAATCGTTCTGTAATATTTGATTGTCTGTGCCGGATGCGGCGTGAATGCCTTATCCGGCCAATAAAATCCTAAAAATTAAATAAGTTGATGTTCTTTc                                                                                                                       >  1:218258/1‑140 (MQ=255)               ttaCGTGCCTGACTTCAAGGGTGGCGATGAAATACGTGGATTAATCGTTCTGTAATATTTGATTGTCTGTGCCGGATGCGGCGTGAATGCCTTATCCGGCCAATAAAATCCTAAAAATTAAATAAGTTGATGTTCTTTCa                                                                                                                      <  2:661228/140‑1 (MQ=255)                 aCGCGCCTGACTTCAAGGGTCGCGATGAAATACGTGGATTAATCGTTCTGTAATATTTGATTGTCTGTGCCGGATGCGGCGTGAATGCCTTATCCGGCCAATAAAATCCTAAAAATTAAATAAGTTGATGTTCTTTCATg                                                                                                                    <  1:3235/140‑1 (MQ=255)                 aCGCGCCTGACTTCAAGGGTCGCGATGAAATACGTGGATTAATCGTTCTGTAATATTTGATTGTCTGTGCCGGATGCGGCGTGAATGCCTTATCCGGCCAATAAAATCCTAAAAATTAAAAAAGTTGATGTTCTTTCATg                                                                                                                    <  2:790741/140‑1 (MQ=255)                   gcgcCTGACTTCAAGGGTCGCGATGAAATACGTGGATTAATCGTTCTGTAATATTTGATTGTCTGTGCCGGATGCGGCGTGAATGCCTTATCCGGCCAATAAAATCCTAAAAATTAAATAAGTTGATGTTCTTTCATGct                                                                                                                  <  2:218258/140‑1 (MQ=255)                    cgcCTGACTTCAAGGGTCGCGATGAAATACGTGGATTAATCGTTCTGTAATATTTGATTGTCTGTGCCGGATGCGGCGTGAATGCCTTATCCGGCCAATAAAATCCTAAAAATTAAATAAGTTGATGTTCTTTCATGctc                                                                                                                 >  1:166988/1‑140 (MQ=255)                    cgcCTGACTTCAAGGGTCGCGATGAAATACGTGGATTAATCGTTCTGTAATATTTGATTGTCTGTGCCGGATGCGGCGTGAATGCCTTATCCGGCCAATAAAATCCTAAAAATTAAATAAGTTGATGTTCTTTCATGctc                                                                                                                 >  1:44979/1‑140 (MQ=255)                    cgcCTGACTTCAAGGGTCGCGATGAAATACGTGGATTAATCGTTCTGTAATATTTGATTGTCTGTGCCGGATGCGGCGTGAATGCCTTATCCGGCCAATAAAATCCTAAAAATTAAATAAGTTGATGTTCTTTCATGctc                                                                                                                 >  2:422043/1‑140 (MQ=255)                                     cgcgATGAAATACGTGGATTAATCGTTCTGTAATATTTGATTGTCTGTGCCGGATGCGGCGTGAATGCCTTATCCGGCCAATAAAATCCTAAAAATTAAATAAGTTGATGTTCTTTCATGct                                                                                                                  >  2:889205/1‑122 (MQ=255)                                     cgcgATGAAATACGTGGATTAATCGTTCTGTAATATTTGATTGTCTGTGCCGGATGCGGCGTGAATGCCTTATCCGGCCAATAAAATCCTAAAAATTAAATAAGTTGATGTTCTTTCATGct                                                                                                                  <  1:889205/122‑1 (MQ=255)                                      gcgATGAAATACGTGGATTAATCGTTCTGTAATATTTGATTGTCTGTGCCGGATGCGGCGTGAATGCCTTATCCGGCCAATAAAATCCTAAAAATTAAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCGTGCctc                                                                                               >  1:139898/1‑140 (MQ=255)                                             aaTACGTGGATTAATCGTTCTGTAATATTTGATTGTCTGTGCCGGATGCGGCGTGAATGCCTTATCCGGCCAATAAAATCCTAAAAATTAAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCGTGCCTCTGGCGGa                                                                                        >  2:545180/1‑140 (MQ=255)                                                 cGTGGATTAATCGTTCTGTAATATTTGATTGTCTGTGCCGGATGCGGCGTGAATGCCTTATCCGGCCAATAAAATCCTAAAAATTAAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCGTGCCTCTGGCGGATGTa                                                                                    <  2:208399/140‑1 (MQ=255)                                                 cGTGGATTAATCGTTCTGTAATATTTGATTGTCTGTGCCGGATGCGGCGTGAATGCCTTATCCGGCCAATAAAATCCTAAAAATTAAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCGTGCCTCTGGCGGATGTa                                                                                    <  2:326732/140‑1 (MQ=255)                                                      ataaaTCGTTCTGTAATATTTGATTGTCTGTGCCGGATGCGGCGTGAATGCCTTATCCGGCCAATAAAATCCTAAAAATTAAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCGTGCCTCTGGCGGATGTACGttt                                                                               >  2:611075/4‑140 (MQ=255)                                                          aTCGTTCTGTAATATTTGATTGTCTGTGCCGGATGCGGCGTGAATGCCTTATCCGGCCAATAAAATCCTAAAAATTAAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCGTGCCTCTGGCGGATGTACGTTTGTCa                                                                           <  2:44979/140‑1 (MQ=255)                                                           tCGTTCTGTAATATTTGATTGTCTGTGCCGGATGCGGCGTGAATGCCTTATCCGGCCAATAAAATCCTAAAAATTAAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCGTGCCTCTGGCGGATGTACGTTTGTCAt                                                                          <  2:139898/140‑1 (MQ=255)                                                           tCGTTCTGTAATATTTGATTGTCTGTGCCGGATGCGGCGTGAATGCCTTATCCGGCCAATAAAATCCTAAAAATTAAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCGTGCCTCTGGCGGATGTACGTTTGTCAt                                                                          <  1:545180/140‑1 (MQ=255)                                                           tCGTTCTGTAATATTTGATTGTCTGTGCCGGATGCGGCGTGAATGCCTTATCCGGCCAATAAAATCCTAAAAATTAAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCGTGCCTCTGGCGGATGTACGTTTGTCAt                                                                          <  2:263683/140‑1 (MQ=255)                                                             gTTCTGTAATATTTGATTGTCTGTGCCGGATGCGGCGTGAATGCCTTATCCGGCCAATAAAATCCTAAAAATTAAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCGTGCCTCTGGCGGATGTACGTTTGTCATGa                                                                        <  1:422043/140‑1 (MQ=255)                                                                      tatTTGATTGTCTGTGCCGGATGCGGCGTGAATGCCTTATCCGGCCAATAAAATCCTAAAAATTAAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCGTGCCTCTGGCGGATGTACGTTTGTCATGAGTCTCActc                                                               >  2:755322/1‑140 (MQ=255)                                                                           gATTGTCTGTGCCGGATGCGGCGTGAATGCCTTATCCGGCCAATAAAATCCTAAAAATTAAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCGTGCCTCTGGCGGATGTACGTTTGTCATGAGTCTCAct                                                                >  1:650174/1‑134 (MQ=255)                                                                           gATTGTCTGTGCCGGATGCGGCGTGAATGCCTTATCCGGCCAATAAAATCCTAAAAATTAAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCGTGCCTCTGGCGGATGTACGTTTGTCATGAGTCTCAct                                                                <  2:650174/134‑1 (MQ=255)                                                                             ttGTCTGTGCCGGATGCGGCGTGAATGCCTTATCCGGCCAATAAAATCCTAAAAATTAAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCGTGCCTCTGGCGGATGTACGTTTGTCATGAGTCTCACTCTGTTGCt                                                        >  1:138414/1‑140 (MQ=255)                                                                                 cTGTGCCGGATGCGGCGTGAATGCCTTATCCGGCCAATAAAATCCTAAAAATTAAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCGTGCCTCTGGCGGATGTACGTTTGTCATGAGTCTCACTCTGTTGCTAAtt                                                    >  1:75088/1‑140 (MQ=255)                                                                                    tgCCGGATGCGGCGTGAATGCCTTATCCGGCCAATAAAATCCTAAAAATTAAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCGTGCCTCTGGCGGATGTACGTTTGTCATGAGTCTCACTCTGTTGCTAATTGcc                                                 >  2:836541/1‑140 (MQ=255)                                                                                       cGGATGCGGCGTGAATGCCTTATCCGGCCAATAAAATCCTAAAAATTAAATAAGTTGATGTTCTTTCATGCTTTTATAAAGGTCGTGCCTCTGGCGGATGTACGTTTGTCATGAGTCTCACTCTGTTGCTAATTGCCGtt                                              <  2:792845/140‑1 (MQ=255)                                                                                          aTGCGGCGTGAATGCCTTATCCGGCCAATAAAATCCTAAAAATTAAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCGTGCCTCTGGCGGATGTACGTTTGTCATGAGTCTCACTCTGTTGCTAATTGCCGTTCGc                                           >  1:203323/1‑140 (MQ=255)                                                                                          aTGCGGCGTGAATGCCTTATCCGGCCAATAAAATCCTAAAAATTAAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCGTGCCTCTGGCGGATGTACGTTTGTCATGAGTCTCACTCTGTTGCTAATTGCCGTTCGc                                           >  2:664461/1‑140 (MQ=255)                                                                                                cgTGAATGCCTTATCCGGCCAATAAAATCCTAAAAATTAAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCGTGCCTCTGGCGGATGTACGTTTGTCATGAGTCTCACTCTGTTGCTAATTGCCGTTCGCTCCTGa                                     >  2:644915/1‑140 (MQ=255)                                                                                                         cTTATCCGGCCAATAAAATCCTAAAAATTAAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCGTGCCTCTGGCGGATGTACGTTTGTCATGAGTCTCACTCTGTTGCTAATTGCCGTTCGCTCCTGAACATCCACt                            <  2:203323/140‑1 (MQ=255)                                                                                                         cTTATCCGGCCAATAAAATCCTAAAAATTAAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCGTGCCTCTGGCGGATGTACGTTTGTCATGAGTCTCACTCTGTTGCTAATTGCCGTTCGCTCCTGAACATCCACt                            <  1:644915/140‑1 (MQ=255)                                                                                                         cTTATCCGGCCAATAAAATCCTAAAAATTAAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCGTGCCTCTGGCGGATGTACGTTTGTCATGAGTCTCACTCTGTTGCTAATTGCCGTTCGCTCCTGAACATCCACt                            <  2:75088/140‑1 (MQ=255)                                                                                                                ggCCAATAAAATCCTAAAAATTAAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCGTGCCTCTGGCGGATGTACGTTTGTCATGAGTCTCACTCTGTTGCTAATTGCCGTTCGCTCCTGAACATCCACTCGATCtt                     >  2:299464/1‑140 (MQ=255)                                                                                                                ggCCAATAAAATCCTAAAAATTAAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCGTGCCTCTGGCGGATGTACGTTTGTCATGAGTCTCACTCTGTTGCTAATTGCCGTTCGCTCCTGAACATCCACTCGATCtt                     >  1:378962/1‑140 (MQ=255)                                                                                                                  ccAATAAAATCCTAAAAATTAAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCGTGCCTCTGGCGGATGTACGTTTGTCATGAGTCTCACTCTGTTGCTAATTGCCGTTCGCTCCTGAACATCCACTCGATCTTCg                   >  2:401337/1‑140 (MQ=255)                                                                                                                  ccAATAAAATCCTAAAAATTAAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCGTGCCTCTGGCGGATGTACGTTTGTCATGAGTCTCACTCTGTTGCTAATTGCCGTTCGCTCCTGAACATCCACTCGATCTTCg                   >  1:698012/1‑140 (MQ=255)                                                                                                                   cAATAAAATCCTAAAAATTAAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCGTGCCTCTGGCGGATGTACGTTTGTCATGAGTCTCACTCTGTTGCTAATTGCCGTTCGCTCCTGAACATCCACTCGATCTTCGc                  >  2:328555/1‑140 (MQ=255)                                                                                                                           tCCTAAAAATTAAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCGTGCCTCTGGCGGATGTACGTTTGTCATGAGTCTCACTCTGTTGCTAATTGCCGTTCGCTCCTGAACATCCACTCGATCTTCGCCTTCTTcc          >  1:477318/1‑140 (MQ=255)                                                                                                                           tCCTAAAAATTAAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCGTGCCTCTGGCGGATGTACGTTTGTCATGAGTCTCACTCTGTTGCTAATTGCCGTTCGCTCCTGAACATCCACTCGATCTTCGCCTTCTTcc          >  1:875604/1‑140 (MQ=255)                                                                                                                             cTAAAAATTAAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCGTGCCTCTGGCGGATGTACGTTTGTCATGAGTCTCACTCTGTTGCTAATTGCCGTTCGCTCCTGAACATCCACTCGATCTTCGCCTTCTTCCgg        <  1:401337/140‑1 (MQ=255)                                                                                                                              tAAAAATTAAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCGTGCCTCTGGCGGATGTACGTTTGTCATGAGTCTCACTCTGTTGCTAATTGCCGTTCGCTCCTGAACATCCACTCGATCTTCGCCTTCTTCCGGt       <  1:664461/140‑1 (MQ=255)                                                                                                                                   aTTAAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCGTGCCTCTGGCGGATGTACGTTTGTCATGAGTCTCACTCTGTTGCTAATTGCCGTTCGCTCCTGAACATCCACTCGATCTTCGCCTTCTTCCGGTTTAtt  <  1:133839/140‑1 (MQ=255)                                                                                                                                      |                                                                                                                                         CGCTTACGAACTTATTACGCGCCTGACTTCAAGGGTCGCGATGAAATACGTGGATTAATCGTTCTGTAATATTTGATTGTCTGTGCCGGATGCGGCGTGAATGCCTTATCCGGCCAATAAAATCCTAAAAATTCAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCGTGCCTCTGGCGGATGTACGTTTGTCATGAGTCTCACTCTGTTGCTAATTGCCGTTCGCTCCTGAACATCCACTCGATCTTCGCCTTCTTCCGGTTTATT  >  NC_000913/4266804‑4267073

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