Predicted mutation
evidence seq id position mutation annotation gene description
RA minE 1,827,842 Δ1 bp coding (1085/2562 nt) mutS → methyl‑directed mismatch repair protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,827,8400C.98.1% 211.8 / ‑1.0 53coding (1083/2562 nt)mutSmethyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base . (23/29);  total (23/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

ACAGCCGGTACTGCGTCAGGTCGGCGACCTGGAACGTATTCTGGCACGTCTGGCTTTACGAACTGCTCGCCCACGCGATCTGGCCCGTATGCGCCACGCTTTCCAGCAACTGCCGGAGCTGCGTGCGCAGTTAGA  >  minE/1827771‑1827905
                                                                     |                                                                 
aCAGCCGGTACTGCGTCAGGTCGGCGACCTGGAACGTATTCTGGCACGTCTGGCTTTACGAACTGCTCGc                                                                   >  1:2887799/1‑70 (MQ=255)
aCAGCCGGTACTGCGTCAGGTCGGCGACCTGGAACGTATTCTGGCACGTCTGGCTTTACGAACTGCTCGc                                                                   >  1:2097595/1‑70 (MQ=255)
 cAGCCGGTACTGCGTCAGGTCGGCGACCTGGAACGTATTCTGGCACGTCTGGCTTTACGAACTGCTCGc                                                                   <  1:2943074/69‑1 (MQ=255)
     cGGTACTGCGTCAGGTCGGCGACCTGGAACGTATTCTGGCACGTCTGGCTTTACGAACTGCTCGCCAcgcg                                                            <  1:1041886/71‑1 (MQ=255)
     cGGTACTGCGTCAGGTCGGCGACCTGGAACGTATTCTGGCACGTCTGGCTTTACGAACTGCTCGCCAcgcg                                                            <  1:1932999/71‑1 (MQ=255)
       gTACTGCGTCAGGTCGGCGACCTGGAACGTATTCTGGCACGTCTGGCTTTACGAACTGCTCGCCACGCGAt                                                          <  1:651444/71‑1 (MQ=255)
        tACTGCGTCAGGTCGGCGACCTGGAACGTATTCTGGCACGTCTGGCTTTACGAACTGCTCGCCACGCGATc                                                         <  1:1225313/71‑1 (MQ=255)
         aCTGCGTCAGGTCGGCGACCTGGAACGTATTCTGGCACGTCTGGCTTTACGAACTGCTCGCCACGCGATc                                                         <  1:3230/70‑1 (MQ=255)
         aCTGCGTCAGGTCGGCGACCTGGAACGTATTCTGGCACGTCTGGCTTTACGAACTGCTCGCCACGCGATCt                                                        >  1:2781206/1‑71 (MQ=255)
          cTGCGTCAGGTCGGCGTCCTGGAACGTATTCTGGCACGTCTGGCTTTACGAACTGCTCGCCACGCGATCTg                                                       <  1:2773051/71‑1 (MQ=255)
          cTGCGTCAGGTCGGCGACCTGGAACGTATTCTGGCACGTCTGGCTTTACGAACTGCTCGCCACGCGATCTg                                                       >  1:716059/1‑71 (MQ=255)
               tCAGGTCGGCGACCTGGAACGTATTCTGGCACGTCTGGCTTTACGAACTGCTCGCCACGCGATCTGGCCCg                                                  >  1:2575562/1‑71 (MQ=255)
               tCAGGTCGGCGACCTGGAACGTATTCTGGCACGTCTGGCTTTACGAACTGCTCGCCACGCGATCTGGCCCg                                                  >  1:182459/1‑71 (MQ=255)
                cAGGTCGGCGACCTGGAACGTATTCTGGCACGTCTGGCTTTACGAACTGCTCGCCACGCGATCTGGCCCg                                                  >  1:1706150/1‑70 (MQ=255)
                 aGGTCGGCGACCTGGAACGTATTCTGGCACGTCTGGCTTTACGAACTGCTCGCCACGCGATCTGGCCCGTa                                                >  1:1546038/1‑71 (MQ=255)
                 aGGTCGGCGACCTGGAACGTATTCTGGCACGTCTGGCTTTACGAACTGCTCGCCACGCGATCTGGCCCGTa                                                >  1:32450/1‑71 (MQ=255)
                 aGGTCGGCGACCTGGAACGTATTCTGGCACGTCTGGCTTTACGAACTGCTCGCCACGCGATCTGGCCCGTa                                                >  1:974963/1‑71 (MQ=255)
                  ggTCGGCGACCTGGAACGTATTCTGGCACGTCTGGCTTTACGAACTGCTCGCCACGCGATCTGGCCCGTAt                                               <  1:2567353/71‑1 (MQ=255)
                     cGGCGACCTGGAACGTATTCTGGCACGTCTGGCTTTACGAACTGCTCGCCACGCGATCt                                                        >  1:1264378/1‑59 (MQ=255)
                     cGGCGACCTGGAACGTATTCTGGCACGTCTGGCTTTACGAACTGCTCGCCACGCGATCTGGCCCGTATgcg                                            >  1:1918158/1‑71 (MQ=255)
                     cGGCGACCTGGAACGTATTCTGGCACGTCTGGCTTTACGAACTGCTCGCCACGCGATCTGGCCCGTATgcg                                            >  1:1571352/1‑71 (MQ=255)
                     cGGCGACCTGGAACGTATTCTGGCACGTCTGGCTTTACGAACTGCTCGCCACGCGATCTGGCCCGTATgcg                                            >  1:604141/1‑71 (MQ=255)
                         gACCTGGAACGTATTCTGGCACGTCTGGCTTTACGAACTGCTCGCCACGCGATCTGGCCCGTATGCGCCAc                                        >  1:102930/1‑71 (MQ=255)
                          aCCTGGAACGTATTCTGGCACGTCTGGCTTTACGAACTGCTCGCCACGCGATCTGGCCCGTATGCGCCACg                                       <  1:730964/71‑1 (MQ=255)
                           catGGAACGTATTCTGGCACGTCTGGCTTTACGAACTGCTCGCCACGCGATCTGGCCCGTATGCGCCACGc                                      >  1:2471490/3‑71 (MQ=255)
                            cTGGAACGTATTCTGGCACGTCTGGCTTTACGAGCTGCTCGCCACGCGATCTGGCCCGTATGCGCCACGCt                                     <  1:1035297/71‑1 (MQ=255)
                            cTGGAACGTATTCTGGCACGTCTGGCTTTACGAACTGCTCGCCACGCGATCTGGCCCGTATGCGCCACGCt                                     <  1:533173/71‑1 (MQ=255)
                            cTGGAACGTATTCTGGCACGTCTGGCTTTACGAACTGCTCGCCACGCGATCTGGCCCGTATGCGCCACGCt                                     <  1:1880766/71‑1 (MQ=255)
                               gAACGTATTCTGGCACGTCTGGCTTTACGAACTGCTCGCCACGCGATCTGGCCCGTATGCGCCACGCTTTc                                  <  1:2479882/71‑1 (MQ=255)
                               gAACGTATTCTGGCACGTCTGGCTTTACGAACTGCTCGCCACGCGATCTGGCCCGTATGCGCCACGCTTTc                                  >  1:1358133/1‑71 (MQ=255)
                               gAACGTATTCTGGCACGTCTGGCTTTACGAACTGCTCGCCACGCGATCTGGCCCGTATGCGCCACGCTTTc                                  <  1:1117271/71‑1 (MQ=255)
                                aaCGTATTCTGGCACGTCTGGCTTTACGAACTGCTCGCCACGCGATCTGGCCCGTATGCGCCACGCTTTc                                  <  1:2838416/70‑1 (MQ=255)
                                 aCGTATTCTGGCACGTCTGGCTTTACGAACTGCTCGCCACGCGATc                                                         >  1:2767782/1‑46 (MQ=255)
                                  cGTATTCTGGCACGTCTGGCTTTACGAACTGCTCGCCACGCGATCTGGCCCGTATGCGCCACGCtt                                    >  1:2403979/1‑66 (MQ=255)
                                   gTATTCTGGCACGTCTGGCTTTACGAACTGCTCGCCACGCGATCTGGCCCGTATGCGCCACGCTTTCCAGc                              <  1:2606384/71‑1 (MQ=255)
                                   gTATTCTGGCACGTCTGGCTTTACGAACTGCTCGCCACGCGATCTGGCCCGTATGCGCCACGCTTTCCAGc                              <  1:1282672/71‑1 (MQ=255)
                                    tATTCTGGCACGTCTGGCTTTACGAACTGCTCGCCACGCGATCTGGCCCGTATGCGCCACGCTTTCCAGCa                             <  1:664009/71‑1 (MQ=255)
                                      ttCTGGCACGTCTGGCTTTACGAACTGCTCGCCACGCGATCTGGCCCGTATGCGCCACGCTTTCCAGCaa                            <  1:1819741/70‑1 (MQ=255)
                                          ggCACGTCTGGCTTTACGAACTGCTCGCCACGCGATCTGGCCCGTATg                                              >  1:1830290/1‑48 (MQ=255)
                                           gCACGTCTGGCTTTACGAACTGCTCGCCACGCGATCTGGCCCGTATGCGCCACGCTTTCCAGCAACTGc                        <  1:2250316/69‑1 (MQ=255)
                                              cGTCTGGCTTTACGAACTGCTCGCCACGCGATCTGGCCCGTATGCGCCACGCTTTCCAGCAACTGCCGGAg                   <  1:239516/71‑1 (MQ=255)
                                                tCTGGCTTTACGAACTGCTCGCCACGCGATCTGGCCCGTATGCGCCACGCTTTCCAGCAACTGCCGGAGCt                 <  1:2711377/71‑1 (MQ=255)
                                                    gCTTTACGAACTGCTCGCCACGCGATCTGGCCCGTATGCGCCACGCTTTCCAGCAACTGCCGGAGCtgcgt             >  1:780706/1‑71 (MQ=255)
                                                    gCTTTACGAACTGCTCGCCACGCGATCTGGCCCGTATGCGCCACGCTTTCCAGCAACTGCCGGAGCtgc               <  1:3048827/69‑1 (MQ=255)
                                                      tttACGAACTGCTCGCCACGCGATCTGGCCCGTATGCGCAACGCTTTCCAGCAACTGCCGGAGCtgcgtg            >  1:775398/1‑70 (MQ=255)
                                                       ttACGAACTGCTCGCCACGCGATCTGGCCCGTATGCGCCACGCTTTCCAGCAACTGCCGGAGCtgcgtg            <  1:1971184/69‑1 (MQ=255)
                                                         aCGAACTGCTCGCCACGCGATCTGGCCCGTATGCGCCACGCTTTCCAGCAACTGCCGGAGCtgcgtgcg          >  1:356316/1‑69 (MQ=255)
                                                         aCGAACTGCTCGCCACGCGATCTGGCCCGTATGCGCCACGCTTTCCAGCAACTGCCGGAGCtgcgtgcg          >  1:397155/1‑69 (MQ=255)
                                                         aCGAACTGCTCGCCACGCGATCTGGCCCGTATGCGCCACGCTTTCCAGCAACTGCCGGAGCTGCGTgcgc         <  1:2512174/70‑1 (MQ=255)
                                                          cGAACTGCTCGCCACGCGATCTGGCCCGTATGCGCCACGCTTTCCAGCAACTGCCGGAGCTGCGTGCGCa        <  1:1185564/70‑1 (MQ=255)
                                                           gAACTGCTCGCCACGCGATCTGGCCCGTATGCGCCACGCTTTCCAGCAACTGCCGGAGCTGCGTGCGCAGt      <  1:2729936/71‑1 (MQ=255)
                                                            aaCTGCTCGCCACGCGATCTGGCCCGTATGCGCCACGCTTTCCAGCAACTGCCGGAGCTGCGTGCGCAGtt     <  1:715991/71‑1 (MQ=255)
                                                             aCTGCTCGCCACGCGATCTGGCCCGTATGCGCCACGCtt                                    <  1:2572966/39‑1 (MQ=38)
                                                             aCTGCTCGCCACGCGATCTGGCCCGTATGCGCCACGCTTTCCAGCAACTg                         >  1:1812230/1‑50 (MQ=39)
                                                             aCTGCTCGCCACGCGATCTGGCCCGTATGCGCCACGCTTTCCAGCAACTGCCGGAGCTGCGTGCGCAGTTa    <  1:153865/71‑1 (MQ=255)
                                                                 cTCGCCCACGCGATCTGGCCCGTATGCGCCACGCTTTCCAGCAACTGCCGGAGCTGCGTGCGCAGTTAGa  <  1:998564/70‑1 (MQ=255)
                                                                     |                                                                 
ACAGCCGGTACTGCGTCAGGTCGGCGACCTGGAACGTATTCTGGCACGTCTGGCTTTACGAACTGCTCGCCCACGCGATCTGGCCCGTATGCGCCACGCTTTCCAGCAACTGCCGGAGCTGCGTGCGCAGTTAGA  >  minE/1827771‑1827905

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