Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,827,932 C→T 34.1% A392V (GCG→GTG)  mutS → methyl‑directed mismatch repair protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,827,9320CT34.1% 55.0 / 34.0 44A392V (GCG→GTG) mutSmethyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base C (7/22);  new base T (6/9);  total (13/31)
Fisher's exact test for biased strand distribution p-value = 3.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCCACGCTTTCCAGCAACTGCCGGAGCTGCGTGCGCAGTTAGAAACTGTCGATAGTGCACCGGTACAGGCGCTACGTGAGAAGATGGGCGAGTTTGCCGAGCTGCGCGATCTGCTGGAGCGAGCAATCATCGACACAC  >  minE/1827862‑1828000
                                                                      |                                                                    
cgcCACGCTTTCCAGCAACTGCCGGAGCTGCGTGCGCAGTAGAAACTGTCGATAGTGCACCGGTACAGgc                                                                      >  1:973004/1‑70 (MQ=255)
  ccACGCTTTCCAGCAACTGCCGGAGCTGCGTGCGCAGTTAGAAACTGTCGATAGTGCACCGGTACAGgcgc                                                                    <  1:1243710/71‑1 (MQ=255)
          tCCAGCAACTGCCGGAGCTGCGTGCGCAGTTAGAAACTGTCGATAGTGCACCGGTACAGGTGCTACGt                                                               >  1:1182078/1‑68 (MQ=255)
                 aCTGCCGGAGCTGCGTGCGCAGTTAGAAACTGTCGATAGTGCACCGGTACAGGCGCTACGTGagaa                                                          <  1:2018227/66‑1 (MQ=255)
                    gCCGGAGCTGCGTGCGCAGTTAGAAACTGTCGATAGTGCACCGGTACAGGCGCTACGTGAGAAGATGGGc                                                   >  1:468670/1‑70 (MQ=255)
                    gCCGGAGCTGCGTGCGCAGTTAGAAACTGTCGATAGTGCACCGGTACAGGCGCTACGTGAGAAGATGGGc                                                   >  1:236606/1‑70 (MQ=255)
                    gCCGGAGCTGCGTGCGCAGTTAGAAACTGTCGATAGTGCACCGGTACAGGCGCTACGTGAGAAGATGGGc                                                   >  1:1564642/1‑70 (MQ=255)
                         aGCTGCGTGCGCAGTTAGAAACTGTCGATAGTGCACCGGTACAGGCGCTACGTGAGAAGATGGGCGAGttt                                             <  1:65188/71‑1 (MQ=255)
                          gCTGCGTGCGCAGTTAGAAACTGTCGATAGTGCACCGGTACAGGCGCTACGTGAGAAGATGGGCGAGtt                                              <  1:1890282/69‑1 (MQ=255)
                               gtgcgCAGTTAGAAACTGTCGATAGTGCACCGGTACAGGCGCTACGTGAGAAGATGGGCGAGTTTGCCGAg                                       <  1:1895208/71‑1 (MQ=255)
                                tgcgCAGTTAGAAACTGTCGATAGTGCACCGGTACAGGCGCTACGTGagaa                                                          <  1:982244/51‑1 (MQ=255)
                                tgcgCAGTTAGAAACTGTCGATAGTGCACCGGTACAGGCGCTACGTGagaa                                                          <  1:405504/51‑1 (MQ=255)
                                tgcgCAGTTAGAAACTGTCGATAGTGCACCGGTACAGGCGCTACGTGagaa                                                          <  1:2023933/51‑1 (MQ=255)
                                tgcgCAGTTAGAAACTGTCGATAGTGCACCGGTACAGGCGCTACGTGagaa                                                          <  1:1874658/51‑1 (MQ=255)
                                tgcgCAGTTAGAAACTGTCGATAGTGCACCGGTACAGGCGCTACGTGAGAAGATGGGCGAGTTTGCCGAg                                       <  1:1216474/70‑1 (MQ=255)
                                       ttAGAAACTGTCGATAGTGCACCGGTACAGGTGCTACGTGAGAAGATGGGc                                                   <  1:762782/51‑1 (MQ=255)
                                       ttAGAAACTGTCGATAGTGCACCGGTACAGGTGCTACGTGAGAAGATGGGc                                                   <  1:208997/51‑1 (MQ=255)
                                       ttAGAAACTGTCGATAGTGCACCGGTACAGGTGCTACGTGAGAAGATGGGCGAGTTTGCCGAGCTGCGCGa                               <  1:903191/71‑1 (MQ=255)
                                       ttAGAAACTGTCGATAGTGCACCGGTACAGGTGCTACGTGAGAAGATGGGCGAGTTTGCCGAGCTGCGCGa                               <  1:1844260/71‑1 (MQ=255)
                                       ttAGAAACTGTCGATAGTGCACCGGTACAGGTGCTACGTGAGAAGATGGGCGAGTTTGCCGAGCTGCGCGa                               <  1:899334/71‑1 (MQ=255)
                                        tAGAAACTGTCGATAGTGCACCGGTACAGGTGCTACGTGAGAAGATGGGCGAGTTTGCCGAGCTGCGCGa                               <  1:1281861/70‑1 (MQ=255)
                                         aGAAACTGTCGATAGTGCACCGGTACAGGTGCTACGTGAGAAGATGGGCGAGTTTGcc                                          >  1:1572571/1‑58 (MQ=255)
                                         aGAAACTGTCGATAGTGCACCGGTACAGGTGCTACGTGAGAAGATGGGCGAGTTTGCCGAGCTGCGCGATc                             >  1:567845/1‑71 (MQ=255)
                                         aGAAACTGTCGATAGTGCACCGGTACAGGTGCTACGTGAGAAGATGGGCGAGTTTGCCGAGCTGCGCGATc                             >  1:1285113/1‑71 (MQ=255)
                                           aaaCTGTCGATAGTGCACCGGTACAGGCGCTACGTGAGAAGATGGGCGAGttt                                             <  1:605073/53‑1 (MQ=255)
                                           aaaCTGTCGATAGTGCACCGGTACAGGCGCTACGTGAGAAGATGGGCGAGttt                                             <  1:1981482/53‑1 (MQ=255)
                                           aaaCTGTCGATAGTGCACCGGTACAGGCGCTACGTGAGAAGATGGGCGAGttt                                             <  1:1984084/53‑1 (MQ=255)
                                           aaaCTGTCGATAGTGCACCGGTACAGGCGATACGTGAGAAGATGGGCGAGttt                                             <  1:1883308/53‑1 (MQ=255)
                                                gTCGATAGTGCACCGGTACAGGTGCTACGTGAGAAGATGGGCGa                                                 <  1:1500058/44‑1 (MQ=255)
                                                    aTAGTGCACCGGTACAGGCGCTACGTGAGAAGATGGGCGAGTTTGCCGAGCTGcgcg                                <  1:1499537/57‑1 (MQ=255)
                                                    aTAGTGCACCGGTACAGGCGCTACGTGAGAAGATGGGCGAGTTTGCCGAGCTGCGCGATCTGCTGgagcga                  <  1:1921347/71‑1 (MQ=255)
                                                     tAGTGCACCGGTACAGGCGCTACGTGAGAAGATGGGCGAGTTTGCCGAGCTGCGCGATCTGCTGgagcgag                 <  1:2111588/71‑1 (MQ=255)
                                                     tAGTGCACCGGTACAGGCGCTACGTGAGAAGATGGGCGAGTTTGCCGAGCTGCGCGATCTGCTGgagcgag                 <  1:527631/71‑1 (MQ=255)
                                                      aGTGCACCGGTACAGGCGCTACGTGAGAAGATGGGCGAGTTTGCCGAGCTGCGCGATCTGCTGgagcgag                 >  1:19658/1‑70 (MQ=255)
                                                        tGCACCGGTACAGGTGCTACGTGAGAAGATGGGCGAGTTTGCCGAGCTGCGCGATCTGCTGGAGCGAGCaa              <  1:1184355/71‑1 (MQ=255)
                                                          cACCGGTACAGGCGCTACGTGAGAAGATGGGCGAGTTTGCCGAGCTGCGCGATCTGCTGGAGCGAGCAatc            <  1:104347/71‑1 (MQ=255)
                                                          cACCGGTACAGGCGCTACGTGAGAAGATGGGCGAGTTTGCCGAGCTGCGCGATCTGCTGGAGCGAGCAat             <  1:1413474/70‑1 (MQ=255)
                                                          cACCGGTACAGGCGCTACGTGAGAAGATGGGCGAGTTTGCCGAGCTGCGCGATCTGCTGGAGCGAGCAat             <  1:1470820/70‑1 (MQ=255)
                                                             cGGTACAGGCGCTACGTGAGAAGATGGGCGAGTTTGCCGAGCTGCGCGATCTGCTGGAGCGAGCAatcat          >  1:939503/1‑70 (MQ=255)
                                                               gTACAGGTGCTACGTGAGAAGATGGGCGAGTTTGCCGAGCTGCGCGATCTGCTGGAGCGAGCAATCATCGa       <  1:1394823/71‑1 (MQ=255)
                                                                tACAGGCGCTACGTGAGAAGTTGGGCGAGTTTGCCGAGCTGCGCGATCTGCTGGAGCGAGCAATCATCGac      <  1:1814259/71‑1 (MQ=255)
                                                                tACAGGCGCTACGTGAGAAGATGGGCGAGTTTGCCGAGCTGCGCGATCTGCTGGAGCGAGCAATCATCGa       >  1:1414088/1‑70 (MQ=255)
                                                                   aGGTGCTACGTGAGAAGATGGGCGAGTTTGCCGAGCTGCGCGATCTGCTGGAGCGAGCAat             >  1:1233684/1‑61 (MQ=255)
                                                                   aGGTGCTACGTGAGAAGATGGGCGAGTTTGCCGAGCTGCGCGATCTGCTGGAGCGAGCAat             >  1:878204/1‑61 (MQ=255)
                                                                    ggCGCTACGTGAGAAGATGGGCGAGTTTGCCGAGCTGCGCGATCTGCTGGAGCGAGCAatcatc         <  1:380357/64‑1 (MQ=255)
                                                                    ggCGCTACGTGAGAAGATGGGCGAGTTTGCCGAGCTGCGCGATCTGCTGGAGCGAGCAATCATCGacacac  >  1:1996649/1‑71 (MQ=255)
                                                                     gcgcTACGTGAGAAGATGGGCGAGTTTGCCGAGCTGCGCGATCTGCTGGAGCGAGCAATCATCGaca     >  1:1218949/1‑67 (MQ=255)
                                                                      |                                                                    
CGCCACGCTTTCCAGCAACTGCCGGAGCTGCGTGCGCAGTTAGAAACTGTCGATAGTGCACCGGTACAGGCGCTACGTGAGAAGATGGGCGAGTTTGCCGAGCTGCGCGATCTGCTGGAGCGAGCAATCATCGACACAC  >  minE/1827862‑1828000

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