Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,678,875 T→C 23.9% A605A (GCA→GCG 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,8750TC23.9% 92.9 / 23.7 75A605A (GCA→GCGfeoBfused ferrous iron transporter, protein B
Reads supporting (aligned to +/- strand):  ref base T (31/26);  new base C (15/3);  total (46/29)
Fisher's exact test for biased strand distribution p-value = 4.99e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTCTTTCAGGCTCTGCCAGGTTTCATCTATCGCACTGAACAGCTCTTCACCGAGGTTAAATTCTGCCGGATTGAACTCTTCGTCCTGAATATTTTCTGCGGTGTAGAGGGTGTTGAGCGTACCCAC  >  minE/2678811‑2678937
                                                                |                                                              
tgtCTTTCAGGCTCTGCCAGGTTTCATCTATCGCACTGAACAGCTCTTCACCGAGGTTAAATTCTGCCg                                                            >  1:2505346/1‑69 (MQ=255)
tgtCTTTCAGGCTCTGCCAGGTTTCATCTATCGCACTGAACAGCTCTTCACCGAGGTTAAATTCTGCCg                                                            >  1:1631564/1‑69 (MQ=255)
 gtCTTTCAGGCTCTGCCAGGTTTCATCTATCGCACTGAACAGCTCTTCACCGAGGTTAAATTCTGCCgg                                                           >  1:908493/1‑69 (MQ=255)
 gtCTTTCAGGCTCTGCCAGGTTTCATCTATCGCACTGAACAGCTCTTCACCGAGGTTAAATTCTGCCgg                                                           >  1:61548/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:3417416/1‑55 (MQ=255)
              tGCCAGGTTTCATCTATCGCACTGAACAGCTCTTCACCGAGGTTAAATTCTGCCg                                                            >  1:656256/1‑55 (MQ=255)
                         aTCTACCGCACTGAACAGCTCTTCACCGAGGTTAAATTCCGCCGGATTg                                                       >  1:3015792/1‑49 (MQ=38)
                         aTCTACCGCACTGAACAGCTCTTCACCGAGGTTAAATTCCGCCGGATTg                                                       >  1:2760117/1‑49 (MQ=38)
                         aTCTACCGCACTGAACAGCTCTTCACCGAGGTTAAATTCCGCCGGATTg                                                       >  1:248948/1‑49 (MQ=38)
                         aTCTACCGCACTGAACAGCTCTTCACCGAGGTTAAATTCCGCCGGATTg                                                       >  1:1271805/1‑49 (MQ=38)
                         aTCTACCGCACTGAACAGCTCTTCACCGAGGTTAAATTCCGCCGGATTg                                                       >  1:2394555/1‑49 (MQ=38)
                         aTCTACCGCACTGAACAGCTCTTCACCGAGGTTAAATTCCGCCGGATTg                                                       >  1:2380382/1‑49 (MQ=38)
                         aTCTACCGCACTGAACAGCTCTTCACCGAGGTTAAATTCCGCCGGATTg                                                       >  1:2554774/1‑49 (MQ=38)
                         aTCTACCGCACTGAACAGCTCTTCACCGAGGTTAAATTCCGCCGGATTg                                                       >  1:3245218/1‑49 (MQ=38)
                         aTCTACCGCACTGAACAGCTCTTCACCGAGGTTAAATTCCGCCGGATTg                                                       >  1:81861/1‑49 (MQ=38)
                         aTCTACCGCACTGAACAGCTCTTCACCGAGGTTAAATTCCGCCGGATTg                                                       >  1:2879922/1‑49 (MQ=38)
                         aTCTACCGCACTGAACAGCTCTTCACCGAGGTTAAATTCCGCCGGATTg                                                       >  1:1443741/1‑49 (MQ=38)
                         aTCTACCGCACTGAACAGCTCTTCACCGAGGTTAAATTCCGCCGGATTg                                                       >  1:3298207/1‑49 (MQ=38)
                         aTCTACCGCACTGAACAGCTCTTCACCGAGGTTAAATTCCGCCGGATTg                                                       >  1:1915493/1‑49 (MQ=38)
                         aTCTACCGCACTGAACAGCTCTTCACCGAGGTTAAATTCCGCCGGATTg                                                       >  1:179694/1‑49 (MQ=38)
                         aTCTACCGCACTGAACAGCTCTTCACAGAGGTTAAATTACGCCGGATTg                                                       >  1:3351751/1‑49 (MQ=255)
                                  aCTGAACAGCTCTTCACCGAGGTTAAATTCTGCCGGATTGAACTCTTCGTCCTGAATATTTTCTGCgg                           <  1:1469646/68‑1 (MQ=255)
                                  aCTGAACAGCTCTTCACCGAGGTTAAATTCTGCCGGATTGAACTCTTCGTCCTGAATATTTTCTGCgg                           <  1:3021121/68‑1 (MQ=255)
                                  aCTGAACAGCTCTTCACCGAGGTTAAATTCTGCCGGATTGAACTCTTCGTCCTGAATATTTGCTGCgg                           <  1:442721/68‑1 (MQ=255)
                                     gAACAGCTCTTCACCGAGGTTAAATTCTGCCGGATTGAACTCTTCGTCCTGAATATTTTCTGCGgtg                         >  1:2783528/1‑67 (MQ=255)
                                     gAACAGCTCTTCACCGAGGTTAAATTCTGCCGGATTGAACTCTTCGTCCTGAATATTTTCTGCGgtg                         >  1:1603186/1‑67 (MQ=255)
                                     gAACAGCTCTTCACCGAGGTTAAATTCTGCCGGATTGAACTCTTCGTCCTGAATATTTTCTGCGGTGTa                       <  1:3187166/69‑1 (MQ=255)
                                     gAACAGCTCTTCACCGAGGTTAAATTCTGCCGGATTGAACTCTTCGTCCTGAATATTTTCTGCGGTGTa                       <  1:2218099/69‑1 (MQ=255)
                                     gAACAGCTCTTCACCGAGGTTAAATTCTGCCGGATTGAACTCTTCGTCCTGAATATTTTCTGCGGTGTa                       <  1:445059/69‑1 (MQ=255)
                                     gAACAGCTCTTCACCGAGGTTAAATTCTGCCGGATTGAACTCTTCGTCCTGAATATTTTCTGCGGTGTa                       <  1:991276/69‑1 (MQ=255)
                                       aCAGCTCTTCACCGAGGTTAAATTCTGCCGGATTGAACTCTTCGTCCTGAATATTTTCTGCGGTGTa                       <  1:435153/67‑1 (MQ=255)
                                         aGCTCTTCACCGAGGTTAAATTCTGCCGGATTGAACTCTTCGTCCTGAATATTTTCTGCGGTGTagag                    <  1:2576217/68‑1 (MQ=255)
                                           ctctTCACCGAGGTTAAATTCTGCCGGATTGAACTCTTCGTCCTGAATATTTTCTGCGGTGTAGAGGgt                 >  1:1425884/1‑69 (MQ=255)
                                           ctctTCACCGAGGTTAAATTCTGCCGGATTGAACTCTTCGTCCTGAATATTTTCTGCGGTGTAGAGGgt                 >  1:2942041/1‑69 (MQ=255)
                                             ctTCACCGAGGTTAAATTCTGCCGGATTGAACTCTTCGTCCTGAata                                     >  1:1724474/1‑47 (MQ=255)
                                             ctTCACCGAGGTTAAATTCTGCCGGATTGAACTCTTCGTCCTGAata                                     >  1:1439063/1‑47 (MQ=255)
                                             ctTCACCGAGGTTAAATTCTGCCGGATTGAACTCTTCGTCCTGAata                                     >  1:2946998/1‑47 (MQ=255)
                                                    gAGGTTAAATTCTGCCGGATTGAACTCTTCGTCCTGAATATTTTCTGCGGTGTAGAGGGTGTTGAGCGt        <  1:1208058/69‑1 (MQ=255)
                                                    gAGGTTAAATTCTGCCGGATTGAACTCTTCGTCCTGAATATTTTCTGCGGTGTAGAGGGTGTTGAGCGt        <  1:1283044/69‑1 (MQ=255)
                                                    gAGGTTAAATTCTGCCGGATTGAACTCTTCGTCCTGAATATTTTCTGCGGTGTAGAGGGTGTTGAGCGt        <  1:3011870/69‑1 (MQ=255)
                                                    gAGGTTAAATTCTGCCGGATTGAACTCTTCGTCCTGAATATTTTCTGCGGTGTAGAGGGTGTTGAGCGt        <  1:1261904/69‑1 (MQ=255)
                                                    gAGGTTAAATTCTGCCGGATTGAACTCTTCGTCCTGAATATTTTCTGCGGTGTAGAGGGTGTTGAGCGt        <  1:2430793/69‑1 (MQ=255)
                                                    gAGGTTAAATTCTGCCGGATTGAACTCTTCGTCCTGAATATTTTCTGCGGTGTAGAGGGTGTTGAGCGt        <  1:2254562/69‑1 (MQ=255)
                                                    gAGGTTAAATTCTGCCGGATTGAACTCTTCGTCCTGAATATTTTCTGCGGTGTAGAGGGTGTTGAGCGt        <  1:1521653/69‑1 (MQ=255)
                                                    gAGGTTAAATTCTGCCGGATTGAACTCTTCGTCCTGAATATTTTCTGCGGTGTAGAGGGTGTTGAGCGt        <  1:1636737/69‑1 (MQ=255)
                                                    gAGGTTAAATTCTGCCGGATTGAACTCTTCGTCCTGAATATTTTCTGCGGTGTAGAGGGTGTTGAGCGt        <  1:199324/69‑1 (MQ=255)
                                                       gTTAAATTCTGCCGGATTGAACTCTTCGTCCTGAATATTTTCTGCGGTGTAGAGGGTGTTGAGCGTAcc     <  1:954986/69‑1 (MQ=255)
                                                       gTTAAATTCTGCCGGATTGAACTCTTCGTCCTGAATATTTTCTGCGGTGTAGAGGGTGTTGAGCGTAcc     <  1:2272675/69‑1 (MQ=255)
                                                       gTTAAATTCTGCCGGATTGAACTCTTCGTCCTGAATATTTTCTGAGGTGTAGAGGGTGTTGAGCGTAcc     <  1:2982356/69‑1 (MQ=255)
                                                        ttAAATTCTGCCGGATTGAACTCTTCGTCCTGAATATTTTCTGCgg                           <  1:2877201/46‑1 (MQ=255)
                                                          aaaTTCTGCCGGATTGAACTCTTCGTCCTGAATATTTTCTGCGGTGTAGAGGGTGTTGAGCGTAcc     >  1:1965891/1‑66 (MQ=255)
                                                          aaaTTCTGCCGGATTGAACTCTTCGTCCTGAATATTTTCTGCGGTGTAGAGGGTGTTGAGCGTAcc     >  1:1693724/1‑66 (MQ=255)
                                                          aaaTTCTGCCGGATTGAACTCTTCGTCCTGAATATTTTCTGCGGTGTAGAGGGTGTTGAGCGTAcc     >  1:3246726/1‑66 (MQ=255)
                                                          aaaTTCTGCCGGATTGAACTCTTCGTCCTGAATATTTTCTGCGGTGTAGAGGGTGTTGAGCGTAcc     >  1:511224/1‑66 (MQ=255)
                                                          aaaTTCTGCCGGATTGAACTCTTCGTCCTGAATATTTTCTGCGGTGTAGAGGGTGTTGAGCGTAcc     >  1:526842/1‑66 (MQ=255)
                                                          aaaTTCTGCCGGATTGAACTCTTCGTCCTGAATATTTTCTGCGGTGTAGAGGGTGTTGAGCGTAcc     >  1:1488287/1‑66 (MQ=255)
                                                          aaaTTCTGCCGGATTGAACTCTTCGTCCTGAATATTTTCTGCGGTGTAGAGGGTGTTGAGCGTAcc     >  1:2475332/1‑66 (MQ=255)
                                                          aaaTTCTGCCGGATTGAACTCTTCGTCCTGAATATTTTCTGCGGTGTAGAGGGTGTTGAGCGTAcc     >  1:1408624/1‑66 (MQ=255)
                                                          aaaTTCTGCCGGATTGAACTCTTCGTCCTGAATATTTTCTGCGGTGTAGAGGGTGTTGAGCGTAcc     >  1:2976213/1‑66 (MQ=255)
                                                          aaaTTCTGCCGGATTGAACTCTTCGTCCTGAATATTTTCTGCGGTGTAGAGGGTGTTGAGCGTAcc     >  1:2686618/1‑66 (MQ=255)
                                                          aaaTTCTGCCGGATTGAACTCTTCGTCCTGAATATTTTCTGCGGTGTAGAGGGTGTTGAGCGTAcc     >  1:2807156/1‑66 (MQ=255)
                                                          aaaTTCTGCCGGATTGAACTCTTCGTCCTGAATATTTTCTGCGGTGTAGAGGGTGTTGAGCGTACCCAc  >  1:3328463/1‑69 (MQ=255)
                                                          aaaTTCTGCCGGATTGAACTCTTCGTCCTGAATATTTTCTGCGGTGTAGAGGGTGTTGAGCGTACCCAc  >  1:1501633/1‑69 (MQ=255)
                                                          aaaTTCCGCCGGATTGAACTCTTCGTCCTGAATATTTTCTGCGGTGTAGAGGGTGTTGAGCGTACCCAc  <  1:2207901/69‑1 (MQ=255)
                                                          aaaTTCCGCCGGATTGAACTCTTCGTCCTGAATATTTTCTGCGGTGTAGAGGGTGTTGAGCGTACCCAc  <  1:550128/69‑1 (MQ=255)
                                                          aaaTTCCGCCGGATTGAACTCTTCGTCCTGAATATTTTCTGCGGTGTAGAGGGTGTTGAGCGTACCCAc  <  1:2555415/69‑1 (MQ=255)
                                                                |                                                              
TGTCTTTCAGGCTCTGCCAGGTTTCATCTATCGCACTGAACAGCTCTTCACCGAGGTTAAATTCTGCCGGATTGAACTCTTCGTCCTGAATATTTTCTGCGGTGTAGAGGGTGTTGAGCGTACCCAC  >  minE/2678811‑2678937

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: