Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,321,385 T→C 40.3% A197A (GCA→GCG cdh ← CDP‑diacylglycerol phosphotidylhydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,321,3850TC40.3% 24.8 / 46.0 67A197A (GCA→GCGcdhCDP‑diacylglycerol phosphotidylhydrolase
Reads supporting (aligned to +/- strand):  ref base T (26/14);  new base C (14/13);  total (40/27)
Fisher's exact test for biased strand distribution p-value = 3.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGCCGTACCATCGCCAGCCCGTAGCGTCCCATATGTTCCCGCGCCTCAGGTACTTCTTCTGCCAGCATCATAAATGGGCTGCGTTGTACCAGTTCGCTTTCCGTTACCCGACGCGCCAGGTATTC  >  minE/2321325‑2321450
                                                            |                                                                 
ctGCCGTACCATCGCCAGCCCGTAGCGTCCCATATGTTCCCGCGCCTCAGGTACTTCTTCTGCCAGcat                                                           <  1:1010635/69‑1 (MQ=255)
ctGCCGTACCATCGCCAGCCCGTAGCGTCCCATATGTTCCCGCGCCTCAGGTACTTCTTCTGCCAGcat                                                           <  1:2497268/69‑1 (MQ=255)
ctGCCGTACCATCGCCAGCCCGTAGCGTCCCATATGTTCCCGCGCCTCAGGTACTTCTTCTGCCAGcat                                                           <  1:1832501/69‑1 (MQ=255)
        ccATCGCCAGCCCGTAGCGTCCCATATGTTCCCGCGCCTCAGGTACTTCTTCTGCCAGCATCATAAAt                                                    >  1:1490559/1‑68 (MQ=255)
        ccATCGCCAGCCCGTAGCGTCCCATATGTTCCCGCGCCTCAGGTACTTCTTCTGCCAGCATCATAAATg                                                   >  1:2674131/1‑69 (MQ=255)
        ccATCGCCAGCCCGTAGCGTCCCATATGTTCCCGCGCCTCAGGTACTTCTTCTGCCAGCATCATAAATg                                                   >  1:487657/1‑69 (MQ=255)
        ccATCGCCAGCCCGTAGCGTCCCATATGTTCCCGCGCCTCAGGTACTTCTTCTGCCAGCATCATAAATg                                                   >  1:420988/1‑69 (MQ=255)
          aTCGCCAGCCCGTAGCGTCCCATATGTTCCCGCGCCTCAGGTACTTCTTCTGCCAGCATCATAAATggg                                                 >  1:633812/1‑69 (MQ=255)
          aTCGCCAGCCCGTAGCGTCCCATATGTTCCCGCGCCTCAGGTACTTCTTCTGCCAGCATCATAAATggg                                                 >  1:1247582/1‑69 (MQ=255)
               cAGCCCGTAGCGTCCCATATGTTCCCGCGCCTCAGGTACTTCTTCTGCCAGCATCATAAATGGGCTGCg                                            >  1:1590753/1‑69 (MQ=255)
               cAGCCCGTAGCGTCCCATATGTTCCCGCGCCTCAGGTACTTCTTCTGCCAGCATCATAAATGGGCTGCg                                            >  1:141419/1‑69 (MQ=255)
                aGCCCGTAGCTTCCCATATGTTCCCGCGCCTCAGGCACTTCTTCCGCCAGCATCATaaa                                                     >  1:2992349/1‑59 (MQ=255)
                aGCCCGTAGCTTCCCATATGTTCCCGCGCCTCAGGCACTTCTTCCGCCAGCATCATAAATGGGCTGCGt                                           >  1:2650792/1‑69 (MQ=255)
                       aGCGTCCCATATGTTCCCGCGCCTCAGGTACTTCTTCTGCCAGCATCATAAATGGGCTGCGtt                                          >  1:885909/1‑63 (MQ=255)
                       aGCGTCCCATATGTTCCCGCGCCTCAGGTACTTCTTCTGCCAGCATCATAAATGGGCTGCGTTg                                         >  1:2673465/1‑64 (MQ=255)
                       aGCGTCCCATATGTTCCCGCGCCTCAGGTACTTCTTCTGCCAGCATCATAAATGGGCTGCGTTg                                         >  1:2378737/1‑64 (MQ=255)
                       aGCGTCCCATATGTTCCCGCGCCTCAGGTACTTCTTCTGCCAGCATCATAAATGGGCTGCGTTg                                         >  1:1082960/1‑64 (MQ=255)
                       aGCGTCCCATATGTTCCCGCGCCTCAGGTACTTCTTCTGCCAGCATCATAAATGGGCTGCGTTg                                         >  1:3087405/1‑64 (MQ=255)
                       aGCGTCCCATATGTTCCCGCGCCTCAGGTACTTCTTCTGCCAGCATCATAAATGGGCTGCGTTg                                         >  1:3193586/1‑64 (MQ=255)
                       aGCGTCCCATATGTTCCCGCGCCTCAGGTACTTCTTCTGCCAGCATCATAAATGGGCTGCGTTg                                         >  1:3219899/1‑64 (MQ=255)
                       aGCGTCCCATATGTTCCCGCGCCTCAGGTACTTCTTCTGCCAGCATCATAAATGGGCTGCGTTg                                         >  1:781128/1‑64 (MQ=255)
                       aGCGTCCAATATGTTCCCGCGCCTCAGGTACTTCTTCTGCCAGCATCATAAATGGGCTGCGTTg                                         >  1:3039529/1‑64 (MQ=255)
                           tCCCATATGTTCCCGCGCCTCAGGCACTTCTTCCGCCAGCATCATAAATGGGCTGCGTTGTAcc                                     >  1:2573139/1‑64 (MQ=255)
                           tCCCATATGTTCCCGCGCCTCAGGCACTTCTTCCGCCAGCATCATAAATGGGCTGCGTTGTACCAGtt                                 >  1:207651/1‑68 (MQ=255)
                           tCCCATATGTTCCCGCGCCTCAGGCACTTCTTCCGCCAGCATCATAAATGGGCTGCGTTGTACCAGTTc                                >  1:39915/1‑69 (MQ=255)
                           tCCCATATGTTCCCGCGCCTCAGGCACTTCTTCCGCCAGCATCATAAATGGGCTGCGTTGTACCAGTTc                                >  1:651453/1‑69 (MQ=255)
                           tCCCATATGTTCCCGCGCCTCAGGCACTTCTTCCGCCAGCATCATAAATGGGCTGCGTTGTACCAGTTc                                >  1:2583885/1‑69 (MQ=255)
                           tCCCATATGTTCCCGCGCCTCAGGCACTTCTTCCGCCAGCATCATAAATGGGCTGCGTTGTACCAGTTc                                >  1:2499849/1‑69 (MQ=255)
                           tCCCATATGTTCCCGCGCCTCAGGCACTTCTTCCGCCAGCATCATAAATGGGCTGCGTTGTACCAGTTc                                >  1:2380685/1‑69 (MQ=255)
                           tCCCATATGTTCCCGCGCCTCAGGCACTTCTTCCGCCAGCATCATAAATGGGCTGCGTTGTACCAGTTc                                >  1:915751/1‑69 (MQ=255)
                           tCCCATATGTTCCCGCGCCTCAGGCACTTCTTCCGCCAGCATCATAAATGGGCTGCGTTGTACCAGTTc                                >  1:2134333/1‑69 (MQ=255)
                           tCCCATATGTTCCCGCGCCTCAGGCACTTCTTCCGCCAGCATCATAAATGGGCTGCGTTGTACCAGTTc                                >  1:1787286/1‑69 (MQ=255)
                           tCCCATATGTTCCCGCGCCTCAAGCACTTCTTCCGCCAGCATCATAAATGGGCTGCGTTGTACCAGTTc                                >  1:159385/1‑69 (MQ=255)
                           tCCAATATGTTCCCGCGCCTCAGGCACTTCTTCCGCCAGCATCATa                                                       >  1:2241052/1‑46 (MQ=255)
                                 atGTTCCCGCGCCTCAGGTACTTCTTCTGCCAGCATCATAAATGGGCTGCGTTGTACCAGTTCGCTTTc                          >  1:807078/1‑69 (MQ=255)
                                  tGTTCCCGCGCCTCAGGTACTTCTTCTGCCAGCATCATa                                                       <  1:2852036/39‑1 (MQ=255)
                                  tGTTCCCGCGCCTCAGGTACTTCTTCTGCCAGCATCATa                                                       <  1:2210830/39‑1 (MQ=255)
                                  tGTTCCCGCGCCTCAGGTACTTCTTCTGCCAGCATCATa                                                       <  1:778200/39‑1 (MQ=255)
                                  tGTTCCCGCGCCTCAGGTACTTCTTCTGCCAGCATCATa                                                       <  1:260452/39‑1 (MQ=255)
                                  tGTTCCCGCGCCTCAGGTACTTCTTCTGCCAGCATCATa                                                       <  1:1218312/39‑1 (MQ=255)
                                  tGTTCCCGCGCCTCAGGTACTTCTTCTGCCAGCATCATa                                                       <  1:1392191/39‑1 (MQ=255)
                                  tGTTCCCGCGCCTCAGGTACTTCTTCTGCCAGCATCATa                                                       <  1:136979/39‑1 (MQ=255)
                                  tGTTCCCGCGCCTCAGGTACTTCTTCTGCCAGCATCATa                                                       <  1:3313722/39‑1 (MQ=255)
                                  tGTTCCCGCGCCTCAGGTACTTCTTCTGCCAGCATCATAAATGGGCTGCGTTGTACCAGTTCGCTTTcc                         >  1:1114926/1‑69 (MQ=255)
                                  tGTTCCCGCGCCTCAGGTACTTCTTCTGCCAGCATCATAAATGGGCTGCGTTGTACCAGTTCGCTTTcc                         >  1:2854855/1‑69 (MQ=255)
                                  tGTTCCCGCGCCTCAGGTACTTCTTCTGCCAGCATCATAAATGGGCTGCGTTGTACCAGTTCGCTTTcc                         >  1:1648344/1‑69 (MQ=255)
                                  tGTTCCCGCGCCTCAGGTACTTCTTCTGCCAGCATCATAAATGGGCTGCGTTGTACCAGTTCGCTTTcc                         >  1:3257549/1‑69 (MQ=255)
                                  tGTTCCCGCGCCTCAGGTACTTCTTCTGCCAGCATCATAAATGGGCTGCGTTGTACCAGTTCGCTTTcc                         >  1:2067098/1‑69 (MQ=255)
                                  tGTTCCCGCGCCTCAGGTACTTCTTCTGCCAGCATCATAAATGGGCTGCGTTGTACCAGTTCGCTTTcc                         >  1:2374702/1‑69 (MQ=255)
                                  tGTTCCCGCGCCTCAGGTACTTCTTCTGCCAGCATCATAAATGGGCTGCGTTGTACCAGTTCGCTTTcc                         >  1:512162/1‑69 (MQ=255)
                                  tGTTCCCGCGCCTCAGGTACTTCTTCTGCCAGCATCATAAATGGGCTGCGTTGTACCAGTTCGCTTTcc                         >  1:1187551/1‑69 (MQ=255)
                                    ttCCCGCGCCTCAGGTACTTCTTCTGCCAGCATCATAAATGGGCTGCGTTGTACCAGTTCGCTTTCCGt                       <  1:2949235/69‑1 (MQ=255)
                                    ttCCCGCGCCTCAGGCACTTCTTCCGCCAGCATCATAAATGGGCTGCGTTGTACCAGTTCGCTTTCCGt                       <  1:1295046/69‑1 (MQ=255)
                                    ttCCCGCGCCTCAGGCACTTCTTCCGCCAGCATCATAAATGGGCTGCGTTGTACCAGTTCGCTTTCCGt                       <  1:3197372/69‑1 (MQ=255)
                                    ttCCCGCGCCTCAGGCACTTCTTCCGCCAGCATCATAAATGGGCTGCGTTGTACCAGTTCGCTTTCCGt                       <  1:3182236/69‑1 (MQ=255)
                                    ttCCCGCGCCTCAGGCACTTCTTCCGCCAGCATCATAAATGGGCTGCGTTGTACCAGTTCGCTTTCCGt                       <  1:3130913/69‑1 (MQ=255)
                                    ttCCCGCGCCTCAGGCACTTCTTCCGCCAGCATCATAAATGGGCTGCGTTGTACCAGTTCGCTTTCCGt                       <  1:1426683/69‑1 (MQ=255)
                                    ttCCCGCGCCTCAGGCACTTCTTCCGCCAGCATCATAAATGGGCTGCGTTGTACCAGTTCGCTTTCCGt                       <  1:1210120/69‑1 (MQ=255)
                                    ttCCCGCGCCTCAGGCACTTCTTCCGCCAGCATCATAAATGGGCTGCGTTGTACCAGTTCGCTTTCCGt                       <  1:801580/69‑1 (MQ=255)
                                    ttCCCGCGCCTCAGGCACTTCTTCCGCCAGCATCATAAATGGGCTGCGTTGTACCAGTTCGCTTTCCGt                       <  1:2944089/69‑1 (MQ=255)
                                    ttCCCGCGCCTCAGGCACTTCTTCCGCCAGCATCATAAATGGGCTGCGTTGTACCAGTTCGCTTTCCGt                       <  1:2304793/69‑1 (MQ=255)
                                    ttCCCGCGCCTCAGGCACTTCTTCCGCCAGCATCATAAATGGGCTGCGTTGTACCAGTTCGCTTTCCGt                       <  1:1665727/69‑1 (MQ=255)
                                    ttCCCGCGCCTCAGGCACTTCTTCCGCCAGCATCATAAATGGGCTGCGTTGTACCAGTTCGCTTTCCGt                       <  1:1996764/69‑1 (MQ=255)
                                    ttCCCGCGCCTCAGGCACTTCTTCCGCCAGCATCATAAATGGGCTGCGTTGTACCAGTTCGCTTTCCGt                       <  1:2502114/69‑1 (MQ=255)
                                    ttCCCGCGCCTCAGGCACTTCTTCCGCCAGCATCATAAATGGGCTGCGTTGTACAAGTTCGCTTTCCGt                       <  1:2136106/69‑1 (MQ=255)
                                       ccGCGCCTCAGGTACTTCTTCTGCCAGCATCATAAATGGGCTGCGTTGTACCAGTTCGCTTTCCGtt                      <  1:1964907/67‑1 (MQ=255)
                                                  gTACTTCTTCTGCCAGCATCATAAATGGGCTGCGTTGTACCAGTTCGCTTTCCGTTACCCGACGCGcc          <  1:2078769/68‑1 (MQ=255)
                                                        cttctGCCAGCATCATAAATGGGCTGCGTTGTACCAGTTCg                               >  1:1202635/1‑41 (MQ=255)
                                                        cttctGCCAGCATCATAAATGGGCTGCGTTGTACCAGTTCg                               >  1:2264634/1‑41 (MQ=255)
                                                        cttctGCCAGCATCATAAATGGGCTGCGTTGTACCAGTTCg                               >  1:535914/1‑41 (MQ=255)
                                                        cttctGCCAGCATCATAAATGGGCTGCGTTGTACCAGTTCg                               >  1:232768/1‑41 (MQ=255)
                                                        cttctGCCAGCATCATAAATGGGCTGCGTTGTACCAGTTCg                               >  1:1236906/1‑41 (MQ=255)
                                                        cttctGCCAGCATCATAAATGGGCTGCGTTGTACCAGTTCg                               >  1:2202107/1‑41 (MQ=255)
                                                        cttctGCCAGCATCATAAATGGGCTGCGTTGTACCAGTTCg                               >  1:2406352/1‑41 (MQ=255)
                                                        cttctGCCAGCATCATAAATGGGCTGCGTTGTACCAGTTCg                               >  1:2487912/1‑41 (MQ=255)
                                                        cttctGCCAGCATCATAAATGGGCTGCGTTGTACCAGTTCg                               >  1:2142916/1‑41 (MQ=255)
                                                        cttctGCCAGCATCATAAATGGGCTGCGTTGTACCAGTTCg                               >  1:1177133/1‑41 (MQ=255)
                                                        cttctGCCAGCATCATAAATGGGCTGCGTTGTACCAGTTCg                               >  1:3129721/1‑41 (MQ=255)
                                                        cttctGCCAGCATCATAAATGGGCTGCGTTGTACCAGTTCg                               >  1:2789392/1‑41 (MQ=255)
                                                        cttctGCCAGCATCATAAATGGGCTGCGTTGTACCAGTTCGCTTTCCGTTACCCGACGCGCCAGGTAtt   >  1:2171689/1‑69 (MQ=255)
                                                        cttctGCCAGCATCATAAATGGGCTGCGTTGTACCAGTTCGCTTTCCGTTACCCGACGCGCCAGGTAtt   >  1:2402911/1‑69 (MQ=255)
                                                        cttctGCCAGCATCATAAATGGGCTGCGTTGTACCAGTTCGCTTTCCGTTACCCGACGCGCCAGGTAtt   >  1:551143/1‑69 (MQ=255)
                                                        cttctGCCAGCATCATAAATGGGCTGCGTTGTACCAGTTCGCTTTCCGTTACCCGACGCGCCAGGTAtt   >  1:2122005/1‑69 (MQ=255)
                                                        cttctGCCAGCATCATAAATGGGCTGCGTTGTACCAGTTCGCTTTCCGTTACCCGACGCGCCAGGTAtt   >  1:1966461/1‑69 (MQ=255)
                                                        cttctGCCAGCATCATAAATGGGCTGCGTTGTACCAGTTCGCTTTCCGTTACCCGACGCGCCAGGTAtt   >  1:2985724/1‑69 (MQ=255)
                                                        cttctGCCAGCATCATAAATGGGCTGCGTTTACCAGTTCGCTTTCCGTTACCCGAcgcg            >  1:2025838/1‑59 (MQ=255)
                                                        cttctGCCAGCATCATAAATGGGCTGCGTTTACCAGTTCGCTTTCCGTTACCCGACGCGCCAGGTATTc  >  1:2022019/1‑69 (MQ=255)
                                                            |                                                                 
CTGCCGTACCATCGCCAGCCCGTAGCGTCCCATATGTTCCCGCGCCTCAGGTACTTCTTCTGCCAGCATCATAAATGGGCTGCGTTGTACCAGTTCGCTTTCCGTTACCCGACGCGCCAGGTATTC  >  minE/2321325‑2321450

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: