Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,778,836 T→C 47.9% intergenic (‑10/‑190) iclR ← / → metH DNA‑binding transcriptional repressor/homocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,778,8360TC47.9% 0.8 / 87.0 96intergenic (‑10/‑190)iclR/metHDNA‑binding transcriptional repressor/homocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent
Reads supporting (aligned to +/- strand):  ref base T (39/11);  new base C (31/15);  total (70/26)
Fisher's exact test for biased strand distribution p-value = 2.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

Only 100 of 122 total aligned reads displayed.

CGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTG  >  minE/2778775‑2778900
                                                             |                                                                
cGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTc                                                            >  1:172724/1‑68 (MQ=255)
cGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTc                                                            >  1:3388658/1‑68 (MQ=255)
    ggCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTa                                                        >  1:3266926/1‑68 (MQ=255)
    ggCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTAt                                                       >  1:3102275/1‑69 (MQ=255)
        aCGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGt                                                    >  1:1161893/1‑68 (MQ=255)
            cggGTTTTCTGCCGCGTTTTGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGaaa                                               <  1:2214773/69‑1 (MQ=255)
            cggGTTTTCTGCCGCGTTTTGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGaaa                                               <  1:2946162/69‑1 (MQ=255)
            cggGTTTTCTGCCGCGTTTTGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGaaa                                               <  1:3029663/69‑1 (MQ=255)
            cggGTTTTCTGCCGCGTTTTGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGaaa                                               <  1:1478114/69‑1 (MQ=255)
            cggGTTTTCTGCCGCGTTTTGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGaaa                                               <  1:2889504/69‑1 (MQ=255)
                ttttCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAAtcat                                           >  1:3402905/1‑69 (MQ=255)
                ttttCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAAtcat                                           >  1:3355568/1‑69 (MQ=255)
                ttttCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAAtcat                                           >  1:739782/1‑69 (MQ=255)
                ttttCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAAtcat                                           >  1:2727181/1‑69 (MQ=255)
                ttttCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAAtcat                                           >  1:1859999/1‑69 (MQ=255)
                         gcgTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATcattttcatttt                                  >  1:3379854/1‑69 (MQ=255)
                         gcgTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATcattttcatttt                                  <  1:3298451/69‑1 (MQ=255)
                         gcgTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATcatttacatttt                                  >  1:78664/1‑69 (MQ=255)
                                    gAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCt                       <  1:1600223/69‑1 (MQ=255)
                                         ggTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTaa                     <  1:926790/66‑1 (MQ=255)
                                         ggTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTaa                     <  1:516005/66‑1 (MQ=255)
                                         ggTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGc                  >  1:2570066/1‑69 (MQ=255)
                                           tGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCaa                <  1:2484115/69‑1 (MQ=255)
                                             cGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAg                         >  1:786392/1‑58 (MQ=255)
                                             cGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAg                         >  1:805351/1‑58 (MQ=255)
                                             cGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAg                         >  1:931254/1‑58 (MQ=255)
                                             cGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAg                         >  1:2090928/1‑58 (MQ=255)
                                               aCCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGt            <  1:388947/69‑1 (MQ=255)
                                               aCCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGt            <  1:491336/69‑1 (MQ=255)
                                               aCCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGt            <  1:1105203/69‑1 (MQ=255)
                                               aCCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGt            <  1:3388204/69‑1 (MQ=255)
                                                ccATGACAGTCTCCTTTTTCTGTATCGTGGAAAtcat                                           >  1:1488152/1‑37 (MQ=255)
                                                ccATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTAtta                             >  1:2941357/1‑50 (MQ=255)
                                                ccATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTAtta                             >  1:3410102/1‑50 (MQ=255)
                                                ccATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTAtta                             >  1:2606724/1‑50 (MQ=255)
                                                ccATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTAtta                             >  1:2466433/1‑50 (MQ=255)
                                                ccATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTAtta                             >  1:3227991/1‑50 (MQ=255)
                                                ccATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTAtta                             >  1:2176657/1‑50 (MQ=255)
                                                ccATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTAtta                             >  1:1980637/1‑50 (MQ=255)
                                                ccATGAAAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTAtta                             >  1:438776/1‑50 (MQ=255)
                                                 cATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATg                   >  1:367640/1‑60 (MQ=255)
                                                 cATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAAt               >  1:529212/1‑64 (MQ=255)
                                                 cATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAAt               >  1:3183215/1‑64 (MQ=255)
                                                 cATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAAt               >  1:1338195/1‑64 (MQ=255)
                                                 cATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAAt               >  1:876880/1‑64 (MQ=255)
                                                 cATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAAt               >  1:78372/1‑64 (MQ=255)
                                                 cATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAAt               >  1:1593337/1‑64 (MQ=255)
                                                 cATGACAGTCTCCTTTTTCTGTATCGTGGAAAtcattttc                                       >  1:1464950/1‑40 (MQ=255)
                                                 cATGACAGTCTCCTTTTTCTGTATCGTGGAAAtcattttc                                       >  1:1643499/1‑40 (MQ=255)
                                                 cATGACAGTCTCCTTTTTCTGTATCGTGGAAAtcattttc                                       >  1:1276280/1‑40 (MQ=255)
                                                 cATGACAGTCTCCTTTTTCTGTATCGTGGAAAtcattttc                                       >  1:1459420/1‑40 (MQ=255)
                                                 cATGACAGTCTCCTTTTTCTGTATCGTGGAAAtcattttc                                       >  1:3133400/1‑40 (MQ=255)
                                                 cATGACAGTCTCCTTTTTCTGTATCGTGGAAAtcattttc                                       >  1:1354082/1‑40 (MQ=255)
                                                 cATGACAGTCTCCTTTTTCTGTATCGTGGAAAtcattttc                                       >  1:1732075/1‑40 (MQ=255)
                                                 cATGACAGTCTCCTTTTTCTGTATCGTGGAAAtcattttc                                       >  1:1117639/1‑40 (MQ=255)
                                                 cATGACAGTCTCCTTTTTCTGTATCGTGGAAAtcattttc                                       >  1:2575579/1‑40 (MQ=255)
                                                 cATGACAGTCTCCTTTTTCTGTATCGTGGAAAtcattttc                                       >  1:2542438/1‑40 (MQ=255)
                                                 cATGACAGTCTCCTTTTTCTGTATCGTGGAAAtcattttc                                       >  1:923748/1‑40 (MQ=255)
                                                 cATGACAGTCTCCTTTTTCTGTATCGTGGAAAtcattttc                                       >  1:2463621/1‑40 (MQ=255)
                                                 cATGACAGTCTCCTTTTTCTGTATCGTGGAAAtcattttc                                       >  1:2352309/1‑40 (MQ=255)
                                                 cATGACAGTCTCCTTTTTCTGTATCGTGGAAAtcatttta                                       >  1:3032190/1‑39 (MQ=39)
                                                     aCAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTactga      >  1:2963506/1‑69 (MQ=255)
                                                     aCAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTactga      >  1:335614/1‑69 (MQ=255)
                                                     aCAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTactga      >  1:3281778/1‑69 (MQ=255)
                                                     aCAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTactga      >  1:1973380/1‑69 (MQ=255)
                                                     aCAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTactga      >  1:685180/1‑69 (MQ=255)
                                                     aCAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTactga      >  1:491728/1‑69 (MQ=255)
                                                     aCAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTactga      >  1:362854/1‑69 (MQ=255)
                                                        gTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTATCAGCGAATGCAATAGTTgctgaac  <  1:734388/69‑2 (MQ=255)
                                                        gTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTATCAGCGAATGCAATAGTTgctgaac  <  1:686187/69‑2 (MQ=255)
                                                        gTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTATCAGCGAATGCAATAGTTgctgaac  <  1:3128876/69‑2 (MQ=255)
                                                        gTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTATCAGCGAATGCAATAGTTgctgaac  <  1:3330516/69‑2 (MQ=255)
                                                        gTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTATCAGCGAATGCAATAGTTgctgaac  <  1:2267749/69‑2 (MQ=255)
                                                        gTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTATCAGCGAATGCAATAGTTgctgaac  <  1:1902011/69‑2 (MQ=255)
                                                        gTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTATCAGCGAATGCAATAGTTgctgaac  <  1:3160179/69‑2 (MQ=255)
                                                        gTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTATCAGCGAATGCAATAGTTgctgaac  <  1:2926381/69‑2 (MQ=255)
                                                        gTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTATCAGCGAATGCAATAGTTgctgaac  <  1:1627075/69‑2 (MQ=255)
                                                          ctctTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGt            >  1:1766261/1‑58 (MQ=255)
                                                          ctctTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGt            >  1:2698569/1‑58 (MQ=255)
                                                          ctctTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGt            >  1:1350939/1‑58 (MQ=255)
                                                          ctctTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTactgaactg  <  1:3345774/68‑1 (MQ=255)
                                                          ctctTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTactgaactg  <  1:270572/68‑1 (MQ=255)
                                                          ctctTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTactgaactg  <  1:266785/68‑1 (MQ=255)
                                                          ctctTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTactgaactg  <  1:260735/68‑1 (MQ=255)
                                                          ctctTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTactgaactg  <  1:91594/68‑1 (MQ=255)
                                                          ctctTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTactgaactg  <  1:3194883/68‑1 (MQ=255)
                                                          ctctTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTactgaactg  <  1:2003215/68‑1 (MQ=255)
                                                          ctctTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTactga      <  1:2080812/64‑1 (MQ=255)
                                                          ctctTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTactga      <  1:1320108/64‑1 (MQ=255)
                                                          ctctTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTactga      <  1:359990/64‑1 (MQ=255)
                                                          ctctTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTactga      <  1:3280720/64‑1 (MQ=255)
                                                          ctctTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTactga      <  1:2613283/64‑1 (MQ=255)
                                                          ctctTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTactga      <  1:2489297/64‑1 (MQ=255)
                                                          ctctTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTactga      <  1:2253817/64‑1 (MQ=255)
                                                             ttttttCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAAt               >  1:3146506/1‑52 (MQ=255)
                                                             ttttttCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAAt               >  1:1151318/1‑52 (MQ=255)
                                                             ttttttCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATa              >  1:1623500/1‑53 (MQ=255)
                                                             ttttttCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATa              >  1:1798204/1‑53 (MQ=255)
                                                             ttttttCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATa              >  1:2545297/1‑53 (MQ=255)
                                                             ttttttCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATa              >  1:1035818/1‑53 (MQ=255)
                                                             |                                                                
CGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTG  >  minE/2778775‑2778900

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: