Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,732,555 A→G 48.9% T69T (ACA→ACG yfhL → predicted 4Fe‑4S cluster‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,732,5550AG48.9% ‑5.2 / 85.3 92T69T (ACA→ACGyfhLpredicted 4Fe‑4S cluster‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (35/12);  new base G (30/15);  total (65/27)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

AACCTGCCAGAAGGTGTGCCCGATCCCCAATACTATTGTGAAAGATCCGGCGCATGTCGAGACAGAAGAACAGTTGTGGGATAAATTTGTGCTGATGCACCACGCGGATAAAATTTAATATTCTA  >  minE/1732492‑1732616
                                                               |                                                             
aaCCTGCCAGAAGGTGTGCCCGATCCCCAATACTATTGTGAAAGATCCGGCGCATGTCGAGACagaag                                                           >  1:2100144/1‑68 (MQ=255)
aaCCTGCCAGAAGGTGTGCCCGATCCCCAATACTATTGTGAAAGATCCGGCGCATGTCGAGACagaag                                                           >  1:2096205/1‑68 (MQ=255)
 aCCTGCCAGAAGGTGTGCCCGATCCCCAATACTATTGTGAAAGATCCGGCGCATGTCGAGACAgaagaa                                                         <  1:3338989/69‑1 (MQ=255)
 aCCTGCCAGAAGGTGTGCCCGATCCCCAATACTATTGTGAAAGATCCGGCGCATGTCGAGACAgaagaa                                                         <  1:1701575/69‑1 (MQ=255)
        aGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAGATCCGCATCATGTCGAGACGGAAGAACAGtt                                                    <  1:3293723/67‑1 (MQ=255)
        aGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAGATCCGCAACATGTCGAGACGGAAGAACAGtt                                                    <  1:2629372/67‑1 (MQ=255)
        aGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAGATCCGCAACATGTCGAGACGGAAGAACAGtt                                                    <  1:2611525/67‑1 (MQ=255)
        aGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAGATCCGCAACATGTCGAGACGGAAGAACAGtt                                                    <  1:746035/67‑1 (MQ=255)
        aGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAGATCCGCAACATGTCGAGACGGAAGAACAGtt                                                    <  1:1828428/67‑1 (MQ=255)
        aGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAGATCCGCAACATGTCGAGACGGAAGAACAGtt                                                    <  1:2410450/67‑1 (MQ=255)
        aGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAGATCCGCAACATGTCGAGACGGAAGAACAGtt                                                    <  1:1919365/67‑1 (MQ=255)
        aGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAGATCCGCAACATGTCGAGACGGAAGAACAGtt                                                    <  1:1952982/67‑1 (MQ=255)
        aGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAGATCCGCAACATGTCGAGACGGAAGAACAGtt                                                    <  1:319714/67‑1 (MQ=255)
        aGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAGATCCGCAACATGTCGAGACGGAAGAACAGtt                                                    <  1:2393860/67‑1 (MQ=255)
        aGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAGATCCGCAACATGTCGAGACGGAAGAACAGtt                                                    <  1:3101962/67‑1 (MQ=255)
        aGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAGATCCGCAACATGTCGAGACGGAAGAACAGtt                                                    <  1:115549/67‑1 (MQ=255)
        aGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAGATCCGCAACATGTCGAGACGGAAGAACAGtt                                                    <  1:1014894/67‑1 (MQ=255)
              tgtgCCCGATCCCCAATACTATTGTGAAAGATCCGGCGCATGTCGAGACAGAAGAACAGTTGTGGGATa                                            <  1:2264571/69‑1 (MQ=255)
              tgtgCCCGATCCCCAATACTATTGTGAAAGATCCGGCGCATGTCGAGACAGAAGAACAGTTGTGGGATa                                            <  1:3154449/69‑1 (MQ=255)
              tgtgCCCGATCCCCAATACTATTGTGAAAGATCCGGCGCATGTCGAGACAGAAGAACAGTTGTGGGATa                                            <  1:122750/69‑1 (MQ=255)
               gtgCCCGATCCCCAATACTATTGTTAAAGATCCGGCGCATGTCGAGACAGAAGAACAGTTGTGGGATa                                            <  1:896779/68‑1 (MQ=255)
                  cccGATCCCCAATACTATTGTGAAAGATCCGGCGCATGTCGAGACAGAAGAACAGTTGTGGGATAAAtt                                        >  1:1296139/1‑69 (MQ=255)
                  cccGATCCCCAATACTATTGTGAAAGATCCGGCGCATGTCGAGACAGAAGAACAGTTGTGGGATAAAtt                                        >  1:502907/1‑69 (MQ=255)
                  cccGATCCCCAATACTATTGTGAAAGATCCGGCGCATGTCGAGACAGAAGAACAGTTGTGGGATAAAtt                                        >  1:3426671/1‑69 (MQ=255)
                  cccGATCCCCAATACTATTGTGAAAGATCCGGCGCATGTCGAGACAGAAGAACAGTTGTGGGATAAAtt                                        >  1:1190798/1‑69 (MQ=255)
                  cccGATCCCCAATACTATTGTGAAAGATACGGCGCATGTCGAGACAGAAGAACAGTTGTGGGATAAAtt                                        >  1:1733942/1‑69 (MQ=255)
                    cGATCCCCAATACTATTGTGAAAGATCCGGCGCATGTCGAGACAGAAGAACAGTTGTGGGATAAATTtg                                      <  1:1442094/69‑1 (MQ=255)
                    cGATCCCCAATACTATTGTGAAAGATCCGGCGCATGTCGAGACAGAAGAACAGTTGTGGGATAAATTtg                                      <  1:2076326/69‑1 (MQ=255)
                    cGATCCCCAATACTATTGTGAAAGATCCGGCGCATGTCGAGACAGAAGAACAGTTGTGGGATAAATTtg                                      <  1:908217/69‑1 (MQ=255)
                    cGATCCCCAATACTATTGTGAAAGATCCGGCGCATGTCGAGACAGAAGAACAGTTGTGGGATAAATTtg                                      <  1:2578891/69‑1 (MQ=255)
                    cGATCCCCAATACTATTGTGAAAGATCCGGCGCATGTCGAGACAGAAGAACAGTTGTGGGATAAATTtg                                      <  1:2908975/69‑1 (MQ=255)
                        ccccAATACTATTGTGAAAGATCCGGCGCATGTCGAGACAGAAGAACAGTTGTGGGATAAATTTGTGCt                                  >  1:2125032/1‑69 (MQ=255)
                        ccccAATACTATTGTGAAAGATCCGGCGCATGTCGAGACAGAAGAACAGTTGTGGGATAAATTTGTGCt                                  >  1:3234049/1‑69 (MQ=255)
                        ccccAATACTATTGTGAAAGATCCGCAACATGTCGAGACGGAAGAACAGtt                                                    <  1:2807411/51‑1 (MQ=255)
                        ccccAATACTATTGTGAAAGATCCGCAACATGTCGAGACGGAAGAACAGtt                                                    <  1:933873/51‑1 (MQ=255)
                                cTATTGTGAAAGATCCGGCGCATGTCGAGACAGAAGAACAGTTGTGGGATa                                            >  1:2680642/1‑51 (MQ=255)
                                cTATTGTGAAAGATCCGGCGCATGTCGAGACAGAAGAACAGTTGTGGGATAAATTTGTGCTGATGcacc                          >  1:403347/1‑69 (MQ=255)
                                cTATTGTGAAAGATCCGGCGCATGTCGAGACAGAAGAACAGTTGTGGGATAAATTTGTGCTGATGcacc                          >  1:2807220/1‑69 (MQ=255)
                                cTATTGTGAAAGATCCGGCGCATGTCGAGACAGAAGAACAGTTGTGGGATAAATTTGTGCTGATGcacc                          >  1:1032801/1‑69 (MQ=255)
                                cTATTGTGAAAGATCCGGCGCATGTCGAGACAGAAGAACAGTTGTGGGATAAATTTGTGCTGATGcacc                          >  1:1344067/1‑69 (MQ=255)
                                cTATTGTGAAAGATCCGGCGCATGTCGAGACAGAAGAACAGTTGTGGGATAAATTTGTGCTGATGcacc                          >  1:1446156/1‑69 (MQ=255)
                                cTATTGTGAAAGATCCGGCGCATGTCGAGACAGAAGAACAGTTGTGGGATAAATTTGTGCTGATGcacc                          >  1:1571012/1‑69 (MQ=255)
                                cTATTGTGAAAGATCCGGCGCATGTCGAGACAGAAGAACAGTTGTGGGATAAATTTGTGCTGATGcacc                          >  1:1192462/1‑69 (MQ=255)
                                cTATTGTGAAAGATCCGGCGCATGTCGAGACAGAAGAACAGTTGTGGGATAAATTTGTGCTGATGcacc                          >  1:1887400/1‑69 (MQ=255)
                                cTATTGTGAAAGATCCGGCGCATGTCGAGACAGAAGAACAGTTGTGGGATAAATTTGTGCTGATGcac                           >  1:2206174/1‑68 (MQ=255)
                                     gtgAAAGATCCGGCGCATGTCGAGACAGAAGAACAGTTGTGGGATAAATTTGTGCTGATGCACCAcgcg                     >  1:885371/1‑69 (MQ=255)
                                             tCCGGCGCATGTCGAGACAGAAGAACAGTTGTGGGATAAATTTGTGCTGATGCACCACGCGGATaaaa              >  1:644290/1‑68 (MQ=255)
                                              ccGGCGCATGTCGAGACAGAAGAACAGTTGTGGGATAAATTTGTGCTGATGCACCACGCGGATaaaa              >  1:2510821/1‑67 (MQ=255)
                                              ccGGCGCATGTCGAGACAGAAGAACAGTTGTGGGATAAATTTGTGCTGATGCACCACGCGGATaaaa              >  1:1276818/1‑67 (MQ=255)
                                              ccGGCGCATGTCGAGACAGAAGAACAGTTGTGGGATAAATTTGTGCTGATGCACCACGCGGATAAAAt             >  1:2526018/1‑68 (MQ=255)
                                              ccGGCGCATGTCGAGACAGAAGAACAGTTGTGGGATAAATTTGTGCTGATGCACCACGCGGATAAAAt             >  1:3159488/1‑68 (MQ=255)
                                              ccGGCGCATGTCGAGACAGAAGAACAGTTGTGGGATAAATTTGTGCTGATGCACCACGCGGATAAAAt             >  1:451710/1‑68 (MQ=255)
                                              ccGGCGCATGTCGAGACAGAAGAACAGTTGTGGGATAAATTTGTGCTGATGCACCACGCGGATAAAAt             >  1:576064/1‑68 (MQ=255)
                                              ccGGCGCATGTCGAGACAGAAGAACAGTTGTGGGATAAATTTGTGCTGATGCACCACGCGGATAAAAt             >  1:44814/1‑68 (MQ=255)
                                              ccGGCGCATGTCGAGACAGAAGAACAGTTGTGGGATAAATTTGTGCTGATGCACCACGCGGATAAAAt             >  1:1281035/1‑68 (MQ=255)
                                                   acATGTCGAGACGGGAGAACAGTTGTGGGATAAAtt                                        >  1:1164294/2‑36 (MQ=37)
                                                   acATGTCGAGACGGAAGAACAGTTGTGGGATAAAtt                                        >  1:3208065/2‑36 (MQ=38)
                                                   acATGTCGAGACGGAAGAACAGTTGTGGGATAAAtt                                        >  1:321777/2‑36 (MQ=38)
                                                   acATGTCGAGACGGAAGAACAGTTGTGGGATAAAtt                                        >  1:407826/2‑36 (MQ=38)
                                                   acATGTCGAGACGGAAGAACAGTTGTGGGATAAAtt                                        >  1:951033/2‑36 (MQ=38)
                                                   acATGTCGAGACGGAAGAACAGTTGTGGGATAAAtt                                        >  1:3248577/2‑36 (MQ=38)
                                                   acATGTCGAGACGGAAGAACAGTTGTGGGATAAAtt                                        >  1:643476/2‑36 (MQ=38)
                                                   acATGTCGAGACGGAAGAACAGTTGTGGGATAAAtt                                        >  1:3375558/2‑36 (MQ=38)
                                                   acATGTCGAGACGGAAGAACAGTTGTGGGATAAAtt                                        >  1:522922/2‑36 (MQ=38)
                                                   acATGTCGAGACGGAAGAACAGTTGTGGGATAAAtt                                        >  1:3421358/2‑36 (MQ=38)
                                                   acATGTCGAGACGGAAGAACAGTTGTGGGATAAAtt                                        >  1:451966/2‑36 (MQ=38)
                                                   acATGTCGAGACGGAAGAACAGTTGTGGGATAAAtt                                        >  1:432104/2‑36 (MQ=38)
                                                   acATGTCGAGACGGAAGAACAGTTGTGGGATAAAtt                                        >  1:3008593/2‑36 (MQ=38)
                                                   acATGTCGAGACGGAAGAACAGTTGTGGGATAAAtt                                        >  1:1087049/2‑36 (MQ=38)
                                                   acATGTCGAGACGGAAGAACAGTTGTGGGATAAAtt                                        >  1:1143772/2‑36 (MQ=38)
                                                   acATGTCGAGACGGAAGAACAGTTGTGGGATAAAtt                                        >  1:1158759/2‑36 (MQ=38)
                                                   acATGTCGAGACGGAAGAACAGTTGTGGGATAAAtt                                        >  1:1276077/2‑36 (MQ=38)
                                                   acATGTCGAGACGGAAGAACAGTTGTGGGATAAAtt                                        >  1:1291762/2‑36 (MQ=38)
                                                   acATGTCGAGACGGAAGAACAGTTGTGGGATAAAtt                                        >  1:1499888/2‑36 (MQ=38)
                                                   acATGTCGAGACGGAAGAACAGTTGTGGGATAAAtt                                        >  1:1523406/2‑36 (MQ=38)
                                                   acATGTCGAGACGGAAGAACAGTTGTGGGATAAAtt                                        >  1:1696093/2‑36 (MQ=38)
                                                   acATGTCGAGACGGAAGAACAGTTGTGGGATAAAtt                                        >  1:2100910/2‑36 (MQ=38)
                                                   acATGTCGAGACGGAAGAACAGTTGTGGGATAAAtt                                        >  1:2451537/2‑36 (MQ=38)
                                                   acATGTCGAGACGGAAGAACAGTTGTGGGATAAAtt                                        >  1:2515433/2‑36 (MQ=38)
                                                   acATGTCGAGACGGAAGAACAGTTGTGGGATAAAtt                                        >  1:3006539/2‑36 (MQ=38)
                                                   acATGTCGAGACGGAAGAACAGTTGTGGGATAAAtt                                        >  1:2918443/2‑36 (MQ=38)
                                                   acATGTCGAGACGGAAGAACAGTTGTGGGATAAAtt                                        >  1:2598753/2‑36 (MQ=38)
                                                   acATGTCGAGACGGAAGAACAGTTGTGGGATAAAtt                                        >  1:2590735/2‑36 (MQ=38)
                                                   acATGTCGAGACGGAAGAACAGTTGTGGGATAAAtt                                        >  1:2588659/2‑36 (MQ=38)
                                                   acATGTCGAGACGGAAGAACAGTTGTGGGATAAAtt                                        >  1:2551505/2‑36 (MQ=38)
                                                        tCGAGACAGAAGAACAGTTGTGGGATAAATTTGTGCTGATGCACCACGCGGATAAAATTTAATATTCTa  >  1:1760908/1‑69 (MQ=255)
                                                        tCGAGACAGAAGAACAGTTGTGGGATAAATTTGTGCTGATGCACCACGCGGATAAAATTTAATATTCTa  <  1:2018145/69‑1 (MQ=255)
                                                          gagaCAGAAGAACAGTTGTGGGATAAATTTGTGCTGATGCACCACGCGGATAAAATTTaa         >  1:1797187/1‑60 (MQ=255)
                                                          gagaCAGAAGAACAGTTGTGGGATAAATTTGTGCTGATGCACCACGCGGATAAAATTTaa         >  1:1738903/1‑60 (MQ=255)
                                                          gagaCAGAAGAACAGTTGTGGGATAAATTTGTGCTGATGCACCACGCGGATAAAATTTaa         >  1:19742/1‑60 (MQ=255)
                                                          gagaCAGAAGAACAGTTGTGGGATAAATTTGTGCTGATGCACCACGCGGATAAAATTTaa         >  1:2845749/1‑60 (MQ=255)
                                                          gagaCAGAAGAACAGTTGTGGGATAAATTTGTGCTGATGCACCACGCGGATAAAATTTaa         >  1:2973244/1‑60 (MQ=255)
                                                          gagaCAGAAGAACAGTTGTGGGATAAATTTGTGCTGATGCACCACGCGGATAAAATTTaa         >  1:355922/1‑60 (MQ=255)
                                                          gagaCAGAAGAACAGTTGTGGGATAAATTTCTGCTGAAGCACCACGCGGAt                  >  1:1914426/1‑51 (MQ=255)
                                                               agaagaACAGTTGTGGGATAAATTTGTGCTGATGCACCACGCGGATAAAAttt           >  1:1213495/1‑53 (MQ=255)
                                                               agaagaACAGTTGTGGGATAAATTTGTGCTGATGCACCACGCGGATAAAAttt           >  1:3045208/1‑53 (MQ=255)
                                                               |                                                             
AACCTGCCAGAAGGTGTGCCCGATCCCCAATACTATTGTGAAAGATCCGGCGCATGTCGAGACAGAAGAACAGTTGTGGGATAAATTTGTGCTGATGCACCACGCGGATAAAATTTAATATTCTA  >  minE/1732492‑1732616

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: