Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,970,351 G→A *230* (TAG→TAA mutH → methyl‑directed mismatch repair protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,970,3510GA100.0% 104.8 / NA 33*230* (TAG→TAAmutHmethyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (17/16);  total (17/16)

CGAATGCGAAAGCGCTTACCGAAGCCATTGGTGCCCGGGGCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAGCCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGATAACCGATCCTAACGCCTGGCTTGCGCTCGGTACGCTGACGCTGCTGGA  >  NC_000913/2970224‑2970484
                                                                                                                               |                                                                                                                                     
cGAATGCGAAAGCGCTTACCGAAGCCATTGGTGCCCGGGGCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTg                                                                                                                            >  2:68823/1‑139 (MQ=255)
          aGCGCTTACCGAAGCCATTGGTGCCCGGGGCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCt                                                                                                                  >  2:1051479/1‑139 (MQ=255)
                 aCCGAAGCCATTGGTGCCCGGGGCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCggg                                                                                                           <  1:990019/139‑1 (MQ=255)
                   cGAAGCCATTGGTGCCCGGGGCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCa                                                                                                         <  1:1051479/139‑1 (MQ=255)
                                  ccGGGGCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTc                                                                                          >  2:691065/1‑139 (MQ=255)
                                      ggCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGc                                                                                                                                >  2:537184/1‑97 (MQ=255)
                                      ggCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGc                                                                                                                                <  1:537184/97‑1 (MQ=255)
                                               aTTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAg                                                                             <  1:278991/139‑1 (MQ=255)
                                                                    ttttATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGc                                                        <  1:691065/139‑1 (MQ=255)
                                                                         tttGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGa                                                                  >  2:172867/1‑124 (MQ=255)
                                                                         tttGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGa                                                                  <  1:172867/124‑1 (MQ=255)
                                                                             aagaagaaTTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCa                                                                              >  2:263189/1‑108 (MQ=255)
                                                                             aagaagaaTTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCa                                                                              <  1:263189/108‑1 (MQ=255)
                                                                                             gTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCa                                                                              <  2:605652/92‑1 (MQ=255)
                                                                                             gTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCa                                                                              >  1:605652/1‑92 (MQ=255)
                                                                                              tGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGtt                                                              >  1:166503/1‑107 (MQ=255)
                                                                                              tGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGtt                                                              <  2:166503/107‑1 (MQ=255)
                                                                                               gCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGATAACCGATCCTAACGCCTGGCt                             >  1:1524770/1‑139 (MQ=255)
                                                                                                cACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGATAACCGATCCTAACGCCTGGCtt                            >  2:1132691/1‑139 (MQ=255)
                                                                                                     ctGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCggg                                                                                                           <  1:1042814/55‑1 (MQ=255)
                                                                                                     ctGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCAGCCGCTTTCAggg                                                                                                           >  2:1042814/1‑55 (MQ=255)
                                                                                                        gCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGcc                                                                                                                     <  1:268464/42‑1 (MQ=255)
                                                                                                        gCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGcc                                                                                                                     >  2:268464/1‑42 (MQ=255)
                                                                                                        gCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAggg                                                                           <  1:108506/84‑1 (MQ=255)
                                                                                                        gCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAggg                                                                           >  2:108506/1‑84 (MQ=255)
                                                                                                        gCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGATAACCGATCCTAACGCCTGGCTTGCGCTCgg                    >  2:650354/1‑139 (MQ=255)
                                                                                                             tCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGATAACCGATCCTAACGCCTGGCTTGCGCTCGGTacgc               >  1:1111882/1‑139 (MQ=255)
                                                                                                                ttttCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCAGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCTTGGATTACCGATCCTAACGCCTGGCTTGCGCTCGGTacgctga            <  1:1132691/139‑1 (MQ=255)
                                                                                                                       aTCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGc                                                        >  1:1330853/1‑88 (MQ=255)
                                                                                                                       aTCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGc                                                        <  2:1330853/88‑1 (MQ=255)
                                                                                                                       aTCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGATa                                                 <  1:771885/95‑1 (MQ=255)
                                                                                                                       aTCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGATa                                                 >  2:771885/1‑95 (MQ=255)
                                                                                                                          cAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGATAACCGATCCTAACGCCTGGCTTGCGCTCGGTACGCTGACGCTGCTGga  <  2:1111882/139‑1 (MQ=255)
                                                                                                                               |                                                                                                                                     
CGAATGCGAAAGCGCTTACCGAAGCCATTGGTGCCCGGGGCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAGCCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGATAACCGATCCTAACGCCTGGCTTGCGCTCGGTACGCTGACGCTGCTGGA  >  NC_000913/2970224‑2970484

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: