Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,221,351 T→C 34.7% I51I (ATT→ATC accC → acetyl‑CoA carboxylase, biotin carboxylase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,221,3510TC34.7% 60.6 / 51.5 90I51I (ATT→ATCaccCacetyl‑CoA carboxylase, biotin carboxylase subunit
Reads supporting (aligned to +/- strand):  ref base T (31/27);  new base C (18/13);  total (50/40)
Fisher's exact test for biased strand distribution p-value = 8.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.41e-01

TCGCTGTGCACTCCAGCGCGGATCGCGATCTAAAACACGTATTACTGGCAGATGAAACGGTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCG  >  minE/2221284‑2221418
                                                                   |                                                                   
tCGCTGTGCACTCCAGCGCGGATCGCTATCTAAAACACGTATTACTGGCAGATGAAACGGTCTGTATTg                                                                    >  1:1677118/1‑69 (MQ=255)
tCGCTGTGCACTCCAGCGCGGATCGCGATCTAAAACACGTATTACTGGCAGATGAAACGGTCTGTATTg                                                                    >  1:2934874/1‑69 (MQ=255)
tCGCTGTGCACTCCAGCGCGGATCGCGATCTAAAACACGTATTACTGGCAGATGAAACGGTCTGTATTg                                                                    >  1:2296701/1‑69 (MQ=255)
tCGCTGTGCACTCCAGCGCGGATCGCGATCTAAAACACGTATTACTGGCAGATGAAACGGTCTGTATTg                                                                    >  1:2307242/1‑69 (MQ=255)
tCGCTGTGCACTCCAGCGCGGATCGCGATCTAAAACACGTATTACTGGCAGATGAAACGGTCTGTATTg                                                                    >  1:317181/1‑69 (MQ=255)
tCGCTGTGCACTCCAGCGCGGATCGCGATCTAAAACACGTATTACTGGCAGATGAAACGGTCTGTATTg                                                                    >  1:3434387/1‑69 (MQ=255)
tCGCTGTGCACTCCAGCGCGGATCGCGATCTAAAACACGTATTACTGGCAGATGAAACGGTCTGTATTg                                                                    >  1:559253/1‑69 (MQ=255)
tCGCTGTGCACTCCAGCGCGGATCGCGATCTAAAACACGTATTACTGGCAGATGAAACGGTCTGTATTg                                                                    >  1:1183261/1‑69 (MQ=255)
tCGCTGTGCACTCCAGCGCGGATCGCGATCTAAAACACGTATTACTGGCAGATGAAACGGTCTGTATTg                                                                    >  1:1161475/1‑69 (MQ=255)
tCGCTGTGCACTCCAGCGCGGATCGCGATCTAAAACACGTATTACTGGCAGATGAAACGGTCTGTATTg                                                                    >  1:1150632/1‑69 (MQ=255)
tCGCTGTGCACTCCAGCGCGGATCGAGATCTAAAACACGTATTACTGGCAGATGAAACGGTCTGTATTg                                                                    >  1:2339418/1‑69 (MQ=255)
  gCTGTGCACTCCAGCGCGGATCGCGATCTAAAACACGTATTACTGGCAGATGAAACGGTCTGTATTgg                                                                   >  1:848267/1‑68 (MQ=255)
  gCTGTGCACTCCAGCGCGGATCGCGATCTAAAACACGTATTACTGGCAGATGAAACGGTCTGTATTGGc                                                                  >  1:677439/1‑69 (MQ=255)
  gCTGTGCACTCCAGCGCGGATCGCGATCTAAAACACGTATTACTGGCAGATGAAACGGTCTGTATTGGc                                                                  >  1:2829975/1‑69 (MQ=255)
  gCTGTGCACTCCAGCGCGGATCGCGATCTAAAACACGTATTACTGGCAGATGAAACGGTCTGTATTGGc                                                                  >  1:2796250/1‑69 (MQ=255)
  gCTGTGCACTCCAGCGCGGATCGCGATCTAAAACACGTATTACTGGCAGATGAAACGGTCTGTATTGGc                                                                  >  1:2677960/1‑69 (MQ=255)
     gtgCACTCCAGCGCGGATCGCGATCTAAAACACGTATTACTGGCAGATGAAACGGTCTGTATTGGCCCt                                                               >  1:1774612/1‑69 (MQ=255)
        cACTCCAGCGCGGATCGCGATCTAAAACACGTATTACTGGCAGATGAAACGGTCTGTATTGGCCCTGc                                                             <  1:2569333/68‑1 (MQ=255)
        cACTCCAGCGCGGATCGCGATCTAAAACACGTATTACTGGCAGATGAAACGGTCTGTATTGGCCCTGCt                                                            >  1:64171/1‑69 (MQ=255)
        cACTCCAGCGCGGATCGCGATCTAAAACACGTATTACTGGCAGATGAAACGGTCTGTATTGGCCCTGCt                                                            >  1:2599368/1‑69 (MQ=255)
        cACTCCAGCGCGGATCGCGATCTAAAACACGTATTACTGGCAGATGAAACGGTCTGTATTGGCCCTGCt                                                            >  1:1959913/1‑69 (MQ=255)
            ccAGCGCGGATCGCGATCTAAAACACGTATTACTGGCAGATGAAACGGTCTGTATTGGCCCTGCTcc                                                          >  1:2918615/1‑67 (MQ=255)
            ccAGCGCGGATCGCGATCTAAAACACGTATTACTGGCAGATGAAACGGTCTGTATTGGCCCTGCTcc                                                          >  1:1992966/1‑67 (MQ=255)
            ccAGCGCGGATCGCGATCTAAAACACGTATTACTGGCAGATGAAACGGTCTGTATTGGCCCTGCTcc                                                          >  1:985686/1‑67 (MQ=255)
            ccAGCGCGGATCGCGATCTAAAACACGTATTACTGGCAGATGAAACGGTCTGTATTGGCCCTGCTcc                                                          >  1:926045/1‑67 (MQ=255)
            ccAGCGCGGATCGCGATCTAAAACACGTATTACTGGCAGATGAAACGGTCTGTATGGGCCCTGCTcc                                                          >  1:1964088/1‑67 (MQ=255)
                       cgcgATCTAAAACACGTATTACTGGCAGATGAAACGGTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGt                                             >  1:1631571/1‑69 (MQ=255)
                       cgcgATCTAAAACACGTATTACTGGCAGATGAAACGGTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGt                                             >  1:1543346/1‑69 (MQ=255)
                       cgcgATCTAAAACACGTATTACTGGCAGATGAAACGGTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGt                                             >  1:1183522/1‑69 (MQ=255)
                       cgcgATCTAAAACACGTATTACTGGCAGATGAAACGGTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGt                                             >  1:1652144/1‑69 (MQ=255)
                       cgcgATCTAAAACACGTATTACTGGCAGATGAAACGGTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGt                                             >  1:1823228/1‑69 (MQ=255)
                             cTAAAACACGTATTACTGGCAGATGAAACGGTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTg                                       >  1:961824/1‑69 (MQ=255)
                                  acacGTATTACTGGCAGATGAAACGGTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTg                                       <  1:1287485/64‑1 (MQ=255)
                                  acacGTATTACTGGCAGATGAAACGGTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTg                                       <  1:2349214/64‑1 (MQ=255)
                                  acacGTATTACTGGCAGATGAAACGGTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTg                                       <  1:888990/64‑1 (MQ=255)
                                  acacGTATTACTGGCAGATGAAACGGTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTg                                       <  1:1648692/64‑1 (MQ=255)
                                  acacGTATTACTGGCAGATGAAACGGTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTg                                       <  1:1945831/64‑1 (MQ=255)
                                  acacGTATTACTGGCAGATGAAACGGTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTg                                       <  1:1926623/64‑1 (MQ=255)
                                  acacGTATTACTGGCAGATGAAACGGTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTg                                       <  1:296425/64‑1 (MQ=255)
                                  acacGTATTACTGGCAGATGAAACGGTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTg                                       <  1:168336/64‑1 (MQ=255)
                                  acacGTATTACTGGCAGATGAAACGGTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTGa                                      <  1:626586/65‑1 (MQ=255)
                                  acacGTATTACTGGCAGATGAAACGGTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTGa                                      <  1:180761/65‑1 (MQ=255)
                                  acacGTATTACTGGCAGATGAAACGGTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTGa                                      <  1:2067206/65‑1 (MQ=255)
                                  acacGTATTACTGGCAGATGAAACGGTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTGa                                      <  1:1756275/65‑1 (MQ=255)
                                  acacGTATTACTGGCAGATCAAACGGTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTGa                                      <  1:1733854/65‑1 (MQ=255)
                                  acacGTATTACTCGCAGATGAAACGGTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTg                                       <  1:2863348/64‑1 (MQ=255)
                                   cacGTATTACTGGCAGATGAAACGGTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTg                                       <  1:3065405/63‑1 (MQ=255)
                                   cacGTATTACTGGCAGATGAAACGGTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTg                                       <  1:124592/63‑1 (MQ=255)
                                   cacGTATTACTGGCAGATGAAACGGTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTg                                       <  1:1143549/63‑1 (MQ=255)
                                           aCTGGCAGATGAAACGGTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCaa                          >  1:1437097/1‑68 (MQ=255)
                                                 aGATGAAACGGTCTGTATCGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAatcatc                     <  1:1454898/67‑1 (MQ=255)
                                                 aGATGAAACGGTCTGTATCGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAatcatc                     <  1:1887488/67‑1 (MQ=255)
                                                 aGATGAAACGGTCTGTATCGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAatcatc                     <  1:2932848/67‑1 (MQ=255)
                                                 aGATGAAACGGTCTGTATCGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAatcatc                     <  1:526114/67‑1 (MQ=255)
                                                 aGATGAAACGGTCTGTATCGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAatcatc                     <  1:3071528/67‑1 (MQ=255)
                                                 aGATGAAACGGTCTGTATCGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAatcatc                     <  1:1634728/67‑1 (MQ=255)
                                                 aGATGAAACGGTCTGTATCGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAatcatc                     <  1:2667253/67‑1 (MQ=255)
                                                 aGATGAAACGGTCTGTATCGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAatcatc                     <  1:354507/67‑1 (MQ=255)
                                                 aGATGAAACGGTCTGTATCGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAatcatc                     <  1:728403/67‑1 (MQ=255)
                                                 aGATGAAACGGTCTGTATCGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAatcatc                     <  1:2519953/67‑1 (MQ=255)
                                                 aGATGAAACGGTCTGTATCGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAatcatc                     <  1:1582556/67‑1 (MQ=255)
                                                 aGATGAAACGGTCTGTATCGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAatcatc                     <  1:2383964/67‑1 (MQ=255)
                                                 aGATGAAACGGTCTGTATCGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAatcatc                     <  1:777585/67‑1 (MQ=255)
                                                   aTGAAACGGTCTGTATCTGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAatca                       >  1:622368/1‑63 (MQ=255)
                                                   aTGAAACGGTCTGTATCGGCCCTGCTCCGTCAGTAAAAAGTTAt                                          >  1:8090/1‑44 (MQ=255)
                                                   aTGAAACGGTCTGTATCGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCaa                          >  1:2905250/1‑60 (MQ=255)
                                                   aTGAAACGGTCTGTATCGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAatca                       >  1:1142388/1‑63 (MQ=255)
                                                   aTGAAACGGTCTGTATCGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAatca                       >  1:1680643/1‑63 (MQ=255)
                                                   aTGAAACGGTCTGTATCGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAatca                       >  1:927342/1‑63 (MQ=255)
                                                   aTGAAACGGTCTGTATCGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAatca                       >  1:513063/1‑63 (MQ=255)
                                                   aTGAAACGGTCTGTATCGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAatca                       >  1:1223614/1‑63 (MQ=255)
                                                   aTGAAACGGTCTGTATCGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAatca                       >  1:1207708/1‑63 (MQ=255)
                                                   aTGAAACGGTCTGTATCGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAatca                       >  1:726532/1‑63 (MQ=255)
                                                   aTGAAACGGTCTGTATCGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAatca                       >  1:736123/1‑63 (MQ=255)
                                                   aTGAAACGGTCTGTATCGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAatca                       >  1:1824975/1‑63 (MQ=255)
                                                   aTGAAACGGTCTGTATCGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAatca                       >  1:1868101/1‑63 (MQ=255)
                                                   aTGAAACGGTCTGTATCGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAatca                       >  1:2728832/1‑63 (MQ=255)
                                                   aTGAAACGGTCTGTATCGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAatca                       >  1:2049337/1‑63 (MQ=255)
                                                   aTGAAACGGTCTGTATCGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAatca                       >  1:2163104/1‑63 (MQ=255)
                                                   aTGAAACGGTCTGTATCGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAatca                       >  1:150211/1‑63 (MQ=255)
                                                   aTGAAACGGTCTGTATCGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAatca                       >  1:2456771/1‑63 (MQ=255)
                                                          ggTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGa          <  1:1449222/69‑1 (MQ=255)
                                                          ggTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGa          <  1:908488/69‑1 (MQ=255)
                                                          ggTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGa          <  1:2676094/69‑1 (MQ=255)
                                                          ggTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGa          <  1:2180120/69‑1 (MQ=255)
                                                           gTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGaa         <  1:1660224/69‑1 (MQ=255)
                                                           gTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGaa         <  1:3270057/69‑1 (MQ=255)
                                                           gTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGaa         <  1:2650231/69‑1 (MQ=255)
                                                           gTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGaa         <  1:138097/69‑1 (MQ=255)
                                                            tCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGaa         <  1:2742084/68‑1 (MQ=255)
                                                                  ttGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCAcc   <  1:162956/68‑1 (MQ=255)
                                                                  ttGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCAcc   >  1:1322940/1‑68 (MQ=255)
                                                                  ttGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCg  >  1:709636/1‑69 (MQ=255)
                                                                  ttGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCg  >  1:1197084/1‑69 (MQ=255)
                                                                  ttGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCg  >  1:1925593/1‑69 (MQ=255)
                                                                  ttGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCg  >  1:55167/1‑69 (MQ=255)
                                                                  ttGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCg  >  1:940701/1‑69 (MQ=255)
                                                                   |                                                                   
TCGCTGTGCACTCCAGCGCGGATCGCGATCTAAAACACGTATTACTGGCAGATGAAACGGTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCG  >  minE/2221284‑2221418

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