Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,679,577 C→G 29.5% G371G (GGG→GGC feoB ← fused ferrous iron transporter, protein B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,679,5770CG29.5% 89.3 / 43.8 71G371G (GGG→GGCfeoBfused ferrous iron transporter, protein B
Reads supporting (aligned to +/- strand):  ref base C (22/28);  new base G (13/8);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 2.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

TTCCCCGGCAAGCCCAGCGCCTGCATCAGACGGTCCATCACAAACGCCGCACGCGCCATATACCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCAATCTGCGGCACCAGTGGCAGCACGGT  >  minE/2679515‑2679645
                                                              |                                                                    
ttCCCCGGCAAGCCCAGCTCCTGCATCAGACGGTCCATCACAAACGCCGCACGCGCCATATACCCGGAg                                                                <  1:3336996/69‑1 (MQ=255)
ttCCCCGGCAAGCCCAGCGCCTGCATCAGACGGTCCATCACAAACGCCGCACGCGCCATATACCCGGAg                                                                <  1:2406075/69‑1 (MQ=255)
ttCCCCGGCAAGCCCAGCGCCTGCATCAGACGGTCCATCACAAACGCCGCACGCGCCATATACCCGGAg                                                                <  1:2740413/69‑1 (MQ=255)
          aGCCCAGCGCCTGCATCAGACGGTCCATCACAAACGCCGCACGCGCCATATACCCGGAGTCCTCaagga                                                      <  1:1027282/69‑1 (MQ=255)
          aGCCCAGCGCCTGCATCAGACGGTCCATCACAAACGCCGCACGCGCCATATACCCGGAGTCCTCaagga                                                      <  1:2931036/69‑1 (MQ=255)
          aGCCCAGCGCCTGCATCAGACGGTCCATCACAAACGCCGCACGCGCCATATACCCGGAGTCCTCaagga                                                      <  1:1486362/69‑1 (MQ=255)
          aGCCCAGCGCCTGCATCAGACGGTCCATCACAAACGCCGCACGCGCCATATACCCGGAGTCCTCaagga                                                      <  1:1614472/69‑1 (MQ=255)
            cccAGCGCCTGCATCAGACGGTCCATCACAAACGCCGCACGCGCCATATACCCGGAGTCCTCaagg                                                       <  1:463966/66‑1 (MQ=255)
             ccAGCGCCTGCATCAGACGGTCCATCACAAACGCCGCACGCGCCATATACCCGGAGTCCTCaaggaagg                                                   <  1:2670718/69‑1 (MQ=255)
             ccAGCGCCTGCATCAGACGGTCCATCACAAACGCCGCACGCGCCATATACCCGGAGTCCTCaaggaagg                                                   <  1:2721548/69‑1 (MQ=255)
             ccAGCGCCTGCATCAGACGGTCCATCACAAACGCCGCACGCGCCATATACCCGGAGTCCTCaaggaagg                                                   <  1:673572/69‑1 (MQ=255)
             ccAGCGCCTGCATCAGACGGTCCATCACAAACGCCGCACGCGCCATATACCCGGAGTCCTCaaggaagg                                                   <  1:677182/69‑1 (MQ=255)
             ccAGCGCCTGCATCAGACGGTCCATCACAAACGCCGCACGCGCCATATACCCGGAGTCCTCaaggaagg                                                   <  1:502410/69‑1 (MQ=255)
             ccAGCGCCTGCATCAGACGGTCCATCACAAACGCCGCACGCGCCATATACCCGGAGTCCTCaaggaagg                                                   <  1:853374/69‑1 (MQ=255)
                           aGACGGTCCATCACAAACGCCGCACGCGCCATATACCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTAc                                     >  1:2484145/1‑69 (MQ=255)
                            gACGGTCCATCACAAACGCCGCACGCGCCATATACCCGGAGTCCTCaaggaagg                                                   <  1:1152932/54‑1 (MQ=255)
                              cGGTCCATCACAAACGCCGCACGCGCCATATACCCGGAGTCCTCAAGGAAGGATAGGAACAGGTAcatc                                  <  1:3387697/69‑1 (MQ=255)
                              cGGTCCATCACAAACGCCGCACGCGCCATATACCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTAcatc                                  <  1:3298094/69‑1 (MQ=255)
                              cGGTCCATCACAAACGCCGCACGCGCCATATACCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTAcatc                                  <  1:1060789/69‑1 (MQ=255)
                              cGGTCCATCACAAACGCCGCACGCGCCATATACCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTAcatc                                  <  1:1130885/69‑1 (MQ=255)
                              cGGTCCATCACAAACGCCGCACGCGCCATATACCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTAcatc                                  <  1:121741/69‑1 (MQ=255)
                              cGGTCCATCACAAACGCCGCACGCGCCATATACCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTAcatc                                  <  1:2908836/69‑1 (MQ=255)
                              cGGTCCATCACAAACGCCGCACGCGCCATATACCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTAcatc                                  <  1:1408680/69‑1 (MQ=255)
                              cGGTCCATCACAAACGCCGCACGCGCCATATACCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTAcatc                                  <  1:1967375/69‑1 (MQ=255)
                              cGGTCCATCACAAACGCCGCACGCGCCATATACCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTAcatc                                  <  1:1510210/69‑1 (MQ=255)
                              cGGTCCATCACAAACGCCGCACGCGCCATATACCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTAcatc                                  <  1:2653602/69‑1 (MQ=255)
                              cGGTCCATCACAAACGCCGCACGCGCCATATACCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTAcatc                                  <  1:302827/69‑1 (MQ=255)
                              cGGTCCATCACAAACGCCGCACGCGCCATATACCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTAcatc                                  <  1:1660907/69‑1 (MQ=255)
                                       acaAACGCCGCACGCGCCATATACCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCAATc                         >  1:711699/1‑69 (MQ=255)
                                       acaAACGCCGCACGCGCCATATACCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCAATc                         >  1:2937850/1‑69 (MQ=255)
                                       acaAACGCCGCACGCGCCATATACCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCAATc                         >  1:261505/1‑69 (MQ=255)
                                       acaAACGCCGCACGCGCCATATACCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCAATc                         >  1:838782/1‑69 (MQ=255)
                                            cgccgcACGCGCCATATACCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCAATCTGCgg                    <  1:42341/69‑1 (MQ=255)
                                                 cACGCGCCATATAGCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCGATCTGCGGCACCa               <  1:65359/69‑1 (MQ=255)
                                                 cACGCGCCATATAGCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCGATCTGCGGCACCa               <  1:551639/69‑1 (MQ=255)
                                                 cACGCGCCATATAGCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCGATCTGCGGCACCa               <  1:295408/69‑1 (MQ=255)
                                                 cACGCGCCATATAGCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCGATCTGCGGCACCa               <  1:415054/69‑1 (MQ=255)
                                                 cACGCGCCATATAGCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCGATCTGCGGCACCa               <  1:2611129/69‑1 (MQ=255)
                                                 cACGCGCCATATAGCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCGATCTGCGGCACCa               <  1:1702344/69‑1 (MQ=255)
                                                 cACGCGCCATATAGCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCGATCTGCGGCACCa               <  1:1867300/69‑1 (MQ=255)
                                                 cACGCGCCATATAGCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCGATCTGCGGCACCa               <  1:2186889/69‑1 (MQ=255)
                                                        cATATAGCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCGATCTg                       >  1:3274847/1‑54 (MQ=255)
                                                        cATATAGCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCGATCTGc                      >  1:709508/1‑55 (MQ=255)
                                                        cATATAGCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCGATCTGc                      >  1:1023979/1‑55 (MQ=255)
                                                        cATATAGCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCGATCTGc                      >  1:2109449/1‑55 (MQ=255)
                                                        cATATAGCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCGATCTGc                      >  1:1963994/1‑55 (MQ=255)
                                                        cATATAGCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCGATCTGc                      >  1:2389807/1‑55 (MQ=255)
                                                        cATATAGCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCGATCTGc                      >  1:3072381/1‑55 (MQ=255)
                                                        cATATAGCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCGATCTGc                      >  1:3076382/1‑55 (MQ=255)
                                                        cATATAGCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCGATCTGc                      >  1:1656895/1‑55 (MQ=255)
                                                        cATATAGCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCGATCTGc                      >  1:3127087/1‑55 (MQ=255)
                                                        cATATAGCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCGATCTGc                      >  1:1588826/1‑55 (MQ=255)
                                                        cATATAGCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCGATCTGc                      >  1:1404514/1‑55 (MQ=255)
                                                        cATATAGCCGGAGTCCTCAAGGAAGGAGAGCAACAGGTACATCATGCCGATCTGc                      >  1:3361197/1‑55 (MQ=255)
                                                            taCCCGGAGTCCTCAAGGAAGGAGAGGAGCAGGTACATCATGCCAATCTGCGGCACCAGTGGCAGCACg    >  1:3134076/1‑69 (MQ=255)
                                                            taCCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATg                               >  1:1953869/1‑42 (MQ=255)
                                                            taCCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCa                            >  1:672854/1‑45 (MQ=255)
                                                            taCCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCa                            >  1:2274723/1‑45 (MQ=255)
                                                            taCCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCAATCTGCGGCAcc                >  1:991609/1‑57 (MQ=255)
                                                            taCCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCAATCTGCGGCAcc                >  1:2167868/1‑57 (MQ=255)
                                                            taCCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCAATCTGCGGCACCAGTGGCAGCACg    >  1:2836649/1‑69 (MQ=255)
                                                            taCCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCAATCTGCGGCACCAGTGGCAGCACg    >  1:271626/1‑69 (MQ=255)
                                                            taCCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCAATCTGCGGCACCAGTGGCAGCACg    >  1:2452904/1‑69 (MQ=255)
                                                            taCCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCAATCTGCGGCACCAGTGGCAGCACg    >  1:2265745/1‑69 (MQ=255)
                                                            taCCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCAATCTGCGGCACCAGTGGCAGCACg    >  1:2260473/1‑69 (MQ=255)
                                                            taCCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCAATCTGCGGCACCAGTGGCAGCACg    >  1:1935720/1‑69 (MQ=255)
                                                            taCCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCAATCTGCGGCACCAGTGGCAGCACg    >  1:1100477/1‑69 (MQ=255)
                                                            taCCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCAATCTGCGGCACCAGTGGCAGCACg    >  1:883486/1‑69 (MQ=255)
                                                            taCCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCAATCTGCGGCACCAGTGGCAGCACg    >  1:922928/1‑69 (MQ=255)
                                                            taCCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCAATCTGCGGCACCAGTGGCAGCACg    >  1:1039580/1‑69 (MQ=255)
                                                            taCCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCAATCTGCGGCACCA‑‑GGCAGCACGgt  >  1:2880284/1‑69 (MQ=255)
                                                              cccGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCAATCTGCGGCACCAGTGGCAGCACGgt  <  1:31126/69‑1 (MQ=255)
                                                              cccGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCAATCTGCGGCACCAGTGGCAGCACGgt  <  1:2933759/69‑1 (MQ=255)
                                                              cccGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCAATCTGCGGCACCAGTGGCAGCACGgt  <  1:2294592/69‑1 (MQ=255)
                                                              cccGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCAATCTGCGGCACCAGTGGCAGCACGgt  <  1:862521/69‑1 (MQ=255)
                                                              |                                                                    
TTCCCCGGCAAGCCCAGCGCCTGCATCAGACGGTCCATCACAAACGCCGCACGCGCCATATACCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCAATCTGCGGCACCAGTGGCAGCACGGT  >  minE/2679515‑2679645

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: