Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,310,934 T→C 52.3% V155V (GTT→GTC hisG → ATP phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,310,9340TC52.3% 7.0 / 77.5 90V155V (GTT→GTChisGATP phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (20/23);  new base C (16/31);  total (36/54)
Fisher's exact test for biased strand distribution p-value = 2.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.24e-01

CTGCTCAAGCGTTATCTCGACCAGAAAGGCATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGAAGTCGCCCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTGGTTTCCACCGGTGCCAC  >  minE/1310869‑1310996
                                                                 |                                                              
cTGCTCTAGCGTTATCTCGACCAGAAAGGCATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGaa                                                             <  1:891934/69‑1 (MQ=255)
cTGCTCAAGCGTTATCTCGACCAGAAAGGCATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGaa                                                             <  1:1059208/69‑1 (MQ=255)
cTGCTCAAGCGTTATCTCGACCAGAAAGGCATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGaa                                                             <  1:2826972/69‑1 (MQ=255)
cTGCTCAAGCGTTATCTCGACCAGAAAGGCATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGaa                                                             <  1:281249/69‑1 (MQ=255)
cTGCTCAAGCGTTATCTCGACCAGAAAGGCATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGaa                                                             <  1:1101443/69‑1 (MQ=255)
cTGCTCAAGCGTTATCTCGACCAGAAAGGCATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGaa                                                             <  1:3082300/69‑1 (MQ=255)
cTGCTCAAGCGTTATCTCGACCAGAAAGGCATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGaa                                                             <  1:3125834/69‑1 (MQ=255)
cTGCTCAAGCGTTATCTCGACCAGAAAGGCATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGaa                                                             <  1:423320/69‑1 (MQ=255)
cTGCTCAAGCGTTATCTCGACCAGAAAGGCATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGaa                                                             <  1:1725404/69‑1 (MQ=255)
cTGCTCAAGCGTTATCTCGACCAGAAAGGCATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGaa                                                             <  1:1451825/69‑1 (MQ=255)
cTGCTCAAGCGTTATCTCGACCAGAAAGGCATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGaa                                                             <  1:1568729/69‑1 (MQ=255)
 tGCTCAAGCGTTATCTCGACCAGAAAGGCATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGaa                                                             <  1:3005853/68‑1 (MQ=255)
     cAAGCGTTATCTCGACCAGAAAGGCATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGaa                                                             <  1:78901/64‑1 (MQ=255)
         cgatATCTCGACCAGAAAGGCATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGAAGTCGcccc                                                     >  1:102532/4‑68 (MQ=255)
         cGTTATCTCGACCAGAAAGGCATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGAAGTCGCCCcg                                                    >  1:3186077/1‑69 (MQ=255)
         cGTTATCTCGACCAGAAAGGCATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGAAGTCGCCCcg                                                    >  1:2299798/1‑69 (MQ=255)
         cGTTATCTCGACCAGAAAGGCATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGAAGTCGCCCcg                                                    >  1:1017028/1‑69 (MQ=255)
         cGTTATCTCGACCAGAAAGGCATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGAAGTCGCCCcg                                                    >  1:3138780/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)
         cGTTATCTCGACCAGAAAGGCATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGAAGTCGCCCcg                                                    >  1:3044762/1‑69 (MQ=255)
          gTTATCTCGACCAGAAAGGCATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGAAGTCGCCCcgc                                                   >  1:1714408/1‑69 (MQ=255)
          gTTATCTCGACCAGAAAGGCATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGAAGTCGCCCcgc                                                   >  1:3380619/1‑69 (MQ=255)
           ttATCTCGACCAGAAAGGCATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGat                                                             <  1:847020/58‑2 (MQ=255)
             aTCTCGACCAGAAAGGCATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGAAGTCGCCCCGCGTg                                                >  1:386483/1‑69 (MQ=255)
             aTCTCGACCAGAAAGGCATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGAAGTCGCCCCGCGTg                                                >  1:2753115/1‑69 (MQ=255)
             aTCTCGACCAGAAAGGCATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGAAGTCGCCCCGCGTg                                                >  1:2806068/1‑69 (MQ=255)
                    ccAGAAAGGCATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGAAGTCGCCCCGCGTGCCgg                                            <  1:2064590/66‑1 (MQ=255)
                    ccAGAAAGGCATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGAAGTCGCCCCGCGTGCCgg                                            <  1:1907997/66‑1 (MQ=255)
                    ccAGAAAGGCATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGAAGTCGCCCCGCGTGCCgg                                            <  1:1373576/66‑1 (MQ=255)
                    ccAGAAAGGCATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGAAGTCGCCCCGCGTGCCgg                                            <  1:2253964/66‑1 (MQ=255)
                    ccAGAAAGGCATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGAAGTCGCCCCGCGTGCCgg                                            <  1:3318877/66‑1 (MQ=255)
                    ccAGAAAGGCATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGAAGTCGCCCCGCGTGCCgg                                            <  1:1492785/66‑1 (MQ=255)
                                 tctTTTAAATCCTGCTTACTGAACGGTTCTGTTGAAGTCGCCCCGCGTGCCGGACTGGCGGATGCGAtt                            >  1:2363788/1‑69 (MQ=255)
                                 tctTTTAAATCCTGCTTACTGAACGGTTCTGTTGAAGTCGCCCCGCGTGCCGGACTGGCGGATGCGAtt                            >  1:693547/1‑69 (MQ=255)
                                 tctTTTAAATCCTGCTTACTGAACGGTTCTGTTGAAGTCGCCCCGCGTGCCGGACTGGCGGATGCGAtt                            >  1:3269528/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:3245139/60‑1 (MQ=255)
                                                    tGAACGGTTCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTgg                  <  1:3296532/60‑1 (MQ=255)
                                                    tGAACGGTTCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTgg                  <  1:3404098/60‑1 (MQ=255)
                                                    tGAACGGTTCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTgg                  <  1:353669/60‑1 (MQ=255)
                                                    tGAACGGTTCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTgg                  <  1:3103589/60‑1 (MQ=255)
                                                    tGAACGGTTCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTgg                  <  1:3006683/60‑1 (MQ=255)
                                                    tGAACGGTTCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTgg                  <  1:397844/60‑1 (MQ=255)
                                                    tGAACGGTTCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTgg                  <  1:539789/60‑1 (MQ=255)
                                                    tGAACGGTTCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTgg                  <  1:573940/60‑1 (MQ=255)
                                                    tGAACGGTTCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTgg                  <  1:95174/60‑1 (MQ=255)
                                                    tGAACGGTTCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTgg                  <  1:656858/60‑1 (MQ=255)
                                                    tGAACGGTTCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTgg                  <  1:978386/60‑1 (MQ=255)
                                                    tGAACGGTTCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTgg                  <  1:2949472/60‑1 (MQ=255)
                                                    tGAACGGTTCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTgg                  <  1:1099197/60‑1 (MQ=255)
                                                    tGAACGGTTCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTgg                  <  1:1100018/60‑1 (MQ=255)
                                                    tGAACGGTTCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTgg                  <  1:1194841/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: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:2352988/60‑1 (MQ=255)
                                                    tGAACGGTTCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTgg                  <  1:3005525/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)
                                                    tGAACGGTTCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTgg                  <  1:2591655/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:1068231/1‑69 (MQ=255)
                                                           ttCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTGGTTTCCACCGGTGCCAc  >  1:1126395/1‑69 (MQ=255)
                                                           ttCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTGGTTTCCACCGGTGCCAc  >  1:823209/1‑69 (MQ=255)
                                                           ttCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTGGTTTCCACCGGTGCCAc  >  1:791238/1‑69 (MQ=255)
                                                           ttCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTGGTTTCCACCGGTGCCAc  >  1:3015913/1‑69 (MQ=255)
                                                           ttCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTGGTTTCCACCGGTGCCAc  >  1:1643613/1‑69 (MQ=255)
                                                           ttCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTGGTTTCCACCGGTGCCAc  >  1:639133/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:2391063/1‑69 (MQ=255)
                                                           ttCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTGGTTTCCACCGGTGCCAc  >  1:3163798/1‑69 (MQ=255)
                                                           ttCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTGGTTTCCACCGGTGCCAc  >  1:2337303/1‑69 (MQ=255)
                                                           ttCTGTCGAAGTTGCGCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTGGTTTCAACCGGTGCCAc  >  1:769078/1‑69 (MQ=255)
                                                                 |                                                              
CTGCTCAAGCGTTATCTCGACCAGAAAGGCATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGAAGTCGCCCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTGGTTTCCACCGGTGCCAC  >  minE/1310869‑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: