Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 597,804 C→T 52.7% A49A (GCC→GCT dmsA → dimethyl sulfoxide reductase, anaerobic, subunit A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE597,8040CT52.7% 1.3 / 60.7 71A49A (GCC→GCTdmsAdimethyl sulfoxide reductase, anaerobic, subunit A
Reads supporting (aligned to +/- strand):  ref base C (9/24);  new base T (21/16);  total (30/41)
Fisher's exact test for biased strand distribution p-value = 1.64e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACATTACCTTTTAGTCGGATTGCGCACGCTGTCGATAGCGCCATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGCCTGTACAGTTAACTGTGGTAGT  >  minE/597762‑597864
                                          |                                                            
aaCATTACCTTTTAGTCGGATTGCGCACGCTGTCGATAGCGCCATTCCAACAAAATCAGACGAAAAGGt                                    >  1:3113800/1‑69 (MQ=255)
aaCATTACCTTTTAGTCGGATTGCGCACGCTGTCGATAGCGCCATTCCAACAAAATCAGACGAAAAGGt                                    <  1:2533997/69‑1 (MQ=255)
aaCATTACCTTTTAGTCGGATTGCGCACGCTGTCGATAGCGCCATTCCAACAAAATCAGACGAAAAGGt                                    <  1:2423035/69‑1 (MQ=255)
aaCATTACCTTTTAGTCGGATTGCGCACGCTGTCGATAGCGCCATTCCAACAAAATCAGACGAAAAGGt                                    >  1:2372025/1‑69 (MQ=255)
      aCCTTTTAGTCGGATTGCGCACGCTGTCGATAGCGCCATTCCAACAAAATCAGACGAAAAGGTTATCTg                              <  1:1676050/69‑1 (MQ=255)
         ttttAGTCGGATTGCGCACGCTGTCGATAGCGCCATTCCAACAAAATCAGACGAAAAGGTTATCTgg                             <  1:1756379/67‑1 (MQ=255)
           ttAGTCGGATTGCGCACGCTGTCGATAGCGCCATTCCAACAAAATCAGACGAAAAGGTTATCTGGAgcg                         >  1:629780/1‑69 (MQ=255)
           ttAGTCGGATTGCGCACGCTGTCGATAGCGCCATTCCAACAAAATCAGACGAAAAGGTTATCTGGAgcg                         >  1:157226/1‑69 (MQ=255)
              gTCGGATTGCGCACGCTGTCGATAGCGCCATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGcc                       <  1:958353/68‑1 (MQ=255)
              gTCGGATTGCGCACGCTGTCGATAGCGCCATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGcc                       <  1:48902/68‑1 (MQ=255)
              gTCGGATTGCGCACGCTGTCGATAGCGCCATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGcc                       <  1:2146763/68‑1 (MQ=255)
              gTCGGATTGCGCACGCTGTCGATAGCGCCATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGcc                       <  1:2034830/68‑1 (MQ=255)
              gTCGGATTGCGCACGCTGTCGATAGCGCCATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGcc                       <  1:118616/68‑1 (MQ=255)
              gTCGGATTGCGCACGCTGTCGATAGCGCCATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGcc                       <  1:122513/68‑1 (MQ=255)
               tCGGATTGCGCACGCTGTCGATAGCGCCATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGcc                       <  1:255233/67‑1 (MQ=255)
                cGGATTGCGCACGCTGTCGATAGCGCCATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGcc                       <  1:1026078/66‑1 (MQ=255)
                cGGATTGCGCACGCTGTCGATAGCGCCATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGcc                       <  1:1251893/66‑1 (MQ=255)
                 ggATTGCGCACGCTGTCGATAGCGCCATTCCAACAAAATCAGACGAAAAGTTTATCTGGAGCGCCTGTa                   <  1:2830454/69‑1 (MQ=255)
                 ggATTGCGCACGCTGTCGATAGCGCCATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGCCTg                     >  1:274796/1‑67 (MQ=255)
                 ggATTGCGCACGCTGTCGATAGCGCCATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGCCTg                     >  1:2517543/1‑67 (MQ=255)
                 ggATTGCGCACGCTGTCGATAGCGCCATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGCCTg                     >  1:500003/1‑67 (MQ=255)
                 ggATTGCGCACGCTGTCGATAGCGCCATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGCCTGTa                   <  1:856925/69‑1 (MQ=255)
                   aTTGCGCACGCTGTCGATAGCGCCATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGcc                       <  1:204630/63‑1 (MQ=255)
                   aTTGCGCACGCTGTCGATAGCGCCATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGCCTGTAc                  <  1:567/68‑1 (MQ=255)
                     tGCGCACGCTGTCGATAGCGCAATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGCCTGTAc                  <  1:1918899/66‑1 (MQ=255)
                         cACGCTGTCGATAGCGCCATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGCCTGTACAGTTAAc           <  1:1052975/69‑1 (MQ=255)
                         cACGCTGTCGATAGCGCCATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGCCTGTACAGTTAAc           <  1:2825391/69‑1 (MQ=255)
                         cACGCTGTCGATAGCGCCATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGCCTGTACAGTTAAc           <  1:874172/69‑1 (MQ=255)
                         cACGCTGTCGATAGCGCCATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGCCTGTACAGTTAAc           <  1:286662/69‑1 (MQ=255)
                         cACGCTGTCGATAGCGCCAATCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGCCTGTACAGTTAAc           <  1:319721/69‑1 (MQ=255)
                               gTCGATAGCGCCATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGCCTGTAc                  <  1:1169745/56‑1 (MQ=255)
                                 cGATAGCGCTATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGCCTGTACa                 >  1:2584807/1‑55 (MQ=255)
                                 cGATAGCGCTATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGCCTGTACAGTTAACtgt        >  1:1023158/1‑64 (MQ=255)
                                 cGATAGCGCTATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGCCTGTACAGTTAACtgt        >  1:517976/1‑64 (MQ=255)
                                 cGATAGCGCTATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGCCTGTACAGTTAACtgt        >  1:835115/1‑64 (MQ=255)
                                 cGATAGCGCTATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGCCTGTACAGTTAACtgt        >  1:2848216/1‑64 (MQ=255)
                                 cGATAGCGCTATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGCCTGTACAGTTAACtgt        >  1:2862849/1‑64 (MQ=255)
                                 cGATAGCGCTATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGCCTGTACAGTTAACtgt        >  1:3037386/1‑64 (MQ=255)
                                 cGATAGCGCTATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGCCTGTACAGTTAACtgt        >  1:3048703/1‑64 (MQ=255)
                                 cGATAGCGCTATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGCCTGTACAGTTAACtgt        >  1:356400/1‑64 (MQ=255)
                                 cGATAGCGCTATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGCCTGTACAGTTAACtgt        >  1:476366/1‑64 (MQ=255)
                                 cGATAGCGCTATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGCCTGTACAGTTAACtgt        >  1:263555/1‑64 (MQ=255)
                                 cGATAGCGCTATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGCCTGTACAGTTAACtgt        >  1:268347/1‑64 (MQ=255)
                                 cGATAGCGCTATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGCCTGTACAGTTAACtgt        >  1:1582571/1‑64 (MQ=255)
                                 cGATAGCGCTATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGCCTGTACAGTTAACtgt        >  1:234549/1‑64 (MQ=255)
                                 cGATAGCGCTATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGCCTGTACAGTTAACtgt        >  1:1585875/1‑64 (MQ=255)
                                 cGATAGCGCTATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGCCTGTACAGTTAACtgt        >  1:2339738/1‑64 (MQ=255)
                                 cGATAGCGCTATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGCCTGTACAGTTAACtgt        >  1:2208690/1‑64 (MQ=255)
                                 cGATAGCGCTATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGCCTGTACAGTTAACtgt        >  1:148881/1‑64 (MQ=255)
                                 cGATAGCGCTATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGCCTGTACAGTTAACtgt        >  1:1707602/1‑64 (MQ=255)
                                 cGATAGCGCTATTCCAACAAAATCAGACGAAAAGGTTATATGGAGCGCCTGTACAGTTAACtgt        >  1:3104685/1‑64 (MQ=255)
                                 cGATAGCGCTAGTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGCCTGTACAGTTAACtgt        >  1:1501181/1‑64 (MQ=255)
                                  gATAGCGCCATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGCCTGTACAGTTAACTGTGGTAGt  >  1:3406264/1‑69 (MQ=255)
                                  gATAGCGCCATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGCCTGTACAGTTAACTGTGGTAGt  >  1:1191187/1‑69 (MQ=255)
                                    tAGCGCTATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGCCTGTTCa                 <  1:1994596/52‑1 (MQ=39)
                                    tAGCGCTATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGCCTGTACa                 <  1:2392480/52‑1 (MQ=255)
                                    tAGCGCTATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGCCTGTACa                 <  1:3283008/52‑1 (MQ=255)
                                    tAGCGCTATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGCCTGTACa                 <  1:943805/52‑1 (MQ=255)
                                    tAGCGCTATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGCCTGTACa                 <  1:1160874/52‑1 (MQ=255)
                                    tAGCGCTATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGCCTGTACa                 <  1:2167991/52‑1 (MQ=255)
                                    tAGCGCTATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGCCTGTACa                 <  1:1758901/52‑1 (MQ=255)
                                    tAGCGCTATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGCCTGTACa                 <  1:558164/52‑1 (MQ=255)
                                    tAGCGCTATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGCCTGTACa                 <  1:532730/52‑1 (MQ=255)
                                    tAGCGCTATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGCCTGTACa                 <  1:503039/52‑1 (MQ=255)
                                    tAGCGCTATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGCCTGTACa                 <  1:2000794/52‑1 (MQ=255)
                                    tAGCGCTATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGCCTGTACa                 <  1:1237932/52‑1 (MQ=255)
                                    tAGCGCTATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGCCTGTACa                 <  1:2048840/52‑1 (MQ=255)
                                    tAGCGCTATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGCCTGTACa                 <  1:445416/52‑1 (MQ=255)
                                    tAGCGCTATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGCCTGTACa                 <  1:2114090/52‑1 (MQ=255)
                                    tAGCGCTATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCCCCTGTACa                 <  1:301011/52‑1 (MQ=255)
                                        gcCATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGCCTGTAc                  <  1:855469/47‑1 (MQ=255)
                                          |                                                            
AACATTACCTTTTAGTCGGATTGCGCACGCTGTCGATAGCGCCATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGCCTGTACAGTTAACTGTGGTAGT  >  minE/597762‑597864

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: