Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,553,364 A→C 49.4% intergenic (‑115/‑103) dsdC ← / → dsdX DNA‑binding transcriptional dual regulator/predicted transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,553,3640AC49.4% 19.4 / 200.8 154intergenic (‑115/‑103)dsdC/dsdXDNA‑binding transcriptional dual regulator/predicted transporter
Reads supporting (aligned to +/- strand):  ref base A (69/9);  new base C (16/60);  total (85/69)
Fisher's exact test for biased strand distribution p-value = 5.09e-18
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.72e-01

Only 100 of 154 total aligned reads displayed.

AGATGAGCTAAAGTGAACCATATCTCAATTCACCTTCATTTTTAGATGTAAATCACTCCATTGATGCAATTTACCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGTTTTTTTGTCC  >  minE/1553305‑1553431
                                                           |                                                                   
agatgagCTAAAGTGAACCATATCTCAATTCACCTTCATTTTTAGATGTAAATCACTCCATTGATGCaa                                                            >  1:2857744/1‑69 (MQ=255)
agatgagCTAAAGTGAACCATATCTCAATTCACCTTCATTTTTAGATGTAAATCACTCCATTGATGCaa                                                            >  1:3301967/1‑69 (MQ=255)
agatgagCTAAAGTGAACCATATCTCAATTCACCTTCATTTTTAGATGTAAATCACTCCATTGATGCaa                                                            >  1:2514393/1‑69 (MQ=255)
agatgagCTAAAGTGAACCATATCTCAATTCACCTTCATTTTTAGATGTAAATCACTCCATTGATGCaa                                                            >  1:1639761/1‑69 (MQ=255)
agatgagCTAAAGTGAACCATATCTCAATTCACCTTCATTTTTAGATGTAAATCACTCCATTGATGCaa                                                            >  1:2053566/1‑69 (MQ=255)
agatgagCTAAAGTGAACCATATCTCAATTCACCTTCATTTTTAGATGTAAATCACTCCATTGATGCaa                                                            >  1:1602684/1‑69 (MQ=255)
agatgagCTAAAGTGAACCATATCTCAATTCACCTTCATTTTTAGATGTAAATCACTCCATTGATGCaa                                                            >  1:2592839/1‑69 (MQ=255)
agatgagCTAAAGTGAACCATATCTCAATTCACCTTCATTTTTAGATGTAAATCACTCCATTGATGCaa                                                            >  1:1222176/1‑69 (MQ=255)
agatgagCTAAAGTGAACCATATCTCAATTCACCTTCATTTTTAGATGTAAATCACTCCATTGATGCaa                                                            >  1:1410221/1‑69 (MQ=255)
   tgagCTAAAGTGAACCATATCTCAATTCACCTTCATTTTTAGATGTAAATCACTCCATTGATGCAAttt                                                         >  1:2074727/1‑69 (MQ=255)
   tgagCTAAAGTGAACCATATCTCAATTCACCTTCATTTTTAGATGTAAATCACTCCATTGATGCAAttt                                                         >  1:3420085/1‑69 (MQ=255)
   tgagCTAAAGTGAACCATATCTCAATTCACCTTCATTTTTAGATGTAAATCACTCCATTGATGCAAttt                                                         >  1:2304907/1‑69 (MQ=255)
   tgagCTAAAGTGAACCATATCTCAATTCACCTTCATTTTTAGATGTAAATCACTCCATTGATGCAAttt                                                         >  1:126620/1‑69 (MQ=255)
   tgagCTAAAGTGAACCATATCTCAATTCACCTTCATTTTTAGATGTAAATCACTCCATTGATGCAAttt                                                         >  1:3154496/1‑69 (MQ=255)
        tAAAGTGAACCATATCTCAATTCACCTTCATTTTTAGATGTAAATCACTCCATTGATGCAAttt                                                         <  1:2343734/64‑1 (MQ=255)
         aaaGTGAACCATATCTCAATTCACCTTCATTTTTAGATGTAAATCACTCCCTTGATGCAATTTAGCTCa                                                   <  1:2100861/69‑1 (MQ=255)
         aaaGTGAACCATATCTCAATTCACCTTCATTTTTAGATGTAAATCACTCCCTTGATGCAATTTAGCTCa                                                   <  1:636096/69‑1 (MQ=255)
         aaaGTGAACCATATCTCAATTCACCTTCATTTTTAGATGTAAATCACTCCCTTGATGCAATTTAGCTCa                                                   <  1:3260198/69‑1 (MQ=255)
         aaaGTGAACCATATCTCAATTCACCTTCATTTTTAGATGTAAATCACTCCCTTGATGCAATTTAGCTCa                                                   <  1:2763993/69‑1 (MQ=255)
         aaaGTGAACCATATCTCAATTCACCTTCATTTTTAGATGTAAATCACTCCCTTGATGCAATTTAGCTCa                                                   <  1:47763/69‑1 (MQ=255)
         aaaGTGAACCATATCTCAATTCACCTTCATTTTTAGATGTAAATCACTCCCTTGATGCAATTTAGCTCa                                                   <  1:2803763/69‑1 (MQ=255)
         aaaGTGAACCATATCTCAATTCACCTTCATTTTTAGATGTAAATCACTCCCTTGATGCAATTTAGCTCa                                                   <  1:951485/69‑1 (MQ=255)
         aaaGTGAACCATATCTCAATTCACCTTCATTTTTAGATGTAAATCACTCCCTTGATGCAATTTAGCTCa                                                   <  1:1069821/69‑1 (MQ=255)
         aaaGTGAACCATATCTCAATTCACCTTCATTTTTAGATGTAAATCACTCCCTTGATGCAATTTAGCTCa                                                   <  1:3262941/69‑1 (MQ=255)
                      tctcAATTCACCTTCATTTTTAGATGTAAATCACTCCCTTGATGCAATTTAGCTCa                                                   <  1:1492768/56‑1 (MQ=255)
                      tctcAATTCACCTTCATTTTTAGATGTAAATCACTCCCTTGATGCAATTTAGCTCa                                                   <  1:2391064/56‑1 (MQ=255)
                      tctcAATTCACCTTCATTTTTAGATGTAAATCACTCCATTGATGCAATTTACCTCATGTGAAAGGCaaa                                      >  1:3062288/1‑69 (MQ=255)
                      tctcAATTCACCTTCATTTTTAGATGTAAATCACTCCATTGATGCAATTTACCTCATGTGAAAGGCaaa                                      >  1:2322632/1‑69 (MQ=255)
                        tcAATTCACCTTCATTTTTAGATGTAAATCACTCCATTGATGCAAtt                                                          >  1:2556496/1‑47 (MQ=255)
                        tcAATTCACCTTCATTTTTAGATGTAAATCACTCCATTGATGCAATTTACCTCATGTGAAAGGCaaa                                      >  1:2121386/1‑67 (MQ=255)
                        tcAATTCACCTTCATTTTTAGATGTAAATCACTCCATTGATGCAATTTACCTCATGTGAAAGGCAAAt                                     >  1:2389676/1‑68 (MQ=255)
                            ttCACCTTCATTTTTAGATGTAAATCACTCCCTTGATGCAATTTAGCTCATGTGAAAGGCAAAttt                                   <  1:2545400/66‑1 (MQ=255)
                            ttCACCTTCATTTTTAGATGTAAATCACTCCCTTGATGCAATTTAGCTCATGTGAAAGGCAAAttt                                   <  1:3034608/66‑1 (MQ=255)
                            ttCACCTTCATTTTTAGATGTAAATCACTCCCTTGATGCAATTTAGCTCATGTGAAAGGCAAAttt                                   <  1:2504551/66‑1 (MQ=255)
                            ttCACCTTCATTTTTAGATGTAAATCACTCCCTTGATGCAATTTAGCTCATGTGAAAGGCAAAttt                                   <  1:2444326/66‑1 (MQ=255)
                            ttCACCTTCATTTTTAGATGTAAATCACTCCCTTGATGCAATTTAGCTCATGTGAAAGGCAAAttt                                   <  1:1929046/66‑1 (MQ=255)
                            ttCACCTTCATTTTTAGATGTAAATCACTCCCTTGATGCAATTTAGCTCATGTGAAAGGCAAAttt                                   <  1:1119119/66‑1 (MQ=255)
                            ttCACCTTCATTTTTAGATGTAAATCACTCCCTTGATGCAATTTAGCTCATGTGAAAGGCAAAttt                                   <  1:809922/66‑1 (MQ=255)
                            ttCACCTTCATTTTTAGATGTAAATCACTCCCTTGATGCAATTTAGCTCATGTGAAAGGCAAAttt                                   <  1:1151124/66‑1 (MQ=255)
                            ttCACCTTCATTTTTAGATGTAAATCACTCCCTTGATGCAATTTAGCTCATGTGAAAGGCAAAttt                                   <  1:1292661/66‑1 (MQ=255)
                            ttCACCTTCATTTTTAGATGTAAATCACTCCCTTGATGCAATTTAGCTCATGTGAAAGGCAAAttt                                   <  1:681269/66‑1 (MQ=255)
                            ttCACCTTCATTTTTAGATGTAAATCACTCCCTTGATGCAATTTAGCTCATGTGAAAGGCAAAttt                                   <  1:552025/66‑1 (MQ=255)
                            ttCACCTTCATTTTTAGATGTAAATCACTCCCTTGATGCAATTTAGCTCATGTGAAAGGCAAAttt                                   <  1:404660/66‑1 (MQ=255)
                            ttCACCTTCATTTTTAGATGTAAATCACTCCCTTGATGCAATTTAGCTCATGTGAAAGGCAAAttt                                   <  1:3301101/66‑1 (MQ=255)
                            ttCACCTTCATTTTTAGATGTAAATCACTCCCTTGATGCAATTTAGCTCATGTGAAAGGCAAAttt                                   <  1:1901096/66‑1 (MQ=255)
                            ttCACCTTCATTTTTAGATGTAAATCACTCCCTTGATGCAATTTAGCTCATGTGAAAGGCAAAttt                                   <  1:2017053/66‑1 (MQ=255)
                              cACCTTCATTTTTAGACGTAAATCACTCCATTGATGCAATTTACCTCATGTGAAAGGCAAATTTTATCg                              <  1:799680/69‑1 (MQ=255)
                                     aTTTTTAGATGTAAATCACTCCCTTGATGCAATTTAGCTCATGTGAAAGGCAAATTTTATCGTTTGTCa                       >  1:600120/1‑69 (MQ=255)
                                     aTTTTTAGATGTAAATCACTCCCTTGATGCAATTTAGCTCATGTGAAAGGCAAATTTTATCGTTTGTCa                       >  1:741435/1‑69 (MQ=255)
                                     aTTTTTAGATGTAAATCACTCCCTTGATGCAATTTAGCTCATGTGAAAGGCAAATTTTATCGTTTGTCa                       >  1:75954/1‑69 (MQ=255)
                                     aTTTTTAGATGTAAATCACTCCCTTGATGCAATTTAGCTCATGTGAAAGGCAAATTTTATCGTTTGTCa                       >  1:96495/1‑69 (MQ=255)
                                     aTTTTTAGATGTAAATCACTCCCTTGATGCAATTTAGCTCATGTGAAAGGCAAATTTTATCGTTTGTCa                       >  1:2433065/1‑69 (MQ=255)
                                     aTTTTTAGATGTAAATCACTCCCTTGATGCAATTTAGCTCATGTGAAAGGCAAATTTTATCGTTTGTCa                       >  1:1035314/1‑69 (MQ=255)
                                     aTTTTTAGATGTAAATCACTCCCTTGATGCAATTTAGCTCATGTGAAAGGCAAATTTTATCGTTTGTCa                       >  1:1745220/1‑69 (MQ=255)
                                     aTTTTTAGATGTAAATCACTCCCTTGATGCAATTTAGCTCATGTGAAAGGCAAATTTTATCGTTTGTCa                       >  1:1662736/1‑69 (MQ=255)
                                     aTTTTTAGATGTAAATCACTCCCTTGATGCAATTTAGCTCATGTGAAAGGCAAATTTTATCGTTTGTCa                       >  1:208868/1‑69 (MQ=255)
                                     aTTTTTAGATGTAAATCACTCCCTTGATGCAATTTAGCTCATGTGAAAGGCAAATTTTATCGTTTGTCa                       >  1:1120073/1‑69 (MQ=255)
                                     aTTTTTAGATGTAAATCACTCCCTTGATGCAATTTAGCTCATGTGAAAGGCAAATTTTATCGTTTGTCa                       >  1:2264939/1‑69 (MQ=255)
                                      tttttAGATGTAAATCACTCCATTGATGCAATTTACCTCATGTGAAAGGCAAATTTTATCGTTTGTCa                       <  1:2374835/68‑1 (MQ=255)
                                      tttttAGATGTAAATCACTCCATTGATGCAATTTACCTCATGTGAAAGGCAAATTTTATCGTTTGTCa                       <  1:2082301/68‑1 (MQ=255)
                                        tttAGATGTAAATCACTCCATTGATGCAATTTACCTCATGTGAAAGGCAAATTTTATCGTTTGTCAg                      <  1:1405881/67‑1 (MQ=255)
                                                tAAATCACTCCATTGATGCAATTTACCTCATGTGAAAGGCAAAt                                     >  1:3010484/1‑44 (MQ=255)
                                                tAAATCACTCCATTGATGCAATTTACCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCgttgt            >  1:1072319/1‑69 (MQ=255)
                                                tAAATCACTCCATTGATGCAATTTACCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCgttgt            >  1:166557/1‑69 (MQ=255)
                                                tAAATCACTCCATTGATGCAATTTACCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCgttgt            >  1:3156003/1‑69 (MQ=255)
                                                tAAATCACTCCATTGATGCAATTTACCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCgttgt            >  1:1457208/1‑69 (MQ=255)
                                                tAAATCACTCCATTGATGCAATTTACCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCgttgt            >  1:1433383/1‑69 (MQ=255)
                                                tAAATCACTCCATTGATGCAATTTACCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCgttgt            >  1:2363570/1‑69 (MQ=255)
                                                tAAATCACTCCATTGATGCAATTTACCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCgttgt            >  1:1213934/1‑69 (MQ=255)
                                                  aaTCACTCCATTGATGCAATTTACCTCATGTGAAAGGCAAAtttt                                  >  1:1008888/1‑45 (MQ=255)
                                                  aaTCACTCCATTGATGCAATTTACCTCATGTGAAAGGCAAATTTTATCg                              >  1:2152268/1‑49 (MQ=255)
                                                  aaTCACTCCATTGATGCAATTTACCTCATGTGAAAGGCAAATTTTATCg                              >  1:25838/1‑49 (MQ=255)
                                                  aaTCACTCCATTGATGCAATTTACCTCATGTGAAAGGCAAATTTTATCg                              >  1:172752/1‑49 (MQ=255)
                                                  aaTCACTCCATTGATGCAATTTACCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGttt          >  1:476889/1‑69 (MQ=255)
                                                  aaTCACTCCATTGATGCAATTTACCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGttt          >  1:370475/1‑69 (MQ=255)
                                                  aaTCACTCCATTGATGCAATTTACCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGttt          >  1:3385586/1‑69 (MQ=255)
                                                  aaTCACTCCATTGATGCAATTTACCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGttt          >  1:3349662/1‑69 (MQ=255)
                                                  aaTCACTCCATTGATGCAATTTACCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGttt          >  1:1792528/1‑69 (MQ=255)
                                                  aaTCACTCCATTGATGCAATTTACCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGttt          >  1:2803110/1‑69 (MQ=255)
                                                  aaTCACTCCATTGATGCAATTTACCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGttt          >  1:3111759/1‑69 (MQ=255)
                                                  aaTCACTCCATTGATGCAATTTACCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGttt          >  1:2007400/1‑69 (MQ=255)
                                                  aaTCACTCCATTGATGCAATTTACCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGttt          >  1:314651/1‑69 (MQ=255)
                                                    tCACTCCCTTGATGCAATTTAGCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGttttt        <  1:1807624/69‑1 (MQ=255)
                                                    tCACTCCCTTGATGCAATTTAGCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGttttt        <  1:2722322/69‑1 (MQ=255)
                                                    tCACTCCCTTGATGCAATTTAGCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGttttt        <  1:920659/69‑1 (MQ=255)
                                                    tCACTCCCTTGATGCAATTTAGCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGttttt        <  1:2173326/69‑1 (MQ=255)
                                                    tCACTCCCTTGATGCAATTTAGCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGttttt        <  1:3224158/69‑1 (MQ=255)
                                                    tCACTCCCTTGATGCAATTTAGCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGttttt        <  1:401386/69‑1 (MQ=255)
                                                    tCACTCCCTTGATGCAATTTAGCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGttttt        <  1:1977095/69‑1 (MQ=255)
                                                    tCACTCCCTTGATGCAATTTAGCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGttttt        <  1:3273330/69‑1 (MQ=255)
                                                     cACTCCATTGATGCAATTTACCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGtttttt       >  1:1162339/1‑69 (MQ=255)
                                                     cACTCCATTGATGCAATTTACCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGtttttt       >  1:3226519/1‑69 (MQ=255)
                                                         ccATTGATGCAATTTACCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGTTTTTTTGt    >  1:2167439/1‑68 (MQ=255)
                                                         ccATTGATGCAATTTACCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGTTTTTTTGt    >  1:806068/1‑68 (MQ=255)
                                                          cATTGATGCAATTTACCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGTTTTTTTGTcc  >  1:148354/1‑69 (MQ=255)
                                                          cATTGATGCAATTTACCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGTTTTTTTGTcc  >  1:570513/1‑69 (MQ=255)
                                                          cATTGATGCAATTTACCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGTTTTTTTGTcc  >  1:3130108/1‑69 (MQ=255)
                                                          cATTGATGCAATTTACCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGTTTTTTTGTcc  >  1:1261157/1‑69 (MQ=255)
                                                          cATTGATGCAATTTACCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGTTTTTTTGTcc  >  1:924637/1‑69 (MQ=255)
                                                          cATTGATGCAATTTACCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGTTTTTTTGTcc  >  1:2408842/1‑69 (MQ=255)
                                                           |                                                                   
AGATGAGCTAAAGTGAACCATATCTCAATTCACCTTCATTTTTAGATGTAAATCACTCCATTGATGCAATTTACCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGTTTTTTTGTCC  >  minE/1553305‑1553431

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: