Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,310,943 C→G 61.8% A158A (GCC→GCG hisG → ATP phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,310,9430CG61.8% 50.8 / 68.3 76A158A (GCC→GCGhisGATP phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (19/10);  new base G (16/31);  total (35/41)
Fisher's exact test for biased strand distribution p-value = 9.67e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.84e-01

CGTTATCTCGACCAGAAAGGCATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGAAGTCGCCCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTGGTTTCCACCGGTGCCAC  >  minE/1310878‑1310996
                                                                 |                                                     
cgatATCTCGACCAGAAAGGCATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGAAGTCGcccc                                                     >  1:102532/4‑68 (MQ=255)
cGTTATCTCGACCAGAAAGGCATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGAAGTCGCCCcg                                                    >  1:1017028/1‑69 (MQ=255)
cGTTATCTCGACCAGAAAGGCATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGAAGTCGCCCcg                                                    >  1:3044762/1‑69 (MQ=255)
cGTTATCTCGACCAGAAAGGCATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGAAGTCGCCCcg                                                    >  1:2299798/1‑69 (MQ=255)
cGTTATCTCGACCAGAAAGGCATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGAAGTCGCCCcg                                                    >  1:3138780/1‑69 (MQ=255)
cGTTATCTCGACCAGAAAGGCATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGAAGTCGCCCcg                                                    >  1:3186077/1‑69 (MQ=255)
cGTTATCTCGACCAGAAAGGCATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGAAGTCGCCCcg                                                    >  1:1219783/1‑69 (MQ=255)
cGTTATCTCGACCAGAAAGGCATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGAAGTCGCCCcg                                                    >  1:1244647/1‑69 (MQ=255)
cGTTATCTCGACCAGAAAGGCATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGAAGTCGCCCcg                                                    >  1:2173536/1‑69 (MQ=255)
 gTTATCTCGACCAGAAAGGCATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGAAGTCGCCCcgc                                                   >  1:3380619/1‑69 (MQ=255)
 gTTATCTCGACCAGAAAGGCATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGAAGTCGCCCcgc                                                   >  1:1714408/1‑69 (MQ=255)
    aTCTCGACCAGAAAGGCATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGAAGTCGCCCCGCGTg                                                >  1:2753115/1‑69 (MQ=255)
    aTCTCGACCAGAAAGGCATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGAAGTCGCCCCGCGTg                                                >  1:2806068/1‑69 (MQ=255)
    aTCTCGACCAGAAAGGCATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGAAGTCGCCCCGCGTg                                                >  1:386483/1‑69 (MQ=255)
           ccAGAAAGGCATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGAAGTCGCCCCGCGTGCCgg                                            <  1:1492785/66‑1 (MQ=255)
           ccAGAAAGGCATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGAAGTCGCCCCGCGTGCCgg                                            <  1:1907997/66‑1 (MQ=255)
           ccAGAAAGGCATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGAAGTCGCCCCGCGTGCCgg                                            <  1:2064590/66‑1 (MQ=255)
           ccAGAAAGGCATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGAAGTCGCCCCGCGTGCCgg                                            <  1:2253964/66‑1 (MQ=255)
           ccAGAAAGGCATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGAAGTCGCCCCGCGTGCCgg                                            <  1:3318877/66‑1 (MQ=255)
           ccAGAAAGGCATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGAAGTCGCCCCGCGTGCCgg                                            <  1:1373576/66‑1 (MQ=255)
                        tctTTTAAATCCTGCTTACTGAACGGTTCTGTTGAAGTCGCCCCGCGTGCCGGACTGGCGGATGCGAtt                            >  1:3269528/1‑69 (MQ=255)
                        tctTTTAAATCCTGCTTACTGAACGGTTCTGTTGAAGTCGCCCCGCGTGCCGGACTGGCGGATGCGAtt                            >  1:2363788/1‑69 (MQ=255)
                        tctTTTAAATCCTGCTTACTGAACGGTTCTGTTGAAGTCGCCCCGCGTGCCGGACTGGCGGATGCGAtt                            >  1:693547/1‑69 (MQ=255)
                        tctTTTAAATCCTGCTTAATGAACGGTTCTGTTGAAGTCGCCCCGCGTGCCGGACTGGCGGATGCGAtt                            >  1:280881/1‑69 (MQ=255)
                        tctTTTAAATCATGCTTACTGAACGGTTCTGTTGAAGTCGCCCCGCGTg                                                >  1:768476/1‑49 (MQ=255)
                            ttAAATCCTGCTTACTGAACGGTTCTGTTGAAGTCGCCCCGCGTGCCGGACTGGCGGATGCGATTTGCg                        >  1:2176903/1‑69 (MQ=255)
                               aaTCCTGCTTACTGAACGGTTCTGTTGAAGTCGCCCCGCGTGCCGGACTGGCGGATGCGATTTGCg                        <  1:2166224/66‑1 (MQ=255)
                                           tGAACGGTTCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTgg                  <  1:539789/60‑1 (MQ=255)
                                           tGAACGGTTCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTgg                  <  1:397844/60‑1 (MQ=255)
                                           tGAACGGTTCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTgg                  <  1:2949472/60‑1 (MQ=255)
                                           tGAACGGTTCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTgg                  <  1:3005525/60‑1 (MQ=255)
                                           tGAACGGTTCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTgg                  <  1:353669/60‑1 (MQ=255)
                                           tGAACGGTTCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTgg                  <  1:3006683/60‑1 (MQ=255)
                                           tGAACGGTTCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTgg                  <  1:95174/60‑1 (MQ=255)
                                           tGAACGGTTCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTgg                  <  1:3103589/60‑1 (MQ=255)
                                           tGAACGGTTCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTgg                  <  1:3404098/60‑1 (MQ=255)
                                           tGAACGGTTCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTgg                  <  1:656858/60‑1 (MQ=255)
                                           tGAACGGTTCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTgg                  <  1:978386/60‑1 (MQ=255)
                                           tGAACGGTTCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTgg                  <  1:3245139/60‑1 (MQ=255)
                                           tGAACGGTTCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTgg                  <  1:573940/60‑1 (MQ=255)
                                           tGAACGGTTCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTgg                  <  1:3296532/60‑1 (MQ=255)
                                           tGAACGGTTCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTgg                  <  1:2591655/60‑1 (MQ=255)
                                           tGAACGGTTCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTgg                  <  1:1644253/60‑1 (MQ=255)
                                           tGAACGGTTCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTgg                  <  1:1764379/60‑1 (MQ=255)
                                           tGAACGGTTCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTgg                  <  1:1822490/60‑1 (MQ=255)
                                           tGAACGGTTCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTgg                  <  1:1194841/60‑1 (MQ=255)
                                           tGAACGGTTCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTgg                  <  1:1100018/60‑1 (MQ=255)
                                           tGAACGGTTCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTgg                  <  1:2201315/60‑1 (MQ=255)
                                           tGAACGGTTCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTgg                  <  1:2248136/60‑1 (MQ=255)
                                           tGAACGGTTCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTgg                  <  1:2248473/60‑1 (MQ=255)
                                           tGAACGGTTCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTgg                  <  1:2260103/60‑1 (MQ=255)
                                           tGAACGGTTCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTgg                  <  1:2266845/60‑1 (MQ=255)
                                           tGAACGGTTCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTgg                  <  1:2352329/60‑1 (MQ=255)
                                           tGAACGGTTCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTgg                  <  1:1099197/60‑1 (MQ=255)
                                           tGAACGGTTCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTgg                  <  1:2352988/60‑1 (MQ=255)
                                           tGAACGGTTCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTgg                  <  1:2471014/60‑1 (MQ=255)
                                           tGAACGGTTCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTgg                  <  1:2504817/60‑1 (MQ=255)
                                           tGAACGGTTCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTgg                  <  1:2584961/60‑1 (MQ=255)
                                               cGGTTCTGTTGAAGTCGCCCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTGGTTTCCACCGGTGc     <  1:2669240/69‑1 (MQ=255)
                                               cGGTTCTGTTGAAGTCGCCCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTGGTTTCCACCGGTGc     <  1:1984866/69‑1 (MQ=255)
                                                  ttCTGTCTAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTGGTTTCCAcc          >  1:3217814/1‑61 (MQ=255)
                                                  ttCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTGGTTTCCACCGGTGCCAc  >  1:791238/1‑69 (MQ=255)
                                                  ttCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTGGTTTCCACCGGTGCCAc  >  1:823209/1‑69 (MQ=255)
                                                  ttCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTGGTTTCCACCGGTGCCAc  >  1:1126395/1‑69 (MQ=255)
                                                  ttCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTGGTTTCCACCGGTGCCAc  >  1:639133/1‑69 (MQ=255)
                                                  ttCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTGGTTTCCACCGGTGCCAc  >  1:1643613/1‑69 (MQ=255)
                                                  ttCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTGGTTTCCACCGGTGCCAc  >  1:1068231/1‑69 (MQ=255)
                                                  ttCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTGGTTTCCACCGGTGCCAc  >  1:520168/1‑69 (MQ=255)
                                                  ttCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTGGTTTCCACCGGTGCCAc  >  1:2048130/1‑69 (MQ=255)
                                                  ttCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTGGTTTCCACCGGTGCCAc  >  1:2076221/1‑69 (MQ=255)
                                                  ttCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTGGTTTCCACCGGTGCCAc  >  1:214077/1‑69 (MQ=255)
                                                  ttCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTGGTTTCCACCGGTGCCAc  >  1:3163798/1‑69 (MQ=255)
                                                  ttCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTGGTTTCCACCGGTGCCAc  >  1:2337303/1‑69 (MQ=255)
                                                  ttCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTGGTTTCCACCGGTGCCAc  >  1:3015913/1‑69 (MQ=255)
                                                  ttCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTGGTTTCCACCGGTGCCAc  >  1:2391063/1‑69 (MQ=255)
                                                  ttCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTGGTTTCAACCGGTGCCAc  >  1:769078/1‑69 (MQ=255)
                                                             tCGCCCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTGGTTTCCACCGGTGc     <  1:3330803/55‑1 (MQ=255)
                                                                 |                                                     
CGTTATCTCGACCAGAAAGGCATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGAAGTCGCCCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTGGTTTCCACCGGTGCCAC  >  minE/1310878‑1310996

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: