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% 146.7 / NA 47*230* (TAG→TAA) mutHmethyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (19/28);  total (19/28)

GCGAATGCGAAAGCGCTTACCGAAGCCATTGGTGCCCGGGGCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAGCCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGATAACCGATCCTAACGCCTGGCTTGCGCTC  >  NC_000913/2970223‑2970464
                                                                                                                                |                                                                                                                 
gcGAATGCGAAAGCGCTTACCGAAGCCATTGGTGCCCGGGGCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCttt                                                                                                         >  2:316534/1‑139 (MQ=255)
 cGAATGCGAAAGCGCTTACCGAAGCCATTGGTGCCCGGGGCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTg                                                                                                        >  2:1382705/1‑139 (MQ=255)
     tGCGAAAGCGCTTACCGAAGCCATTGGTGCCCGGGGCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCt                                                                                                    <  1:369078/139‑1 (MQ=255)
              gcTTACCGAAGTCATTGGTGCCCGGGGAGAACTGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTc                                                                                           <  2:350741/139‑1 (MQ=255)
              gcTTACCGAAGCCATTGGTGCCCGGGGCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTc                                                                                           <  2:887180/139‑1 (MQ=255)
              gcTTACCGAAGCCATAGGTGCCCGGGGCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCACTTTc                                                                                           <  1:1552783/139‑1 (MQ=255)
                  aCCGAAGCCATTGGTGCCCGGGGCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCggg                                                                                       <  2:917352/139‑1 (MQ=255)
                  aCCGAAGCCATTGGTGCCCGGGGCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCggg                                                                                       <  2:1561581/139‑1 (MQ=255)
                   ccGAAGCCATTGGTGCCCGGGGCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGc                                                                                      <  1:1313612/139‑1 (MQ=255)
                    cGAAGCCATTGGTGCCCGGGGCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCa                                                                                     <  1:1494459/139‑1 (MQ=255)
                     gAAGCCATTGGTGCCCGGGGCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCat                                                                                    <  1:1507883/139‑1 (MQ=255)
                     gAAGCCATTGGTGCCCGGGGCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCat                                                                                    <  1:316534/139‑1 (MQ=255)
                       aGCCATTGGTGCCCGGGGCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCatat                                                                                  >  2:897420/1‑139 (MQ=255)
                       aGCCATTGGTGCCCGGGGCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCa                                                                                     >  2:487350/1‑136 (MQ=255)
                       aGCCATTGGTGCCCGGGGCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCa                                                                                     <  1:487350/136‑1 (MQ=255)
                            ttGGTGCCCGGGGCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTa                                                                             >  2:1578445/1‑139 (MQ=255)
                                   ccGGGGCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTc                                                                      <  1:951359/139‑1 (MQ=255)
                                    cGGGGCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCa                                                                     >  2:1267606/1‑139 (MQ=255)
                                       ggCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCAttt                                                                  >  2:679719/1‑139 (MQ=255)
                                               taTTCTGGCGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGGGCACTACTGGCCCGTAATTTTCAGATCCAGTAACCATCGTTTTGACCTGCCGCTTTCCGGGCATATAATTACCTCTTCATTTTTTTGGCa                                                          <  2:411811/138‑1 (MQ=255)
                                                aTTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAg                                                         <  1:1578445/139‑1 (MQ=255)
                                                aTTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAg                                                         <  1:897420/139‑1 (MQ=255)
                                                                              aagaagaaTTTCATCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTc                                                                                           <  1:1305866/75‑1 (MQ=255)
                                                                              aagaagaaTTTCATCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTc                                                                                           >  2:1305866/1‑75 (MQ=255)
                                                                              aagaagaaTTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTc                                                                                           >  2:1157606/1‑75 (MQ=255)
                                                                              aagaagaaTTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTc                                                                                           <  1:1157606/75‑1 (MQ=255)
                                                                              aagaagaaTTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCg                                                                                         <  1:1284173/77‑1 (MQ=255)
                                                                              aagaagaaTTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCg                                                                                         >  2:1284173/1‑77 (MQ=255)
                                                                              aagaagaaTTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGATAAc                           >  2:797553/1‑139 (MQ=255)
                                                                                                      ctGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAggg                                                       >  1:473292/1‑87 (MQ=255)
                                                                                                      ctGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAggg                                                       <  2:473292/87‑1 (MQ=255)
                                                                                                      ctGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGATAACCGATCCTAACGCCTGGCTTGCGCt   <  1:270574/139‑1 (MQ=255)
                                                                                                       tGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGATAACCGATCCTAACGCCTGGCTTGCGCTc  <  1:1382705/139‑1 (MQ=255)
                                                                                                         gCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTc                                                                                           >  2:1331064/1‑48 (MQ=255)
                                                                                                         gCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTc                                                                                           <  1:521001/48‑1 (MQ=255)
                                                                                                         gCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTc                                                                                           <  1:1331064/48‑1 (MQ=255)
                                                                                                         gCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTc                                                                                           >  2:521001/1‑48 (MQ=255)
                                                                                                         gCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGc                                                                                      <  1:178557/53‑1 (MQ=255)
                                                                                                         gCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGc                                                                                      >  2:178557/1‑53 (MQ=255)
                                                                                                         gCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCa                                                          <  1:1419378/81‑1 (MQ=255)
                                                                                                         gCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCa                                                          >  2:1419378/1‑81 (MQ=255)
                                                                                                         gCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATCTTTTTGGCAggg                                                       <  1:274075/84‑1 (MQ=255)
                                                                                                         gCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATCTTTTTGGCAggg                                                       >  2:274075/1‑84 (MQ=255)
                                                                                                              tCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGc                                    >  1:1566407/1‑98 (MQ=255)
                                                                                                              tCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGc                                    <  2:1566407/98‑1 (MQ=255)
                                                                                                                             gTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTa                                         >  1:810189/1‑78 (MQ=255)
                                                                                                                             gTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTa                                         <  2:810189/78‑1 (MQ=255)
                                                                                                                                |                                                                                                                 
GCGAATGCGAAAGCGCTTACCGAAGCCATTGGTGCCCGGGGCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAGCCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGATAACCGATCCTAACGCCTGGCTTGCGCTC  >  NC_000913/2970223‑2970464

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