Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,397,557	G→T	G49V (GGT→GTT)	mutL →	methyl‑directed mismatch repair protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,397,557	0	G	T	100.0%	240.0 / NA	71	G49V (GGT→GTT)	mutL	methyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (35/36);  total (35/36)

GGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGC  >  NC_000913/4397460‑4397653                                                                                                   |                                                                                                 ggtgaggtgGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTg                                                                                                <  1:685077/100‑1 (MQ=255) ggtgaggtgGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTg                                                                                                >  2:919752/1‑100 (MQ=255)  gtgaggtgGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTgg                                                                                               >  2:1338989/1‑100 (MQ=255)   tgaggtgGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTggg                                                                                              <  2:538085/100‑1 (MQ=255)     aggtgGTCGAGCGACCTGCGTCGGTAGTCAAAGAATTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGc                                                                                            <  2:468711/100‑1 (MQ=255)       gtggtCGCGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACTCGTTTGGGCg                                                                                           <  1:566231/99‑1 (MQ=255)       gtggtCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGc                                                                                            <  2:1641907/98‑1 (MQ=255)          gtCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAAc                                                                                       <  1:669501/100‑1 (MQ=255)            cGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACtt                                                                                     >  1:867040/1‑100 (MQ=255)              aGCGACCTGCGTCTGTAGTCAAAGAACTAGTGGAAAACATCTTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTa                                                                                    <  2:541581/99‑1 (MQ=255)                cGACCTGCGTCGGTAGTCAAAGAACGAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATc                                                                                  <  1:888798/99‑1 (MQ=255)                 gACCTGCGTCGGTAGTCAAAGAATTAGTGGAAAACAGGCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATc                                                                                  <  2:144573/98‑1 (MQ=255)                 gACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCg                                                                                >  2:1043038/1‑100 (MQ=255)                 gACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCg                                                                                >  2:878444/1‑100 (MQ=255)                          cGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGa                                                                       >  2:229379/1‑100 (MQ=255)                          cGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGa                                                                       >  1:779971/1‑100 (MQ=255)                           ggTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTg                                                                        <  1:693376/98‑1 (MQ=255)                              aGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAAc                                                                   >  2:1125608/1‑100 (MQ=255)                                 cAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGc                                                                >  2:1443520/1‑100 (MQ=255)                                 cAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGCGCGAAACTTAACCGCATTCGTGATAACgg                                                                 >  1:1582802/1‑99 (MQ=255)                                 cAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑AAGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGc                                                                >  1:1538309/1‑100 (MQ=255)                                     gAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCg                                                            <  2:334691/100‑1 (MQ=255)                                      aaCTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCgg                                                           <  2:142246/100‑1 (MQ=255)                                       aCTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGt                                                          <  2:710606/100‑1 (MQ=255)                                        cTAGTGGAAAACACCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTCGTGCGAAACTTATCTGCATTCGAAATAACGGCTGctg                                                           >  2:1667670/1‑96 (MQ=255)                                           gTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCa                                                      >  2:980691/1‑100 (MQ=255)                                             ggAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGGTGcacttgaggt                                                    >  2:1092058/1‑91 (MQ=255)                                             ggAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCaaa                                                    >  2:991614/1‑100 (MQ=255)                                                aaaCAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCaaaaaa                                                 >  1:1758192/1‑100 (MQ=255)                                                    aGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATg                                             <  2:867046/100‑1 (MQ=255)                                                    aGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATg                                             <  2:1368645/100‑1 (MQ=255)                                                     gCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGa                                            <  2:249349/100‑1 (MQ=255)                                                     gCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAAAGGCTGCGGTATCAAAAAAGATGa                                            >  2:1331861/1‑100 (MQ=255)                                                       cTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgc                                          >  1:977997/1‑100 (MQ=255)                                                        tCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAATATGAgct                                         <  2:956185/100‑1 (MQ=255)                                                        tCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgct                                         <  2:669792/100‑1 (MQ=255)                                                                ggTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctgg                                 >  2:874915/1‑100 (MQ=255)                                                                ggTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATGCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctgg                                 >  1:1157604/1‑100 (MQ=255)                                                                gccccaaCGCGT‑TCGATATCAATATCGATACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctg                                  >  2:96915/7‑99 (MQ=25)                                                                 gTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggc                                <  1:1426825/100‑1 (MQ=255)                                                                 ccccaaCGCGT‑TCGATATCAATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcg                               >  1:1094585/6‑100 (MQ=255)                                                                  tGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggc                                <  1:441868/99‑1 (MQ=255)                                                                      aCGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcgctgg                           <  1:229378/100‑1 (MQ=255)                                                                       cgcgTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcgctggc                          <  2:200846/100‑1 (MQ=255)                                                                       cgcgTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcgctggc                          <  2:989199/100‑1 (MQ=255)                                                                          gTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTc                        >  2:1497566/1‑99 (MQ=255)                                                                          cgttCGATATCAATATCGATACGCGTCGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTc                        >  2:817324/4‑100 (MQ=37)                                                                           tATCGATATTGATATCGA‑ACTCGTTGGGTCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCg                       <  2:287670/99‑1 (MQ=255)                                                                            aTCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTc                     >  2:199380/1‑100 (MQ=255)                                                                             tCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCa                    >  2:318767/1‑100 (MQ=255)                                                                               gATATTGATATCGA‑ACGCGTTGGGGCGATACTTATCCGCATTCGTGATAACGGCTGCGGTATCAATAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCAt                   <  1:1091362/99‑1 (MQ=255)                                                                               gATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATg                  <  2:1035189/100‑1 (MQ=255)                                                                                atatTGATATCGA‑ACGCGTTGGTGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCCTGc                 >  1:489342/1‑100 (MQ=255)                                                                                atatTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATg                  <  2:948394/99‑1 (MQ=255)                                                                                 tatTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGcc                >  1:619756/1‑100 (MQ=255)                                                                                 tatTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGcc                <  2:530797/100‑1 (MQ=255)                                                                                 tatTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGcc                <  2:826601/100‑1 (MQ=255)                                                                                    tGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGccacc             <  1:919746/100‑1 (MQ=255)                                                                                    tGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATGCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGccacc             <  2:418531/100‑1 (MQ=255)                                                                                     gATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGccacca            >  1:671691/1‑100 (MQ=255)                                                                                       tatCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGt          >  2:1252527/1‑100 (MQ=255)                                                                                        atCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTa         >  1:1669685/1‑100 (MQ=255)                                                                                         tCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTTATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTa         <  2:708718/99‑1 (MQ=255)                                                                                         tCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTaa        >  1:1667710/1‑100 (MQ=255)                                                                                         tCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTaa        >  2:389569/1‑100 (MQ=255)                                                                                          cGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTaaa       <  1:4569/100‑1 (MQ=255)                                                                                          cGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACCGCTGCGGTATCAAAATAGAAGAGCTGGCGCTGGAGCTGGCTCGTCATGCCAACAGTa         >  1:444372/1‑98 (MQ=255)                                                                                           gA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTaaaa      >  1:39168/1‑100 (MQ=255)                                                                                           gA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTaaaa      <  1:874909/100‑1 (MQ=255)                                                                                               cgcgTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGAATAAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCg   <  1:502919/100‑1 (MQ=255)                                                                                                gcgTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGc  <  2:218129/100‑1 (MQ=255)                                                                                                   |                                                                                                 GGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGC  >  NC_000913/4397460‑4397653

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