Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,063,788 T→A 11.6% intergenic (+138/+24) sodB → / ← ydhP superoxide dismutase, Fe/predicted transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,063,7880TA11.6% 192.1 / 18.8 78intergenic (+138/+24)sodB/ydhPsuperoxide dismutase, Fe/predicted transporter
Reads supporting (aligned to +/- strand):  ref base T (26/42);  new base A (3/6);  total (30/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.15e-01

AGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGA  >  minE/1063718‑1063856
                                                                      |                                                                    
aGGTCGGATAAGGAACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGt                                                                      >  1:2617983/1‑71 (MQ=255)
 ggTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGtt                                                                     <  1:4126829/71‑1 (MQ=255)
   tCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTAt                                                                   >  1:477898/1‑71 (MQ=255)
    cGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGAAACGCAGAAGGTTATc                                                                  >  1:512098/1‑71 (MQ=255)
     ggATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATcc                                                                 <  1:1608416/71‑1 (MQ=255)
     ggATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATcc                                                                 <  1:3646879/71‑1 (MQ=255)
      gATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCt                                                                >  1:2234043/1‑71 (MQ=255)
       aTAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCtt                                                               >  1:887475/1‑71 (MQ=255)
       aTAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCtt                                                               >  1:1967622/1‑71 (MQ=255)
        tAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTc                                                              >  1:2754802/1‑71 (MQ=255)
         aaGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCtt                                                               >  1:2108447/1‑69 (MQ=255)
         aaGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCt                                                             >  1:2698783/1‑71 (MQ=255)
               cTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTACCCTTCTGCGttt                                                       <  1:1647735/71‑1 (MQ=255)
                 cgccgcATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGttttt                                                     >  1:2115702/1‑71 (MQ=255)
                 cgccgcATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGttttt                                                     >  1:3606099/1‑71 (MQ=255)
                 cgccgcATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGttttt                                                     >  1:653407/1‑71 (MQ=255)
                 cgccgcATCCGGCAAGATAAATCGCACGTTTCAGCAACTGTAACGCAGAAGGCTATCCTTCTGCGTTTTTg                                                    >  1:2125879/1‑71 (MQ=39)
                  gccgcATCCGGAAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTg                                                    >  1:3421303/1‑71 (MQ=255)
                     gcATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGttt                                                 <  1:170757/71‑1 (MQ=255)
                       aTCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTaa                                               <  1:4003008/71‑1 (MQ=255)
                          cGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAAtt                                             <  1:592307/70‑1 (MQ=255)
                          cGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAAt                                              >  1:1566966/1‑69 (MQ=255)
                             cAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGATAACCtt                                                               <  1:520621/49‑1 (MQ=38)
                             cAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGATAACCtt                                                               <  1:665775/49‑1 (MQ=38)
                             cAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGATAACCtt                                                               <  1:950027/49‑1 (MQ=38)
                             cAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGATAACCtt                                                               <  1:1010487/49‑1 (MQ=38)
                             cAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGATAACCtt                                                               <  1:1070316/49‑1 (MQ=38)
                               aGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGtt                                                  <  1:2177105/60‑1 (MQ=255)
                               aGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGtt                                                  <  1:1482440/60‑1 (MQ=255)
                               aGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGtt                                                  <  1:423835/60‑1 (MQ=255)
                               aGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGATAACCtt                                                               <  1:3376079/47‑1 (MQ=37)
                               aGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGATAACCTTCt                                                             >  1:2332126/1‑49 (MQ=37)
                               aGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGATAACCTTCt                                                             >  1:1750660/1‑49 (MQ=37)
                               aGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGATAACCTTCt                                                             >  1:1522982/1‑49 (MQ=37)
                                gATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGt                                       <  1:486713/70‑1 (MQ=255)
                                   aaaTCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAg                                    >  1:1477351/1‑70 (MQ=255)
                                    aaTCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTTTCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCa                                  <  1:1941930/71‑1 (MQ=255)
                                    aaTCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCa                                  <  1:2525638/71‑1 (MQ=255)
                                      tCGCACGTTTTCAGCAACTGTAACGCAGAAGGTTa                                                                    <  1:590531/35‑1 (MQ=255)
                                      tCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTa                                                                    <  1:1596328/35‑1 (MQ=255)
                                      tCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTa                                                                    <  1:1591212/35‑1 (MQ=255)
                                      tCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTa                                                                    <  1:3149232/35‑1 (MQ=255)
                                      tCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTa                                                                    <  1:3900919/35‑1 (MQ=255)
                                      tCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTa                                                                    <  1:3659482/35‑1 (MQ=255)
                                      tCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTa                                                                    <  1:2888532/35‑1 (MQ=255)
                                      tCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTa                                                                    <  1:271323/35‑1 (MQ=255)
                                      tCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTa                                                                    <  1:570798/35‑1 (MQ=255)
                                      tCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTa                                                                    <  1:2582395/35‑1 (MQ=255)
                                      tCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTa                                                                    <  1:2575252/35‑1 (MQ=255)
                                      tCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTa                                                                    <  1:1990061/35‑1 (MQ=255)
                                       cGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGttt                                                       <  1:1006311/47‑1 (MQ=255)
                                       cGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGc                                   <  1:3913394/67‑1 (MQ=255)
                                       cGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACg                               >  1:2152988/1‑71 (MQ=255)
                                         cACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCa                             <  1:2386859/71‑1 (MQ=255)
                                         cACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCa                             <  1:4044740/71‑1 (MQ=255)
                                         cACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCa                             <  1:427142/71‑1 (MQ=255)
                                          aCGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGtt                                                        <  1:255517/43‑1 (MQ=255)
                                            gTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCa                             <  1:1218940/68‑1 (MQ=255)
                                              tGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTg                        >  1:2937436/1‑71 (MQ=255)
                                              tGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTg                        >  1:3064568/1‑71 (MQ=255)
                                               gTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGt                       >  1:1863353/1‑71 (MQ=255)
                                               gTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGt                       >  1:3915736/1‑71 (MQ=255)
                                                tCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGtt                      >  1:3754742/1‑71 (MQ=255)
                                                  aGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTc                    <  1:2713146/71‑1 (MQ=255)
                                                    cAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAg                  >  1:4139487/1‑71 (MQ=255)
                                                    cAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAg                  >  1:2398013/1‑71 (MQ=255)
                                                     aaCTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCaa                                 >  1:1795687/1‑55 (MQ=255)
                                                     aaCTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCa                                  >  1:3685501/1‑54 (MQ=255)
                                                         gTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAg                  >  1:630544/1‑66 (MQ=255)
                                                         gTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAg                  >  1:2959892/1‑66 (MQ=255)
                                                          tAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGgttgtt            <  1:1852419/71‑1 (MQ=255)
                                                          tAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGgttgtt            <  1:1908488/71‑1 (MQ=255)
                                                          tAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGgttgtt            <  1:1620879/71‑1 (MQ=255)
                                                          tAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGgttg              <  1:545053/69‑1 (MQ=255)
                                                          tAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGgttg              <  1:4253078/69‑1 (MQ=255)
                                                          tAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGgttg              <  1:1629615/69‑1 (MQ=255)
                                                          tAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGgttg              <  1:3893902/69‑1 (MQ=255)
                                                           aaCGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGgttgtt            <  1:820361/70‑1 (MQ=255)
                                                                  aaGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCt    <  1:128056/71‑1 (MQ=255)
                                                                    ggTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGa  <  1:379993/71‑1 (MQ=255)
                                                                      |                                                                    
AGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGA  >  minE/1063718‑1063856

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: