Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,857,878 T→C R262R (CGT→CGC mutS → methyl‑directed mismatch repair protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,857,8780TC97.3% 124.3 / ‑3.9 37R262R (CGT→CGCmutSmethyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (20/16);  total (20/17)
Fisher's exact test for biased strand distribution p-value = 4.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01

GTTGAATCTGCAATTTGGGACCCGCGATCTGGTCGGTTTTGGCGTCGAGAACGCGCCGCGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAGTATGCGAAAGATACCCAACGTACGACTCTGCCGCATATTCGTTCCATCACCATGGAACGTGAGCAGGACAGCATCATTATGGATGCCGCGACGCGTCGTAATCTGGAAATCACCCAGAACCTGGCGGGTGGTGCGGAAAATACGCTGGCTTCTGTGCTCGACTGCACCGTCACGCCGATGGGCAGCCGTATGCTGAAACGCTG  >  NC_000913/2857728‑2858021
                                                                                                                                                      |                                                                                                                                               
gTTGAATCTGCAATTTGGGACCCGCGATCTGGTCGGTTTTGGCGTCGAGAACGCGCCGCGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAGTATGCGAAAGATACCCAACGTACGACTCTGCCGCATATTCGTTCCATCACCATGGAACGt                                                                                                                                                 <  1:580739/151‑1 (MQ=255)
   gAATCTGCAATTTGGGACCCGCGATCTGGTCGGTTTTGGCGTCGAGAACGCGCCGCGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAGTATGCGAAAGATACCCAACGTACGACTCTGCCGCATATTCGTTCCATCACCATGGAACGCGAg                                                                                                                                              <  2:446470/151‑1 (MQ=255)
         gCAATTTGGGACCCGCGATCTGGTCGGTTTTGGCGTCGAGAACGCGCCGCGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAGTATGCGAAAGATACCCAACGTACGACTCTGCCGCATATTCGTTCCATCACCATGGAACGTGAGCAGGa                                                                                                                                         <  2:455755/150‑1 (MQ=255)
         gCAATTTGGGACCCGCGATCTGGTCGGTTTTGGCGTCGAGAACGCGCCGCGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAGTATGCGAAAGATACCCAACGTACGACTCTGCCGCATATTCGTTCCATCACCATGGAACGCGAGCAgg                                                                                                                                          <  1:140558/149‑1 (MQ=255)
                          aTCTGGTCGGTTTTGGCGTCGAGAACGCGCCGCGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAGTATGCGAAAGATACCCAACGTACGACTCTGCCGCATATTCGTTCCATCACCATGGAACGCGAGCAGGACAGCATCATTATGGATGc                                                                                                                       >  2:35998/1‑151 (MQ=255)
                              ggTCGGTTTTGGCGTCGAGAACGCGCCGCGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAGTATGCGAAAGATACCCAACGTACGACTCTGCCGCATATTCGTTCCATCACCATGGAACGCGAGCAGGACAGCATCATTATGGATGCcgcg                                                                                                                   <  2:390115/151‑1 (MQ=255)
                                                 aaCGCGCCGCGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAGTATGCGAAAGATACCCAACGTACGACTCTGCCGCATATTCGTTCCATCACCATGGAACGCGAGCAGGACAGCATCATTATGGATGCCGCGACGCGTCGTAATCTGGaa                                                                                                 >  2:457082/1‑150 (MQ=255)
                                                            cgGACTTTGTGCTGCCGGTTGTCTGTTGCAGTATGCGAAAGATACCCAACGTACGACTCTGCCGCATATTCGTTCCATCACCATGGAACGCGAGCAGGACAGCATCATTATGGATGCCGCGACGCGTCGTAATCTGGAAATCACCCAGAAc                                                                                     >  2:68149/1‑151 (MQ=255)
                                                                   tgtgCTGCCGGTTGTCTGTTGCAGTATGCGAAAGATACCCAACGTACGACTCTGCCGCATATTCGTTCCATCACCATGGAACGCGAGCAGGACAGCATCATTATGGATGCCGCGACGCGTCGTAATCTGGAAATCACCCAGAACCTGGCg                                                                               <  2:546353/150‑1 (MQ=255)
                                                                      gctgcCGGTTGTCTGTTGCAGTATGCGAAAGATACCCAACGTACGACTCTGCCGCATATTCGTTCCATCACCATGGAACGCGAGCAGGACAGCATCATTATGGATGCCGCGACGCGTCGTAATCTGGAAATCACCCAGAACCTGGCGggtg                                                                           >  2:505151/1‑151 (MQ=255)
                                                                      gctgcCGGTTGTCTGTTGCAGTATGCGAAAGATACCCAACGTACGACTCTGCCGCATATTCGTTCCATCACCATGGAACGCGAGCAGGACAGCATCATTATGGATGCCGCGACGCGTCGTAATCTGGAAATCACCCAGAACCTGGCGggt                                                                            <  1:603540/150‑1 (MQ=255)
                                                                      gctgcCGGTTGTCTGTTGCAGTATGCGAAAGATACCCAACGTACGACTCTGCCGCATATTCGTTCCATCACCATGGAACGCGAGCAGCACAGCATCATTATGGATGCCGCGACGCGTCGTAATCTGGAAATCACCCAGAACCTGGCGggtg                                                                           >  2:549598/1‑151 (MQ=255)
                                                                       ctgcCGGTTGTCTGTTGCAGTATGCGAAAGATACCCAACGTACGACTCTGCCGCATATTCGTTCCATCACCATGGAACGCGAGCAGGACAGCATCATTATGGATGCCGCGACGCGTCGTAATCTGGAAATCACCCAGAACCTGGCGggtgg                                                                          >  1:28643/1‑151 (MQ=255)
                                                                           cGGTTGTCTGTTGCAGTATGCGAAAGATACCCAACGTACGACTCTGCCGCATATTCGTTCCATCACCATGGAACGCGAGCAGGACAGCATCATTATGGATGCCGCGACGCGTCGTAATCTGGAAATCACCCAGAACCTGGCGGGTGGTGCg                                                                      <  2:149236/151‑1 (MQ=255)
                                                                            ggTTGTCTGTTGCAGTATGCGAAAGATACCCAACGTACGACTCTGCCGCATATTCGTTCCATCACCATGGAACGCGAGCAGGACAGCATCATTATGGATGCCGCGACGCGTCGTAATCTGGAAATCACCCAGAACCTGGCGGGTGGTGCgg                                                                     >  2:535983/1‑151 (MQ=255)
                                                                             gTTGTCTGTTGCAGTATGCGAAAGATACCCAACGTACGACTCTGCCGCATATTCGTTCCATCACCATGGAACGCGAGCAGGACAGCATCATTATGGATGCCGCGACGCGTCGTAATCTGGAAATCACCCAGAACCTGGCGGGGGGTGCgg                                                                     >  1:465448/1‑150 (MQ=255)
                                                                                 tCTGTTGCAGTATGCGAAAGATACCCAACGTACGACTCTGCCGCATATTCGTTCCATCACCATGGAACGCGAGCAGGACAGCATCATTATGGATGCCGCGACGCGTCGTAATCTGGAAATCACCCAGAACCTGGCGGGTGGTGCGGAAAAt                                                                >  1:118638/1‑151 (MQ=255)
                                                                                    gTTGCAGTATGCGAAAGATACCCAACGTACGACTCTGCCGCATATTCGTTCCATCACCATGGAACGCGAGCAGGACAGCATCATTATGGATGCCGCGACGCGTCGTAATCTGGAAATCACCCAGAACCTGGCGGGTGGTGCGGAAAATAc                                                              <  1:446034/150‑1 (MQ=255)
                                                                                                 aaaGATACCCAACGTACGACTCTGCCGCATATTCGTTCCATCACCATGGAACGCGAGCAGGACAGCATCATTATGGATGCCGCGACGCGTCGTAATCTGGAAATCACCCAGAACCTGGCGGGTGGTGCGGAAAATACGCTGGCTTCTGTGc                                                <  1:29491/151‑1 (MQ=255)
                                                                                                     aTACCCAACGTACGACTCTGCCGCATATTCGTTCCATCACCATGGAACGCGAGCAGGACAGCATCATTATGGATGCCGCGACGCGTCGTAATCTGGAAATCACCCAGAACCTGGCGGGTGGTGCGGAAAATACGCTGGCTTCTGTGCTCg                                             >  2:530568/1‑150 (MQ=255)
                                                                                                        cccAACGTACGACTCTGCCGCATATTCGTTCCATCACCATGGAACGCGAGCAGGACAGCATCATTATGGATGCCGCGACGCGTCGTAATCTGGAAATCACCCAGAACCTGGCGGGTGGTGCGGAAAATACGCTGGCTTCTGTGCTCGACTg                                         >  1:248213/1‑151 (MQ=255)
                                                                                                         ccAACGTACGACTCTGCCGCATATTCGTTCCATCACCATGGAACGCGAGCAGGACAGCATCATTATGGATGCCGCGACGCGTCGTAATCTGGAAATCACCCAGAACCTGGCGGGTGGTGCGGAAAATACGCTGGCTTCTGTGCTCGACTGc                                        <  1:12753/151‑1 (MQ=255)
                                                                                                           aaCGTACGACTCTGCCGCATATTCGTTCCATCACCATGGAACGCGAGCAGGACAGCATCATTATGGATGCCGCGACGCGTCGTAATCTGGAAATCACCCAGAACCTGGCGGGTGGTGCGGAAAATACGCTGGCTTCTGTGCTCGACTGCa                                       >  1:545875/1‑150 (MQ=255)
                                                                                                              gTACGACTCTGCCGCATATTCGTTCCATCACCATGGAACGCGAGCAGGACAGCATCATTATGGATGCCGCGACGCGTCGTAATCTGGAAATCACCCAGAACCTGGCGGGTGGTGCGGAAAATACGCTGGCTTCTGTGCTCGACTGCACCGt                                   >  1:166840/1‑151 (MQ=255)
                                                                                                                 cGACTCTGCCGCATATTCGTTCCATCACCATGGAACGCGAGCAGGACAGCATCATTATGGATGCCGCGACGCGTCGTAATCTGGAAATCACCCAGAACCTGGCGGGTGGTGCGGAAAATACGCTGGCTTCTGTGCTCGACTGCACCGTCAc                                <  1:293675/151‑1 (MQ=255)
                                                                                                                  gACTCTGCCGCATATTCGTTCCATCACCATGGAACGCGAGCAGGACAGCATCATTATGGATGCCGCGACGCGTCGTAATCTGGAAATCACCCAGAACCTGGCGGGTGGTGCGGAAAATACGCTGGCTTCTGTGCTCGACTGCACCGTCACg                               <  2:450865/151‑1 (MQ=255)
                                                                                                                        gCCGCATATTCGTTCCATCACCATGGAACGCGAGCAGGACAGCATCATTATGGATGCCGCGACGCGTCGTAATCTGGAAATCACCCAGAACCTGGCGGGTGGTGCGGAAAATACGCTGGCTTCTGTGCTCGACTGCACCGTCACGCCGATg                         >  1:133815/1‑151 (MQ=255)
                                                                                                                           gCATATTCGTTCCATCACCATGGAACGCGAGCAGGACAGCATCATTATGGATGCCGCGACGCGTCGTAATCTGGAAATCACCCAGAACCTGGCGGGTGGTGCGGAAAATACGCTGGCTTCTGTGCTCGACTGCACCGTCACGCCGATGGGc                      <  1:457090/151‑1 (MQ=255)
                                                                                                                            cATATTCGTTCCATCACCATGGAACGCGAGCAGGACAGCATCATTATGGATGCCGCGACGCGTCGTAATCTGGAAATCACCCAGAACCTGGCGGGTGGTGCGGAAAATACGCTGGCTTCTGTGCTCGACTGCACCGTCACGCCGATGGGc                      >  1:304977/1‑150 (MQ=255)
                                                                                                                              tatTCGTTCCATCACCATGGAACGCGAGCAGGACAGCATCATTATGGATGCCGCGACGCGTCGTAATCTGGAAATCACCCAGAACCTGGCGGGTGGTGCGGAAAATACGCTGGCTTCTGTGCTCGACTGCACCGTCACGCCGATGGGCAGc                   >  2:49910/1‑151 (MQ=255)
                                                                                                                               atTCGTTCCATCACCATGGAACGCGAGCAGGACAGCATCATTATGGATGCCGCGACGCGTCGTAATCTGGAAATCACCCAGAACCTGGCGGGTGGTGCGGAAAATACGCTGGCTTCTGTGCTCGACTGCACCGTCACGCCGATGGGCAGcc                  >  1:370121/1‑151 (MQ=255)
                                                                                                                                  cGTTCCATCACCATGGAACGCGAGCAGGACAGCATCATTATGGATGCCGCGACGCGTCGTAATCTGGAAATCACCCAGAACCTGGCGGGTGGTGCGGAAAATACGCTGGCTTCTGTGCTCGACTGCACCGTCACGCCGATGGGCAGCCGt                <  1:68655/150‑1 (MQ=255)
                                                                                                                                   gTTCCATCACCATGGAACGCGAGCAGGACAGCATCATTATGGATGCCGCGACGCGTCGTAATCTGGAAATCACCCAGAACCTGGCGGGTGGTGCGGAAAATACGCTGGCTTCTGTGCTCGACTGCACCGTCACGCCGATGGGCAGCCGTAt              >  2:458809/1‑151 (MQ=255)
                                                                                                                                    ttCCATCACCATGGAACGCGAGCAGGACAGCATCATTATGGATGCCGCGACGCGTCGTAATCTGGAAATCACCCAGAACCTGGCGGGTGGTGCGGAAAATACGCTGGCTTCTGTGCTCGACTGCACCGTCACGCCGATGGGCAGCCGTATg             >  2:516027/1‑151 (MQ=255)
                                                                                                                                     tCCATCACCATGGAACGCGAGCAGGACAGCATCATTATGGATGCCGCGACGCGTCGTAATCTGGAAATCACCCAGAACCTGGCGGGTGGTGCGGAAAATACGCTGGCTTCTGTGCTCGACTGCACCGTCACGCCGATGGGCAGCCGTATg             <  2:58740/150‑1 (MQ=255)
                                                                                                                                      ccATCACCATGGAACGCGAGCAGGACAGCATCATTATGGATGCCGCGTCGCGTCGTAATCTGGAAATCACCCAGAACCTGGCGGGTGATGCGGAAAATACGCTGGCTTCTGTGCTCGACTGCACCGTCACGCCGATGGGCAGCCGTATGCt           <  1:565238/151‑1 (MQ=255)
                                                                                                                                              aTGGAACGCGAGCAGGACAGCATCATTATGGATGCCGCGACGCGTCGTAATCTGGAAATCACCCAGAACCTGGCGGGTGGTGCGGAAAATACGCTGGCTTCTGTGCTCGACTGCACCGTCACGCCGATGGGCAGCCGTATGCTGAAACgct   <  2:548171/151‑1 (MQ=255)
                                                                                                                                               tGGAACGCGAGCAGGACAGCATCATTATGGATGCCGCGACGCGTCGTAATCTGGAAATCACCCAGAACCTGGCGGGTGGTGCGGAAAATACGCTGGCTTCTGTGCTCGACTGCACCGTCACGCCGATGGGCAGCCGTATGCTGAAACgctg  >  1:598958/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                               
GTTGAATCTGCAATTTGGGACCCGCGATCTGGTCGGTTTTGGCGTCGAGAACGCGCCGCGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAGTATGCGAAAGATACCCAACGTACGACTCTGCCGCATATTCGTTCCATCACCATGGAACGTGAGCAGGACAGCATCATTATGGATGCCGCGACGCGTCGTAATCTGGAAATCACCCAGAACCTGGCGGGTGGTGCGGAAAATACGCTGGCTTCTGTGCTCGACTGCACCGTCACGCCGATGGGCAGCCGTATGCTGAAACGCTG  >  NC_000913/2857728‑2858021

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