Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 452,550 G→A 34.3% R172H (CGC→CAC)  phr → deoxyribodipyrimidine photolyase, FAD‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE452,5500GA34.3% 46.8 / 33.3 64R172H (CGC→CAC) phrdeoxyribodipyrimidine photolyase, FAD‑binding
Reads supporting (aligned to +/- strand):  ref base G (22/20);  new base A (15/7);  total (37/27)
Fisher's exact test for biased strand distribution p-value = 2.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAG  >  minE/452483‑452616
                                                                   |                                                                  
tCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTgcg                                                                    >  1:1129515/1‑68 (MQ=255)
tCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTgcg                                                                    >  1:1145737/1‑68 (MQ=255)
tCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTgcg                                                                    >  1:945079/1‑68 (MQ=255)
tCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTgcg                                                                    >  1:701079/1‑68 (MQ=255)
tCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTgcg                                                                    >  1:1204091/1‑68 (MQ=255)
tCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTgcg                                                                    >  1:415312/1‑68 (MQ=255)
tCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTgcg                                                                    >  1:2187536/1‑68 (MQ=255)
tCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTgcg                                                                    >  1:2178370/1‑68 (MQ=255)
tCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTgcg                                                                    >  1:1928257/1‑68 (MQ=255)
tCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGACGCCGGAGTGCGTCGCTgcg                                                                    >  1:998948/1‑68 (MQ=255)
  tttACGCCTTTTAAGAATGCCTGGCTGAAACGGGTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCGcc                                                                  >  1:2938102/1‑68 (MQ=255)
  tttACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCGcc                                                                  >  1:1123235/1‑68 (MQ=255)
  tttACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCGcc                                                                  >  1:490187/1‑68 (MQ=255)
  tttACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCGcc                                                                  >  1:445826/1‑68 (MQ=255)
  tttACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCGcc                                                                  >  1:337611/1‑68 (MQ=255)
  tttACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCGcc                                                                  >  1:134946/1‑68 (MQ=255)
       gCCTTTTAAAAACGCATGGTTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCACCAAAAg                                                             >  1:2717854/1‑68 (MQ=255)
       gCCTTTTAAAAACGCATGGTTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCACCAAAAGt                                                            >  1:181935/1‑69 (MQ=255)
       gCCTTTTAAAAACGCATGGTTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCACCAAAAGt                                                            >  1:2603392/1‑69 (MQ=255)
       gCCTTTTAAAAACGCATGGTTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCACCAAAAGt                                                            >  1:2605418/1‑69 (MQ=255)
       gCCTTTTAAAAACGCATGGTTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCACCAAAAGt                                                            >  1:265014/1‑69 (MQ=255)
       gCCTTTTAAAAACGCATGGTTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCACCAAAAGt                                                            >  1:447106/1‑69 (MQ=255)
       gCCTTTTAAAAACGCATGGTTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCACCAAAAGt                                                            >  1:857962/1‑69 (MQ=255)
       gCCTTTTAAAAACGCATGGTTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCACCAAAAGt                                                            >  1:2808036/1‑69 (MQ=255)
       gCCTTTTAAAAACGCATGGTTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCACCAAAAGt                                                            >  1:850829/1‑69 (MQ=255)
       gCCTTTTAAAAACGCATGGTTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCACCAAAAGt                                                            >  1:524756/1‑69 (MQ=255)
       gCCTTTTAAAAACGCATGGTTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCACCAAAAGt                                                            >  1:3177408/1‑69 (MQ=255)
       gCCTTTTAAAAACGCATGGTTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCACCAAAAGt                                                            >  1:1557723/1‑69 (MQ=255)
       gCCTTTTAAAAACGCATGGTTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCACCAAAAGt                                                            >  1:985546/1‑69 (MQ=255)
       gCCTTTTAAAAACGCATGGTTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCACCAAAAGt                                                            >  1:468708/1‑69 (MQ=255)
       gCCTTTTAAAAACGCATGGTTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCACCAAAAGt                                                            >  1:489261/1‑69 (MQ=255)
                    gCCTGGCTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGa                                               <  1:2918476/69‑1 (MQ=255)
                       tGGCTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATc                                             <  1:60730/68‑1 (MQ=255)
                       tGGCTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATc                                             <  1:329972/68‑1 (MQ=255)
                       tGGCTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATc                                             <  1:608943/68‑1 (MQ=255)
                       tGGCTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATc                                             <  1:215564/68‑1 (MQ=255)
                       tGGCTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATc                                             <  1:2194774/68‑1 (MQ=255)
                       tGGCTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATc                                             <  1:881025/68‑1 (MQ=255)
                        ggCTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGa                                           <  1:1303224/69‑1 (MQ=255)
                        ggCTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGa                                           <  1:94977/69‑1 (MQ=255)
                        ggCTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGa                                           <  1:1185179/69‑1 (MQ=255)
                        ggCTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGa                                           <  1:1171606/69‑1 (MQ=255)
                        ggCTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGa                                           <  1:2296814/69‑1 (MQ=255)
                        ggCTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGa                                           <  1:24416/69‑1 (MQ=255)
                        ggCTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGa                                           <  1:2671803/69‑1 (MQ=255)
                        ggCTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGa                                           <  1:415627/69‑1 (MQ=255)
                        ggCTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGa                                           <  1:2560149/69‑1 (MQ=255)
                        ggCTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGAGGATCGa                                           <  1:2097339/69‑1 (MQ=255)
                            gAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATaga                                       <  1:870012/69‑1 (MQ=255)
                            gAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATaga                                       <  1:1468952/69‑1 (MQ=255)
                            gAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATaga                                       <  1:2519815/69‑1 (MQ=255)
                                              gATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGccatccat                        >  1:521613/1‑66 (MQ=255)
                                              gATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGccatccat                        >  1:612876/1‑66 (MQ=255)
                                              gATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGccatccat                        >  1:2737384/1‑66 (MQ=255)
                                              gATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGccatccat                        >  1:2660126/1‑66 (MQ=255)
                                              gATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGccatccat                        >  1:2542580/1‑66 (MQ=255)
                                              gATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGccatccat                        >  1:2519245/1‑66 (MQ=255)
                                              gATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGccatccat                        >  1:2031065/1‑66 (MQ=255)
                                                  ccGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTg                 >  1:798798/1‑69 (MQ=255)
                                                  ccGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTg                 >  1:1023626/1‑69 (MQ=255)
                                                  ccGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTg                 >  1:171342/1‑69 (MQ=255)
                                                  ccGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTg                 >  1:316003/1‑69 (MQ=255)
                                                  ccGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTg                 >  1:2482986/1‑69 (MQ=255)
                                                  ccGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTg                 >  1:2247162/1‑69 (MQ=255)
                                                  ccGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTg                 >  1:2243116/1‑69 (MQ=255)
                                                  ccGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTg                 >  1:2131850/1‑69 (MQ=255)
                                                         gCGTCGCTGCACCAAAAGTTCGTAGTAGCGGATCGATAGATCCCGCGCCATCCATTACGCTGAATTATc          <  1:347568/69‑1 (MQ=255)
                                                         gCGTCGCTGCACCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCGCGCCATCCATTACGCTGAATTATc          <  1:1356270/69‑1 (MQ=255)
                                                         gCGTCGCTGCACCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCGCGCCATCCATTACGCTGAATTATc          <  1:3322340/69‑1 (MQ=255)
                                                         gCGTCGCTGCACCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCGCGCCATCCATTACGCTGAATTATc          <  1:3016770/69‑1 (MQ=255)
                                                         gCGTCGCTGCACCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCGCGCCATCCATTACGCTGAATTATc          <  1:861349/69‑1 (MQ=255)
                                                         gCGTCGCTGCACCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCGCGCCATCCATTACGCTGAATTATc          <  1:2411340/69‑1 (MQ=255)
                                                         gCGTCGCTGCACCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCGCGCCATCCATTACGCTGAATTATc          <  1:1967639/69‑1 (MQ=255)
                                                          cGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATcc         >  1:2829707/1‑69 (MQ=255)
                                                                 gcgcCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAg  <  1:1255017/69‑1 (MQ=255)
                                                                 gcgcCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAg  <  1:3282142/69‑1 (MQ=255)
                                                                 gcgcCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAg  <  1:1499769/69‑1 (MQ=255)
                                                                   |                                                                  
TCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAG  >  minE/452483‑452616

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: