Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,678,841 T→C 36.6% I617V (ATA→GTA)  feoB ← fused ferrous iron transporter, protein B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,678,8410TC36.6% 59.6 / 63.0 101I617V (ATA→GTA) feoBfused ferrous iron transporter, protein B
Reads supporting (aligned to +/- strand):  ref base T (46/18);  new base C (15/22);  total (61/40)
Fisher's exact test for biased strand distribution p-value = 2.94e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 101 total aligned reads displayed.

TGGGGTTCATCAGTACGCTAAGGCTGAAGGTGTCTTTCAGGCTCTGCCAGGTTTCATCTATCGCACTGAACAGCTCTTCACCGAGGTTAAATTCTGCCGGATTGAACTCT  >  minE/2678781‑2678890
                                                            |                                                 
tGGGGTTCATCAGTACGCTAAGGCTGAAGGTGTCTTTCAGGCTCTGCCAGGTTTCATCTATCGCACTg                                            <  1:301548/68‑1 (MQ=255)
tGGGGTTCATCAGTACGCTAAGGCTGAAGGTGTCTTTCAGGCTCTGCCAGGTTTCATCTATCGCACTg                                            <  1:3191415/68‑1 (MQ=255)
tGGGGTTCATCAGTACGCTAAGGCTGAAGGTGTCTTTCAGGCTCTGCCAGGTTTCATCTATCGCACTg                                            <  1:2301060/68‑1 (MQ=255)
tGGGGTTCATCAGTACGCTAAGGCTGAAGGTGTCTTTCAGGCTCTGCCAGGTTTCATCTATCGCACTg                                            <  1:1936204/68‑1 (MQ=255)
tGGGGTTCATCAGTACGCTAAGGCTGAAGGTGTCTTTCAGGCTCTGCCAGGTTTCATCTATCGCACTg                                            <  1:1746934/68‑1 (MQ=255)
tGGGGTTCATCAGTACGCTAAGGCTGAAGGTGTCTTTCAGGCTCTGCCAGGTTTCATCTATCGCACTg                                            <  1:1633604/68‑1 (MQ=255)
tGGGGTTCATCAGTACGCTAAGGCTGAAGGTGTCTTTCAGGCTCTGCCAGGTTTCATCTATCGCACTg                                            <  1:2442985/68‑1 (MQ=255)
tGGGGTTCATCAGTACGCTAAGGCTGAAGGTGTCTTTCAGGCTCTGCCAGGTTTCATCTATCGCACTg                                            <  1:2484842/68‑1 (MQ=255)
tGGGGTTCATCAGTACGCTAAGGCTGAAGGTGTCTTTCAGGCTCTGCCAGGTTTCATCTATCGCACTg                                            <  1:296655/68‑1 (MQ=255)
  gggTTCATCAGTACGCTAAGGCTGAAGGTGTCTTTCAGGCTCTGCCAGGTTTCATCTATCGc                                                <  1:3280973/62‑1 (MQ=255)
  gggTTCATCAGTACGCTAAGGCTGAAGGTGTCTTTCAGGCTCTGCCAGGTTTCATCTATCGc                                                <  1:203746/62‑1 (MQ=255)
  gggTTCATCAGTACGCTAAGGCTGAAGGTGTCTTTCAGGCTCTGCCAGGTTTCATCTATCGc                                                <  1:2978867/62‑1 (MQ=255)
  gggTTCATCAGTACGCTAAGGCTGAAGGTGTCTTTCAGGCTCTGCCAGGTTTCATCTATCGc                                                <  1:1549897/62‑1 (MQ=255)
  gggTTCATCAGTACGCTAAGGCTGAAGGTGTCTTTCAGGCTCTGCCAGGTTTCATCTATCGc                                                <  1:1270424/62‑1 (MQ=255)
  gggTTCATCAGTACGCTAAGGCTGAAGGTGTCTTTCAGGCTCTGCCAGGTTTCATCTATCGCACTGAAc                                         >  1:3133439/1‑69 (MQ=255)
  gggTTCATCAGTACGCTAAGGCTGAAGGTGTCTTTCAGGCTCTGCCAGGTTTCATCTATCGCACTGAAc                                         >  1:2672679/1‑69 (MQ=255)
  gggTTCATCAGTACGCTAAGGCTGAAGGTGTCTTTCAGGCTCTGCCAGGTTTCATCTATCGCACTGAAc                                         >  1:370723/1‑69 (MQ=255)
  gggTTCATCAGTACGCTAAGGCTGAAGGTGTCTTTCAGGCTCTGCCAGGTTTCATCTATCGCACTGAAc                                         >  1:2326199/1‑69 (MQ=255)
              aCGCTAAGGCTGAAGGTGTCTTTCAGGCTCTGCCAGGTTTCATCTATCg                                                 >  1:1696461/1‑49 (MQ=255)
              aCGCTAAGGCTGAAGGTGTCTTTCAGGCTCTGCCAGGTTTCATCTATCg                                                 >  1:1198385/1‑49 (MQ=255)
              aCGCTAAGGCTGAAGGTGTCTTTCAGGCTCTGCCAGGTTTCATCTATCg                                                 >  1:3433288/1‑49 (MQ=255)
              aCGCTAAGGCTGAAGGTGTCTTTCAGGCTCTGCCAGGTTTCATCTATCg                                                 >  1:490123/1‑49 (MQ=255)
              aCGCTAAGGCTGAAGGTGTCTTTCAGGCTCTGCCAGGTTTCATCTATCg                                                 >  1:2665934/1‑49 (MQ=255)
              aCGCTAAGGCTGAAGGTGTCTTTCAGGCTCTGCCAGGTTTCATCTATCg                                                 >  1:1686786/1‑49 (MQ=255)
              aCGCTAAGGCTGAAGGTGTCTTTCAGGCTCTGCCAGGTTTCATCTATCg                                                 >  1:2112838/1‑49 (MQ=255)
              aCGCTAAGGCTGAAGGTGTCTTTCAGGCTCTGCCAGGTTTCATCTATCg                                                 >  1:761466/1‑49 (MQ=255)
              aCGCTAAGGCTGAAGGTGTCTTTCAGGCTCTGCCAGGTTTCATCTATCg                                                 >  1:1521033/1‑49 (MQ=255)
              aCGCTAAGGCTGAAGGTGTCTTTCAGGCTCTGCCAGGTTTCATCTATCg                                                 >  1:1282098/1‑49 (MQ=255)
               cGCTAAGGCTGAAGGTGTCTTTCAGGCTCTGCCAGGTTTCATCTATCGCACTGAACAGCTCTTCACCGa                            >  1:1810509/1‑69 (MQ=255)
               cGCTAAGGCTGAAGGTGTCTTTCAGGCTCTGCCAGGTTTCATCTATCGCACTGAACAGCTCTTCACCGa                            >  1:1586638/1‑69 (MQ=255)
               cGCTAAGGCTGAAGGTGTCTTTCAGGCTCTGCCAGGTTTCATCTATCGCACTGAACAGCTCTTCACCGa                            >  1:2925204/1‑69 (MQ=255)
               cGCTAAGGCTGAAGGTGTCTTTCAGGCTCTGCCAGGTTTCATCTATCGCACTGAACAGATCTTCACCGa                            >  1:562020/1‑69 (MQ=255)
                      gCTGAAGGTGTCTTTCAGGCTCTGCCAGGTTTCATCTATCGCACTGAACAGCTCTTCACCGAGGTTaaa                     >  1:1021824/1‑69 (MQ=255)
                      gCTGAAGGTGTCTTTCAGGCTCTGCCAGGTTTCATCTATCGCACTGAACAGCTCTTCACCGAGGTTaa                      >  1:238936/1‑68 (MQ=255)
                      gCTGAAGGTGTCTTTCAGGCTCTGCCAGGTTTCATCTACCGCACTGACCAGCTCTTCACCGAGGTTaaa                     <  1:274098/69‑1 (MQ=255)
                      gCTGAAGGTGTCTTTCAGGCTCTGCCAGGTTTCATCTACCGCACTGAACAGCTCTTCACCGAGGTTaaa                     <  1:2202948/69‑1 (MQ=255)
                      gCTGAAGGTGTCTTTCAGGCTCTGCCAGGTTTCATCTACCGCACTGAACAGCTCTTCACCGAGGTTaaa                     <  1:952855/69‑1 (MQ=255)
                      gCTGAAGGTGTCTTTCAGGCTCTGCCAGGTTTCATCTACCGCACTGAACAGCTCTTCACCGAGGTTaaa                     <  1:1260587/69‑1 (MQ=255)
                      gCTGAAGGTGTCTTTCAGGCTCTGCCAGGTTTCATCTACCGCACTGAACAGCTCTTCACCGAGGTTaaa                     <  1:2611976/69‑1 (MQ=255)
                      gCTGAAGGTGTCTTTCAGGCTCTGCCAGGTTTCATCTACCGCACTGAACAGCTCTTCACCGAGGTTaaa                     <  1:2857383/69‑1 (MQ=255)
                      gCTGAAGGTGTCTTTCAGGCTCTGCCAGGTTTCATCTACCGCACTGAACAGCTCTTCACCGAGGTTaaa                     <  1:947456/69‑1 (MQ=255)
                      gCTGAAGGTGTCTTTCAGGCTCTGCCAGGTTTCATCTACCGCACTGAACAGCTCTTCACCGAGGTTaaa                     <  1:778645/69‑1 (MQ=255)
                      gCTGAAGGTGTCTTTCAGGCTCTGCCAGGTTTCATCTACCGCACTGAACAGCTCTTCACCGAGGTTaaa                     <  1:3044981/69‑1 (MQ=255)
                      gCTGAAGGTGTCTTTCAGGCTCTGCCAGGTTTCATCTACCGCACTGAACAGCTCTTCACCGAGGTTa                       <  1:2351611/67‑1 (MQ=255)
                      gCTGAAGGTGTCTTTCAGGCTCTGCCAGGTTTCATCTACCGCACTGAACAGCTCTTCACCGAGGTTa                       <  1:926340/67‑1 (MQ=255)
                      gCTGAAGGTGTCTTTCAGGCTCTGCCAGGTTTCATCTACCGCACTGAACAGCTCTTCACCGAGGTTa                       <  1:2441556/67‑1 (MQ=255)
                      gCTGAAGGTGTCTTTCAGGCTCTGCCAGGTTTCATCTACCGCACTGAACAGCTCTTCACCGAGGTTa                       <  1:45841/67‑1 (MQ=255)
                      gCTGAAGGTGTCTTTCAGGCTCTGCCAGGTTTCATCTACCGCACTGAACAGCTCTTCACCGAGGTTa                       <  1:952479/67‑1 (MQ=255)
                      gCTGAAGGTGTCTTTCAGGCTCTGCCAGGTTTCATCTACCGCACTGAACAGCTCTTCACCGAGGTTa                       <  1:1271446/67‑1 (MQ=255)
                      gCTGAAGGTGTCTTTCAGGCTCTGCCAGGTTTCATCTACCGCACTGAACAGCTCTTCACCGAGGTTa                       <  1:1271198/67‑1 (MQ=255)
                      gCTGAAGGTGTCTTTCAGGCTCTGCCAGGTTTCATCTACCGCACTGAACAGCTCTTCACCGAGGTTa                       <  1:1984077/67‑1 (MQ=255)
                      gCTGAAGGTGTCTTTCAGGCTCTGCCAGGTTTCATCTACCGCACTGAACAGCTCTTCACCGAGGTTa                       <  1:2665898/67‑1 (MQ=255)
                      gCTGAAGGTGTCTTTCAGGCTCTGCCAGGTTTCATCTACCGCACTGAACAGCTCTTCACCGAGGTTa                       <  1:3202328/67‑1 (MQ=255)
                      gCTGAAGGTGTCTTTCAGGCTCTGCCAGGTTTCATCTACCGCACTGAACAGCTCTTCACCGAGGTTa                       <  1:3144509/67‑1 (MQ=255)
                       cTGAAGGTGTCTTTCAGGCTCTGCCAGGTTTCATCTACCGCACTGAACAGCTCTTCACCGAGGTTa                       <  1:2598075/66‑1 (MQ=255)
                              tgtCTTTCAGGCTCTGCCAGGTTTCATCTATCGCACTGAACAGCTCTTCACCGAGGTTAAATTCTGCCg             >  1:1631564/1‑69 (MQ=255)
                              tgtCTTTCAGGCTCTGCCAGGTTTCATCTATCGCACTGAACAGCTCTTCACCGAGGTTAAATTCTGCCg             >  1:2505346/1‑69 (MQ=255)
                               gtCTTTCAGGCTCTGCCAGGTTTCATCTATCGCACTGAACAGCTCTTCACCGAGGTTAAAt                    >  1:1066391/1‑61 (MQ=255)
                               gtCTTTCAGGCTCTGCCAGGTTTCATCTATCGCACTGAACAGCTCTTCACCGAGGTTAAAt                    >  1:2722494/1‑61 (MQ=255)
                               gtCTTTCAGGCTCTGCCAGGTTTCATCTATCGCACTGAACAGCTCTTCACCGAGGTTAAAt                    >  1:1699019/1‑61 (MQ=255)
                               gtCTTTCAGGCTCTGCCAGGTTTCATCTATCGCACTGAACAGCTCTTCACCGAGGTTAAAt                    >  1:2174496/1‑61 (MQ=255)
                               gtCTTTCAGGCTCTGCCAGGTTTCATCTATCGCACTGAACAGCTCTTCACCGAGGTTAAAt                    >  1:1446761/1‑61 (MQ=255)
                               gtCTTTCAGGCTCTGCCAGGTTTCATCTATCGCACTGAACAGCTCTTCACCGAGGTTAAATTCTGCCgg            >  1:61548/1‑69 (MQ=255)
                               gtCTTTCAGGCTCTGCCAGGTTTCATCTATCGCACTGAACAGCTCTTCACCGAGGTTAAATTCTGCCgg            >  1:908493/1‑69 (MQ=255)
                               gtCTTTCAGGCTCTGCCAGGTTTCATCTATCGCACTGAAAAGCTCTTCACCGAGGTTAAATTCTGCCgg            >  1:1608703/1‑69 (MQ=255)
                                  tttCAGGCTCTGCCAGGTTTCATCTATCGCACTGAACAGCTCTTCACCGAGGTTAAATTCTGCCGGAt          >  1:829156/1‑68 (MQ=255)
                                  tttCAGGCTCTGCCAGGTTTCATCTATCGCACTGAACAGCTCTTCACCGAGGTTAAATTCTGCCGGAt          >  1:1840923/1‑68 (MQ=255)
                                  tttCAGGCTCTGCCAGGTTTCATCTATCGCACTGAACAGCTCTTCACCGAGGTTAAATTCTGCCGGAt          >  1:257038/1‑68 (MQ=255)
                                      aGGCTCTGCCAGGTTTCATCTATCGCACTGAACAGCTCTTCACCGAGGTTAAATTCTGCCGGATTGaa      <  1:2252989/68‑1 (MQ=255)
                                      aGGCTCTGCCAGGTTTCATCTATCGCACTGAACAGCTCTTCACCGAGGTTAAATTCTGCCGGATTGaa      <  1:775809/68‑1 (MQ=255)
                                         ctctGCCAGGTTTCATCTATCGCACTGAACAGCTCTTCACCGAGGTTAAATTCTGCCGGATTGAActct  >  1:289461/1‑69 (MQ=255)
                                            tGCCAGGTTTCATCTATCGCACTGAACAGCTCTTCACCGAGGTTAAATTCTGcag             <  1:3079173/55‑3 (MQ=255)
                                            tGCCAGGTTTCATCTATCGCACTGAACAGCTCTTCACCGAGGTTAAATTCTGcag             <  1:3264546/55‑3 (MQ=255)
                                            tGCCAGGTTTCATCTATCGCACTGAACAGCTCTTCACCGAGGTTAAATTCTGCCg             >  1:656256/1‑55 (MQ=255)
                                            tGCCAGGTTTCATCTATCGCACTGAACAGCTCTTCACCGAGGTTAAATTCTGCCg             >  1:3417416/1‑55 (MQ=255)
                                                     tCATCTATCGCACTGAACAGCTCTTCACCGAGGTTa                       >  1:827595/1‑36 (MQ=255)
                                                     tCATCTATCGCACTGAACAGCTCTTCACCGAGGTTa                       >  1:1817447/1‑36 (MQ=255)
                                                     tCATCTATCGCACTGAACAGCTCTTCACCGAGGTTa                       >  1:3359838/1‑36 (MQ=255)
                                                     tCATCTATCGCACTGAACAGCTCTTCACCGAGGTTa                       >  1:2099529/1‑36 (MQ=255)
                                                     tCATCTATCGCACTGAACAGCTCTTCACCGAGGTTa                       >  1:2128209/1‑36 (MQ=255)
                                                     tCATCTATCGCACTGAACAGCTCTTCACCGAGGTTa                       >  1:2178593/1‑36 (MQ=255)
                                                     tCATCTATCGCACTGAACAGCTCTTCACCGAGGTTa                       >  1:3058115/1‑36 (MQ=255)
                                                     tCATCTATCGCACTGAACAGCTCTTCACCGAGGTTa                       >  1:2278388/1‑36 (MQ=255)
                                                     tCATCTATCGCACTGAACAGCTCTTCACCGAGGTTa                       >  1:1180163/1‑36 (MQ=255)
                                                     tCATCTATCGCACTGAACAGCTCTTAACCGAGGTTa                       >  1:1500445/1‑36 (MQ=255)
                                                       aTCTACCGCACTGAACAGCTCTTCACCGAGGTTAAATTCCGCCGGATTg        >  1:81861/1‑49 (MQ=38)
                                                       aTCTACCGCACTGAACAGCTCTTCACCGAGGTTAAATTCCGCCGGATTg        >  1:2760117/1‑49 (MQ=38)
                                                       aTCTACCGCACTGAACAGCTCTTCACCGAGGTTAAATTCCGCCGGATTg        >  1:1443741/1‑49 (MQ=38)
                                                       aTCTACCGCACTGAACAGCTCTTCACCGAGGTTAAATTCCGCCGGATTg        >  1:1271805/1‑49 (MQ=38)
                                                       aTCTACCGCACTGAACAGCTCTTCACCGAGGTTAAATTCCGCCGGATTg        >  1:179694/1‑49 (MQ=38)
                                                       aTCTACCGCACTGAACAGCTCTTCACCGAGGTTAAATTCCGCCGGATTg        >  1:1915493/1‑49 (MQ=38)
                                                       aTCTACCGCACTGAACAGCTCTTCACCGAGGTTAAATTCCGCCGGATTg        >  1:3298207/1‑49 (MQ=38)
                                                       aTCTACCGCACTGAACAGCTCTTCACCGAGGTTAAATTCCGCCGGATTg        >  1:3245218/1‑49 (MQ=38)
                                                       aTCTACCGCACTGAACAGCTCTTCACCGAGGTTAAATTCCGCCGGATTg        >  1:3015792/1‑49 (MQ=38)
                                                       aTCTACCGCACTGAACAGCTCTTCACCGAGGTTAAATTCCGCCGGATTg        >  1:2380382/1‑49 (MQ=38)
                                                       aTCTACCGCACTGAACAGCTCTTCACCGAGGTTAAATTCCGCCGGATTg        >  1:2394555/1‑49 (MQ=38)
                                                       aTCTACCGCACTGAACAGCTCTTCACCGAGGTTAAATTCCGCCGGATTg        >  1:248948/1‑49 (MQ=38)
                                                       aTCTACCGCACTGAACAGCTCTTCACCGAGGTTAAATTCCGCCGGATTg        >  1:2554774/1‑49 (MQ=38)
                                                       aTCTACCGCACTGAACAGCTCTTCACCGAGGTTAAATTCCGCCGGATTg        >  1:2879922/1‑49 (MQ=38)
                                                       aTCTACCGCACTGAACAGCTCTTCACAGAGGTTAAATTACGCCGGATTg        >  1:3351751/1‑49 (MQ=255)
                                                            |                                                 
TGGGGTTCATCAGTACGCTAAGGCTGAAGGTGTCTTTCAGGCTCTGCCAGGTTTCATCTATCGCACTGAACAGCTCTTCACCGAGGTTAAATTCTGCCGGATTGAACTCT  >  minE/2678781‑2678890

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