Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 771,789 G→C 73.6% A33G (GCC→GGC)  hflD ← predicted lysogenization regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE771,7890GC73.6% 105.7 / 39.0 72A33G (GCC→GGC) hflDpredicted lysogenization regulator
Reads supporting (aligned to +/- strand):  ref base G (15/4);  new base C (25/28);  total (40/32)
Fisher's exact test for biased strand distribution p-value = 2.98e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

CCGCCAGCGTCGAGCTGGGGTTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCGGCATCACAATGCCCCTGGTGAGCGAGTTGTTGCACCAGGCGTGCCGACTGACAAATA  >  minE/771727‑771845
                                                              |                                                        
ccGCCAGCGTCGAGCTGGGGTTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCGGCATcac                                                    <  1:2301421/69‑1 (MQ=255)
ccGCCAGCGTCGAGCTGGGGTTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCGGCATcac                                                    <  1:620509/69‑1 (MQ=255)
             gCTGGGGTTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCGGc                                                         <  1:2787063/51‑1 (MQ=255)
             gCTGGGGTTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCGGc                                                         <  1:1377802/51‑1 (MQ=255)
                      cATATCGATAATACTGTTGAGTGAGACGTGTAGCGCATCGCCATCACAATGCCCTTGGTGAGc                                    >  1:1736013/1‑63 (MQ=255)
                      cATATCGATAATACTGTTGAGTGAGACGTGTAGCGCATCGCCATCACAATGCCCCt                                           >  1:1749207/1‑56 (MQ=255)
                      cATATCGATAATACTGTTGAGTGAGACGTGTAGCGCATCGCCATCACAATGCCCCTGGTGAg                                     >  1:135050/1‑62 (MQ=255)
                      cATATCGATAATACTGTTGAGTGAGACGTGTAGCGCATCGCCATCACAATGCCCCTGGTGAGc                                    >  1:3107672/1‑63 (MQ=255)
                      cATATCGATAATACTGTTGAGTGAGACGTGTAGCGCATCGCCATCACAATGCCCCTGGTGAGc                                    >  1:234068/1‑63 (MQ=255)
                      cATATCGATAATACTGTTGAGTGAGACGTGTAGCGCATCGCCATCACAATGCCCCTGGTGAGc                                    >  1:2405538/1‑63 (MQ=255)
                      cATATCGATAATACTGTTGAGTGAGACGTGTAGCGCATCGCCATCACAATGCCCCTGGTGAGc                                    >  1:2423222/1‑63 (MQ=255)
                      cATATCGATAATACTGTTGAGTGAGACGTGTAGCGCATCGCCATCACAATGCCCCTGGTGAGc                                    >  1:1002380/1‑63 (MQ=255)
                      cATATCGATAATACTGTTGAGTGAGACGTGTAGCGCATCGCCATCACAATGCCCCTGGTGAGc                                    >  1:2952230/1‑63 (MQ=255)
                      cATATCGATAATACTGTTGAGTGAGACGTGTAGCGCATCGCCATCACAATGCCCCTGGTGAGc                                    >  1:2760514/1‑63 (MQ=255)
                      cATATCGATAATACTGTTGAGTGAGACGTGTAGCGCATCGCCATCACAATGCCCCTGGTGAGc                                    >  1:3274721/1‑63 (MQ=255)
                      cATATCGATAATACTGTTGAGTGAGACGTGTAGCGCATCGCCATCACAATGCCCCTGGTGAGc                                    >  1:3337125/1‑63 (MQ=255)
                      cATATCGATAATACTGTTGAGTGAGACGTGTAGCGCATCGCCATCACAATGCCCCTGGTGAGc                                    >  1:3353593/1‑63 (MQ=255)
                      cATATCGATAATACTGTTGAGTGAGACGTGTAGCGCATCGCCATCACAATGCCCCTGGTGAGc                                    >  1:391869/1‑63 (MQ=255)
                      cATATCGATAATACTGTTGAGTGAGACGTGTAGCGCATCGCCATCACAATGCCCCTGGTGAGc                                    >  1:512751/1‑63 (MQ=255)
                      cATATCGATAATACTGTTGAGTGAGACGTGTAGCGCATCGCCATCACAATGCCCCTGGTGAGc                                    >  1:58239/1‑63 (MQ=255)
                      cATATCGATAATACTGTTGAGTGAGACGTGTAGCGCATCGCCATCACAATGCCCCTGGTGAGc                                    >  1:929439/1‑63 (MQ=255)
                      cATATCGATAATACTGTTGAGTGAGACGTGTAGCGCATCGCCATCACAATGCCCCTGGTGAGc                                    >  1:980777/1‑63 (MQ=255)
                      cATATCGATAATACTGTTGAGTGAGACGTGTAGCGCATCGCCATCACAATGCCCCTGGTGAGc                                    >  1:1876396/1‑63 (MQ=255)
                      cATATCGATAATACTGTTGAGTGAGACGTGTAGCGCATCGCCATCACAATGCCCCTGGTGAGc                                    >  1:1125419/1‑63 (MQ=255)
                      cATATCGATAATACTGTTGAGTGAGACGTGTAGCGCATCGCCATCACAATGCCCCTGGTGAGc                                    >  1:1153984/1‑63 (MQ=255)
                      cATATCGATAATACTGTTGAGTGAGACGTGTAGCGCATCGCCATCACAATGCCCCTGGTGAGc                                    >  1:130145/1‑63 (MQ=255)
                      cATATCGATAATACTGTTGAGTGAGACGTGTAGCGCATCGCCATCACAATGCCCCTGGTGAGc                                    >  1:1307113/1‑63 (MQ=255)
                      cATATCGATAATACTGTTGAGTGAGACGTGTAGCGCATCGCCATCACAATGCCCCTGCTGAGc                                    >  1:1431157/1‑63 (MQ=255)
                      cATATCGATAATACTGTTGAGTGAGACGTGTAGCGCATCGCCATCACAATGCCCATGGTGAGc                                    >  1:3023800/1‑63 (MQ=255)
                           cAATAATACTGTTGAGTGAGACGTGTAGCGCATCGGCATCACAATGCCCCTGGTGAGCGAGttgttg                           >  1:1608655/1‑67 (MQ=255)
                           cAATAATACTGTTGAGTGAGACGTGTAGCGCATCGGCATCACAATGCCCCTGGTGAGCGAGTTGTTGc                          >  1:1059923/1‑68 (MQ=255)
                               aataCTGTTGAGTGAGACGTGTAGCGCATCGCCATCACAATGCCCCTGGTGAGCGAGttgttg                           <  1:1946235/63‑1 (MQ=255)
                               aataCTGTTGAGTGAGACGTGTAGCGCATCGCCATCACAATGCCCCTGGTGAGCGAGttgttg                           <  1:492641/63‑1 (MQ=255)
                               aataCTGTTGAGTGAGACGTGTAGCGCATCGCCATCACAATGCCCCTGGTGAGCGAGttgttg                           <  1:470607/63‑1 (MQ=255)
                               aataCTGTTGAGTGAGACGTGTAGCGCATCGCCATCACAATGCCCCTGGTGAGCGAGttgttg                           <  1:439517/63‑1 (MQ=255)
                               aataCTGTTGAGTGAGACGTGTAGCGCATCGCCATCACAATGCCCCTGGTGAGCGAGttgttg                           <  1:214207/63‑1 (MQ=255)
                               aataCTGTTGAGTGAGACGTGTAGCGCATCGCCATCACAATGCCCCTGGTGAGCGAGttgttg                           <  1:3386159/63‑1 (MQ=255)
                               aataCTGTTGAGTGAGACGTGTAGCGCATCGCCATCACAATGCCCCTGGTGAGCGAGttgttg                           <  1:1361693/63‑1 (MQ=255)
                               aataCTGTTGAGTGAGACGTGTAGCGCATCGCCATCACAATGCCCCTGGTGAGCGAGttgttg                           <  1:337789/63‑1 (MQ=255)
                               aataCTGTTGAGTGAGACGTGTAGCGCATCGCCATCACAATGCCCCTGGTGAGCGAGttgttg                           <  1:220296/63‑1 (MQ=255)
                               aataCTGTTGAGTGAGACGTGTAGCGCATCGCCATCACAATGCCCCTGGTGAGCGAGttgttg                           <  1:2220759/63‑1 (MQ=255)
                               aataCTGTTGAGTGAGACGTGTAGCGCATCGCCATCACAATGCCCCTGGTGAGCGAGttgttg                           <  1:767501/63‑1 (MQ=255)
                               aataCTGTTGAGTGAGACGTGTAGCGCATCGCCATCACAATGCCCCTGGTGAGCGAGttgttg                           <  1:135819/63‑1 (MQ=255)
                               aataCTGTTGAGTGAGACGTGTAGCGCATCGCCATCACAATGCCCCTGGTGAGCGAGttgttg                           <  1:3139306/63‑1 (MQ=255)
                               aataCTGTTGAGTGAGACGTGTAGCGCATCGCCATCACAATGCCCCTGGTGAGCGAGttgttg                           <  1:3136849/63‑1 (MQ=255)
                               aataCTGTTGAGTGAGACGTGTAGCGCATCGCCATCACAATGCCCCTGGTGAGCGAGttgttg                           <  1:769905/63‑1 (MQ=255)
                               aataCTGTTGAGTGAGACGTGTAGCGCATCGCCATCACAATGCCCCTGGTGAGCGAGttgttg                           <  1:3106153/63‑1 (MQ=255)
                               aataCTGTTGAGTGAGACGTGTAGCGCATCGCCATCACAATGCCCCTGGTGAGCGAGttgttg                           <  1:1051518/63‑1 (MQ=255)
                               aataCTGTTGAGTGAGACGTGTAGCGCATCGCCATCACAATGCCCCTGGTGAGCGAGttgttg                           <  1:889123/63‑1 (MQ=255)
                               aataCTGTTGAGTGAGACGTGTAGCGCATCGCCATCACAATGCCCCTGGTGAGCGAGttgttg                           <  1:2497705/63‑1 (MQ=255)
                               aataCTGTTGAGTGAGACGTGTAGCGCATCGCCATCACAATGCCCCTGGTGAGCGAGttgttg                           <  1:2583964/63‑1 (MQ=255)
                               aataCTGTTGAGTGAGACGTGTAGCGCATCGCCATCACAATGCCCCTGGTGAGCGAGttgttg                           <  1:2687503/63‑1 (MQ=255)
                               aataCTGTTGAGTGAGACGTGTAGCGCATCGCCATCACAATGCCCCTGGTGAGCGAGttgttg                           <  1:1611481/63‑1 (MQ=255)
                               aataCTGTTGAGTGAGACGTGTAGCGCATCGCCATCACAATGCCCCTGGTGAGCGAGttgttg                           <  1:2818087/63‑1 (MQ=255)
                               aataCTGTTGAGTGAGACGTGTAGCGCATCGCCATCACAATGCCCCTGGTGAGCGAGttgttg                           <  1:812860/63‑1 (MQ=255)
                               aataCTGTTGAGTGAGACGTGTAGCGCATCGCCATCACAATGCCCCTGGTGAGCGAGttgttg                           <  1:2898983/63‑1 (MQ=255)
                               aataCTGTTGAGTGAGACGTGTAGCGCATCGCCATCACAATGCCCCTGGTGAGCGAGttgttg                           <  1:3065821/63‑1 (MQ=255)
                               aataCTGTTGAGTGAGACGTGTAGCGCATCGCCATCACAATGCCCCTGGTGAGCGAGttgttg                           <  1:3098164/63‑1 (MQ=255)
                                  aCTGTTGAGTGAGACGTGTAGCGCATCGGCATCACAATGCCCCTGGTGAGCGAGTTGTTGCACCAGGc                   >  1:1862840/1‑68 (MQ=255)
                                  aCTGTTGAGTGAGACGTGTAGCGCATCGGCATCACAATGCCCCTGGTGAGCGAGTTGTTGCACCAGGc                   >  1:1851574/1‑68 (MQ=255)
                                  aCTGTTGAGTGAGACGTGTAGCGCATCGGCATCACAATGCCCCTGGTGAGCGAGTTGTTGCACCAGGCg                  >  1:1020040/1‑69 (MQ=255)
                                  aCTGTTGAGTGAGACGTGTAGCGCATCGGCATCACAATGCCCCTGGTGAGCGAGTTGTTGCACCAGGCg                  >  1:1902296/1‑69 (MQ=255)
                                  aCTGTTGAGTGAGACGTGTAGCGCATCGGCATCACAATGCCCCTGGTGAGCGAGTTGTTGCACCAGGCg                  >  1:692321/1‑69 (MQ=255)
                                  aCTGTTGAGTGAGACGTGTAGCGCATCGGCATCACAATGCCCCTGGTGAGCGAGTTGTTGCACCAGGCg                  >  1:611496/1‑69 (MQ=255)
                                  aCTGTTGAGTGAGACGTGTAGCGCATCGGCATCACAATGCCCCTGGTGAGCGAGTTGTTGCACCAGGCg                  >  1:3298613/1‑69 (MQ=255)
                                  aCTGTTGAGTGAGACGTGTAGCGCATCGGCATCACAATGCCCCTGGTGAGCGAGTTGTTGCACCAGGCg                  >  1:2882147/1‑69 (MQ=255)
                                  aCTGTTGAGTGAGACGTGTAGCGCATCGGCATCACAATGCACCTGGTGAGCGAGTTGTTGCACCAGGCg                  >  1:1799775/1‑69 (MQ=255)
                                         agtgagACGTGTAGCGCATCGCCATCACAATGCCCCTGGTGAGCGAGttgttg                           <  1:1081009/53‑1 (MQ=255)
                                           tgagACGTGTAGCGCATCGGCATCACAATGCCCCTGGTGCGCGAGTTGTTGCACCAGGCGTGCCGACt          >  1:504031/1‑68 (MQ=255)
                                           tgagACGTGTAGCGCATCGGCATCACAATGCCCCTGGTGAGCGAGTTGTTGCACCAGGCGTGCCGACt          >  1:312459/1‑68 (MQ=255)
                                           tgagACGTGTAGCGCATCGGCATCACAATGCCCCTGGTGAGCGAGTTGTTGCACCAGGCGTGCCGACt          >  1:2238260/1‑68 (MQ=255)
                                                  tgtAGCGCATCGGCATCACAATGCCCCTGGTGAGCGAGTTGTTGCACCAGGCGTGCCGACTGACAAATa  >  1:1650378/1‑69 (MQ=255)
                                                              |                                                        
CCGCCAGCGTCGAGCTGGGGTTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCGGCATCACAATGCCCCTGGTGAGCGAGTTGTTGCACCAGGCGTGCCGACTGACAAATA  >  minE/771727‑771845

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: