Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,223,383 T→C 35.1% F165F (TTT→TTC panF → pantothenate:sodium symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,223,3830TC35.1% 32.5 / 29.5 54F165F (TTT→TTCpanFpantothenate:sodium symporter
Reads supporting (aligned to +/- strand):  ref base T (15/20);  new base C (8/11);  total (23/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTATCGGCGGTGCGCGCCTGCTGGAAACCGCGGCGGGTATTCCTTATGAAACCGGGCTGCTGATTTTTGGTATCAGCATTGCGTTATATACCGCCTTTGGTGGCTTTCGCGCCAGCGTGCTGAAC  >  minE/2223315‑2223440
                                                                    |                                                         
tttATCGGCGGTGCGCGCCTGCTGGAAACCGCGGCGGGTATTCCTTATGAAACCGGGCTGCTGAttttt                                                           >  1:2640657/1‑69 (MQ=255)
  tATCGGCGGTGCGCGCCTGCTGGAAACCGCGGCGGGTATTCCTTATGAAACCGGGCTGCTGATTTTTgg                                                         >  1:83460/1‑69 (MQ=255)
  tATCGGCGGTGCGCGCCTGCTGGAAACCGCGGCGGGTATTCCTTATGAAACCGGGCTGCTGATTTTTgg                                                         >  1:2007219/1‑69 (MQ=255)
  tATCGGCGGTGCGCGCCTGCTGGAAACCGCGGCGGGTATTCCTTATGAAACCGGGCTGCTGATTTTTgg                                                         >  1:2668449/1‑69 (MQ=255)
  tATCGGCGGTGCGCGCCTGCTGGAAACCGCGGCGGGTATTCCTTATGAAACCGGGCTGCTGATTTTTgg                                                         >  1:1833303/1‑69 (MQ=255)
  tATCGGCGGTGCGCGCCTGCTGGAAACCGCGGCGGGTATTCCTTATGAAACCGGGCTGCTGATTTTTgg                                                         >  1:3148200/1‑69 (MQ=255)
  tATCGGCGGTGCGCGCCTGCTGGAAACCGCGGCGGGTATTCCTTATGAAACCGGGCTGCTGATTTTTgg                                                         >  1:1351253/1‑69 (MQ=255)
  tATCGGCGGTGCGCGCCTGCTGGAAACCGCGGCGGGTATTCCTTATGAAACCGGGCTGCTGATTTTTgg                                                         >  1:767905/1‑69 (MQ=255)
         ggTGCGCGCCTGCTGGAAACCGCGGCGGGTATTCCTTCTGAAACCGGGCTGCTGATTTTTGGTATCAGc                                                  <  1:402326/69‑1 (MQ=255)
         ggTGCGCGCCTGCTGGAAACCGCGGCGGGTATTCCTTATGAAACCGGGCTGCTGATTTTTGGTATCAGc                                                  <  1:368979/69‑1 (MQ=255)
         ggTGCGCGCCTGCTGGAAACCGCGGCGGGTATTCCTTATGAAACCGGGCTGCTGATTTTTGGTATCAGc                                                  <  1:1909241/69‑1 (MQ=255)
         ggTGCGCGCCTGCTGGAAACCGCGGCGGGTATTCCTTATGAAACCGGGCTGCTGATTTTTGGTATCAGc                                                  <  1:1824287/69‑1 (MQ=255)
         ggTGCGCGCCTGCTGGAAACCGCGGCGGGTATTCCTTATGAAACCGGGCTGCTGATTTTTGGTATCAGc                                                  <  1:3125819/69‑1 (MQ=255)
         ggTGCGCGCCTGCTGGAAACCGCGGCGGGTATTCCTTATGAAACCGGGCTGCTGATTTTTGGTATCAGc                                                  <  1:1636142/69‑1 (MQ=255)
           tGCGCGCCTGCTGGAAACCGCGGCGGGTATTCCTTATGAAACCGGGCTGCTGATTTTTGGTATCAGCAt                                                <  1:356599/69‑1 (MQ=255)
           tGCGCGCCTGCTGGAAACCGCGGCGGGTATTCCTTATGAAACCGGGCTGCTGATTTTTGGTATCAGCAt                                                <  1:2623227/69‑1 (MQ=255)
            gcgcgcCTGCTGGAAACCGCGGCGGGTATTCCTTATGAAACCGGGCTGCTGATTTTTGGTATCAGCAtt                                               >  1:2154299/1‑69 (MQ=255)
                gcCTGCTGGAAACCGCGGCGGGTATTCCTTATGAAACCGGGCTGCTGATTTTTGGTATCAGCATTGCGt                                           <  1:1812441/69‑1 (MQ=255)
                    gctgGAAACCGCGGCGGGTATTCCTTATGAAACCGGGCTGCTGATTTTTGGTATCAGCAt                                                <  1:1112710/60‑1 (MQ=255)
                        gAAACCGCGGCGGGTATTCCTTATGAAACCGGGCTGCTGATTTTTGGTATCAGc                                                  <  1:1005968/54‑1 (MQ=255)
                        gAAACCGCGGCGGGTATTCCTTATGAAACCGGGCTGCTGATTTTTGGTATCAGCATTGCGTTATATAcc                                   >  1:3342216/1‑69 (MQ=255)
                               cggcggGTATTCCTTATGAAACCGGGCTGCTGATTTTCGGTATCAGCATTGCGTTGTATACCGCCTTtg                            <  1:3416784/69‑1 (MQ=255)
                               cggcggGTATTCCTTATGAAACCGGGCTGCTGATTTTCGGTATCAGCATTGCGTTGTATACCGCCTTtg                            <  1:540768/69‑1 (MQ=255)
                               cggcggGTATTCCTTATGAAACCGGGCTGCTGATTTTCGGTATCAGCATTGCGTTGTATACCGCCTTtg                            <  1:376917/69‑1 (MQ=255)
                               cggcggGTATTCCTTATGAAACCGGGCTGCTGATTTTCGGTATCAGCATTGCGTTGTATACCGCCTTtg                            <  1:2487215/69‑1 (MQ=255)
                               cggcggGTATTCCTTATGAAACCGGGCTGCTGATTTTCGGTATCAGCATTGCGTTGTATACCGCCTTtg                            <  1:2554256/69‑1 (MQ=255)
                               cggcggGTATTCCTTATGAAACCGGGCTGCTGATTTTCGGTATCAGCATTGCGTTGTATACCGCCTTtg                            <  1:1049931/69‑1 (MQ=255)
                               cggcggGTATTCCTTATGAAACCGGGCTGCTGATTTTCGGTATCAGCATTGCGTTGTATACCGCCTTtg                            <  1:116908/69‑1 (MQ=255)
                               cggcggGTATTCCTTATGAAACCGGGCTGCTGATTTTCGGTATCAGCATTGCGTTGTATACCGCCTTtg                            <  1:1744701/69‑1 (MQ=255)
                               cggcggGTATTCCTTATGAAACCGGGCTGCTGATTTTCGGTATCAGCATTGCGTTGTATACCGCCTTtg                            <  1:1289121/69‑1 (MQ=255)
                               cggcggGTATTCCTTATGAAACCGGGCTGCTGATTTTCGGTATCAGCATTGCGTTGTATACCGCCTTtg                            <  1:1356164/69‑1 (MQ=255)
                               cggcggGTATTCCTTATGAAACCGGGCTGCTGATTTTCGGTATCAGCATTGCGTTGTATACCGCCTTtg                            <  1:1624875/69‑1 (MQ=255)
                                   gggTATTCCTTATGAAACCGGGCTGCTGATTTTTGGTATCAGCATTGCGTTATATACCGCCTTtggtgg                        >  1:1054765/1‑69 (MQ=255)
                                          ccTTATGAAACCGGGCTGCTGATTTTTGGTATCAGCATTGCGTTATATACCGCCTTTGGTGGCTTTcgc                 >  1:798873/1‑69 (MQ=255)
                                          ccTTATGAAACCGGGCTGCTGATTTTTGGTATCAGCATTGCGTTATATACCGCCTTTGGTGGCTTTcgc                 >  1:2000214/1‑69 (MQ=255)
                                          ccTTATGAAACCGGGCTGCTGATTTTTGGTATCAGCATTGCGTTATATACCGCCTTTGGTGGCTTTcgc                 >  1:1942267/1‑69 (MQ=255)
                                          ccTTATGAAACCGGGCTGCTGATTTTTGGTATCAGCATTGCGTTATATACCGCCTTTGGTGGCTTTcgc                 >  1:3249285/1‑69 (MQ=255)
                                          ccTTATGAAACCGGGCTGCTGATTTTTGGTATCAGCATTGCGTTATATACCGCCTTTGGTGGCTTTcgc                 >  1:150459/1‑69 (MQ=255)
                                           cTTATGAAACCGGGCTGCTGATTTTTGGTATCAGc                                                  <  1:2538542/35‑1 (MQ=255)
                                             tATGAAGCCGGGCTGCTGATTTTCGGTATCAGCATTGCGTTGTATACCGCCttt                             >  1:436529/1‑54 (MQ=255)
                                             tATGAAACCGGGCTGCTGATTTTCGGTATCAGCATTGCGTTGTATACCGCCttt                             >  1:169972/1‑54 (MQ=255)
                                             tATGAAACCGGGCTGCTGATTTTCGGTATCAGCATTGCGTTGTATACCGCCttt                             >  1:127454/1‑54 (MQ=255)
                                             tATGAAACCGGGCTGCTGATTTTCGGTATCAGCATTGCGTTGTATACCGCCttt                             >  1:2409259/1‑54 (MQ=255)
                                             tATGAAACCGGGCTGCTGATTTTCGGTATCAGCATTGCGTTGTATACCGCCttt                             >  1:1816637/1‑54 (MQ=255)
                                             tATGAAACCGGGCTGCTGATTTTCGGTATCAGCATTGCGTTGTATACCGCCttt                             >  1:2860052/1‑54 (MQ=255)
                                             tATGAAACCGGGCTGCTGATTTTCGGTATCAGCATTGCGTTGTATACCGCCttt                             >  1:2337765/1‑54 (MQ=255)
                                             tATGAAAACGGGCTGCTGATTTTCGGTATCAGCATTGCGTTGTATAc                                    >  1:797007/1‑47 (MQ=38)
                                                gAAACCGGGCTGCTGATTTTTGGTATCAGCATTGCGTTATATACCGCCTTTGTTGGCTTTCGCGCCAg            <  1:3178958/68‑1 (MQ=255)
                                                gAAACCGGGCTGCTGATTTTTGGTATCAGCATTGCGTTATATACCGCCTTTGGTGGCTTTCGCGCCAg            <  1:673971/68‑1 (MQ=255)
                                                      gggCTGCTGATTTTTGGTATCAGCATTGCGTTATATACCGCCTTTGGTGGCTTTCGCGCCAGCGTGCTg     <  1:3241884/69‑1 (MQ=255)
                                                      gggCTGCTGATTTTTGGTATCAGCATTGCGTTATATACCGCCTTTGGTGGCTTTCGCGCCAGCGTGCTg     <  1:553877/69‑1 (MQ=255)
                                                       ggCTGCTGATTTTTGGTATCAGCATTGCGTTATATACCGCCTTTGGTGGCTTTCGCGCCAGCGTGCTg     <  1:563470/68‑1 (MQ=255)
                                                       ggCTGCTGATTTTTGGTATCAGCATTGCGTTATATACCGCCTTTGGTGGCTTTCGCGCCAGCGTGCTg     <  1:2266077/68‑1 (MQ=255)
                                                         ctgctgATTTTTGGTATCAGCATTGCGTTATATACCGCCTTTGGTGGCTTTCGCGCCAGCGTGCTGAAc  <  1:3037812/69‑1 (MQ=255)
                                                         ctgctgATTTTTGGTATCAGCATTGCGTTATATACCGCCTTTGGTGGCTTTCGCGCCAGCGTGCTGAAc  <  1:2165985/69‑1 (MQ=255)
                                                                    |                                                         
TTTATCGGCGGTGCGCGCCTGCTGGAAACCGCGGCGGGTATTCCTTATGAAACCGGGCTGCTGATTTTTGGTATCAGCATTGCGTTATATACCGCCTTTGGTGGCTTTCGCGCCAGCGTGCTGAAC  >  minE/2223315‑2223440

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: