Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,845,301 C→T 56.3% D232N (GAC→AAC)  cysD ← sulfate adenylyltransferase, subunit 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,845,3010CT56.3% 18.1 / 69.2 87D232N (GAC→AAC) cysDsulfate adenylyltransferase, subunit 2
Reads supporting (aligned to +/- strand):  ref base C (9/29);  new base T (42/7);  total (51/36)
Fisher's exact test for biased strand distribution p-value = 4.34e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 112 total aligned reads displayed.

CAGCGGCCAGCAGCCCAGCGTACGGAAACGCACCATCCGTTTTTTAATCACTTCGCCCGGTTGCAGGTCGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCGAGA  >  minE/1845233‑1845368
                                                                    |                                                                   
cAGCGGCCAGCTGCCCAGCGTACGGAAACGCACCATCCGTTTTTTAATCACTTCGCCCGGTTGCAGGTc                                                                     <  1:2099748/69‑1 (MQ=255)
cAGCGGCCAGCAGCCCAGCGTACGGAAACGCACCATCCGTTTTTTAATCACTTCGCCCGGTTGCAGGTc                                                                     <  1:388127/69‑1 (MQ=255)
cAGCGGCCAGCAGCCCAGCGTACGGAAACGCACCATCCGTTTTTTAATCACTTCGCCCGGTTGCAGGTc                                                                     >  1:1125582/1‑69 (MQ=255)
cAGCGGCCAGCAGCCCAGCGTACGGAAACGCACCATCCGTTTTTTAATCACTTCGCCCGGTTGCAGGTc                                                                     <  1:2896096/69‑1 (MQ=255)
cAGCGGCCAGCAGCCCAGCGTACGGAAACGCACCATCCGTTTTTTAATCACTTCGCCCGGTTGCAGGTc                                                                     >  1:466569/1‑69 (MQ=255)
cAGCGGCCAGCAGCCCAGCGTACGGAAACGCACCATCCGTTTTTTAATCACTTCGCCCGGTTGCAGGTc                                                                     <  1:1392534/69‑1 (MQ=255)
cAGCGGCCAGCAGCCCAGCGTACGGAAACGCACCATCCGTTTTTTAATCACTTCGCCCGGTTGCAGGTc                                                                     <  1:2103721/69‑1 (MQ=255)
cAGCGGCCAGCAGCCCAGCGTACGGAAACGCACCATCCGTTTTTTAATCACTTCGCCCGGTTGCAGGTc                                                                     <  1:2703886/69‑1 (MQ=255)
cAGCGGCCAGCAGCCCAGCGTACGGAAACGCACCATCCGTTTTTTAATCACTTCGCCCGGTTGCAGGTc                                                                     <  1:3113942/69‑1 (MQ=255)
                   gTACGGAAACGCACCATCCGTTTTTTAATCACTTCGCCCGGTTGCAGGTCGATACGGTTGTCatcaatc                                                  >  1:377943/1‑69 (MQ=255)
                   gTACGGAAACGCACCATCCGTTTTTTAATCACTTCGCCCGGTTGCAGGTCGATACGGTTGTCatcaatc                                                  >  1:122899/1‑69 (MQ=255)
                   gTACGGAAACGCACCATCCGTTTTTTAATCACTTCGCCCGGTTGCAGGTCGATACGGTTGTCatcaat                                                   >  1:521187/1‑68 (MQ=255)
                        gAAACGCACCATCCGTTTTTTAATCACTTCGCCCGGTTGCAGGTCGATACGGTTGTCATCAatcatcat                                             <  1:975768/69‑1 (MQ=255)
                        gAAACGCACCATCCGTTTTTTAATCACTTCGCCCGGTTGCAGGTCGATACGGTTGTCATCAatcatcat                                             <  1:223156/69‑1 (MQ=255)
                           aCGCACCATCCGTTTTTTAATCACTTCGCCCGGTTGCAGGTCGATACGGTTGTCATCAATCATCATCaa                                          <  1:2691315/69‑1 (MQ=255)
                                 cATCCGTTTTTTAATCACTTCGCCCGGTTGCAGGTCGATACGGTTGTCATCAatcatcat                                             <  1:1132807/60‑1 (MQ=255)
                                     cGTTTTTTAATCACTTCGCCCGGTTGCAGGTCGATACGGTTGTCATCAATCATCATCAACATACCGTcg                                >  1:1602721/1‑69 (MQ=255)
                                     cGTTTTTTAATCACTTCGCCCGGTTGCAGGTCGATACGGTTGTCATCAATCATCATCAACATACCGTcg                                >  1:897425/1‑69 (MQ=255)
                                     cGTTTTTTAATCACTTCGCCCGGTTGCAGGTCGATACGGTTGTCATCAATCATCATCAACATACCGTcg                                >  1:2704760/1‑69 (MQ=255)
                                     cGTTTTTTAATCACTTCGCCCGGTTGCAGGTCGATACGGTTGTCATCAATCATCATCAACATACCGTcg                                >  1:2858576/1‑69 (MQ=255)
                                       ttttttAATCACTTCGCCCGATTGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTcgcg                              <  1:2895425/69‑1 (MQ=255)
                                       ttttttAATCACTTCGCCCGATTGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTcgcg                              <  1:749442/69‑1 (MQ=255)
                                       ttttttAATCACTTCGCCCGATTGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTcgcg                              <  1:19629/69‑1 (MQ=255)
                                       ttttttAATCACTTCGCCCGATTGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTcgcg                              <  1:2088270/69‑1 (MQ=255)
                                       ttttttAATCACTTCGCCCGATTGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTcgcg                              <  1:2886547/69‑1 (MQ=255)
                                        tttttAATCACTTCGCCCGATTGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTcgcg                              <  1:3100896/68‑1 (MQ=255)
                                               tCACTTCGCCCGGTTGCAGGTCGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGa                       <  1:858823/68‑1 (MQ=255)
                                               tCACTTCGCCCGGTTGCAGGTCGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGa                       <  1:388012/68‑1 (MQ=255)
                                               tCACTTCGCCCGGTTGCAGGTCGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGa                       <  1:3202641/68‑1 (MQ=255)
                                                cACTTCGCCCGGTTGCAGGTCGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGa                       <  1:1289153/67‑1 (MQ=255)
                                                     cGCCCGATTGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTcgcg                              <  1:2146372/55‑1 (MQ=255)
                                                       cccGCTTGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTCCAGAACCGGACg               >  1:1335153/1‑68 (MQ=255)
                                                       cccGATTGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTCCAGAACCGGAc                >  1:1788984/1‑67 (MQ=255)
                                                       cccGATTGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTCCAGAACCGGACg               >  1:2310490/1‑68 (MQ=255)
                                                       cccGATTGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTCCAGAACCGGACg               >  1:2411541/1‑68 (MQ=255)
                                                       cccGATTGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTCCAGAACCGGACg               >  1:553793/1‑68 (MQ=255)
                                                       cccGATTGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTCCAGAACCGGACg               >  1:2615273/1‑68 (MQ=255)
                                                       cccGATTGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTCCAGAACCGGACg               >  1:2858479/1‑68 (MQ=255)
                                                       cccGATTGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTCCAGAACCGGACg               >  1:220109/1‑68 (MQ=255)
                                                       cccGATTGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTCCAGAACCGGACg               >  1:635553/1‑68 (MQ=255)
                                                       cccGATTGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTCCAGAACCGGACg               >  1:1965156/1‑68 (MQ=255)
                                                       cccGATTGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTCCAGAACCGGACg               >  1:833126/1‑68 (MQ=255)
                                                       cccGATTGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTCCAGAACCGGACg               >  1:846047/1‑68 (MQ=255)
                                                       cccGATTGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTCCAGAACCGGACg               >  1:1521937/1‑68 (MQ=255)
                                                       cccGATTGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTCCAGAACCGGACg               >  1:1339143/1‑68 (MQ=255)
                                                       cccGATTGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTCCAGAACCGGACg               >  1:987563/1‑68 (MQ=255)
                                                       cccGATTGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTCCAGAACCGGACg               >  1:1264233/1‑68 (MQ=255)
                                                       cccGATTGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTCCAGAACCGGACg               >  1:120868/1‑68 (MQ=255)
                                                       cccGATTGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTCCAGAACCGGACg               >  1:101849/1‑68 (MQ=255)
                                                       cccGATTGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTACAGAACCGGACg               >  1:1519777/1‑68 (MQ=255)
                                                          ggTTGCAGGTCGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTc            >  1:793767/1‑68 (MQ=255)
                                                          gattGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGtt                            >  1:2774254/3‑52 (MQ=255)
                                                          gattGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTc                           >  1:2548974/3‑53 (MQ=255)
                                                          gattGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTCCAGAACCGGACGtt             >  1:990389/3‑67 (MQ=255)
                                                          gattGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTCCAGAACCGGACGtt             >  1:2072251/3‑67 (MQ=255)
                                                          gattGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTCCAGAACCGGACGtt             >  1:2630485/3‑67 (MQ=255)
                                                          gattGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTCCAGAACCGGACGtt             >  1:578148/3‑67 (MQ=255)
                                                          gattGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTCCAGAACCGGACGtt             >  1:2297671/3‑67 (MQ=255)
                                                          gattGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTCCAGAACCGGACGtt             >  1:1972634/3‑67 (MQ=255)
                                                          gattGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTCCAGAACCGGACGtt             >  1:3428392/3‑67 (MQ=255)
                                                          gattGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTCCAGAACCGGACGtt             >  1:1866992/3‑67 (MQ=255)
                                                          gattGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTCCAGAACCGGACGtt             >  1:1659684/3‑67 (MQ=255)
                                                          gattGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTCCAGAACCGGACGtt             >  1:1709409/3‑67 (MQ=255)
                                                          gattGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTCCAGAACCGGACGtt             >  1:3391723/3‑67 (MQ=255)
                                                          gattGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTCCAGAACCGGACGt              >  1:3228848/3‑66 (MQ=255)
                                                          gattGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTCCAGAACCGGACGTTc            >  1:432895/3‑68 (MQ=255)
                                                          gattGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTAGCGTTCCAGAACCGGAc                >  1:2370415/3‑64 (MQ=255)
                                                                 ggTCGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCga    <  1:1816165/69‑1 (MQ=255)
                                                                 ggTCGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCga    <  1:841110/69‑1 (MQ=255)
                                                                 ggTCGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCga    <  1:1911265/69‑1 (MQ=255)
                                                                 ggTCGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCga    <  1:864866/69‑1 (MQ=255)
                                                                 ggTCGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCga    <  1:878464/69‑1 (MQ=255)
                                                                 ggTCGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCga    <  1:886641/69‑1 (MQ=255)
                                                                 ggTCGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCga    <  1:1408448/69‑1 (MQ=255)
                                                                 ggTCGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCga    <  1:1336236/69‑1 (MQ=255)
                                                                 ggTCGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCga    <  1:447119/69‑1 (MQ=255)
                                                                 ggTCGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCga    <  1:1010045/69‑1 (MQ=255)
                                                                 ggTCGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCga    <  1:2814777/69‑1 (MQ=255)
                                                                 ggTCGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCga    <  1:274074/69‑1 (MQ=255)
                                                                 ggTCGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCga    <  1:294564/69‑1 (MQ=255)
                                                                 ggTCGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCga    <  1:314461/69‑1 (MQ=255)
                                                                 ggTCGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCga    <  1:2295882/69‑1 (MQ=255)
                                                                 ggTCGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCga    <  1:373128/69‑1 (MQ=255)
                                                                 ggTCGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCga    <  1:3404124/69‑1 (MQ=255)
                                                                  gTCGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCgag   <  1:3161041/69‑1 (MQ=255)
                                                                   tCGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCgaga  <  1:3202175/69‑1 (MQ=255)
                                                                    cGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGCACCGGACGTTCCGCAGCgaga  >  1:733206/1‑68 (MQ=255)
                                                                    cGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACg               <  1:79632/55‑1 (MQ=255)
                                                                    cGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACg               <  1:3393690/55‑1 (MQ=255)
                                                                    cGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACg               <  1:2142513/55‑1 (MQ=255)
                                                                    cGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCgaga  >  1:414398/1‑68 (MQ=255)
                                                                    cGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCgaga  >  1:2302711/1‑68 (MQ=255)
                                                                    cGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCgaga  >  1:2002357/1‑68 (MQ=255)
                                                                    cGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCgaga  >  1:2295556/1‑68 (MQ=255)
                                                                    cGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCgaga  >  1:2439147/1‑68 (MQ=255)
                                                                    cGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCgaga  >  1:3042495/1‑68 (MQ=255)
                                                                    cGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCgaga  >  1:2963415/1‑68 (MQ=255)
                                                                    cGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCgaga  >  1:409437/1‑68 (MQ=255)
                                                                    cGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCgaga  >  1:412582/1‑68 (MQ=255)
                                                                    cGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCgag   >  1:1297572/1‑67 (MQ=255)
                                                                    |                                                                   
CAGCGGCCAGCAGCCCAGCGTACGGAAACGCACCATCCGTTTTTTAATCACTTCGCCCGGTTGCAGGTCGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCGAGA  >  minE/1845233‑1845368

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: