Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,063,791 T→A 12.5% intergenic (+141/+21) 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,7910TA12.5% 192.4 / 14.3 73intergenic (+141/+21)sodB/ydhPsuperoxide dismutase, Fe/predicted transporter
Reads supporting (aligned to +/- strand):  ref base T (34/29);  new base A (3/6);  total (37/36)
Fisher's exact test for biased strand distribution p-value = 3.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATA  >  minE/1063721‑1063860
                                                                      |                                                                     
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:1522982/1‑49 (MQ=37)
                            aGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGATAACCTTCt                                                                 >  1:1750660/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)
                                    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:3064568/1‑71 (MQ=255)
                                           tGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTg                            >  1:2937436/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:1620879/71‑1 (MQ=255)
                                                       tAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGgttgtt                <  1:1852419/71‑1 (MQ=255)
                                                       tAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGgttgtt                <  1:1908488/71‑1 (MQ=255)
                                                       tAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGgttg                  <  1:545053/69‑1 (MQ=255)
                                                       tAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGgttg                  <  1:3893902/69‑1 (MQ=255)
                                                       tAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGgttg                  <  1:1629615/69‑1 (MQ=255)
                                                       tAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGgttg                  <  1:4253078/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)
                                                                    tATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACa    >  1:3957188/1‑70 (MQ=255)
                                                                    tATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACa    >  1:524136/1‑70 (MQ=255)
                                                                     aTCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACAt   >  1:40132/1‑70 (MQ=255)
                                                                     aTCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACAt   >  1:722605/1‑70 (MQ=255)
                                                                     aTCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACAt   >  1:727941/1‑70 (MQ=255)
                                                                     aTCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACAt   >  1:1345275/1‑70 (MQ=255)
                                                                     aTCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACAt   >  1:2725324/1‑70 (MQ=255)
                                                                     aTCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACAt   >  1:286761/1‑70 (MQ=255)
                                                                      tCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATa  <  1:2123951/70‑1 (MQ=255)
                                                                      tCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATa  <  1:4030315/70‑1 (MQ=255)
                                                                      tCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATa  <  1:3902751/70‑1 (MQ=255)
                                                                      tCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATa  <  1:3484673/70‑1 (MQ=255)
                                                                      |                                                                     
TCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATA  >  minE/1063721‑1063860

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: