Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,678,500 C→T 24.6% A730A (GCG→GCA 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,5000CT24.6% 123.1 / 32.6 97A730A (GCG→GCAfeoBfused ferrous iron transporter, protein B
Reads supporting (aligned to +/- strand):  ref base C (40/33);  new base T (10/14);  total (50/47)
Fisher's exact test for biased strand distribution p-value = 3.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGATATCCACCCGGCTACGCGCGCGGCGCAGCAGACCGATAACCACGATGTTAAACAGGATAACCGCCAGAATGCACACCAGGCTGTAAGTTGGATGCTGACTGTAGCTGGCGACCTGATAGAACAATG  >  minE/2678436‑2678564
                                                                |                                                                
cGACATCTACCCGGCTACGCGCGCGGCGCAGCAGACCGATAACCACGATGTTAAACAGGATAACTGcc                                                               >  1:787601/1‑68 (MQ=255)
cGACATCTACCCGGCTACGCGCGCGGCGCAGCAGACCGATAACCACGATGTTAAACAGGATAACTGcc                                                               >  1:310191/1‑68 (MQ=255)
cGACATCTACCCGGCTACGCGCGCGGCGCAGCAGACCGATAACCACGATGTTAAACAGGATAACTGcc                                                               >  1:2017241/1‑68 (MQ=255)
cGACATCTACCCGGCTACGCGCGCGGCGCAGCAGACCGATAACCACGATGTTAAACAGGATAACTGcc                                                               >  1:3044112/1‑68 (MQ=255)
cGACATCTACCCGGCTACGCGCGCGGCGCAGCAGACCGATAACCACGATGTTAAACAGGATAACTGcc                                                               >  1:526199/1‑68 (MQ=255)
cGACATCTACCCGGCTACGCGCGCGGCGCAGCAGACCGATAACCACGATGTTAAACAGGATAACTGcc                                                               >  1:1353612/1‑68 (MQ=255)
cGACATCTACCCGGCTACGCGCGCGGCGCAGCAGACCGATAACCACGATGTTAAACAGGATAACTGcc                                                               >  1:2759899/1‑68 (MQ=255)
cGACATCTACCCGGCTACGCGCGCGGCGCAGCAGACCGATAACCACGATGTTAAACAGGATAACTGcc                                                               >  1:2091660/1‑68 (MQ=255)
cGACATCTACCCGGCTACGCGCGCGGCGCAGCAGACCGATAACCACGATGTTAAACAGGATAACTGcc                                                               >  1:1435405/1‑68 (MQ=255)
cGACATCTACCCGGCTACGCGCGCGGCGCAGCAGACCGATAACCACGATGTTAAACAGGATAACTGcc                                                               >  1:2857587/1‑68 (MQ=255)
  atatCCACCCGGCTACGCGCGCGGCGCAGCAGACCGATAACCACGATGTTAAACAGGATAACCGCCAg                                                             <  1:3046060/68‑1 (MQ=255)
  atatCCACCCGGCTACGCGCGCGGCGCAGCAGACCGATAACCACGATGTTAAACAGGATAACCGCCAg                                                             <  1:3234908/68‑1 (MQ=255)
  atatCCACCCGGCTACGCGCGCGGCGCAGCAGACCGATAACCACGATGTTAAACAGGATAACCGCCAg                                                             <  1:592509/68‑1 (MQ=255)
  atatCCACCCGGCTACGCGCGCGGCGCAGCAGACCGATAACCACGATGTTAAACAGGATAACCGCCAg                                                             <  1:1460038/68‑1 (MQ=255)
  atatCCACCCGGCTACGCGCGCGGCGCAGCAGACCGATAACCACGATGTTAAACAGGATAACCGCCAg                                                             <  1:2178971/68‑1 (MQ=255)
  atatCCACCCGGCTACGCGCGCGGCGCAGCAGACCGATAACCACGATGTTAAACAGGATAACCGCCAg                                                             <  1:765625/68‑1 (MQ=255)
   tatCCACCCGGCTACGCGCGCGGCGCAGCAGACCGATAACCACGATGTTAAACAGGATAACCGCCAg                                                             <  1:425169/67‑1 (MQ=255)
   tatCCACCCGGCTACGCGCGCGGCGCAGCAGACCGATAACCACGATGTTAAACAGGATAACCGCCAGaa                                                           <  1:1387149/69‑1 (MQ=255)
   tatCCACCCGGCTACGCGCGCGGCGCAGCAGACCGATAACCACGATGTTAAACAGGATAACCGCCAGaa                                                           <  1:1728701/69‑1 (MQ=255)
   tatCCACCCGGCTACGCGCGCGGCGCAGCAGACCGATAACCACGATGTTAAACAGGATAACCGCCAGaa                                                           <  1:2031386/69‑1 (MQ=255)
   tatCCACCCGGCTACGCGCGCGGCGCAGCAGACCGATAACCACGATGTTAAACAGGATAACCGCCAGaa                                                           <  1:1855724/69‑1 (MQ=255)
   tatCCACCCGGCTACGCGCGCGGCGCAGCAGACCGATAACCACGATGTTAAACAGGATAACCGCCAGaa                                                           <  1:12264/69‑1 (MQ=255)
   tatCCACCCGGCTACGCGCGCGGCGCAGCAGACCGATAACCACGATGTTAAACAGGATAACCGCCAGaa                                                           <  1:3325712/69‑1 (MQ=255)
     tCCACCCGGCTACGCGCGCGGCGCAGCAGACCGATAACCACGATGTTAAACAGGATAACCGCCAGaa                                                           >  1:2686028/1‑67 (MQ=255)
     tCCACCCGGCTACGCGCGCGGCGCAGCAGACCGATAACCACGATGTTAAACAGGATAACCGCCAGAATg                                                         >  1:2772111/1‑69 (MQ=255)
     tCCACCCGGCTACGCGCGCGGCGCAGCAGACCGATAACCACGATGTTAAACAGGATAACCGCCAGAATg                                                         <  1:2946656/69‑1 (MQ=255)
     tCCACCCGGCTACGCGCGCGGCGCAGCAGACCGATAACCACGATGTTAAACAGGATAACCGCCAGAATg                                                         <  1:2277041/69‑1 (MQ=255)
     tCCACCCGGCTACGCGCGCGGCGCAGCAGACCGATAACCACGATGTTAAACAGGATAACCGCCAGAATg                                                         >  1:2450124/1‑69 (MQ=255)
     tCCACCCGGCTACGCGCGCGGCGCAGCAGACCGATAACCACGATGTTAAACAGGATAACCGCCAGAATg                                                         <  1:3358108/69‑1 (MQ=255)
     tCCACCCGGCTACGCGCGCGGCGCAGCAGACCGATAACCACGATGTTAAACAGGATAACCGCCAGAATg                                                         >  1:397564/1‑69 (MQ=255)
     tCCACCCGGCTACGCGCGCGGCGCAGCAGACCGATAACCACGATGTTAAACAGGATAACCGCCAGAATg                                                         <  1:1662097/69‑1 (MQ=255)
     tCCACCCGGCTACGCGCGCGGCGCAGCAGACCGATAACCACGATGTTAAACAGGATAACCGCCAGAATg                                                         <  1:62209/69‑1 (MQ=255)
     tCCACCCGGCTACGCGCGCGGCGCAGCAGACCGATAACCACGATGTTAAACAGGATAACCGCCAGAATg                                                         <  1:66266/69‑1 (MQ=255)
     tCCACCCGGCTACGCGCGCGGCGCAGCAGACCGATAACCACGATGTTAAACAGGATAACCGCCAGAATg                                                         >  1:711813/1‑69 (MQ=255)
     tCCACCCGGCTACGCGCGCGGCGCAGCAGACCGATAACCACGATGTTAAACAGGATAACCGCCAGAATg                                                         >  1:712697/1‑69 (MQ=255)
       cACCCGGCTACGCGCGCGGCGCAGCAGACCGATAACCACGATGTTAAACAGGATAACCGCCAGAATg                                                         >  1:3431299/1‑67 (MQ=255)
       cACCCGGCTACGCGCGCGGCGCAGCAGACCGATAACCACGATGTTAAACAGGATAACCGCCAGAATg                                                         >  1:592668/1‑67 (MQ=255)
       cACCCGGCTACGCGCGCGGCGCAGCAGACCGATAACCACGATGTTAAACAGGATAACCGCCAGAATg                                                         >  1:606953/1‑67 (MQ=255)
       cACCCGGCTACGCGCGCGGCGCAGCAGACCGATAACCACGATGTTAAACAGGATAACCGCCAGAATg                                                         >  1:625449/1‑67 (MQ=255)
       cACCCGGCTACGCGCGCGGCGCAGCAGACCGATAACCACGATGTTAAACAGGATAACCGCCAGAATg                                                         >  1:1014907/1‑67 (MQ=255)
       cACCCGGCTACGCGCGCGGCGCAGCAGACCGATAACCACGATGTTAAACAGGATAACCGCCAGAATg                                                         >  1:2383130/1‑67 (MQ=255)
       cACCCGGCTACGCGCGCGGCGCAGCAGACCGATAACCACGATGTTAAACAGGATAACCGCCAGAATg                                                         >  1:1001716/1‑67 (MQ=255)
           cGGCTACGCGCGCGGCGCAGCAGACCGATAACCACGATGTTAAACAGGATAACCGCCAg                                                             <  1:2090669/59‑1 (MQ=255)
               tACGCGCGCGGCGCAGCAGACCGATAACCACGATGTTAAACAGGATAACCGCCAGAATGCACACCAGGc                                               >  1:1444779/1‑69 (MQ=255)
               tACGCGCGCGGCGCAGCAGACCGATAACCACGATGTTAAACAGGATAACCGCCAGAATGCACACCAGGc                                               >  1:2651750/1‑69 (MQ=255)
               tACGCGCGCGGCGCAGCAGACCGATAACCACGATGTTAAACAGGATAACCGCCAGAATGCACACCAGGc                                               >  1:480383/1‑69 (MQ=255)
               tACGCGCGCGGCGCAGCAGACCGATAACCACGATGTTAAACAGGATAACCGCCAGAATGCACACCAGGc                                               >  1:2820070/1‑69 (MQ=255)
                        ggcgCAGCAGACCGATAACCACGATGTTAAACAGGATAACCGCCAGAATGCACACCAGGCTGTAAGTTg                                      <  1:207627/69‑1 (MQ=255)
                        ggcgCAGCAGACCGATAACCACGATGTTAAACAGGATAACCGCCAGAATGCACACCAGGCTGTAAGTTg                                      <  1:45773/69‑1 (MQ=255)
                        ggcgCAGCAGACCGATAACCACGATGTTAAACAGGATAACCGCCAGAATGCACACCAGGCTGTAAGTTg                                      <  1:2016341/69‑1 (MQ=255)
                             agcagACCGATAACCACGATGTTAAACAGGATAACTGCCAGAATGCACACCAGGCTGTAAGTTGATGCt                                <  1:3204598/69‑1 (MQ=255)
                             agcagACCGATAACCACGATGTTAAACAGGATAACCGCCAg                                                             <  1:765122/41‑1 (MQ=255)
                              gcagACCGATAACCACGATGTTAAACAGGATAACTGCCAGAATGCACACCAGGCTGTAAGTTGGATGCt                                <  1:3337050/69‑1 (MQ=255)
                              gcagACCGATAACCACGATGTTAAACAGGATAACTGCCAGAATGCACACCAGGCTGTAAGTTGGATGCt                                <  1:1364844/69‑1 (MQ=255)
                              gcagACCGATAACCACGATGTTAAACAGGATAACTGCCAGAATGCACACCAGGCTGTAAGTTGGATGCt                                <  1:3341143/69‑1 (MQ=255)
                              gcagACCGATAACCACGATGTTAAACAGGATAACTGCCAGAATGCACACCAGGCTGTAAGTTGGATGCt                                <  1:1859654/69‑1 (MQ=255)
                              gcagACCGATAACCACGATGTTAAACAGGATAACTGCCAGAATGCACACCAGGCTGTAAGTTGGATGCt                                <  1:2035843/69‑1 (MQ=255)
                              gcagACCGATAACCACGATGTTAAACAGGATAACTGCCAGAATGCACACCAGGCTGTAAGTTGGATGCt                                <  1:1762730/69‑1 (MQ=255)
                              gcagACCGATAACCACGATGTTAAACAGGATAACTGCCAGAATGCACACCAGGCTGTAAGTTGGATGCt                                <  1:2181736/69‑1 (MQ=255)
                              gcagACCGATAACCACGATGTTAAACAGGATAACTGCCAGAATGCACACCAGGCTGTAAGTTGGATGCt                                <  1:1613115/69‑1 (MQ=255)
                              gcagACCGATAACCACGATGTTAAACAGGATAACTGCCAGAATGCACACCAGGCTGTAAGTTGGATGCt                                <  1:2619549/69‑1 (MQ=255)
                              gcagACCGATAACCACGATGTTAAACAGGATAACTGCCAGAATGCACACCAGGCTGTAAGTTGGATGCt                                <  1:1248869/69‑1 (MQ=255)
                               cagACCGATAACCACGATGTTAAACAGGATAACTGCCAGAATGCACACCAGGCTGTAAGTTGGATGCt                                <  1:1328906/68‑1 (MQ=255)
                               cagACCGATAACCACGATGTTAAACAGGATAACTGCCAGAATGCACACCAGGCTGTAAGTTGGATGCt                                <  1:2062881/68‑1 (MQ=255)
                                   ccGATAACCACGTTGTTAAACAGGATAACCGCCAGAATGCACACCAGGCTGTAAGTTGGATGCTGAc                             >  1:3121245/1‑67 (MQ=255)
                                   ccGATAACCACGATGTTAAACAGGATAACCGCCAGAATGCACACCAGGCTGTAAg                                         >  1:1697641/1‑55 (MQ=255)
                                   ccGATAACCACGATGTTAAACAGGATAACCGCCAGAATGCACACCAGGCTGTAAGTTGGATGCTGAc                             >  1:888508/1‑67 (MQ=255)
                                   ccGATAACCACGATGTTAAACAGGATAACCGCCAGAATGCACACCAGGCTGTAAGTTGGATGCTGAc                             >  1:3424653/1‑67 (MQ=255)
                                   ccGATAACCACGATGTTAAACAGGATAACCGCCAGAATGCACACCAGGCTGTAAGTTGGATGCTGAc                             >  1:35517/1‑67 (MQ=255)
                                   ccGATAACCACGATGTTAAACAGGATAACCGCCAGAATGCACACCAGGCTGTAAGTTGGATGCTGAc                             >  1:489124/1‑67 (MQ=255)
                                           cACGATGTTAAACAGGATAACTGCCAGAATGCACAGCAGGCTGTAAGTTGGATGCt                                <  1:1220756/56‑1 (MQ=255)
                                                tGTTAAACAGGATAACCGCCAGAATGCACACCAGGCTGTAAGTTGGATGCTGACTGTAGCTGGCGACCt              <  1:187255/69‑1 (MQ=255)
                                                tGTTAAACAGGATAACCGCCAGAATGCACACCAGGCTGTAAGTTGGATGCTGACTGTAGCTGGCGACCt              <  1:183273/69‑1 (MQ=255)
                                                tGTTAAACAGGATAACCGCCAGAATGCACACCAGGCTGTAAGTTGGATGCTGACTGTAGCTGGCGACCt              <  1:400092/69‑1 (MQ=255)
                                                tGTTAAACAGGATAACCGCCAGAATGCACACCAGGCTGTAAGTTGGATGCTGACTGTAGCTGGCGACCt              <  1:467412/69‑1 (MQ=255)
                                                   tAAACAGGATAACCGCCAGAATGCACACCAGGCTGTAAGTTGGATGCTGACTGTAGCTGGCGACCt              >  1:841736/1‑66 (MQ=255)
                                                   tAAACAGGAAAACCGCCAGAATGCACACCAGGCTGTAAGTTGGATGCTGACTGTAGCTGGCGACCt              >  1:826654/1‑66 (MQ=255)
                                                          gATAACCGCCAGAATGCACACCAGGCTGTAAGTTGGATGCTGACTGTAGCTGGCGACCTGATAGAACaa    >  1:2152451/1‑69 (MQ=255)
                                                          gATAACCGCCAGAATGCACACCAGGCTGTAAGTTGGATGCTGACTGTAGCTGGCGACCTGATAGAACaa    >  1:1890672/1‑69 (MQ=255)
                                                          gATAACCGCCAGAATGCACACCAGGCTGTAAGTTGGATGCTGACTGTAGCTGGCGACCTGATAGAACaa    >  1:1478948/1‑69 (MQ=255)
                                                           aTAACCGCCAGAATGCACACCAGGCTGTAAGTTGGATGCTGACTGTAGCTGGc                   >  1:690486/1‑53 (MQ=255)
                                                           aTAACCGCCAGAATGCACACCAGGCTGTAAGTTGGATGCTGACTGTAGCTGGCGACCTGATAGAACAAt   >  1:2428344/1‑69 (MQ=255)
                                                           aTAACCGCCAGAATGCACACCAGGCTGTAAGTTGGATGCTGACTGTAGCTGGCGACCTGATAGAACAAt   >  1:808356/1‑69 (MQ=255)
                                                           aTAACCGCCAGAATGCACACCAGGCTGTAAGTTGGATGCTGACTGTAGCTGGCGACCTGATAGAACAAt   >  1:1145819/1‑69 (MQ=255)
                                                           aTAACCGCCAGAATGCACACCAGGCTGTAAGTTGGATGCTGACTGTAGCTGGCGACCTGATAGAACAAt   >  1:720247/1‑69 (MQ=255)
                                                           aTAACCGCCAGAATGCACACCAGGCTGTAAGTTGGATGCTGACTGTAGCTGGCGACCTGATAGAACAAt   >  1:135597/1‑69 (MQ=255)
                                                           aTAACCGCCAGAATGCACACCAGGCTGTAAGTTGGATGCTGACTGTAGCTGGCGACCTGATAGAACAAt   >  1:1439895/1‑69 (MQ=255)
                                                           aTAACCGCCAGAATGCACACCAGGCTGTAAGTTGGATGCTGACTGTAGCTGGCGACCTGATAGAACAAt   >  1:2901531/1‑69 (MQ=255)
                                                           aTAACCGCCAGAATGCACACCAGGCTGTAAGTTGGATGCTGACTGTAGCTGGCGACCTGATAGAACAAt   >  1:2787245/1‑69 (MQ=255)
                                                           aTAACCGCCAGAATGCACACCAGGCTGTAAGTTGGATGCTGACTGTAGCTGGCGACCTGATAGAACAAt   >  1:2627123/1‑69 (MQ=255)
                                                           aTAACCGCCAGAATGCACACCAGGCTGTAAGTTGGATGCTGACTGTAGCTGGCGACCTGATAGAACAAt   >  1:2602629/1‑69 (MQ=255)
                                                            tAACCGCCAGAATGCACACCAGGCTGTAAGTTGGATGCt                                <  1:3398017/39‑1 (MQ=255)
                                                            tAACCGCCAGAATGCACACCAGGCTGTAAGTTGGATGCt                                <  1:3394218/39‑1 (MQ=255)
                                                            tAACCGCCAGAATGCACACCAGGCTGTAAGTTGGATGCt                                <  1:2549553/39‑1 (MQ=255)
                                                            tAACCGCCAGAATGCACACCAGGCTGTAAGTTGGATGCt                                <  1:938964/39‑1 (MQ=255)
                                                             aaCCGCCAGAATGCACACCAGGCTGTAAGTTGGATGCt                                <  1:1945186/38‑1 (MQ=255)
                                                             aaCCGCCAGAATGCACACCAGGCTGTAAGTTGGATGCTGACTGTAGCTGGCGACCTGATAGAACAATg  >  1:2311263/1‑68 (MQ=255)
                                                                |                                                                
CGATATCCACCCGGCTACGCGCGCGGCGCAGCAGACCGATAACCACGATGTTAAACAGGATAACCGCCAGAATGCACACCAGGCTGTAAGTTGGATGCTGACTGTAGCTGGCGACCTGATAGAACAATG  >  minE/2678436‑2678564

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: