Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,609,702 C→T 23.8% E151E (GAG→GAA fre ← flavin reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,609,7020CT23.8% 80.6 / 22.6 42E151E (GAG→GAAfreflavin reductase
Reads supporting (aligned to +/- strand):  ref base C (17/15);  new base T (6/4);  total (23/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01

AGCCCGCTTCCGGTTGTTCAACCACCGGCACCACTTGCAGACCAGGATGCTTCAACGAAAGCGCCTCAAGCTCGCAGAGATCATACAGATGCTGCTCTTCACGCCCGCCCCAGTAAATGGTGATATCACGGTTTGGGT  >  W3110S.gb/3609632‑3609769
                                                                      |                                                                   
aGCCCGCTTCCGGTTGTTCAACCACCGGCACCACTTGCAGACCAGGATGCTTCAACGAAAGCGCCTCAAGc                                                                     <  1:4280318/71‑1 (MQ=255)
   ccGCTTCCGGTTGTTCAACCACCGGCACCACTTGCAGACCAGGATGCTTCAACGAAAGCGCCTCAAGCtag                                                                  >  1:4248671/1‑69 (MQ=255)
   ccGCTTCCGGTTGTTCAACCACCGGCACCACTTGCAGACCAGGATGCTTCAACGAAAGCGCCTCAAGCtag                                                                  >  1:4243039/1‑69 (MQ=255)
         ccGGTTGTTCAACCACCGGCACCACTTGCAGACCAGGATGCTTCAACGAAAGCGCCTCAAGCTCGCaga                                                              <  1:1781378/69‑1 (MQ=255)
         ccGGTTGTTCAACCACCGGCACCACTTGCAGACCAGGATGCTTCAACGAAAGCGCCTCAAGCTCGCa                                                                >  1:444065/1‑67 (MQ=255)
               gttCAACCACCGGCACCACTTGCAGACCAGGATGCTTCAACGAAAGCGCCTCAAGCTCGCAGAGATCATa                                                       <  1:4815300/70‑1 (MQ=255)
               gttCAACCACCGGCACCACTTGCAGACCAGGATGCTTCAACGAAAGCGACTCAAGCTCGCAGAGATCATa                                                       >  1:256542/1‑70 (MQ=255)
                   aaCCACCGGCACCACTTGCAGACCAGGATGCTTCAACGAAAGCGCCTCAAGCTCGCAGAGATCATACAGa                                                   <  1:2085256/70‑1 (MQ=255)
                   aaCCACCGGCACCACTTGCAGACCAGGATGCTTCAACGAAAGCGCCTCAAGCTCGCAGAGATCATACAGAt                                                  <  1:2434960/71‑1 (MQ=255)
                   aaCCACCGGCACCACTTGCAGACCAGGATGCTTCAACGAAAGCGACTCAAGCTCGCAGAGATCATACAg                                                    <  1:4098493/69‑1 (MQ=255)
                     ccaccGGCACCACTTGCAGACCAGGATGCTTCAACGAAAGCGCCTCAAGTTCGCAGAGATCATACAGAtgc                                                <  1:1420348/71‑1 (MQ=255)
                       accGGCACCACTTGCAGACCAGGATGCTTCAACGAAAGCGCCTCAAGCTCGCAGAGATCATACAGAtgctg                                              >  1:3597781/1‑71 (MQ=255)
                            caccacTTGCAGACCAGGATGCTTCAACGAAAGCGCCTCAAGCTCGCAGAGATCATACAGATGCTGCTCtt                                         <  1:2094371/71‑1 (MQ=255)
                              ccacTTGCAGACCAGGATGCTTCAACGAAAGCGCCTCAAGCTCGCa                                                                <  1:2653734/46‑1 (MQ=255)
                                    gCAGACCAGGATGCTTCAACGAAAGCGCCTCAAGCTCGCAGAGATCATACAGATGCTGCTCTTCAcgcccg                                 >  1:3424127/1‑71 (MQ=255)
                                          cAGGATGCTTCAACGAAAGCGCCTCAAGCTCGCAGAGATCATACAGATGCTGCTCTTCACGCCCGCCCCAg                           >  1:3812327/1‑71 (MQ=255)
                                            ggATGCTTCAACGAAAGCGCCTCAAGTTCGCAGAGATCATACAGATGCTGCTCTTCACGCCCGCCCCAGTa                         <  1:2358862/71‑1 (MQ=255)
                                            ggATGCTTCAACGAAAGCGCCTCAAGCTCGCAGAGATCATa                                                       >  1:276200/1‑41 (MQ=255)
                                            ggATGCTTCAACGAAAGCGCCTCAAGCTCGCAGAGATCATa                                                       >  1:1252822/1‑41 (MQ=255)
                                            ggATGCTTCAACGAAAGCGCCTCAAGCTCGCAGAGATCATACAGATGCTGCTCTTc                                        >  1:2566968/1‑56 (MQ=255)
                                            ggATGCTTCAACGAAAGCGCCTCAAGCTCGCAGAGATCATACAGATGCTGCTCTTCa                                       >  1:1715897/1‑57 (MQ=255)
                                              aTGCTTCAACGAAAGCGCCTCAAGCTCGCAGAGATCATACAGATGCTGctct                                          >  1:611940/1‑52 (MQ=255)
                                              aTGCTTCAACGAAAGCGCCTCAAGCTCGCAGAGATCATACAGATGCTGCTCTTCACGCCCGCCCCAGTaaa                       <  1:2581742/71‑1 (MQ=255)
                                               tGCTTCAACGAAAGCGCCTCAAGCTCGCAGAGATCATACAGATGCTGCTCTTCACGCCCGCCCCAGTAAAt                      <  1:2279388/71‑1 (MQ=255)
                                                 cTTCAACGAAAGCGCCTCAAGTTCGCAGAGATCATACAGATGCTGCTCTTCACGCCCGCCCCAGTAAATgg                    >  1:3348283/1‑71 (MQ=255)
                                                 cTTCAACGAAAGCGCCTCAAGCTCGCAGAGATCATACAGATGCTGCTCTTCAcgcccgc                                >  1:2179151/1‑59 (MQ=255)
                                                    cAACGAAAGCGCCTCAAGTTCGCAGAGATCATACAGATGCTGCTCTTCACGCCCGCCCCAGTAAATGGTGa                 >  1:267433/1‑71 (MQ=255)
                                                    cAACGAAAGCGCCTCAAGTTCGCAGAGATCATACAGATGCTGCTCTTCACGCCCGCCCCAGTAAATGGTGa                 >  1:1000485/1‑71 (MQ=255)
                                                    cAACGAAAGCGCCTCAAGCTCGCAGAGATCATACAGATGCTGCTCTTCACGCCCGCCCCAGTAAATGGt                   <  1:4498261/69‑1 (MQ=255)
                                                     aaCGAAAGCGCCTCAAGCTCGCAGAGATCATACAGATGCTGCTCTTCACGCCCGCCCCAGTAAATGGTGat                <  1:173465/71‑1 (MQ=255)
                                                      aCGAAAGCGCCTCAAGCTCGCAGAGATCATACAGATGCTGCTCTTCACGCCCGCCCCa                            >  1:1440772/1‑58 (MQ=255)
                                                      aCGAAAGCGCCTCAAGCTCGCAGAGATCATACAGATGATGCTCTTc                                        >  1:3164891/1‑46 (MQ=255)
                                                        gAAAGCGCCTCAAGCTCGCAGAGATCATACAGATGCTGCTCTTCACGCCCGCCCCAGTAAATGGTGatat              >  1:1517527/1‑70 (MQ=255)
                                                            gcgcCTCAAGCTCGCAGAGATCATACAGATGCTGCTCTTCACGCCCGCCCCAGTAAATGGt                   <  1:3914594/61‑1 (MQ=255)
                                                              gcCTCAAGTTCGCAGAGATCATACAGATGCTGCTCTTCACGCCCGCCCCAGTAAAt                      >  1:4571728/1‑56 (MQ=255)
                                                              gcCTCAAGTTCGCAGAGATCATACAGATGCTGCTCTTCACGCCCGCCCCAGTAAAt                      >  1:478049/1‑56 (MQ=255)
                                                              gcCTCAAGTTCGCAGAGATCATACAGATGCTGCTCTTCACGCCCGCCCCAGTAAAt                      >  1:790163/1‑56 (MQ=255)
                                                                 tCAAGTTCGCAGAGATCATACAGATGCTGCTCTTCACGCCCGCCCCAGTaa                        <  1:3610260/51‑1 (MQ=255)
                                                                 tCAAGTTCGCAGAGATCATACAGATGCTGCTCTTCACGCCCGCCCCAGTaa                        <  1:1748529/51‑1 (MQ=255)
                                                                   aaGCTCGCAGAGATCATACAGATGCTGCTCTTCACGCCCGCCCCAGTAAATg                     >  1:2665416/1‑52 (MQ=255)
                                                                   aaGCTCGCAGAGATCATACAGATGCTGCTCTTCACGCCCGCCCCAGTAAATGGTGATATCACGGTTTGGGt  <  1:3589373/71‑1 (MQ=255)
                                                                   aaGCTCGCAGAGATCATACAGATGCTGCTCTTCACGCCCGCCCCAGTAAATGGTGATATCACGGTTTGGGt  <  1:2562670/71‑1 (MQ=255)
                                                                   aaGCTCGCAGAGATCATACAGATGCTGCTCTTCACGCCCGCCCCAGTAAATGGTGATATCACGGTTTGGGt  <  1:4737813/71‑1 (MQ=255)
                                                                      |                                                                   
AGCCCGCTTCCGGTTGTTCAACCACCGGCACCACTTGCAGACCAGGATGCTTCAACGAAAGCGCCTCAAGCTCGCAGAGATCATACAGATGCTGCTCTTCACGCCCGCCCCAGTAAATGGTGATATCACGGTTTGGGT  >  W3110S.gb/3609632‑3609769

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: