Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,845,199 G→T 58.8% P266T (CCG→ACG)  cysD ← sulfate adenylyltransferase, subunit 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,845,1990GT58.8% 0.2 / 106.9 68P266T (CCG→ACG) cysDsulfate adenylyltransferase, subunit 2
Reads supporting (aligned to +/- strand):  ref base G (13/15);  new base T (29/11);  total (42/26)
Fisher's exact test for biased strand distribution p-value = 4.27e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.31e-01

CGCGGTCAATCACGCGGCCCTGACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCCAGCAGCCCAGCGTACGGAAACGCAC  >  minE/1845133‑1845265
                                                                  |                                                                  
cgcgGTCAATCACGCGGCCCTGACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCAgt                                                                >  1:3995339/1‑71 (MQ=255)
cgcgGTCAATCACGCGGCCCTGACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCAgt                                                                >  1:1736184/1‑71 (MQ=255)
       aaTCACGCGGCCCTGACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTtgtg                                                         <  1:3993707/71‑1 (MQ=255)
         tCACGCGGCCCTGACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTtg                                                           >  1:1865864/1‑67 (MQ=255)
         tCACGCGGCCCTGACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCa                                                       >  1:4055338/1‑71 (MQ=255)
             gcgGCCCTGACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTg                                                   <  1:651525/71‑1 (MQ=255)
             gcgGCCCTGACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTg                                                   <  1:537581/71‑1 (MQ=255)
             gcgGCCCTGACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTg                                                   <  1:3843774/71‑1 (MQ=255)
                     gACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCATTTGACTCcacc                                           >  1:1509455/1‑71 (MQ=255)
                     gACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCATTTGACTCcacc                                           >  1:1608494/1‑71 (MQ=255)
                     gACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCATTTGACTCcacc                                           >  1:3611259/1‑71 (MQ=255)
                     gACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCATTTGACTCcac                                            >  1:1567048/1‑70 (MQ=255)
                     gACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCATTTGACTCcac                                            >  1:3978301/1‑70 (MQ=255)
                     gACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCATTTGACTCcac                                            >  1:2790119/1‑70 (MQ=255)
                     gACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCATTTGACTCca                                             >  1:3341279/1‑69 (MQ=255)
                     gACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCATTTGACTCca                                             >  1:1177856/1‑69 (MQ=255)
                     gACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCcacc                                           >  1:848908/1‑71 (MQ=255)
                     gACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCcac                                            >  1:1732141/1‑70 (MQ=255)
                      aCGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCATTTGACTCcaccg                                          >  1:3531633/1‑71 (MQ=255)
                       cGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCcaccgc                                         <  1:3953593/71‑1 (MQ=255)
                         ttCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCAgtgt                                                              <  1:2954971/48‑1 (MQ=255)
                         ttCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCAgtgt                                                              <  1:3424284/48‑1 (MQ=255)
                          tCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCATTTGACTCcaccgcacc                                      >  1:1621804/1‑71 (MQ=255)
                           cACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCcaccgcacc                                      <  1:2412911/70‑1 (MQ=255)
                              tggtggtggAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCAttt                                                    <  1:2952906/53‑1 (MQ=255)
                              tggtggtggAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCAttt                                                    <  1:2249635/53‑1 (MQ=255)
                              tggtggtggAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCAttt                                                    <  1:3479560/53‑1 (MQ=255)
                              tggtggtggAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCAttt                                                    <  1:2607283/53‑1 (MQ=255)
                              tggtggtggAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCAttt                                                    <  1:695844/53‑1 (MQ=255)
                              tggtggtggAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCATTTGACTCCACCGCACCGGt                                   >  1:3778542/1‑70 (MQ=255)
                              tggtggtggAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCAttt                                                    <  1:1123/53‑1 (MQ=255)
                                gtggtggAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAg                                >  1:2354294/1‑71 (MQ=255)
                                gtggtggAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAg                                >  1:2666561/1‑71 (MQ=255)
                                gtggtggAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAg                                >  1:1783013/1‑71 (MQ=255)
                                   gtggAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCgg                             >  1:567843/1‑71 (MQ=255)
                                   gtggAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCgg                             >  1:2268850/1‑71 (MQ=255)
                                    tggAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGc                            <  1:3006529/71‑1 (MQ=255)
                                    tggAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGc                            <  1:206057/71‑1 (MQ=255)
                                    tggAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGc                            <  1:2670776/71‑1 (MQ=255)
                                     ggAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCa                                                       <  1:207666/43‑1 (MQ=255)
                                     ggAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGc                            >  1:1423812/1‑70 (MQ=255)
                                     ggAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGc                            >  1:942381/1‑70 (MQ=255)
                                         aCCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCcag                         >  1:1456385/1‑69 (MQ=255)
                                         aCCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCcag                         >  1:1252290/1‑69 (MQ=255)
                                            aGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCcagcaggc                    <  1:2180930/71‑3 (MQ=255)
                                            aGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCCAGCAGcc                    <  1:3532325/71‑1 (MQ=255)
                                               aTCTCTTCGATGATTTCCGTCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCCAGCAGCCCAg                 >  1:988109/1‑71 (MQ=255)
                                               aTCTCTTCGATGATTTCCGTCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCCAGCAGCCCAg                 >  1:2159113/1‑71 (MQ=255)
                                               aTCTCTTCGATGATTTCCGTCAGTGTTTGTACATTTGACTCCACCGCACCGGTCAGCGGCCAGCAGCCCAg                 >  1:4118037/1‑71 (MQ=255)
                                                tctcTTCGATGATTTCCGTCTGTGTTTGTGCATTTGa                                                  >  1:3315762/1‑37 (MQ=37)
                                                tctcTTCGATGATTTCCGTCAGTGTTTGTGCATTTGa                                                  >  1:2599210/1‑37 (MQ=255)
                                                tctcTTCGATGATTTCCGTCAGTGTTTGTGCATTTGa                                                  >  1:3226304/1‑37 (MQ=255)
                                                tctcTTCGATGATTTCCGTCAGTGTTTGTGCATTTGa                                                  >  1:1650191/1‑37 (MQ=255)
                                                   ctTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCcaccg                                          >  1:3150209/1‑42 (MQ=255)
                                                    ttCGATGATTTCCGTCAGTGTTTGTGCATTTGACTCCACCGCACCGGTAAGCGGCCAGGAGCCCAGCg               <  1:2257988/68‑1 (MQ=255)
                                                    ttCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCCAGCAGCCCAGCGt              <  1:3558351/69‑1 (MQ=255)
                                                    ttCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCCAGCAGCCCAGCGt              <  1:111917/69‑1 (MQ=255)
                                                    ttCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCCAGCAGCCCAGCGt              <  1:1999958/69‑1 (MQ=255)
                                                    ttCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCCAGCAGCCCAGCGt              <  1:553640/69‑1 (MQ=255)
                                                    ttCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGAGGCCAGCAGCCCAGCGt              <  1:1524598/69‑1 (MQ=255)
                                                      cGATGATTTCCGTCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCCAGCAGCCCAg                 >  1:1163989/1‑64 (MQ=255)
                                                      cGATGATTTCCGTCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCCAGCAGCCCAGCGTACg           >  1:2748069/1‑70 (MQ=255)
                                                      cGATGATTTCCGTCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCCAGCAGCCCAGCGTACg           >  1:3408517/1‑70 (MQ=255)
                                                      cGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCca                          >  1:1161504/1‑55 (MQ=255)
                                                              tCCGTCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCcagca                       >  1:1033795/1‑50 (MQ=255)
                                                              tCCGTCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCcagca                       >  1:2000910/1‑50 (MQ=255)
                                                              tCCGTCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCcagca                       >  1:2402332/1‑50 (MQ=255)
                                                               ccGTCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCCAGCAGCCCAGCGTACGGAAACGCAc  >  1:1871381/1‑70 (MQ=255)
                                                                  |                                                                  
CGCGGTCAATCACGCGGCCCTGACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCCAGCAGCCCAGCGTACGGAAACGCAC  >  minE/1845133‑1845265

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: