Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,701,272	C→T	Q106* (CAG→TAG)	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,701,272	0	C	T	100.0%	74.9 / NA	23	Q106* (CAG→TAG)	malY	negative regulator of MalT activity/cystathionine beta‑lyase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (15/8);  total (15/8)

CACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCCTTCTGTCATCTATATGGTTTCAGAACTGATTCGTCAGTGGTCTGAAACAGGTGAAGGCGTGGTGATCCACACACCCGCCTATGACGCATTTTACAAGGCCATTGAAGGTAACCAGCGCACAGTAATGCCCGTTGCTTTAGAGAAGCAGGCTGATGGTTGGTTTTGCGATATGGGCAAGTTGGAAGCCGTGTTGGCGAAACCAGAATGTAAAATTATGCTCCTGTGTAGCCCACAGAATCCTACCGGGAAAGTGTGGAC  >  NC_000913/1701034‑1701495                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                cACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCCTTCTGTCATCTATATGGTTTCAGAACTGATTCGTTAGTGGTCTGaa                                                                                                                                                                                                                      >  1:306420/1‑250 (MQ=255)    tgctgACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCCTTCTGTCATCTATATGGTTTCAGAACTGATTCGTTAGTGGTCTGAAAc                                                                                                                                                                                                                    >  1:215008/1‑249 (MQ=255)     gctgACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCCTTCTGTCATCTATATGGTTTCAGAACTGATTCGTTAGTGGTCTGAAACa                                                                                                                                                                                                                   <  2:330351/249‑1 (MQ=255)         aCCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCCTTCTGTCATCTATATGGTTTCAGAACTGATTCGTTAGTGTTCTGAAACAGGTGaa                                                                                                                                                                                                             <  1:319482/251‑1 (MQ=255)            tGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGATGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCCTTCTGTCATCTATATGGTTTCAGAACTGATTCGTTAGTGGTCTGAAACAGGTGAAGGc                                                                                                                                                                                                          >  1:405645/1‑251 (MQ=255)              ttACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGCGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAACCGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCCTTCTGTCATCTATATGGTTTCCGAACTGATTCGTTAGTGGTCTGAAACAGGTGAAGGCgg                                                                                                                                                                                                        >  2:383321/1‑250 (MQ=255)                                                  gCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCCTTCTGTCATCTATATGGTTTCAGAACTGATTCGTTAGTGGTCTGAAACAGGTGAAGGCGTGGTGATCCACACACCCGCCTATGACGCATTTTACa                                                                                                                                                                     >  1:64939/1‑250 (MQ=255)                                                            ttATCGGGGCGTTGAATCAGCGCCTGATGCCCGGCGTATTGGGCTACAGCCGCTGAAAAAACGATGATTTTCTCGCGCCTATGCCCCACTGTTTTCCCCCCCAGCATTCCCCCCCCATCGATTCTAAGACGGTGGTGTATGGCCCTTCTGTCATCTATAGGTTTTCAGAACTGATTCGTTAGTGGTCTGAAACAGGTGAGGGCGTGGTGACCCACACCCCCGCCTATGACGCATTTTACAAGGCCATTGa                                                                                                                                                           <  1:259190/250‑1 (MQ=255)                                                                ggAGGCGCTGAACCAGCGCCTGATGCGCGGCGTTTTTGGCTACAGCCGCTGAAAAAACGATGAGTTTCCCGCGGCTATTGCCCTCTGTTTTTCCCCCCAGCATTACCCCGCCATCGTTTCTCAGACGGTGGTGTATGGCCCTTCTGTCTTCTATATGTTTTAAGAACTGATTCGTTAGTGGTCTAAAACAGGTGAAGGCGTGGTGACCCACACACCCGCCTATGACGCATTTTACAAGGCCATTGAAGGTa                                                                                                                                                      <  2:64939/250‑1 (MQ=255)                                                                                        gccggggtgattgggCTACAGCCGTTGGAAAAACGATGAGTTTCCCGGGGTTATTGCCCGCTGTTTTTCCCCCCAGCATTCCACCGCCATGGATTCTCAGACGGTGGTGTAGGGCCCTTCTTTCATCTATAGGGTTTCAGAACTGATTCGTTAGTGGTCGGAAACGGGGGAAGGCGTGGTGACCCCCACACCCGCCTATGACGCATTTTACAAGGCCATTGAAGGTAACCAGCGCACAGTAAGGCCCGTTg                                                                                                                              <  2:390468/238‑1 (MQ=255)                                                                                                  tttGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCCTTCTGTCATCTATATGGTTTCAGAACTGATTCGTTAGTGGTCTGAAACAGGTGAAGGCGTGGTGATCCACACACCCGCCTATGACGCATTTTACAAGGCCATTGAAGGTAACCAGCGCACAGTAATGCCCGTTGCTTTAgaga                                                                                                                     >  1:383475/1‑250 (MQ=255)                                                                                                   ttGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCCTTCTGTCATCTATATGGTTTCAGAACTGATTCGTTAGTGGTCTGAAACAGGTGAAGGCGTGGTGATCCACACACCCGCCTATGACGCATTTTACAAGGCCATTGAAGGTAACCAGCGCACAGTAATGCCCGTTGCTTTAgaga                                                                                                                     >  2:240945/1‑249 (MQ=255)                                                                                                                          gATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCCTTCTGTCATCTATATGGTTTCAGAACTGATTCGTTAGTGGTCTGAAACAGGTGAAGGCGTGGTGATCCACACACCCGCCTATGACGCATTTTACAAGGCCATTGAAGGTAACCAGCGCACAGTAATGCCCGTTGCTTTAGAGAAGCAGGCTGATGGTTGGTTTTGCg                                                                                             >  2:76272/1‑250 (MQ=255)                                                                                                                                             ttGCCCACTGGTTTTCCACACAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCCTTCTGTCATCTATATGGTTTCAGAACTGATTCGTTAGTGGTCTGAAACAGGTGAAGGCGTGGTGATCCACACACCCGCCTATGACGCATTTTACAAGGCCATTGAAGGTAACCAGCGCACAGTAATGCCCGTTGCTTTAGAGAAGCAGGCTGATGGTTGGTTTTGCGATATGGGCAAGTTGGAAGcc                                                                         >  1:117272/1‑251 (MQ=255)                                                                                                                                                   aCTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCCTTCTGTCATCTATATGGTTTCAGAACTGATTCGTTAGTGGTCTGAAACAGGTGAAGGCGTGGTGATCCACACACCCGCCTATGACGCATTTTACAAGGCCATTGAAGGTAACCAGCGCACAGTAATGCCCGTTGCTTTAGAGAAGCAGGCTGATGGTTGGTTTTGCGATATGGGCAAGTTGGAAGCCGTGTTg                                                                   >  2:184775/1‑251 (MQ=255)                                                                                                                                                        ttttCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCCTTCTGTCATCTATATGGTTTCAGAACTGATTCGTTAGTGGTCTGAAACAGGTGAAGGCGTGGTGATCCACACACCCGCCTATGACGCATTTTACAAGGCCATTGAAGGTAACCAGCGCACAGTAATGCCCGTTGCTTTAGAGAAGCAGGCTGATGGTTGGTTTTGCGATATGGGCAAGTTGGAAGCCGTGTTGGCGa                                                               >  1:262540/1‑250 (MQ=255)                                                                                                                                                            ccACCCAGCATTCCCCCGCCATCGTTTCTCAGACGGTGTTGTATGCCCTTTCTTTCATCTATATGTTTTCAAACCTGTTTCGTTAGTGGTCTGAAACAGGTGAAGGCGTGGTGTTCCACCCCCCCGCCTATGACGCATTTTACAAGGCCATTGAAGGTAACCAGCGCACAGTAATGCCCGTTGCTTTAGAGAAGCAGGCTGATGGTTGGTTTTGCGATATGGGCAAGTTGGAAGCCGTGTTGGCGAAACCa                                                          <  1:302917/251‑1 (MQ=255)                                                                                                                                                             cccccAGCATTCCCCCGCCATCGATTCTCAGACGGTGGTGTAGGGCCCTTCTGTCATCTATATGGTTTCAGAACTGATTCGTTAGTGGTCTGAAACAGGTGAAGGCGGGGTGACCCACCCCCCCGCCTATGACGCATTTTACAAGGCCATTGAAGGTAACCAgc                                                                                                                                                <  1:404114/162‑1 (MQ=255)                                                                                                                                                             cACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCCTTCTGTCATCTATATGGTTTCAGAACTGATTCGTTAGTGGTCTGAAACAGGTGAAGGCGTGGTGATCCACACACCCGCCTATGACGCATTTTACAAGGCCATTGAAGGTAACCAgc                                                                                                                                                >  2:404114/1‑164 (MQ=255)                                                                                                                                                                             gCCATCGATTCTCAGACGGTGGTGTATGGCCCTTCTGTCATCTATATGGTTTCAGAACTGATTCGTTAGTGGTCTGAAACAGGTGAAGGCGTGGTGATCCACACACCCGCCTATGACGCATTTTACAAGGCCATTGAAGGTAACCAGCGCACAGTAATGCCCGTTGCTTTAGAGAAGCAGGCTGATGGTTGGTTTTGCGATATGGGCAAGTTGGAAGCCGTGTTGGCGAAACCAGAATGTAAAATTATGCt                                         >  1:88556/1‑251 (MQ=255)                                                                                                                                                                                                  gtgtATGGCCCGTCTGTCATCTATATGGTTCAAGAACTGATTCGTTAGGGGTCTGAAACAGGTGAAGGGGTGGTGATCCACACACCCGCCTATGACGCATTTTACAAGGCCATTGAAGGTAACCAGCGCACAGTAATGCCCGTTGCTTTAGAGAAGCAGGCTGATGGTTGGTTTTGCGATATGGGCAAGTTGGAAGCCGTGTTGGCGAAACCAGAATGTAAAATTATGCTCCTGTGTAGCCCACAGAATcc                    <  1:205414/251‑1 (MQ=255)                                                                                                                                                                                                        ggCCCTTCTGTCATCTATATGGTTTCAGAACTGATTCGTTAGTGGTCTGAAACAGGTGAAGGCGTGGTGATCCACACACCCGCCTATGACGCATTTTACAAGGCCATTGAAGGTAACCAGCGCACAGTAATGCCCGTTGCTTTAGAGAAGCAGGCTGATGGTTGGTTTTGCGATATGGGCAAGTTGGAAGCCGTGTTGGCGAAACCAGAATGTAAAATTATGCTCCTGTGTAGCCCACAGAATCCTACCgg              >  2:299306/1‑251 (MQ=255)                                                                                                                                                                                                                    aTCTATATGGTTTCAGAACTGATTCGTTAGTGGTCTGAAACAGGTGAAGGCGTGGTGATCCACACACCCGCCTATGACGCATTTTACAAGGCCATTGAAGGTAACCAGCGCACAGTAATGCCCGTTGCTTTAGAGAAGCAGGCTGATGGTTGGTTTTGCGATATGGGCAAGTTGGAAGCCGTGTTGGCGAAACCAGAATGTAAAATTATGCTCCTGTGTAGCCCACAGAATCCTACCGGGAAAGTGTGGAc  >  2:182041/1‑251 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                CACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCCTTCTGTCATCTATATGGTTTCAGAACTGATTCGTCAGTGGTCTGAAACAGGTGAAGGCGTGGTGATCCACACACCCGCCTATGACGCATTTTACAAGGCCATTGAAGGTAACCAGCGCACAGTAATGCCCGTTGCTTTAGAGAAGCAGGCTGATGGTTGGTTTTGCGATATGGGCAAGTTGGAAGCCGTGTTGGCGAAACCAGAATGTAAAATTATGCTCCTGTGTAGCCCACAGAATCCTACCGGGAAAGTGTGGAC  >  NC_000913/1701034‑1701495

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