Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,702,016	T→C	Y354H (TAT→CAT)	malY →	negative regulator of MalT activity/cystathionine beta‑lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,702,016	0	T	C	100.0%	85.0 / NA	26	Y354H (TAT→CAT)	malY	negative regulator of MalT activity/cystathionine beta‑lyase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (15/11);  total (15/11)

TCCCACAATCCACTTATCTGGCATGGCTTGATTTACGTCCGTTGAATATTGACGACAACGCGTTGCAAAAAGCACTTATCGAACAAGAAAAAGTCGCGATCATGCCGGGGTATACCTACGGTGAAGAAGGTCGTGGTTTTGTCCGTCTCAATGCCGGCTGCCCACGTTCGAAACTGGAAAAAGGTGTGGCTGGATTAATTAACGCCATCCGCGCTGTTCGTTAAC  >  NC_000913/1701906‑1702130                                                                                                               |                                                                                                                   tCCCACAATCCACTTATCTGGCATGGCTTGATTTACGTCCGTTGAATATTGACGACAACGCGTTGCAAAAAGCACTTATCGAACAAGAAAAAGTCGCGATCATGCCGGGGCATACCTACGGTGAAGAAGGTCg                                                                                              <  1:278707/133‑1 (MQ=255)          ccACTTATCTGGCATGGCTTGATTTACGTCCGTTGAATATTGACGACAACGCGTTGCAAAAAGCACTTATCGAACAAGAAAAAGTCGCGATCATGCCGGGGCATACCTACGGTGAAGAAGGTCGTGGTTTTGt                                                                                     >  1:44461/1‑133 (MQ=255)           cACTTATCTGGCATGGCTTGATTTACGTCCGTTGAATATTGACGACAACGCGTTGCAAAAAGCACTTATCGAACAAGAAAAAGTCGCGATCATGCCGGGGCATACCTACGGTGAAGAAGGTCGTGGTTTTGTc                                                                                    <  2:65169/133‑1 (MQ=255)               tATCTGGCATGGCTTGATTTACGTCCGTTGAATATTGACGACAACGCGTTGCAAAAAGCACTTATCGAACAAGAAAAAGTCGCGATCATGCCGGGGCATACCTACGGTGAAGAAGGTCGTGGTTTTGTCCGtc                                                                                <  2:67100/133‑1 (MQ=255)                     gCATGGCTTGATTTACGTCCGTTGAATATTGACGACAACGCGTTGCAAAAAGCACTTATCGAACAAGAAAAAGTCGCGATCATGCCGGGGCATACCTACGGGAAAGAAGGGCGGTGTGTTGTCCCTCTCAATg                                                                          >  1:195748/1‑133 (MQ=255)                        tGGCTTGATTTACGTCCGTTGAATATTGACGACAACGCGTTGCAAAAAGCACTTATCGAACAAGAAAAAGTCGCGATCATGCCGGGGCATACCTACGGTGAAGAAGGTCGTGGTTTTGTCCGTCTCAATGCCg                                                                       >  1:65438/1‑133 (MQ=255)                           cTTGATTTACGTCCGTTGAATATTGACGACAACGCGTTGCAAAAAGCACTTATCGAACAAGAAAAAGTCGCGATCATGCCGGGGCATACCTACGGTGAAGAAGGTCGTGGTTTTGTCCGTCTCAATGCCGGCt                                                                    >  1:349549/1‑133 (MQ=255)                            ttGATTTACGTCCGTTGAATATTGACGACAACGCGTTGCAAAAAGCACTTATCGAACAAGAAAAAGTCGCGATCATGCCGGGGCATACCTACGGTGAAGAAGGTCGTGGTTTTGTCCGTCTCAATGCCGGCTg                                                                   >  1:96473/1‑133 (MQ=255)                               aTTTACGTCCGTTGACTATTGACGGCAAGGCGTTGCAAAAATCACTTATCGAACAAGAAAAAGTCGCGATCATGCCGGGGCATACCTACGGTGAAGAAGGTCGTGGTTTTGTCCGTCTCAATGCCGGCTGccc                                                                <  2:207708/133‑1 (MQ=255)                               aTTTACGTCCGTTGAATATTGACGACAACGCGTTGCAAAAAGCACTTATCGAACAAGAAAAAGTCGCGATCATGCCGGGGCATACCTACGGTGAAGAAGGTCGTGGTTTTGTCCGTCTCAATGCCGGCTGccc                                                                >  1:288245/1‑133 (MQ=255)                                     gTCCGTTGAATATTGACGACAACGCGTTGCAAAAAGCACTTATCGAACAAGAAAAAGTCGCGATCATGCCGGGGCATACCTACGGTGAAGAAGGTCGTGGGTTTTTCCTTCTCAATGGCGGCGGCCCACGTTc                                                          >  2:167606/1‑133 (MQ=255)                                       ccGTTGAATATTGACGACAACGCGTTGCAAAAAGCACTTATCGAACAAGAAAAAGTCGCGATCATGCCGGGGCATACCTACGGTGAAGAAGGTCGTGGTTTTGTCCGTCTCAATGCCGGCTGCCCACGTTCGa                                                        <  2:112117/133‑1 (MQ=255)                                         gTTGAATATTGACGACAACGCGTTGCAAAAAGCACTTATCGAACAAGAAAAAGTCGCGATCATGCCGGGGCATACCTACGGTGAAGAAGGTCGTGGTTTTGTCCGTCTCAATGCCGGCTGCCCACGTTCGaaa                                                      <  2:290721/133‑1 (MQ=255)                                                    acgacAACGCGTTGCAAAAAGCACTTATCGAACAAGAAAAAGTCGCGATCATGCCGGGGCATACCTACGGTGAAGAAGGTCGTGGTTTTGTTCGTCTCAAGGCCCGCTGCCCACGTTGGAAATTGGAAAAAgg                                           >  2:148021/1‑133 (MQ=255)                                                    acgacAACGCGTTGCAAAAAGCACTTATCGAACAAGAAAAAGTCGCGATCATGCCGGGGCATACCTACGGTGAAGAAGGTCGTGGTTTTGTCCGTCTCAATGCCGGCTGCCCACGTTCGAAACTGGAAAAAgg                                           >  1:400507/1‑133 (MQ=255)                                                       acAACGCGTTGCAAAAAGCACTTATCGAACAAGAAAAAGTCGCGATCATGCCGGGGCATACCTACGGTGAAGAAGGTCGTGGTTTTGTCCGTCTCAATGCCGGCTGCCCACGTTCGAAACTGGAAAAAGgtgt                                        >  1:300630/1‑133 (MQ=255)                                                         aaCGCGTTGCAAAAAGCACTTATCGAACAAGAAAAAGCCGCGATCATGCCGGGGAAGACCAACGGTAAAGAAGGTCGTGGTTTTGTCCGTCTCAATGCCGGCTGCCCACGTTCGAAGCTGGAAAAGGGTGTgg                                      <  2:365365/133‑1 (MQ=255)                                                           cgcgTTGCAAAAAGCACTTATCGAACAAGAAAAAGTCGCGATCATGCCGGGGCATACCTACGGTGAAGAAGGTCGTGGTTTTGTCCGTCTCGATGCCGGCGGCCGACGATCGAAACTGGAAAAAGGTGTGGCt                                    >  2:132403/1‑133 (MQ=255)                                                           cgcgTTGCAAAAAGCACTTATCGAACAAGAAAAAGTCGCGATCATGCCGGGGCATACCTACGGTGAAGAAGGTCGTGGTTTTGTCCGTCTCAATGCCGGCTGCCCACGTTCGAAACTGGAAAAAGGTGTGGCt                                    >  2:507304/1‑133 (MQ=255)                                                               ttGCAAAACGCACTTCTCTAACAAGAAAAAATCGCGATCATGCCGGGGCATACCTACGGTGAAGAAGGTCGTGGTTTTGTCCGTCTCAATGCCGGCTGCCCACGTTCGAAACTGGAAAAAGGTGTGGCTGGat                                <  1:148021/133‑1 (MQ=255)                                                                      aaGCACTTATCGAACAAGAAAAAGTCGCGATCATGCCGGGGCATACCTACGGTGAAGAAGGTCGTGGTTTTGTCCGTCTCAATGCCGGCTGCCCACGTTCGAAACTGGAAAAAGGTGTGGCTGGAttaattaa                         <  2:65438/133‑1 (MQ=255)                                                                            ttATCGAACAAGAAAAAGTCGCGATCATGCCGGGGCATACCTACGGTGAAGAAGGTCGTGGTTTTGTCCGTCTCAATGCCGGCTGCCCACGTTCGAAACTGGAAAAAGGTGTGGCTGGATTAATTAACGCCAt                   >  1:446533/1‑133 (MQ=255)                                                                            ttATCGAACAAGAAAAAGTCGCGATCATGCCGGGGCATACCTACGGTGAAGAAGGTCGTGGTTTTGTCCGTCTCAATGCCGGCTGCCCACGTTCGAAACTGGAAAAAGGTGTGGCTGGATTAATTAACGCCAt                   >  2:366413/1‑133 (MQ=255)                                                                              aTCGAACAAGAAAAAGTCGCGATCATGCCGGGGCATACCTACGGTGAAGAAGGCCGTGGTTTTGTCCGTCTCAATGCCGGCTGCCCACGTTCGAAACTGGAAAAAGGTGTGGCTGGATTAATTAACGCCATcc                 <  2:300630/133‑1 (MQ=255)                                                                               tCGAACAAGAAAAAGTCGCGATCATGCCGGGGCATACCTACGGTGAAGAAGGTCGTGGTTTTGTCCGTCTCAATGCCGGCTGCCCACGTTCGAAACTGGAAAAAGGTGTGGCTGGATTAATTAACGCCATCcg                <  1:357578/133‑1 (MQ=255)                                                                                    cAAGAAAAAGTCGCGATCATGCCGGGGCATACCTACGGTGAAGAAGGTCGTGGTTTTGTCCGTCTCAATGCCGGCTGCCCACGTTCGAAACTGGAAAAAGGTGTGGCTGGATTAATTAACGCCATCCGCGCTg           >  1:24480/1‑133 (MQ=255)                                                                                             gTCGCGATCATGCCGGGGCATACCTACGGTGAAGAAGGTCGTGGTTTTGTCCGTCTCAATGCCGGCTGCCCACGTTCGAAACTGGAAAAAGGTGTGGCTGGATTAATTAACGCCATCCGCGCTGTTCGTTAAc  <  2:24480/133‑1 (MQ=255)                                                                                                               |                                                                                                                   TCCCACAATCCACTTATCTGGCATGGCTTGATTTACGTCCGTTGAATATTGACGACAACGCGTTGCAAAAAGCACTTATCGAACAAGAAAAAGTCGCGATCATGCCGGGGTATACCTACGGTGAAGAAGGTCGTGGTTTTGTCCGTCTCAATGCCGGCTGCCCACGTTCGAAACTGGAAAAAGGTGTGGCTGGATTAATTAACGCCATCCGCGCTGTTCGTTAAC  >  NC_000913/1701906‑1702130

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: ‑