Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,319,007 C→T 44.3% K295K (AAG→AAA ugd ← UDP‑glucose 6‑dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,319,0070CT44.3% 19.4 / 53.0 72K295K (AAG→AAAugdUDP‑glucose 6‑dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (38/2);  new base T (17/15);  total (55/17)
Fisher's exact test for biased strand distribution p-value = 4.09e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGAAGTTATCTGAACCGCTCTTCATAATCAGACGATAAATACCCACCACTTGCGGCTTGCGTGACAAAATGGCATCGGCAATAAAATCTTTACGCGTGCGGTTAGCATCGACAATTG  >  minE/1318951‑1319068
                                                        |                                                             
cGGAAGTTTTCTGAACCGCTCTTCATAATCAGACGATAAATACCCACCACTTGCGGCTTGCGTGACaaa                                                   >  1:1809874/1‑69 (MQ=255)
cGGAAGTTATCTGAACCGCTCTTCATAATCAGACGATAAATACCCACCACTTGCGGCTTGCGTGACaaa                                                   >  1:2532542/1‑69 (MQ=255)
  gAAGTTATCTGAACCGCTCTTCATAATCAGACGATAAATACCCACCACTTGCGGCTTGCGTGACAAAAt                                                 >  1:1204903/1‑69 (MQ=255)
  gAAGTTATCTGAACCGCTCTTCATAATCAGACGATAAATACCCACCACTTGCGGCTTGCGTGACAAAAt                                                 >  1:1253334/1‑69 (MQ=255)
  gAAGTTATCTGAACCGCTCTTCATAATCAGACGATAAATACCCACCACTTGCGGCTTGCGTGACAAAAt                                                 >  1:1389569/1‑69 (MQ=255)
  gAAGTTATCTGAACCGCTCTTCATAATCAGACGATAAATACCCACCACTTGCGGCTTGCGTGACAAAAt                                                 >  1:2057176/1‑69 (MQ=255)
      ttATCTGAACCGCTCTTCATAATCAGACGATAAATACCCACCACTTGCGGCTTGCGTGACAAAATGGc                                              >  1:2277379/1‑68 (MQ=255)
      ttATCTGAACCGCTCTTCATAATCAGACGATAAATACCCACCACTTGCGGCTTGCGTGACAAAATGGCa                                             >  1:964035/1‑69 (MQ=255)
      ttATCTGAACCGCTCTTCATAATCAGACGATAAATACCCACCACTTGCGGCTTGCGTGACAAAATGGCa                                             >  1:3219266/1‑69 (MQ=255)
       tATCTGAACCGCTCTTCATAATCAGACGATAAATACCCACCACTTGCGGCTTGCGTGACAAAAt                                                 >  1:913983/1‑64 (MQ=255)
       tATCTGAACCGCTCTTCATAATCAGACGATAAATACCCACCACTTGCGGCTTGCGTGACAAAATgg                                               >  1:2241952/1‑66 (MQ=255)
       tATCTGAACCGCTCTTCATAATCAGACGATAAATACCCACCACTTGCGGCTTGCGTGACAAAATgg                                               >  1:64685/1‑66 (MQ=255)
       tATCTGAACCGCTCTTCATAATCAGACGATAAATACCCACCACTTGCGGCTTGCGTGACAAAATgg                                               >  1:3435642/1‑66 (MQ=255)
       tATCTGAACCGCTCTTCATAATCAGACGATAAATACCCACCACTTGCGGCTTGCGTGACAAAATgg                                               >  1:3223094/1‑66 (MQ=255)
       tATCTGAACCGCTCTTCATAATCAGACGATAAATACCCACCACTTGCGGCTTGCGTGACAAAATgg                                               >  1:1267224/1‑66 (MQ=255)
       tATCTGAACCGCTCTTCATAATCAGACGATAAATACCCACCACTTGCGGCTTGCGTGACAAAATgg                                               >  1:1198224/1‑66 (MQ=255)
       tATCTGAACCGCTCTTCATAATCAGACGATAAATACCCACCACTTGCGGCTTGCGTGACAAAATgg                                               >  1:3146801/1‑66 (MQ=255)
       tATCTGAACCGCTCTTCATAATCAGACGATAAATACCCACCACTTGCGGCTTGCGTGACAAAATgg                                               >  1:1034172/1‑66 (MQ=255)
           tGAACCGCTCTTCATAATCAGACGATAAATACCCACCACTTGCGGCTTGCGTGACAAAATGGCATCGGc                                        <  1:526115/69‑1 (MQ=255)
           tGAACCGCTCTTCATAATCAGACGATAAATACCCACCACTTGCGGCTTGCGTGACAAAATGGCATCGGc                                        >  1:3167990/1‑69 (MQ=255)
           tGAACCGCTCTTCATAATCAGACGATAAATACCCACCACTTGCGGCTTGCGTGACAAAATGGCATCGGc                                        >  1:1797911/1‑69 (MQ=255)
                             cAGACGATAAATACCCACCACTTGCGGTTTGCGTGACAAAATGGCATCGGCAATAAAATCTTTAtgcg                       <  1:1770985/68‑1 (MQ=255)
                             cAGACGATAAATACCCACCACTTGCGGTTTGCGTGACAAAATGGCATCGGCAATAAAATCTTTAcgcg                       <  1:934476/68‑1 (MQ=255)
                             cAGACGATAAATACCCACCACTTGCGGTTTGCGTGACAAAATGGCATCGGCAATAAAATCTTTAcgcg                       <  1:2850058/68‑1 (MQ=255)
                             cAGACGATAAATACCCACCACTTGCGGTTTGCGTGACAAAATGGCATCGGCAATAAAATCTTTAcgcg                       <  1:3049459/68‑1 (MQ=255)
                             cAGACGATAAATACCCACCACTTGCGGTTTGCGTGACAAAATGGCATCGGCAATAAAATCTTTAcgcg                       <  1:2581781/68‑1 (MQ=255)
                             cAGACGATAAATACCCACCACTTGCGGTTTGCGTGACAAAATGGCATCGGCAATAAAATCTTTAcgcg                       <  1:1002/68‑1 (MQ=255)
                             cAGACGATAAATACCCACCACTTGCGGTTTGCGTGACAAAATGGCATCGGCAATAAAATCTTTAcgcg                       <  1:1059489/68‑1 (MQ=255)
                             cAGACGATAAATACCCACCACTTGCGGTTTGCGTGACAAAATGGCATCGGCAATAAAATCTTTAcgcg                       <  1:24362/68‑1 (MQ=255)
                             cAGACGATAAATACCCACCACTTGCGGTTTGCGTGACAAAATGGCATCGGCAATAAAATCTTTAcgcg                       <  1:1066018/68‑1 (MQ=255)
                             cAGACGATAAATACCCACCACTTGCGGTTTGCGTGACAAAATGGCATCGGCAATAAAATCTTTAcgcg                       <  1:1257695/68‑1 (MQ=255)
                             cAGACGATAAATACCCACCACTTGCGGTTTGCGTGACAAAATGGCATCGGCAATAAAATCTTTAcgcg                       <  1:1352553/68‑1 (MQ=255)
                             cAGACGATAAATACCCACCACTTGCGGTTTGCGTGACAAAATGGCATCGGCAATAAAATCTTTAcgcg                       <  1:523333/68‑1 (MQ=255)
                             cAGACGATAAATACCCACCACTTGCGGTTTGCGTGACAAAATGGCATCGGCAATAAAATCTTTAcgcg                       <  1:1575434/68‑1 (MQ=255)
                             cAGACGATAAATACCCACCACTTGCGGTTTGCGTGACAAAATGGCATCGGCAATAAAATCTTTAcgcg                       <  1:820094/68‑1 (MQ=255)
                              aGACGATAAATACCCACCACTTGCGGTTTGCGTGACAAAATGGCATCGGCAATAAAATCTTTAcgcg                       <  1:1516198/67‑1 (MQ=255)
                                  gATAAATACCCACCACTTGCGGTTTGCGTGACAAAATGGCATCGGCAAt                                     >  1:917819/1‑49 (MQ=255)
                                  gATAAATACCCACCACTTGCGGTTTGCGTGACAAAATGGCATCGGCAATa                                    >  1:581864/1‑50 (MQ=255)
                                  gATAAATACCCACCACTTGCGGTTTGCGTGACAAAATGGCATCGGCAATa                                    >  1:3287359/1‑50 (MQ=255)
                                  gATAAATACCCACCACTTGCGGTTTGCGTGACAAAATGGCATCGGCAATa                                    >  1:954395/1‑50 (MQ=255)
                                  gATAAATACCCACCACTTGCGGTTTGCGTGACAAAATGGCATCGGCAATa                                    >  1:3317254/1‑50 (MQ=255)
                                  gATAAATACCCACCACTTGCGGTTTGCGTGACAAAATGGCATCGGCAATa                                    >  1:417273/1‑50 (MQ=255)
                                  gATAAATACCCACCACTTGCGGTTTGCGTGACAAAATGGCATCGGCAATa                                    >  1:3134045/1‑50 (MQ=255)
                                  gATAAATACCCACCACTTGCGGTTTGCGTGACAAAATGGCATCGGCAATa                                    >  1:3119448/1‑50 (MQ=255)
                                  gATAAATACCCACCACTTGCGGTTTGCGTGACAAAATGGCATCGGCAATa                                    >  1:987109/1‑50 (MQ=255)
                                  gATAAATACCCACCACTTGCGGTTTGCGTGACAAAATGGCATCGGCAATa                                    >  1:2960064/1‑50 (MQ=255)
                                  gATAAATACCCACCACTTGCGGTTTGCGTGACAAAATGGCATCGGCAATa                                    >  1:1187204/1‑50 (MQ=255)
                                  gATAAATACCCACCACTTGCGGTTTGCGTGACAAAATGGCATCGGCAATa                                    >  1:2666912/1‑50 (MQ=255)
                                  gATAAATACCCACCACTTGCGGTTTGCGTGACAAAATGGCATCGGCAATa                                    >  1:2283899/1‑50 (MQ=255)
                                  gATAAATACCCACCACTTGCGGTTTGCGTGACAAAATGGCATCGGCAATa                                    >  1:2214077/1‑50 (MQ=255)
                                  gATAAATACCCACCACTTGCGGTTTGCGTGACAAAATGGCATCGGCAATa                                    >  1:182781/1‑50 (MQ=255)
                                  gATAAATACCCACCACTTGCGGTTTGCGTGACAAAATGGCATCGGCAATa                                    >  1:174552/1‑50 (MQ=255)
                                  gATAAATACCAACCACTTGCGGTTTGCGTGACAAAATGGCATCGGCAATa                                    >  1:1827491/1‑50 (MQ=255)
                                    tAAATACCCACCACTTGCGGCTTGCGTGACAAAATGGCATCGGCAATAAAATCTTTACGCGTGCGGtt                >  1:2557472/1‑68 (MQ=255)
                                    tAAATACCCACCACTTGCGGCTTGCGTGACAAAATGGCATCGGCAATAAAATCTTTACGCGTGCGGtt                >  1:1319999/1‑68 (MQ=255)
                                    tAAATACCCACCACTTGCGGCTTGCGTGACAAAATGGCATCGGCAATAAAATCTTTACGCGTGCGGtt                >  1:1577153/1‑68 (MQ=255)
                                     aaaTACCCACCACTTGCGGCTTGCGTGACAAAATGGCATCGGCAATAAAATCTTTACGCGTGCGGTTa               >  1:1217933/1‑68 (MQ=255)
                                     aaaTACCCACCACTTGCGGCTTGCGTGACAAAATGGCATCGGCAATAAAATCTTTACGCGTGCGGTTAg              >  1:2814756/1‑69 (MQ=255)
                                     aaaTACCCACCACTTGCGGCTTGCGTGACAAAATGGCATCGGCAATAAAATCTTTACGCGTGCGGTTAg              >  1:1349754/1‑69 (MQ=255)
                                     aaaTACCCACCACTTGCGGCTTGCGTGACAAAATGGCATCGGCAATAAAATCTTTACGCGTGCGGTTAg              >  1:1436871/1‑69 (MQ=255)
                                     aaaTACCCACCACTTGCGGCTTGCGTGACAAAATGGCATCGGCAATAAAATCTTTACGCGTGCGGTTAg              >  1:556437/1‑69 (MQ=255)
                                     aaaTACCCACCACTTGCGGCTTGCGTGACAAAATGGCATCGGCAATAAAATCTTTACGCGTGCGGTTAg              >  1:2029337/1‑69 (MQ=255)
                                     aaaTACCCACCACTTGCGGCTTGCGTGACAAAATGGCATCGGCAATAAAATCTTTACGCGTGCGGTTAg              >  1:2237582/1‑69 (MQ=255)
                                     aaaTACCCACCACTTGCGGCTTGCGTGACAAAATGGCATCGGCAATAAAATCTTTACGAGt                      >  1:3103384/1‑61 (MQ=255)
                                     aaaTAACCACCACTTGCGGCTTGCGTGACAAAATGGCATCGGCAATAAAATCTTTACGCGTGCGGTTa               >  1:1037163/1‑68 (MQ=255)
                                            caccacTTGCGGCTTGCGTGACAAAATGGCATCGGCAATAAAATCTTTACGCGTGCGGTTAGCATCGAc       <  1:1890772/69‑1 (MQ=255)
                                               cacTTGCGGCTTGCGTGACAAAATGGCATCGGCAATAAAATCTTTACGCGTGCGGTTAGCATCGACAAt    >  1:1686322/1‑69 (MQ=255)
                                                 cTTGCGGCTTGCGTGACAAAATGGCATCGGCAATAAAATCTTTACGCGTGCGGTTAGCATCGACAAtt   >  1:794956/1‑68 (MQ=255)
                                                 cTTGCGGCTTGCGTGACAAAATGGCATCGGCAATAAAATCTTTACGCGTGCGGTTAGCATCGACAATTg  >  1:1697390/1‑69 (MQ=255)
                                                 cTTGCGGCTTGCGTGACAAAATGGCATCGGCAATAAAATCTTTACGCGTGCGGTTAGCATCGACAATTg  >  1:1892051/1‑69 (MQ=255)
                                                 cTTGCGGCTTGCGTGACAAAATGGCATCGGCAATAAAATCTTTACGCGTGCGGTTAGCATCGACAATTg  >  1:920398/1‑69 (MQ=255)
                                                 cTTGCGGCTTGCGTGACAAAATGGCATCGGCAATAAAATCTTTACGCGTGCGGTTAGCATCGACAATTg  >  1:2547421/1‑69 (MQ=255)
                                                        |                                                             
CGGAAGTTATCTGAACCGCTCTTCATAATCAGACGATAAATACCCACCACTTGCGGCTTGCGTGACAAAATGGCATCGGCAATAAAATCTTTACGCGTGCGGTTAGCATCGACAATTG  >  minE/1318951‑1319068

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: