Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,680,432 G→T 26.2% A86A (GCC→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,680,4320GT26.2% 243.3 / 74.7 152A86A (GCC→GCAfeoBfused ferrous iron transporter, protein B
Reads supporting (aligned to +/- strand):  ref base G (45/67);  new base T (28/12);  total (73/79)
Fisher's exact test for biased strand distribution p-value = 1.60e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 153 total aligned reads displayed.

CAGTTGTAGCGTCAGGTACAGGTTACGCTCAAGGTTAGACGCATCCACCACGTTAATCAGCAGGTCGGCGTCGCCACTCAAAATGTAGTGACAGGCGATTTGCTCATCGAGCGAGGTCTGCGATGAGATGGTG  >  minE/2680366‑2680498
                                                                  |                                                                  
cagTTGTAGCGTCAGGTACAGGTTACGCTCAAGGTTAGACGCATCCACCACGTTAATCAGCAGGTCGGc                                                                  <  1:1798534/69‑1 (MQ=255)
cagTTGTAGCGTCAGGTACAGGTTACGCTCAAGGTTAGACGCATCCACCACGTTAATCAGCAGGTCGGc                                                                  <  1:610114/69‑1 (MQ=255)
   ttGTAGCGTCAGGTACAGGTTACGCTCAAGGTTAGACGCATCCACCACGTTAATCAGCAGGTCGGCGTc                                                               <  1:2788042/69‑1 (MQ=255)
     gTAGCGTCAGGTACAGGTTACGCTCAAGGTTAGACGCATCCACCACGTTAATCAGCAGGTCGGCGTCGc                                                             <  1:874354/69‑1 (MQ=255)
     gTAGCGTCAGGTACAGGTTACGCTCAAGGTTAGACGCATCCACCACGTTAATCAGCAGGTCGGCGTCGc                                                             <  1:1881564/69‑1 (MQ=255)
     gTAGCGTCAGGTACAGGTTACGCTCAAGGTTAAACGCATCCACCACGTTAATCAGCAGGTCGGCGTCGc                                                             <  1:2893190/69‑1 (MQ=255)
           tctggTACAGGTTACGCTCAAGGTTAGACGCATCCACCACGTTAATCAGCAGGTCTGCGTCGCCACTCa                                                       <  1:1422963/66‑1 (MQ=255)
             aGGTACAGGTTACGCTCAAGGTTAGACGCATCCACCACGTTAATCAGCAGGTCGGCGTCGCCACTCaaa                                                     <  1:1852199/69‑1 (MQ=255)
             aGGTACAGGTTACGCTCAAGGTTAGACGCATCCACCACGTTAATCAGCAGGTCGGCGTCGCCACTCaaa                                                     <  1:2656181/69‑1 (MQ=255)
                 aCATGTTACGCTCAAGGTTAGACGCATCCACCACGTTAATCAGCAGGTCTGCGTCGCCACTCa                                                       <  1:322838/63‑1 (MQ=255)
                 aCAGGTTACGCTCAAGGTTAGACGCATCCACCACGTTAATCAGCAGGTCTGCGTCGCCACTCa                                                       <  1:1759307/63‑1 (MQ=255)
                 aCAGGTTACGCTCAAGGTTAGACGCATCCACCACGTTAATCAGCAGGTCTGCGTCGCCACTCa                                                       <  1:1538168/63‑1 (MQ=255)
                 aCAGGTTACGCTCAAGGTTAGACGCATCCACCACGTTAATCAGCAGGTCTGCGTCGCCACTCa                                                       <  1:2146030/63‑1 (MQ=255)
                 aCAGGTTACGCTCAAGGTTAGACGCATCCACCACGTTAATCAGCAGGTCTGCGTCGCCACTCa                                                       <  1:2599265/63‑1 (MQ=255)
                 aCAGGTTACGCTCAAGGTTAGACGCATCCACCACGTTAATCAGCAGGTCTGCGTCGCCACTCa                                                       <  1:2740001/63‑1 (MQ=255)
                 aCAGGTTACGCTCAAGGTTAGACGCATCCACCACGTTAATCAGCAGGTCTGCGTCGCCACTCa                                                       <  1:679040/63‑1 (MQ=255)
                 aCAGGTTACGCTCAAGGTTAGACGCATCCACCACGTTAATCAGCAGGTCGGCGTCGCCACTCAAAATGt                                                 >  1:539619/1‑69 (MQ=255)
                 aCAGGTTACGCTCAAGGTTAGACGCATCCACCACGTTAATCAGCAGGTCGGCGTCGCCACTCAAAATGt                                                 >  1:1016565/1‑69 (MQ=255)
                  cAGGTTACGCTCAAGGTTAGACGCATCCACCACGTTAATCAGCAGGTCGGCGTCGCCACTCAAAATGTa                                                <  1:2977340/69‑1 (MQ=255)
                  cAGGTTACGCTCAAGGTTAGACGCATCCACCACGTTAATCAGCAGGTCGGCGTCGCCACTCAAAATGTa                                                >  1:1485110/1‑69 (MQ=255)
                  cAGGTTACGCTCAAGGTTAGACGCATCCACCACGTTAATCAGCAGGTCGGCGTCGCCACTCAAAATGTa                                                <  1:2739102/69‑1 (MQ=255)
                  cAGGTTACGCTCAAGGTTAGACGCATCCACCACGTTAATCAGCAGGTCGGCGTCGCCACTCAAAATGTa                                                <  1:145505/69‑1 (MQ=255)
                  cAGGTTACGCTCAAGGTTAGACGCATCCACCACGTTAATCAGCAGGTCGGCGTCGCCACTCAAAATGTa                                                <  1:2847287/69‑1 (MQ=255)
                  cAGGTTACGCTCAAGGTTAGACGCATCCACCACGTTAATCAGCAGGTCGGCGTCGCCACTCAAAATGTa                                                >  1:661836/1‑69 (MQ=255)
                  cAGGTTACGCTCAAGGTTAGACGCATCCACCACGTTAATCAGCAGGTCGGCGTCGCCACTCAAAATGTa                                                <  1:670096/69‑1 (MQ=255)
                  cAGGTTACGCTCAAGGTTAGACGCATCCACCACGTTAATCAGCAGGTCGGCGTCGCCACTCAAAATGTa                                                <  1:694681/69‑1 (MQ=255)
                  cAGGTTACGCTCAAGGTTAGACGCATCCACCACGTTAATCAGCAGGTCGGCGTCGCCACTCAAAATGTa                                                >  1:1510963/1‑69 (MQ=255)
                  cAGGTTACGCTCAAGGTTAGACGCATCCACCACGTTAATCAGCAGGTCGGCGTCGCCACTCAAAATGTa                                                <  1:2063021/69‑1 (MQ=255)
                  cAGGTTACGCTCAAGGTTAGACGCATCCACCACGTTAATCAGCAGGTCGGCGTCGCCACTCAAAATGTa                                                <  1:2868045/69‑1 (MQ=255)
                  cAGGTTACGCTCAAGGTTAGACGCATCCACCACGTTAATCAGCAGGTCGGCGTCGCCACTCAAAATGTa                                                <  1:2885226/69‑1 (MQ=255)
                  cAGGTTACGCTCAAGGTTAGACGCATCCACCACGTTAATCAGCAGGTCGGCGTCGCCACTCAAAATGTa                                                <  1:3068334/69‑1 (MQ=255)
                  cAGGTTACGCTCAAGGTTAGACGCATCCACCACGTTAATCAGCAGGTCGGCGTCGCCACTCAAAATGTa                                                <  1:819143/69‑1 (MQ=255)
                  cAGGTTACGCTCAAGGTTAGACGCATCCACCACGTTAATCAGCAGGTCGGCGTCGCCACTCAAAATGTa                                                <  1:2946813/69‑1 (MQ=255)
                       tACGCTCAAGGTTAGACGCATCCACCACGTTAATCAGCAGGTCGGCGTCGCCACTCAAAATGTAGTGa                                            >  1:2513439/1‑68 (MQ=255)
                       tACGCTCAAGGTTAGACGCATCCACCACGTTAATCAGCAGGTCGGCGTCGCCACTCAAAATGTAGTGAc                                           >  1:2557809/1‑69 (MQ=255)
                       tACGCTCAAGGTTAGACGCATCCACCACGTTAATCAGCAGGTCGGCGTCGCCACTCAAAATGTAGTGAc                                           >  1:2440639/1‑69 (MQ=255)
                       tACGCTCAAGGTTAGACGCATCCACCACGTTAATCAGCAGGTCGGCGTCGCCACTCAAAATGTAGTGAc                                           >  1:3314713/1‑69 (MQ=255)
                       tACGCTCAAGGTTAGACGCATCCACCACGTTAATCAGCAGGTCGGCGTCGCCACTCAAAATGTAGTGAc                                           >  1:2402441/1‑69 (MQ=255)
                       tACGCTCAAGGTTAGACGCATCCACCACGTTAATCAGCAGGTCGGCGTCGCCACTCAAAATGTAGTGAc                                           >  1:2401971/1‑69 (MQ=255)
                       tACGCTCAAGGTTAGACGCATCCACCACGTTAATCAGCAGGTCGGCGTCGCCACTCAAAATGTAGTGAc                                           >  1:1186558/1‑69 (MQ=255)
                       tACGCTCAAGGTTAGACGCATCCACCACGTTAATCAGCAGGTCGGCGTCGCCACTCAAAATGTAGTGAc                                           >  1:1482971/1‑69 (MQ=255)
                        aCGCTCATGGTTAGACGCATCCACCACGTTAATCAGCAGGTCGGCGTCGCCACTCAAAATGTAGTGACa                                          >  1:709007/1‑69 (MQ=255)
                        aCGCTCAAGGTTAGACGCATCCACCACGTTAATCAGCAGGTCGGCGTCGCCACTCATAATGTAGTGACa                                          >  1:1628605/1‑69 (MQ=255)
                        aCGCTCAAGGTTAGACGCATCCACCACGTTAATCAGCAGGTCGGCGTCGCCACTCAAAATGTAGTGAc                                           >  1:3398261/1‑68 (MQ=255)
                        aCGCTCAAGGTTAGACGCATCCACCACGTTAATCAGCAGGTCGGCGTCGCCACTCAAAATGTAGTGACa                                          >  1:3285110/1‑69 (MQ=255)
                        aCGCTCAAGGTTAGACGCATCCACCACGTTAATCAGCAGGTCGGCGTCGCCACTCAAAATGTAGTGACa                                          >  1:776771/1‑69 (MQ=255)
                        aCGCTCAAGGTTAGACGCATCCACCACGTTAATCAGCAGGTCGGCGTCGCCACTCAAAATGTAGTGACa                                          >  1:2672308/1‑69 (MQ=255)
                        aCGCTCAAGGTTAGACGCATCCACCACGTTAATCAGCAGGTCGGCGTCGCCACTCAAAATGTAGTGACa                                          >  1:2834179/1‑69 (MQ=255)
                        aCGCTCAAGGTTAGACGCATCCACCACGTTAATCAGCAGGTCGGCGTCGCCACTCAAAATGTAGTGACa                                          >  1:1570034/1‑69 (MQ=255)
                        aCGCTCAAGGTTAGACGCATCCACCACGTTAATCAGCAGGTCGGCGTCGCCACTCAAAATATAGTGAc                                           >  1:2197024/1‑68 (MQ=255)
                                      aCGCATCCACCACGTTAATCAGCAGGTCGGCGTCGCCACTCAAAATGTAGTGACAGGCGATTTGCTCa                             <  1:1359473/68‑1 (MQ=255)
                                      aCGCATCCACCACGTTAATCAGCAGGTCGGCGTCGCCACTCAAAATGTAGTGACAGGCGATTTGCTCa                             <  1:2476547/68‑1 (MQ=255)
                                      aCGCATCCACCACGTTAATCAGCAGGTCGGCGTCGCCACTCAAAATGTAGTGACAGGCGATTTGCTCa                             <  1:1335466/68‑1 (MQ=255)
                                      aCGCATCCACCACGTTAATCAGCAGGTCGGCGTCGCCACTCAAAATGTAGTGACAGGCGATTTGCTCa                             <  1:1971995/68‑1 (MQ=255)
                                      aCGCATCCACCACGTTAATCAGCAGGTCGGCGTCGCCACTCAAAATGTAGTGACAGGCGATTTGCTCa                             <  1:1558305/68‑1 (MQ=255)
                                       cGCATCCACCACGTTAATCAGCAGGTCTGCGTCGCCACTCa                                                       <  1:642775/41‑1 (MQ=255)
                                        gCATCCACCACGTTAATCAGCAGGTCTGCGTCGCCACTCAAAATGTAGTGACAGGCGATTTGCTc                              >  1:1470052/1‑65 (MQ=255)
                                        gCATCCACCACGTTAATCAGCAGGTCTGCGTCGCCACTCAAAATGTAGTGACAGGCGATTTGCTc                              >  1:2957028/1‑65 (MQ=255)
                                        gCATCCACCACGTTAATCAGCAGGTCTGCGTCGCCACTCAAAATGTAGTGACAGGCGATTTGCTCa                             >  1:3307232/1‑66 (MQ=255)
                                        gCATCCACCACGTTAATCAGCAGGTCTGCGTCGCCACTCAAAATGTAGTGACAGGCGATTTGCTCa                             >  1:3341619/1‑66 (MQ=255)
                                        gCATCCACCACGTTAATCAGCAGGTCTGCGTCGCCACTCAAAATGTAGTGACAGGCGATTTGCTCa                             >  1:986631/1‑66 (MQ=255)
                                        gCATCCACCACGTTAATCAGCAGGTCTGCGTCGCCACTCAAAATGTAGTGACAGGCGATTTGCTCa                             >  1:412124/1‑66 (MQ=255)
                                        gCATCCACCACGTTAATCAGCAGGTCTGCGTCGCCACTCAAAATGTAGTGACAGGCGATTTGCTCa                             >  1:782414/1‑66 (MQ=255)
                                        gCATCCACCACGTTAATCAGCAGGTCTGCGTCGCCACTCAAAATGTAGTGACAGGCGATTTGCTCa                             >  1:3127321/1‑66 (MQ=255)
                                        gCATCCACCACGTTAATCAGCAGGTCTGCGTCGCCACTCAAAATGTAGTGACAGGCGATTTGCTCa                             >  1:803734/1‑66 (MQ=255)
                                        gCATCCACCACGTTAATCAGCAGGTCTGCGTCGCCACTCAAAATGTAGTGACAGGCGATTTGCTCa                             >  1:92590/1‑66 (MQ=255)
                                        gCATCCACCACGTTAATCAGCAGGTCTGCGTCGCCACTCAAAATGTAGTGACAGGCGATTTGCTCa                             >  1:1415248/1‑66 (MQ=255)
                                        gCATCCACCACGTTAATCAGCAGGTCTGCGTCGCCACTCAAAATGTAGTGACAGGCGATTTGCTCa                             >  1:2375666/1‑66 (MQ=255)
                                        gCATCCACCACGTTAATCAGCAGGTCTGCGTCGCCACTCAAAATGTAGTGACAGGCGATTTGCTCa                             >  1:2372437/1‑66 (MQ=255)
                                        gCATCCACCACGTTAATCAGCAGGTCTGCGTCGCCACTCAAAATGTAGTGACAGGCGATTTGCTCa                             >  1:2293828/1‑66 (MQ=255)
                                        gCATCCACCACGTTAATCAGCAGGTCTGCGTCGCCACTCAAAATGTAGTGACAGGCGATTTGCTCa                             >  1:2179657/1‑66 (MQ=255)
                                        gCATCCACCACGTTAATCAGCAGGTCTGCGTCGCCACTCAAAATGTAGTGACAGGCGATTTGCTCa                             >  1:2613646/1‑66 (MQ=255)
                                        gCATCAACCACGTTAATCAGCAGGTCTGCGTCGCCACTCAAAATGTAGTGACAGGCGATTTGCTCa                             >  1:2005592/1‑66 (MQ=255)
                                             caccacGTTAATCAGCAGGTCGGCGTCGCCACTCAAAATGTAGTGACAGGCGATTTGCTCATcgagcga                     >  1:2172813/1‑69 (MQ=255)
                                             caccacGTTAATCAGCAGGTCGGCGTCGCCACTCAAAATGTAGTGACAGGCGATTTGCTCATcgagcga                     >  1:1145457/1‑69 (MQ=255)
                                             caccacGTTAATCAGCAGGTCGGCGTCGCCACTCAAAATGTAGTGACAGGCGATTTGCTCATcgagcga                     >  1:3031253/1‑69 (MQ=255)
                                             caccacGTTAATCAGCAGGTCGGCGTCGCCACTCAAAATGTAGTGACAGGCGATTTGCTCATcgagcga                     >  1:3089616/1‑69 (MQ=255)
                                                  cGTTAATCAGCAGGTCGGCGTCGCCACTCAAAATGTAGTGACAGGCg                                      <  1:1911613/47‑1 (MQ=255)
                                                  cGTTAATCAGCAGGTCGGCGTCGCCACTCAAAATGTAGTGACAGGCg                                      <  1:2165460/47‑1 (MQ=255)
                                                  cGTTAATCAGCAGGTCGGCGTCGCCACTCAAAATGTAGTGACAGGCg                                      <  1:1838435/47‑1 (MQ=255)
                                                  cGTTAATCAGCAGGTCGGCGTCGCCACTCAAAATGTAGTGACAGGCg                                      <  1:2341791/47‑1 (MQ=255)
                                                  cGTTAATCAGCAGGTCGGCGTCGCCACTCAAAATGTAGTGACAGGCg                                      <  1:406873/47‑1 (MQ=255)
                                                    ttAATCAGCAGGTCGGCGTCGCCACTCAAAATGTAGTGACAGGCGATTTGCTCATCGAGCGAGGTCTGc              <  1:1756140/69‑1 (MQ=255)
                                                    ttAATCAGCAGGTCGGCGTCGCCACTCAAAATGTAGTGACAGGCGATTTGCTCATCGAGCGAGGTCTGc              <  1:3289193/69‑1 (MQ=255)
                                                    ttAATCAGCAGGTCGGCGTCGCCACTCAAAATGTAGTGACAGGCGATTTGCTCATCGAGCGAGGTCTGc              >  1:2889340/1‑69 (MQ=255)
                                                    ttAATCAGCAGGTCGGCGTCGCCACTCAAAATGTAGTGACAGGCGATTTGCTCATCGAGCGAGGTCTGc              <  1:2675887/69‑1 (MQ=255)
                                                    ttAATCAGCAGGTCGGCGTCGCCACTCAAAATGTAGTGACAGGCGATTTGCTCATCGAGCGAGGTCTGc              <  1:2536691/69‑1 (MQ=255)
                                                    ttAATCAGCAGGTCGGCGTCGCCACTCAAAATGTAGTGACAGGCGATTTGCTCATCGAGCGAGGTCTGc              <  1:348494/69‑1 (MQ=255)
                                                      aaTCAGCAGGTCGGCGTCGCCACTCAAAATGTAGTGACAGGCGATTTGCTCATCGAGCGAGGTCTGCg             >  1:1790069/1‑68 (MQ=255)
                                                                tCGGCGTCGCCACTCAAAATGTAGTGACAGGCGATTTGCTCAt                            >  1:948862/1‑43 (MQ=255)
                                                                tCGGCGTCGCCACTCAAAATGTAGTGACAGGCGATTTGCTCATCGAGCGAGGTCTGCGATGAGAtggtg  >  1:352632/1‑69 (MQ=255)
                                                                tCGGCGTCGCCACTCAAAATGTAGTGACAGGCGATTTGCTCATCGAGCGAGGTCTGCGATGAGAtggtg  >  1:2542143/1‑69 (MQ=255)
                                                                tCGGCGTCGCCACTCAAAATGTAGTGACAGGCGATTTGCTCATCGAGCGAGGTCTGCGATGAGAtggtg  >  1:3574/1‑69 (MQ=255)
                                                                tCGGCGTCGCCACTCAAAATGTAGTGACAGGCGATTTGCTCATCGAGCGAGGTCTGCGATGAGAtggtg  >  1:2820284/1‑69 (MQ=255)
                                                                tCGGCGTCGCCACTCAAAATGTAGTGACAGGCGATTTGCTCATCGAGCGAGGTCTGCGATGAGAtggtg  >  1:84984/1‑69 (MQ=255)
                                                                tCGGCGTCGCCACTCAAAATGTAGTGACAGGCGATTTGCTCATCGAGCGAGGTCTGCGATGAGAtggtg  >  1:1015633/1‑69 (MQ=255)
                                                                tCGGCGTCGCCACTCAAAATGTAGTGACAGGCGATTTGCTCATCGAGCGAGGTCTGCGATGAGAtggtg  >  1:2439137/1‑69 (MQ=255)
                                                                 cGGCGTCGCCACTCAAAATGTAGTGACAGGCGATTTGCTCATcgagcg                      <  1:256631/48‑1 (MQ=255)
                                                                 cGGCGTCGCCACTCAAAATGTAGTGACAGGCGATTTGCTCATcgagcg                      <  1:541536/48‑1 (MQ=255)
                                                                 cGGCGTCGCCACTCAAAATGTAGTGACAGGCGATTTGCTCATcgagcg                      <  1:1880598/48‑1 (MQ=255)
                                                                  |                                                                  
CAGTTGTAGCGTCAGGTACAGGTTACGCTCAAGGTTAGACGCATCCACCACGTTAATCAGCAGGTCGGCGTCGCCACTCAAAATGTAGTGACAGGCGATTTGCTCATCGAGCGAGGTCTGCGATGAGATGGTG  >  minE/2680366‑2680498

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: