Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,051,975 A→C 26.5% A124A (GCT→GCG sodC ← superoxide dismutase, Cu, Zn

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,051,9750AC26.5% 104.5 / 46.5 83A124A (GCT→GCGsodCsuperoxide dismutase, Cu, Zn
Reads supporting (aligned to +/- strand):  ref base A (25/36);  new base C (8/14);  total (33/50)
Fisher's exact test for biased strand distribution p-value = 8.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGG  >  minE/1051918‑1052037
                                                         |                                                              
cgcTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGt                                                      <  1:3121010/68‑1 (MQ=255)
cgcTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGt                                                      <  1:2338668/68‑1 (MQ=255)
cgcTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGt                                                      <  1:2385897/68‑1 (MQ=255)
cgcTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGt                                                      <  1:503978/68‑1 (MQ=255)
cgcTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGt                                                      <  1:3421540/68‑1 (MQ=255)
cgcTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGt                                                      <  1:2062288/68‑1 (MQ=255)
cgcTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGt                                                      <  1:2036149/68‑1 (MQ=255)
cgcTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGt                                                      <  1:2531139/68‑1 (MQ=255)
cgcTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGt                                                      <  1:1911090/68‑1 (MQ=255)
cgcTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGt                                                      <  1:2650736/68‑1 (MQ=255)
cgcTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGt                                                      <  1:1564300/68‑1 (MQ=255)
cgcTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGt                                                      <  1:1496333/68‑1 (MQ=255)
cgcTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGt                                                      <  1:815599/68‑1 (MQ=255)
cgcTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGt                                                      <  1:796161/68‑1 (MQ=255)
cgcTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGt                                                      <  1:1257591/68‑1 (MQ=255)
   tttGTCTTTGACTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTCGCCTTGCCGTCatt                                                  <  1:1725280/69‑1 (MQ=255)
   tttGTCTTTGACTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTCGCCTTGCCGTCatt                                                  <  1:2761798/69‑1 (MQ=255)
   tttGTCTTTGACTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTCGCCTTGCCGTCatt                                                  <  1:610412/69‑1 (MQ=255)
   tttGTCTTTGACTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTCGCCTTGCCGTCatt                                                  <  1:2610942/69‑1 (MQ=255)
   tttGTCTTTGACTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTCGCCTTGCCGTCatt                                                  <  1:1814832/69‑1 (MQ=255)
   tttGTCTTTGACTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTCGCCTTGCCGTCatt                                                  <  1:723576/69‑1 (MQ=255)
   tttGTCTTTGACTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTCGCCTTGCCGTCatt                                                  <  1:1717766/69‑1 (MQ=255)
   tttGTCTTTGACTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTCGCCTTGCCGTCatt                                                  <  1:1195800/69‑1 (MQ=255)
   tttGTCTTTGACTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTCGCCTTGCCGTCatt                                                  <  1:1052654/69‑1 (MQ=255)
   tttGTCTTTGACTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTCGCCTTGCCGTCatt                                                  <  1:465934/69‑1 (MQ=255)
   tttGTCTTTGACTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTCGCCTTGCCGTCatt                                                  <  1:2790474/69‑1 (MQ=255)
   tttGTCTTTGACTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTCGCCTTGCCGTCatt                                                  <  1:1423938/69‑1 (MQ=255)
    ttGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCatta                                                 >  1:655789/1‑69 (MQ=255)
    ttGTCTTTGACTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTCGCCTTGCCGTCatt                                                  <  1:2859549/68‑1 (MQ=255)
                      aGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGt                                                      <  1:511459/46‑1 (MQ=255)
                        tGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTgcaggcag                             <  1:78885/69‑1 (MQ=255)
                        tGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTgcaggcag                             <  1:2018190/69‑1 (MQ=255)
                        tGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTgcaggcag                             <  1:1213350/69‑1 (MQ=255)
                        tGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTgcaggcag                             <  1:1710060/69‑1 (MQ=255)
                        tGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTgcaggcag                             <  1:879414/69‑1 (MQ=255)
                          aTTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGAt                           <  1:1825047/69‑1 (MQ=255)
                          aTTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGAt                           <  1:1807332/69‑1 (MQ=255)
                             tCAGACGAGGCGCGATGACGGCATCGGTAGCTTTg                                                          >  1:1606623/1‑35 (MQ=255)
                             tCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTgac                                            >  1:2601056/1‑49 (MQ=255)
                             tCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTgac                                            >  1:3371443/1‑49 (MQ=255)
                             tCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTgac                                            >  1:2476357/1‑49 (MQ=255)
                             tCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTgac                                            >  1:2708189/1‑49 (MQ=255)
                             tCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTgac                                            >  1:2712116/1‑49 (MQ=255)
                             tCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTgac                                            >  1:3057257/1‑49 (MQ=255)
                             tCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTgac                                            >  1:2256628/1‑49 (MQ=255)
                             tCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTgac                                            >  1:795795/1‑49 (MQ=255)
                             tCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTgac                                            >  1:1663197/1‑49 (MQ=255)
                             tCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTgac                                            >  1:1705076/1‑49 (MQ=255)
                             tCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTgac                                            >  1:1865927/1‑49 (MQ=255)
                             tCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTgac                                            >  1:535955/1‑49 (MQ=255)
                             tCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTgac                                            >  1:1960221/1‑49 (MQ=255)
                             tCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTgac                                            >  1:204182/1‑49 (MQ=255)
                             tCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTga                                             >  1:2062642/1‑48 (MQ=255)
                             tCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTgca                                  >  1:1959020/1‑59 (MQ=255)
                             tCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTgca                                  >  1:1404001/1‑59 (MQ=255)
                             tCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTgca                                  >  1:915776/1‑59 (MQ=255)
                             tCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTgca                                  >  1:3085868/1‑59 (MQ=255)
                             tCAGACGAGGCGCCATGACGGCATCGGTAGCTTTGCCGTCATTATTgac                                            >  1:1443583/1‑49 (MQ=255)
                               aGACGAGGCGCGATGACGGCATCGGTCGCCTTGCCGTCatt                                                  <  1:3350893/41‑1 (MQ=38)
                                 aCGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTTCAGGCAGATCGCCTa                     <  1:1849029/68‑1 (MQ=255)
                                 aCGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTa                     <  1:41662/68‑1 (MQ=255)
                                 aCGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTa                     <  1:2287629/68‑1 (MQ=255)
                                 aCGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTa                     <  1:2932557/68‑1 (MQ=255)
                                 aCGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTa                     <  1:2963912/68‑1 (MQ=255)
                                 aCGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTa                     <  1:2753016/68‑1 (MQ=255)
                                 aCGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTa                     <  1:1366610/68‑1 (MQ=255)
                                 aCGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTa                     <  1:179715/68‑1 (MQ=255)
                                 aCGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTa                     <  1:2002611/68‑1 (MQ=255)
                                 aCGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTa                     <  1:180250/68‑1 (MQ=255)
                                 aCGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTa                     <  1:1848738/68‑1 (MQ=255)
                                 aCGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTa                     <  1:1809792/68‑1 (MQ=255)
                                           aTGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTa                     <  1:758030/58‑1 (MQ=255)
                                                ggCATCGGTCGCCTTGCCGTCATTATTGACTACCa                                       >  1:2474577/1‑35 (MQ=21)
                                                ggCATCGGTCGCCTTGCCGTCATTATTGACTACCAGTgca                                  >  1:103408/1‑40 (MQ=25)
                                                ggCATCGGTCGCCTTGCCGTCATTATTGACTACCAGTgca                                  >  1:88191/1‑40 (MQ=25)
                                                ggCATCGGTCGCCTTGCCGTCATTATTGACTACCAGTgca                                  >  1:1057259/1‑40 (MQ=25)
                                                ggCATCGGTCGCCTTGCCGTCATTATTGACTACCAGTgca                                  >  1:1880384/1‑40 (MQ=25)
                                                ggCATCGGTCGCCTTGCCGTCATTATTGACTACCAGTgca                                  >  1:1030232/1‑40 (MQ=25)
                                                ggCATCGGTCGCCTTGCCGTCATTATTGACTACCAGTgca                                  >  1:2333446/1‑40 (MQ=25)
                                                ggCATCGGTCGCCTTGCCGTCATTATTGACTACCAGTgca                                  >  1:2893577/1‑40 (MQ=25)
                                                   aTCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTgg  >  1:217521/1‑69 (MQ=255)
                                                   aTCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTgg  >  1:1277583/1‑69 (MQ=255)
                                                   aTCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTgg  >  1:3390090/1‑69 (MQ=255)
                                                         aGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGc           >  1:1142765/1‑54 (MQ=255)
                                                         aGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGc           >  1:2547158/1‑54 (MQ=255)
                                                         |                                                              
CGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGG  >  minE/1051918‑1052037

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: