Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,701,088	(C)5→4	coding (132/1173 nt)	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,084	0	C	.	100.0%	187.2 / NA	45	coding (128/1173 nt)	malY	negative regulator of MalT activity/cystathionine beta‑lyase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (27/18);  total (27/18)

AGGATAAGATATGTTCGATTTTTCAAAGGTCGTGGATCGTCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTC  >  NC_000913/1700947‑1701216                                                                                                                                          |                                                                                                                                     aGGATAAGATATGTTCGATTTTTCAAAGGTCGTGGATCGTCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGccc                                                                                                                                    <  2:665983/140‑1 (MQ=255) aGGATAAGATATGTTCGATTTTTCAAAGGTCGTGGATCGTCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGccc                                                                                                                                    >  2:190291/1‑140 (MQ=255) aGGATAAGATATGTTCGATTTTTCAAAGGTCGTGGATCGTCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGccc                                                                                                                                    >  2:172103/1‑140 (MQ=255)            tGTTCGATTTTTCAAAGGTCGTGGATCGTCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTG‑CCCCTGCATTATCg                                                                                                                        >  2:497668/1‑140 (MQ=255)             gTTCGATTTTTCAAAGGTCGTGGATCGTCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTG‑CCCCTGCATTATCGa                                                                                                                       <  2:141184/140‑1 (MQ=255)              ttCGATTTTTCAAAGGTCGTGGATCGTCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTG‑CCCCTGCATTATCGAg                                                                                                                      >  1:846651/1‑140 (MQ=255)               tCGATTTTTCAAAGGTCGTGGATCGTCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTG‑CCCCTGCATTATCGAgg                                                                                                                     <  1:946444/140‑1 (MQ=255)                 gATTTTTCAAAGGTCGTGGATCGTCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTG‑CCCCTGCATTATCGAGgcg                                                                                                                   <  1:312927/140‑1 (MQ=255)                  aTTTTTCAAAGGTCGTGGATCGTCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTG‑CCCCTGCATTATCGAGgcgc                                                                                                                  >  1:431432/1‑140 (MQ=255)                       tCAAAGGTCGTGGATCGTCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTG‑CCCCTGCATTATCGAGGCGCTGAAt                                                                                                             <  2:522317/140‑1 (MQ=255)                           aGGTCGTGGATCGTCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTG‑CCCCTGCATTATCGAGGCGCTGAATCAgc                                                                                                         <  1:190291/140‑1 (MQ=255)                             gTCGTGGATCGTCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTG‑CCCCTGCAATATCGAGGCGCTGAATCAgcgc                                                                                                       <  1:172103/140‑1 (MQ=255)                                gTGGATCGTCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTG‑CCCCTGCATTATCGAGGCGCTGAATCAGCGCCTg                                                                                                    >  2:83982/1‑140 (MQ=255)                                 tGGATCGTCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTG‑CCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGa                                                                                                   >  1:898342/1‑140 (MQ=255)                                  ggATCGTCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTG‑CCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGAt                                                                                                  >  2:1111984/1‑140 (MQ=255)                                    aTCGTCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCTCTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTG‑CCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGc                                                                                                <  2:712822/140‑1 (MQ=255)                                     tCGTCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTG‑CCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCa                                                                                               <  1:192436/140‑1 (MQ=255)                                       gTCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTG‑CCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACg                                                                                             <  1:497668/140‑1 (MQ=255)                                                   ggTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTG‑CCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCt                                                                                 >  1:187815/1‑140 (MQ=255)                                                    gtgtACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTG‑CCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTa                                                                                <  1:1111984/140‑1 (MQ=255)                                                       tACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTG‑CCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAg                                                                             >  2:106877/1‑140 (MQ=255)                                                              tGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTG‑CCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTgg                                                                      >  2:130089/1‑140 (MQ=255)                                                              tGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTG‑CCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTgg                                                                      >  2:667747/1‑140 (MQ=255)                                                              tGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTG‑CCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTgg                                                                      >  2:666765/1‑140 (MQ=255)                                                              tGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTG‑CCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTgg                                                                      >  1:28663/1‑140 (MQ=255)                                                              tGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTG‑CCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTgg                                                                      >  1:740105/1‑140 (MQ=255)                                                              tGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTG‑CCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTgg                                                                      >  1:768005/1‑140 (MQ=255)                                                                        tCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTG‑CCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATg                                                            >  2:1134271/1‑140 (MQ=255)                                                                                    tCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTG‑CCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCgg                                                <  2:740105/140‑1 (MQ=255)                                                                                        cACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTG‑CCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTAt                                            >  1:731159/1‑140 (MQ=255)                                                                                             ctgACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTG‑CCCCTGCATTATCGAGGCGCTGAATCAGCGCCTTATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGccc                                       >  2:1126683/1‑140 (MQ=255)                                                                                             ctgACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTG‑CCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGccc                                       >  1:343816/1‑140 (MQ=255)                                                                                             ctgACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTG‑CCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGccc                                       >  2:186790/1‑140 (MQ=255)                                                                                             ctgACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTG‑CCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGccc                                       >  2:389515/1‑140 (MQ=255)                                                                                                 ccTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTG‑CCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTg                                   <  1:666765/140‑1 (MQ=255)                                                                                                     ttACCGTTCACGATTTCAGACATGGATTTTGCCACTG‑CCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGcc                                        >  1:384341/1‑131 (MQ=255)                                                                                                     ttACCGTTCACGATTTCAGACATGGATTTTGCCACTG‑CCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGcc                                        <  2:384341/131‑1 (MQ=255)                                                                                                     ttACCGTTCACGATTTCAGACATGGATTTTGCCACTG‑CCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGttt                               >  1:610297/1‑140 (MQ=255)                                                                                                     ttACCGTTCACGATTTCAGACATGGATTTTGCCACTG‑CCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGttt                               >  2:337622/1‑140 (MQ=255)                                                                                                       aCCGTTCACGATTTCAGACATGGATTTTGCCACTG‑CCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTc                             <  1:389515/140‑1 (MQ=255)                                                                                                           ttCACGATTTCAGACATGGATTTTGCCACTG‑CCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCAcc                         >  2:880563/1‑140 (MQ=255)                                                                                                             cACGATTTCAGACATGGATTTTGCCACTG‑CCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGtttt                              <  1:727844/133‑1 (MQ=255)                                                                                                             cACGATTTCAGACATGGATTTTGCCACTG‑CCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGtttt                              >  2:727844/1‑133 (MQ=255)                                                                                                                  tttCAGACATGGATTTTGCCACTG‑CCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCAtt                  >  1:945616/1‑140 (MQ=255)                                                                                                                       gACATGGATTTTGCCACTG‑CCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACAcc             <  1:880563/140‑1 (MQ=255)                                                                                                                         cATGGATTTTGCCACTG‑CCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGc           <  2:945616/140‑1 (MQ=255)                                                                                                                            ggATTTTGCCACTG‑CCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCAt        <  1:878719/140‑1 (MQ=255)                                                                                                                                  tGCCACTG‑CCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATtc  >  1:368325/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                                     AGGATAAGATATGTTCGATTTTTCAAAGGTCGTGGATCGTCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTC  >  NC_000913/1700947‑1701216

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