Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,975,554	A→G	T137A (ACC→GCC)	wecE →	dTDP‑4‑dehydro‑6‑deoxy‑D‑glucose transaminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,975,554	0	A	G	100.0%	149.7 / NA	43	T137A (ACC→GCC)	wecE	dTDP‑4‑dehydro‑6‑deoxy‑D‑glucose transaminase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (19/24);  total (19/24)

GTGGCGCAAAAATCGTTTTTGTGGATGTTCGCCCGGACACCATGAACATCGACGAAACGCTGATTGAAGCGGCGATCACCGACAAAACGCGCGTTATCGTGCCGGTCCATTACGCGGGTGTGGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCCGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAA  >  NC_000913/3975417‑3975688                                                                                                                                          |                                                                                                                                       gTGGCGCAAAAATCGTTTTTGTGGATGTTCGCCCGGACACCATGAACATCGACGAAACGCTGATTGAAGCGGCGATCACCGACAAAACGCGCGTTATCGTGCCGGTCCATTACGCGGGTGTGGCCTGCGAAATGGACGcc                                                                                                                                      >  2:485575/1‑140 (MQ=255)  tGGCGCAAAAATCGTTTTTGTGGATGTTCGCCCGGACACCATGAACATCGACGAAACGCTGATTGAAGCGGCGATCACCGACAAAACGCGCGTTATCGTGCCGGTCCATTACGCGGGTGTGGCCTGCGAAATGGACGCCa                                                                                                                                     >  2:755991/1‑140 (MQ=255)  tGGCGCAAAAATCGTTTTTGTGGATGTTCGCCCGGACACCATGAACATCGACGAAACGCTGATTGAAGCGGCGATCACCGACAAAACGCGCGTTATCGTGCCGGTCCATTACGCGGGTGTGGCCTGCGAAATGGACGCCa                                                                                                                                     >  2:61067/1‑140 (MQ=255)     cgcAAAAATCGTTTTTGTGGATGTTCGCCCGGACACCATGAACATCGACGAAACGCTGATTGAAGCGGCGATCACCGACAAAACGCGCGTTATCGTGCCGGTCCATTACGCGGGTGTGGCCTGCGAAATGGACGCCATTa                                                                                                                                  >  1:281882/1‑140 (MQ=255)          aaaTCGTTTTTGTGGATGTTCGCCCGGACACCATGAACATCGACGAAACGCTGATTGAAGCGGCGATCACCGACAAAACGCGCGTTATCGTGCCGGTCCATTACGCGGGTGTGGCCTGCGAAATGGACGCCATTATGGCg                                                                                                                             >  2:208351/1‑140 (MQ=255)            aTCGTTTTTGTGGATGTTCGCCCGGACACCATGAACATCGACGAAACGCTGATTGAAGCGGCGATCACCGACAAAACGCGCGTTATCGTGCCGGTCCATTACGCGGGTGTGGCCTGCGAAATGGACGCCATTATGGCGtt                                                                                                                           >  2:334574/1‑140 (MQ=255)                  tttGTGGATGTTCGCCCGGACACCATGAACATCGACGAAACGCTGATTGAAGCGGCGATCACCGACAAAACGCGCGTTATCGTGCCGGTCCATTACGCGGGTGTGGCCTGCGAAATGGACGCCATTATGGCGTTGGCGaa                                                                                                                     >  2:659097/1‑140 (MQ=255)                    tgtgGATGTTCGCCCGGACACCATGAACATCGACGAAACGCTGATTGAAGCGGCGATCACCGACAAAACGCGCGTTATCGTGCCGGTCCATTACGCGGGTGTGGCCTGCGAAATGGACGCCATTATGGCGTTGGCGaaaa                                                                                                                   <  2:281882/140‑1 (MQ=255)                            ttCGCCCGGACACCATGAACATCGACGAAACGCTGATTGAAGCGGCGATCACCGACAAAACGCGCGTTATCGTGCCGGTCCATTACGCGGGTGTGGCCTGCGAAATGGACGCCATTATGGCGTTGGCGAAAAAGCATAAt                                                                                                           <  1:831275/140‑1 (MQ=255)                                cccGGACACCATGAACATCGACGAAACGCTGATTGAAGCGGCGATCACCGACAAAACGCGCGTTATCGTGCCGGTCCATTACGCGGGTGTGGCCTGCGAAATGGACGCCATTATGGCGTTGGCGAAAAAGCATAAtttgt                                                                                                       <  1:63723/140‑1 (MQ=255)                                cccGGACACCATGAACATCGACGAAACGCTGATTGAAGCGGCGATCACCGACAAAACGCGCGTTATCGTGCCGGTCCATTACGCGGGTGTGGCCTGCGAAATGGACGCCATTATGGCGTTGGCGAAAAAGCATAAtttgt                                                                                                       <  2:839305/140‑1 (MQ=255)                                 ccGGACACCATGAACATCGACGAAACGCTGATTGAAGCGGCGATCACCGACAAAACGCGCGTTATCGTGCCGGTACATTACGCGGGTGTGGCCTGCGAAATGGACGCCATTATGGCGTTGGCGAAAAAGCATAAtttgtt                                                                                                      <  1:860574/140‑1 (MQ=255)                                  cGGACACCATGAACATCGACGAAACGCTGATTGAAGCGGCGATCACCGACAAAACGCGCGTTATCGTGCCGGTCCATTACGCGGGTGTGGCCTGCGAAATGGACGCCATTATGGCGTTGGCGAAAAAGCATAAtttgttt                                                                                                     >  1:446133/1‑140 (MQ=255)                                   ggACACCATGAACATCGACGAAACGCTGATTGAAGCGGCGATCACCGACAAAACGCGCGTTATCGTGCCGGTCCATTACGCGGGTGTGGCCTGCGAAATGGACGCCATTATGGCGTTGGCGAAAAAGCATAAtttgtttg                                                                                                    <  1:659097/140‑1 (MQ=255)                                          aTGAACATCGACGAAACGCTGATTGAAGCGGCGATCACCGACAAAACGCGCGTTATCGTGCCGGTCCATTACGCGGGTGTGGCCTGCGAAATGGACGCCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTaga                                                                                             <  1:587634/140‑1 (MQ=255)                                           tGAACATCGACGAAACGCTGATTGAAGCGGCGATCACCGACAAAACGCGCGTTATCGTGCCGGTCCATTACGCGGGTGTGGCCTGCGAAATGGACGCCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTagaa                                                                                            >  2:400554/1‑140 (MQ=255)                                               cATCGACGAAACGCTGATTGAAGCGGCAATCACCGACAAAACGCGCGTTATCGTGCCGGTCCATTACGCGGGTGTGGCCTGCGAAATGGACGCCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATg                                                                                        >  1:883156/1‑140 (MQ=255)                                                   gacgaAACGCTGATTGAAGCGGCGATCACCGACAAAACGCGCGTTATCGTGCCGGTCCATTACGCGGGTGTGGCCTGCGAAATGGACGCCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTaga                                                                                             <  2:340425/131‑1 (MQ=255)                                                   gacgaAACGCTGATTGAAGCGGCGATCACCGACAAAACGCGCGTTATCGTGCCGGTCCATTACGCGGGTGTGGCCTGCGAAATGGACGCCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTaga                                                                                             >  1:340425/1‑131 (MQ=255)                                                   gacgaAACGCTGATTGAAGCGGCGATCACCGACAAAACGCGCGTTATCGTGCCGGTCCATTACGCGGGTGTGGCCTGCGAAATGGACGCCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCCGc                                                                                    >  2:47773/1‑140 (MQ=255)                                                          cGCTGATTGAAGCGGCGATCACCGACAAAACGCGCGTTATCGTGCCGGTCCATTACGCGGGTGTGGCCTGCGAAATGGACGCCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCCGCTCAGGGc                                                                             <  1:61067/140‑1 (MQ=255)                                                          cGCTGATTGAAGCGGCGATCACCGACAAAACGCGCGTTATCGTGCCGGTCCATTACGCGGGTGTGGCCTGCGAAATGGACGCCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCCGCTCAGGGc                                                                             <  1:755991/140‑1 (MQ=255)                                                                 tGAAGCGGCGATCACCGACAAAACGCGCGTTATCGTGCCGGTCCATTACGC‑GGTGTGGCCTGCGAAATGGACGCCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCCGCTCAGGGCGTGATGTc                                                                     >  1:559366/1‑140 (MQ=255)                                                                    aGCGGCGATCACCGACAAAACGCGCGTTATCGTGCCGGTCCATTACGCGGGTGTGGCCTGCGAAATGGACGCCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCCGCTCAGGGCGTGATGTCCa                                                                   <  2:517332/140‑1 (MQ=255)                                                                    aGCGGCGATCACCGACAAAACGCGCGTTATCGTGCCGGTCCATTACGCGGGTGTGGCCTGCGAAATGGACGCCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCCGCTCAGGGCGTGATGTCCa                                                                   <  1:208351/140‑1 (MQ=255)                                                                           aTCACCGACAAAACGCGCGTTATCGTGCCGGTCCATTACGCGGGTGTGGCCTGCGAAATGGACGCCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCCGCTCAGGGCGTGATGTCCACTTACaa                                                            >  1:122625/1‑140 (MQ=255)                                                                           aTCACCGACAAAACGCGCGTTATCGTGCCGGTCCATTACGCGGGTGTGGCCTGCGAAATGGACGCCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCCGCTCAGGGCGTGATGTCCACTTACaa                                                            <  1:389160/140‑1 (MQ=255)                                                                                   cAAAACGCGCGTTATCGTGCCGGTCCATTACGCGGGTGTGGCCTGCGAAATGGACGCCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCCGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTg                                                    <  1:415555/140‑1 (MQ=255)                                                                                     aaaCGCGCGTTATCGTGCCGGTCCATTACGCGGGTGTGGCCTGCGAAATGGACGCCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCCGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCa                                                  <  1:485575/140‑1 (MQ=255)                                                                                        cgcgcgTTATCGTGCCGGTCCATTACGCGGGTGTGGCCTGCGAAATGGACGCCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCCGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTg                                               <  1:47773/140‑1 (MQ=255)                                                                                          cgcgTTATCGTGCCGGTCCATTACGCGGGTGTGGCCTGCGAAATGGACGCCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCCGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTggg                                             <  1:286841/140‑1 (MQ=255)                                                                                          cgcgTTATCGTGCCGGTCCATTACGCGGGTGTGGCCTGCGAAATGGACGCCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCCGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTggg                                             <  1:77285/140‑1 (MQ=255)                                                                                            cgTTATCGTGCCGGTCCATTACGCGGGTGTGGCCTGCGAAATGGACGCCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCCGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGaa                                           >  1:286388/1‑140 (MQ=255)                                                                                               tATCGTGCCGGTCCATTACGCGGGTGTGGCCTGCGAAATGGACGCCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCCGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCa                                        <  1:645008/140‑1 (MQ=255)                                                                                                   gTGCCGGTCCATTACGC‑GGTGTGGCCTGCGAAATGGACGCCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCCGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGt                                   <  2:559366/140‑1 (MQ=255)                                                                                                    tGCCGGTCCATTACGCGGGTGTGGCCTGCGAAATGGACGCCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCCGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGt                                   <  1:400554/140‑1 (MQ=255)                                                                                                               tACGCGGGTGTGGCCTGCGAAATGGACGCCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCCGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGgctg                        >  2:359649/1‑140 (MQ=255)                                                                                                                      gtgtGGCCTGCGAAATGGACGCCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCCGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGc                 <  2:539703/140‑1 (MQ=255)                                                                                                                       tgtgGCCTGCGAAATGGACGCCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCCGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCt                <  1:763025/140‑1 (MQ=255)                                                                                                                           gCCTGCGAAATGGACGCCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCCGCCCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTg                            <  1:603455/124‑1 (MQ=255)                                                                                                                           gCCTGCGAAATGGACGCCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCCGCCCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTg                            >  2:603455/1‑124 (MQ=255)                                                                                                                                     tGGACGCCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCCGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCaaa  >  2:400922/1‑140 (MQ=255)                                                                                                                                     tGGACGCCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCCGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCaaa  >  1:978885/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                                       GTGGCGCAAAAATCGTTTTTGTGGATGTTCGCCCGGACACCATGAACATCGACGAAACGCTGATTGAAGCGGCGATCACCGACAAAACGCGCGTTATCGTGCCGGTCCATTACGCGGGTGTGGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCCGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAA  >  NC_000913/3975417‑3975688

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