Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,397,867	C→T	Y152Y (TAC→TAT)	mutL →	DNA mismatch repair protein MutL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,397,867	0	C	T	100.0%	193.9 / NA	54	Y152Y (TAC→TAT)	mutL	DNA mismatch repair protein MutL
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (29/25);  total (29/25)

CCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTACAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCG  >  NC_000913/4397732‑4398001                                                                                                                                        |                                                                                                                                       ccTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTataaca                                                                                                                                    <  1:327819/140‑1 (MQ=255)           aCTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTATAACACCCcggcgcg                                                                                                                          >  2:1049021/1‑140 (MQ=255)                aCGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTATAACACCCCGGCGCGGCGCa                                                                                                                     >  1:743569/1‑140 (MQ=255)                     ccGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTATAACACCCCGGCGCGGCGCAAATTc                                                                                                                >  1:531364/1‑140 (MQ=255)                       gCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCCGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTATAACACCCCGGCGCGGCACAAATTCCt                                                                                                              >  1:860029/1‑140 (MQ=255)                            aCAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTATAACACCCCGGCGCGGCGCAAATTCCTGCGCa                                                                                                         <  2:809298/140‑1 (MQ=255)                                cagGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTATAACACCCCGGCGCGGCGCAAATTCCTGCGCACCga                                                                                                     >  2:4541/1‑140 (MQ=255)                                               gCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTATAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTaa                                                                                      <  1:1049021/140‑1 (MQ=255)                                               gCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTATAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTaa                                                                                      <  2:93531/140‑1 (MQ=255)                                               gCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTATAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTaa                                                                                      >  2:754055/1‑140 (MQ=255)                                               gCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTATAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTaa                                                                                      >  1:358657/1‑140 (MQ=255)                                                 ttATGCGGAAGGCGCCGATA‑GAACCTTACCGGAAAACCGGCCGGGCCATCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTATAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACca                                                                                   <  1:946353/139‑1 (MQ=255)                                                 cTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTATAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCg                                                                                             >  1:656797/1‑131 (MQ=255)                                                 cTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTATAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCg                                                                                             <  2:656797/131‑1 (MQ=255)                                                   aTGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTATAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACcac                                                                                  <  1:8475/140‑1 (MQ=255)                                                    tgtcGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTATAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACca                                                                                   >  1:880413/4‑138 (MQ=255)                                                    tgtcGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTATAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACca                                                                                   <  2:880413/135‑1 (MQ=255)                                                    tGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTATAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACcaca                                                                                 <  2:370577/140‑1 (MQ=255)                                                    tGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTATAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACcaca                                                                                 <  2:496149/140‑1 (MQ=255)                                                    tGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTATAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACcaca                                                                                 >  1:1107217/1‑140 (MQ=255)                                                    tGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTATAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACcaca                                                                                 <  2:743569/140‑1 (MQ=255)                                                    tGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTATAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACcaca                                                                                 <  1:4541/140‑1 (MQ=255)                                                     gCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTATAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACAt                                                                                >  2:1064643/1‑140 (MQ=255)                                                        gAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTATAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATtga                                                                             <  2:531364/140‑1 (MQ=255)                                                         aaGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTATAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATtgat                                                                            >  1:176257/1‑140 (MQ=255)                                                         aaGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTATAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATtgat                                                                            >  1:884449/1‑140 (MQ=255)                                                          aGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTATAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATtgatg                                                                           >  2:648224/1‑140 (MQ=255)                                                          aGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTATAACACCCCGGCGCGGCACAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATtgatg                                                                           >  2:45810/1‑140 (MQ=255)                                                              cgcgATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTATAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACca                                                                                   <  1:910759/128‑1 (MQ=255)                                                              cgcgATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTATAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACca                                                                                   >  2:910759/1‑128 (MQ=255)                                                                cgATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTATAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATTGATGAGatca                                                                     >  1:359163/1‑140 (MQ=255)                                                                   tatGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTATAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATTGATGAGATCATcc                                                                  <  2:241737/140‑1 (MQ=255)                                                                             aCGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTATAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATTGATGAGATCATCCGCCGCATTgc                                                        >  2:1003631/1‑140 (MQ=255)                                                                             aCGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTATAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATTGATGAGATCATCCGCCGCATTgc                                                        >  2:787357/1‑140 (MQ=255)                                                                             aCGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTATAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATTGATGAGATCATCCGCCGCATTgc                                                        >  2:510810/1‑140 (MQ=255)                                                                             aCGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTATAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATTGATGAGATCATCCGCCGCATTgc                                                        >  2:403643/1‑140 (MQ=255)                                                                              cGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTATAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATTGATGAGATCATCCGCCGCATTgcg                                                       <  1:1003631/140‑1 (MQ=255)                                                                                          cggcgCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTATAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTc                                           <  2:358657/140‑1 (MQ=255)                                                                                           ggcgCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTATAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCg                                          <  2:176257/140‑1 (MQ=255)                                                                                               cATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTATAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGt                                      <  2:359163/140‑1 (MQ=255)                                                                                                         ggACGACGCTGGAGGTGCTGGATCTGTTCTATAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAAc                            >  1:359124/1‑140 (MQ=255)                                                                                                         ggACGACGCTGGAGGTGCTGGATCTGTTCTATAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAAc                            <  1:754055/140‑1 (MQ=255)                                                                                                         ggACGACGCTGGAGGTGCTGGATCTGTTCTATAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAAc                            <  1:787357/140‑1 (MQ=255)                                                                                                          gacgacGCTGGAGGTGCTGGATCTGTTCTATAACACCCCGGCGCGGCACAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAAcc                           <  2:860029/140‑1 (MQ=255)                                                                                                           acgacgCTGGAGGTGCTGGATCTGTTCTATAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCt                          >  1:30428/1‑140 (MQ=255)                                                                                                                  tGGAGGTGCTGGATCTGTTCTATAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCAt                   <  1:1064643/140‑1 (MQ=255)                                                                                                                         gCTGGATCTGTTCTATAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAAc                            >  2:135007/1‑124 (MQ=255)                                                                                                                         gCTGGATCTGTTCTATAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAAc                            <  1:135007/124‑1 (MQ=255)                                                                                                                         gCTGGATCTGTTCTATAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTa            >  1:792020/1‑140 (MQ=255)                                                                                                                          cTGGATCTGTTCTATAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTaa           >  2:163043/1‑140 (MQ=255)                                                                                                                          cTGGATCTGTTCTATAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTaa           >  2:505581/1‑140 (MQ=255)                                                                                                                               tCTGTTCTATAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATtg      <  1:510810/140‑1 (MQ=255)                                                                                                                               tCTGTTCTATAACACCCCGGCGCGGCACAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATtg      <  1:45810/140‑1 (MQ=255)                                                                                                                                  gTTCTATAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGc   >  1:340235/1‑140 (MQ=255)                                                                                                                                   ttCTATAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCg  >  2:246101/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                       CCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTACAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCG  >  NC_000913/4397732‑4398001

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