Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,892,019:1 +GC 100% intergenic (‑437/+54) gudP ← / ← ygdL predicted D‑glucarate transporter/conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,892,0191.G97.9% 159.6 / ‑3.4 47intergenic (‑437/+54)gudP/ygdLpredicted D‑glucarate transporter/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (21/25);  total (22/25)
Fisher's exact test for biased strand distribution p-value = 4.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
*minE1,892,0192.C97.9% 158.6 / ‑3.7 47intergenic (‑437/+54)gudP/ygdLpredicted D‑glucarate transporter/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (21/25);  total (22/25)
Fisher's exact test for biased strand distribution p-value = 4.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CTGGCATTTGCCTGAAAAGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCA‑‑TAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACC  >  minE/1891951‑1892078
                                                                      ||                                                           
cTGGCTTTTGCCTGAAAAGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCag                                                              <  1:3162540/70‑2 (MQ=255)
            tGAAAAGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAAc                                                 >  1:3022377/1‑71 (MQ=255)
            tGAAAAGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAAc                                                 >  1:3547098/1‑71 (MQ=255)
             gAAAAGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAAcc                                                >  1:135972/1‑71 (MQ=255)
              aaaaGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAAcc                                                <  1:876678/70‑1 (MQ=255)
              aaaaGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCg                                               >  1:336944/1‑71 (MQ=255)
               aaaGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCgg                                              <  1:3582697/71‑1 (MQ=255)
               aaaGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCgg                                              >  1:3092738/1‑71 (MQ=255)
                aaGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCggg                                             >  1:4124870/1‑71 (MQ=255)
                aaGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCggg                                             >  1:1467700/1‑71 (MQ=255)
                aaGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCggg                                             >  1:2025014/1‑71 (MQ=255)
                aaGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCggg                                             >  1:3404617/1‑71 (MQ=255)
                aaGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCggg                                             >  1:1203488/1‑71 (MQ=255)
                aaGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCgg                                              <  1:3750739/70‑1 (MQ=255)
                       gattcAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATa                                      <  1:1502913/71‑1 (MQ=255)
                       gattcAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATa                                      <  1:1775677/71‑1 (MQ=255)
                        attcAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATa                                      <  1:873904/70‑1 (MQ=255)
                           cAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGa                                   <  1:387865/70‑1 (MQ=255)
                           cAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGa                                   <  1:1152587/70‑1 (MQ=255)
                            aGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGa                                                   <  1:717044/53‑1 (MQ=255)
                            aGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGa                                                   <  1:3512737/53‑1 (MQ=255)
                            aGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGa                                                   <  1:2647261/53‑1 (MQ=255)
                            aGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGa                                                   <  1:2845215/53‑1 (MQ=255)
                            aGCACCCCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGa                                                   <  1:2282030/53‑1 (MQ=255)
                              cACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGtt                               >  1:36675/1‑71 (MQ=255)
                                ccGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGaa                                                  >  1:1064891/1‑50 (MQ=255)
                                ccGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGaa                                                  >  1:2592759/1‑50 (MQ=255)
                                  gCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTcgcg                           >  1:566275/1‑71 (MQ=255)
                                    tAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGt                          <  1:3734112/70‑1 (MQ=255)
                                     aaaaCGACATTTACCGCTCGCTGAACATATAGC‑‑TAAGCACAGAACCg                                               >  1:1243749/1‑47 (MQ=255)
                                     aaaaCGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCg                                               >  1:467656/1‑48 (MQ=255)
                                      aaaCGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGTATAAGACGTTCGCGtc                         <  1:1413482/69‑1 (MQ=255)
                                       aaCGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGtc                         <  1:3072304/68‑1 (MQ=255)
                                          gACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCgg                                         <  1:1357234/49‑1 (MQ=255)
                                            cATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCg                  <  1:2517524/70‑1 (MQ=255)
                                              tttACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATaa                                     <  1:2044171/49‑1 (MQ=255)
                                              tttACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATaa                                     <  1:1816954/49‑1 (MQ=255)
                                              tttACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATaa                                     <  1:4079426/49‑1 (MQ=255)
                                              tttACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATaa                                     <  1:3086344/49‑1 (MQ=255)
                                                tACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATaa                                     >  1:1440787/1‑47 (MQ=255)
                                                tACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATaa                                     >  1:2562933/1‑47 (MQ=255)
                                                tACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATaa                                     >  1:2063704/1‑47 (MQ=255)
                                                 aCCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTg             <  1:2150612/70‑1 (MQ=255)
                                                   cgctcgctGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTcgc                            <  1:2637422/53‑1 (MQ=255)
                                                   cgctcgctGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGAtt          >  1:3970011/1‑71 (MQ=255)
                                                   cgctcgctGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTAGCGtct                        >  1:2518707/1‑57 (MQ=255)
                                                        gctGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGttt      <  1:1240862/70‑1 (MQ=38)
                                                           gAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAAcc  >  1:2494156/1‑71 (MQ=255)
                                                                      ||                                                           
CTGGCATTTGCCTGAAAAGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCA‑‑TAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACC  >  minE/1891951‑1892078

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: