Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,398,693	C→T	L428F (CTT→TTT)	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,398,693	0	C	T	100.0%	118.9 / NA	34	L428F (CTT→TTT)	mutL	DNA mismatch repair protein MutL
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (12/22);  total (12/22)

TCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGCCAGTGGCAGAACGTTGGCTGCGTCAGGCACAATT  >  NC_000913/4398566‑4398826                                                                                                                                |                                                                                                                                      tcCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCATTTGCGGCGa                                                                                                                              >  2:760/1‑137 (MQ=255) tcCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCATTTGCGGCGa                                                                                                                              <  1:760/137‑1 (MQ=255)   cGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCATTTGCGGCGAACAGt                                                                                                                         >  2:575508/1‑140 (MQ=255)   cGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCATTTGCGGCGAACAGt                                                                                                                         >  2:366801/1‑140 (MQ=255)               ggCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCATTTGCGGCGAACAGTCAGAGTTTTGGt                                                                                                             <  1:767945/140‑1 (MQ=255)               ggCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCATTTGCGGCGAACAGTCAGAGTTTTGGt                                                                                                             <  1:366801/140‑1 (MQ=255)               ggCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCATTTGCGGCGAACAGTCAGAGTTTTGGt                                                                                                             <  1:575508/140‑1 (MQ=255)                gCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCATTTGCGGCGAACAGTCAGAGTTTTGGTc                                                                                                            <  1:1138463/140‑1 (MQ=255)                  cGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCATTTGCGGCGAACAGTCAGAGTTTTGGTCgg                                                                                                          <  1:1170291/140‑1 (MQ=255)                  cGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCATTTGCGGCGAACAGTCAGAGTTTTGGTCgg                                                                                                          <  2:184485/140‑1 (MQ=255)                  cGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCATTTGCGGCGAACAGTCAGAGTTTTGGTCgg                                                                                                          <  2:1176081/140‑1 (MQ=255)                     aTGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCATTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTa                                                                                                       >  1:260633/1‑140 (MQ=255)                     aTGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCATTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTa                                                                                                       <  2:1042189/140‑1 (MQ=255)                     aTGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCATTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTa                                                                                                       <  2:673143/140‑1 (MQ=255)                         gcAGCCAGGCTACCAGAAACAGCAAGGTGAAGCGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCATTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGa                                                                                                   <  1:1044000/140‑1 (MQ=255)                          cAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCATTTGCGGCGAACAGTCAGAGTTTTGGTCgg                                                                                                          >  2:935272/1‑132 (MQ=255)                          cAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCATTTGCGGCGAACAGTCAGAGTTTTGGTCgg                                                                                                          <  1:935272/132‑1 (MQ=255)                           aGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCATTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACt                                                                                                 <  1:400421/140‑1 (MQ=255)                           aGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCATTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACt                                                                                                 >  2:907514/1‑140 (MQ=255)                               aGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCATTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCg                                                                                             <  1:232336/140‑1 (MQ=255)                                            cAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCATTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACtg                                                                                <  1:907514/140‑1 (MQ=255)                                            cAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCATTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACtg                                                                                <  2:260633/140‑1 (MQ=255)                                               cAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCATTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGc                                                                             >  2:131968/1‑140 (MQ=255)                                                          gtATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCATTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAgc                                                                  <  2:208353/140‑1 (MQ=255)                                                             tCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCATTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAgcgc                                                                <  1:55906/139‑1 (MQ=255)                                                             tCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCATTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAgcgc                                                                >  2:55906/1‑139 (MQ=255)                                                             tCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCATTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGcgcg                                                               >  1:207579/1‑140 (MQ=255)                                                                ccAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCATTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACg                                                            >  1:216327/1‑140 (MQ=255)                                                                         gCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCATTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACAttt                                                   <  2:216327/140‑1 (MQ=255)                                                                                    cgcCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCATTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATcc                                        >  2:1210554/1‑140 (MQ=255)                                                                                        gATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCATTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGc                                    <  1:1210554/140‑1 (MQ=255)                                                                                        gATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCATTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGc                                    <  1:885660/140‑1 (MQ=255)                                                                                                                     aGGAACCTGCATTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGCCAGTGGCAGAACGTTGGCTGCGTCAGGca       <  2:207579/140‑1 (MQ=255)                                                                                                                          ccTGCATTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCCGGAGCGCGACGGCAACATTTCACTTTTATCCTTGCCAGTGGCAGAACGTTGGCTGCGTCAGGCACAAtt  >  1:232063/1‑140 (MQ=255)                                                                                                                                |                                                                                                                                      TCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGCCAGTGGCAGAACGTTGGCTGCGTCAGGCACAATT  >  NC_000913/4398566‑4398826

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