Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,660,974 C→T 24.8% A217A (GCG→GCA yfgF ← predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,660,9740CT24.8% 168.4 / 47.2 132A217A (GCG→GCAyfgFpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (53/46);  new base T (14/19);  total (67/65)
Fisher's exact test for biased strand distribution p-value = 3.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.74e-01

Only 100 of 132 total aligned reads displayed.

TTTTCATTTAACGGCATGCATAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGTTTTAATTGCGAATAATAGCTA  >  minE/1660915‑1661034
                                                           |                                                            
ttttCATTTAACGGCATGCATAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGATTGCGa                                                          <  1:3050008/64‑1 (MQ=255)
ttttCATTTAACGGCATGCATAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGATTGCGa                                                          <  1:1712178/64‑1 (MQ=255)
ttttCATTTAACGGCATGCATAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGATTGCGa                                                          <  1:1810507/64‑1 (MQ=255)
ttttCATTTAACGGCATGCATAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGATTGCGa                                                          <  1:739885/64‑1 (MQ=255)
ttttCATTTAACGGCATGCATAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGATTGCGa                                                          <  1:741080/64‑1 (MQ=255)
ttttCATTTAACGGCATGCATAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGATTGCGa                                                          <  1:1825925/64‑1 (MQ=255)
ttttCATTTAACGGCATGCATAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGATTGCGa                                                          <  1:598092/64‑1 (MQ=255)
ttttCATTTAACGGCATGCATAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGATTGCGa                                                          <  1:496725/64‑1 (MQ=255)
ttttCATTTAACGGCATGCATAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGATTGCGa                                                          <  1:413019/64‑1 (MQ=255)
ttttCATTTAACGGCATGCATAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGATTGCGa                                                          <  1:358847/64‑1 (MQ=255)
ttttCATTTAACGGCATGCATAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGATTGCGa                                                          <  1:347161/64‑1 (MQ=255)
ttttCATTTAACGGCATGCATAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGATTGCGa                                                          <  1:2706014/64‑1 (MQ=255)
ttttCATTTAACGGCATGCATAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGATTGCGa                                                          <  1:2767580/64‑1 (MQ=255)
            ggCATGCATAACAGCAATAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGAc                                          <  1:1667603/68‑1 (MQ=255)
                tgcatAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTgg                                     <  1:2281853/69‑1 (MQ=255)
                tgcatAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTgg                                     <  1:402822/69‑1 (MQ=255)
                tgcatAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTgg                                     <  1:3353533/69‑1 (MQ=255)
                tgcatAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTgg                                     <  1:1805614/69‑1 (MQ=255)
                tgcatAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTgg                                     <  1:1017782/69‑1 (MQ=255)
                tgcatAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTgg                                     <  1:1636422/69‑1 (MQ=255)
                tgcatAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTgg                                     <  1:1144348/69‑1 (MQ=255)
                tgcatAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTgg                                     <  1:1440342/69‑1 (MQ=255)
                  catAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCtttttt                                                <  1:248219/56‑1 (MQ=255)
                  catAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCtttttt                                                <  1:2372782/56‑1 (MQ=255)
                  catAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCtttttt                                                <  1:2937113/56‑1 (MQ=255)
                  catAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCtttttt                                                <  1:2237673/56‑1 (MQ=255)
                  catAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGc                                    >  1:675683/1‑68 (MQ=255)
                  catAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGc                                    >  1:2838105/1‑68 (MQ=255)
                  catAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGc                                    >  1:1125128/1‑68 (MQ=255)
                  catAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGc                                    >  1:253517/1‑68 (MQ=255)
                  catAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGc                                    >  1:1209412/1‑68 (MQ=255)
                  catAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGc                                    >  1:2296602/1‑68 (MQ=255)
                  catAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGc                                    >  1:2884247/1‑68 (MQ=255)
                  catAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGc                                    >  1:2151139/1‑68 (MQ=255)
                  catAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGc                                    >  1:1275840/1‑68 (MQ=255)
                  catAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGc                                    >  1:2008165/1‑68 (MQ=255)
                  catAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGc                                    >  1:774824/1‑68 (MQ=255)
                  catAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGc                                    >  1:611472/1‑68 (MQ=255)
                  catAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGc                                    >  1:3108693/1‑68 (MQ=255)
                   atAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGATTGCGa                                                          <  1:673590/45‑1 (MQ=255)
                     aaCAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCa                                <  1:2881877/69‑1 (MQ=255)
                     aaCAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCa                                <  1:151686/69‑1 (MQ=255)
                     aaCAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCa                                <  1:1736476/69‑1 (MQ=255)
                      aCAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCttttt                                                 >  1:2822347/1‑51 (MQ=255)
                      aCAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCttttt                                                 >  1:2083692/1‑51 (MQ=255)
                      aCAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCttttt                                                 >  1:1008702/1‑51 (MQ=255)
                      aCAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCttttt                                                 >  1:1030814/1‑51 (MQ=255)
                      aCAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCaa                               >  1:1611819/1‑69 (MQ=255)
                      aCAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCaa                               >  1:2879874/1‑69 (MQ=255)
                      aCAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCa                                >  1:2415306/1‑68 (MQ=255)
                      aCAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCa                                <  1:865453/68‑1 (MQ=255)
                      aCAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCa                                <  1:2909898/68‑1 (MQ=255)
                         gCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGCTGACTTTGGCATCAACct                            >  1:3161847/1‑69 (MQ=255)
                                tAAAGCACCTAATGCCAGTAGCCAGATTGCGAACTCTTTTTTGGTAACCTTGGCATCAACCTGCTGtt                      >  1:618787/1‑68 (MQ=255)
                                tAAAGCACCTAATGCCAGTAGCCAGATTGCGAACTCTTTTTTGGTAACCTTGGCATCAACCTGCTGtt                      >  1:3058766/1‑68 (MQ=255)
                                tAAAGCACCTAATGCCAGTAGCCAGATTGCGAACTCTTTTTTGGTAACCTTGGCATCAACCTGCTGtt                      >  1:1233808/1‑68 (MQ=255)
                                tAAAGCACCTAATGCCAGAAGCCAGATTGCGAACTCTTTTTTGGt                                             >  1:572495/1‑45 (MQ=255)
                                 aaaGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGtttt                    <  1:2369356/69‑1 (MQ=255)
                                 aaaGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGtttt                    <  1:3204566/69‑1 (MQ=255)
                                 aaaGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGtttt                    <  1:1073442/69‑1 (MQ=255)
                                 aaaGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGtttt                    <  1:488819/69‑1 (MQ=255)
                                 aaaGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGtttt                    <  1:1985655/69‑1 (MQ=255)
                                 aaaGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGtttt                    <  1:1394112/69‑1 (MQ=255)
                                 aaaGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGtttt                    <  1:698414/69‑1 (MQ=255)
                                 aaaGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGtttt                    <  1:1674663/69‑1 (MQ=255)
                                  aaGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGtttt                    <  1:3347195/68‑1 (MQ=255)
                                          aaTGCCAGTAGCCAGATTGCGAACTCTTTTTTGGTAACCTTGGCATCAACctgct                         >  1:3421934/1‑55 (MQ=255)
                                          aaTGCCAGTAGCCAGATTGCGAACTCTTTTTTGGTAACCTTGGCATCAACctgct                         >  1:97728/1‑55 (MQ=255)
                                          aaTGCCAGTAGCCAGATTGCGAACTCTTTTTTGGTAACCTTGGCATCAACctgct                         >  1:2808050/1‑55 (MQ=255)
                                          aaTGCCAGTAGCCAGATTGCGAACTCTTTTTTGGTAACCTTGGCATCAACctgct                         >  1:1891365/1‑55 (MQ=255)
                                          aaTGCCAGTAGCCAGATTGCGAACTCTTTTTTGGTAACCTTGGCATCAACctgct                         >  1:276299/1‑55 (MQ=255)
                                          aaTGCCAGTAGCCAGATTGCGAACTCTTTTTTGGTAACCTTGGCATCAACctgct                         >  1:2346295/1‑55 (MQ=255)
                                          aaTGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGTTTTAATTGCg             >  1:1435460/1‑67 (MQ=255)
                                          aaTGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGTTTTAATTGCGaa           <  1:1511451/69‑1 (MQ=255)
                                          aaTGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGTTTTAATTGCGaa           <  1:2830383/69‑1 (MQ=255)
                                          aaTGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGTTTTAATTGCGaa           <  1:3276285/69‑1 (MQ=255)
                                          aaTGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGTTTTAATTGCGaa           <  1:2111107/69‑1 (MQ=255)
                                          aaTGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGTTTTAATTGCGa            >  1:2672859/1‑68 (MQ=255)
                                          aaTGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGTTTTAATTGCGa            >  1:3430600/1‑68 (MQ=255)
                                           aTGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGTTTTAATTGCGaa           <  1:425764/68‑1 (MQ=255)
                                               cAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCa                                   >  1:2954506/1‑40 (MQ=255)
                                               cAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACctgatg                        >  1:3160898/1‑51 (MQ=255)
                                               cAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGTTTTAATTGCGaataat       >  1:770735/1‑68 (MQ=255)
                                               cAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGTTTTAATTGCGAataata      >  1:1029704/1‑69 (MQ=255)
                                               cAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGTTTTAATTGCGAataata      >  1:1707436/1‑69 (MQ=255)
                                               cAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGTTTTAATTGCGAataata      >  1:189861/1‑69 (MQ=255)
                                                 gTAGCCAGAGCGAGAACTCTTTTTTGGTGACTTTgg                                     <  1:184842/36‑1 (MQ=255)
                                                    gCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATc                                 >  1:1999638/1‑37 (MQ=255)
                                                    gCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCa                                >  1:3212930/1‑38 (MQ=255)
                                                    gCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCa                                >  1:2481012/1‑38 (MQ=255)
                                                    gCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCa                                >  1:2266237/1‑38 (MQ=255)
                                                    gCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCa                                >  1:1548447/1‑38 (MQ=255)
                                                    gCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCa                                >  1:2685427/1‑38 (MQ=255)
                                                    gCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCa                                <  1:3036051/38‑1 (MQ=255)
                                                    gCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCa                                >  1:14264/1‑38 (MQ=255)
                                                    gCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCa                                >  1:1247457/1‑38 (MQ=255)
                                                    gCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGTTTTAATTGCGAATAATAGCTa  >  1:3307676/1‑68 (MQ=255)
                                                    gCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGTTTTAATTGCGAATAATAGCTa  >  1:2337709/1‑68 (MQ=255)
                                                    gCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGTTTTAATTGCGAATAATAGCTa  >  1:1342782/1‑68 (MQ=255)
                                                    gCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGTTTTAATTGCGAATAATAGCTa  >  1:2662690/1‑68 (MQ=255)
                                                           |                                                            
TTTTCATTTAACGGCATGCATAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGTTTTAATTGCGAATAATAGCTA  >  minE/1660915‑1661034

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: