Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,403,261 A→G 39.5% A368A (GCA→GCG rtn → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,403,2610AG39.5% 17.8 / 30.3 48A368A (GCA→GCGrtnconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (27/2);  new base G (9/10);  total (36/12)
Fisher's exact test for biased strand distribution p-value = 5.71e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAGGCGTCAAATTTGGTATTAACATTGCGCCGGACCATCTGCACAGCGAAAGCTTTAAAGCAGATATCCAGAAACTGCTCACTTCCCTGCCCGCACACCATTTCCAGAT  >  minE/1403200‑1403308
                                                             |                                               
tAGGCGTCAAATTTGGTATTAACATTGCGCCGGACCATCTGCACAGCGAAAGCTTTAAAGCAGata                                             >  1:584142/1‑66 (MQ=255)
tAGGCGTCAAATTTGGTATTAACATTGCGCCGGACCATCTGCACAGCGAAAGCTTTAAAGCAGata                                             >  1:2893130/1‑66 (MQ=255)
tAGGCGTCAAATTTGGTATTAACATTGCGCCGGACCATCTGCACAGCGAAAGCTTTAAAGCAGata                                             >  1:2114244/1‑66 (MQ=255)
 aGGCGTCAAATTTGGTATTAAGATTGCGCCGGACCATCTGCACAGCGAAAGCTTTAAAGCAGATATcc                                          >  1:2294002/1‑68 (MQ=255)
 aGGCGTCAAATTTGGTATTAACATTGCGCCGGACCATCTGCACAGCGAAAGCTTTAAAGCAGATATcc                                          >  1:2297790/1‑68 (MQ=255)
 aGGCGTCAAATTTGGTATTAACATTGCGCCGGACCATCTGCACAGCGAAAGCTTTAAAGCAGATATcc                                          >  1:2782946/1‑68 (MQ=255)
 aGGCGTCAAATTTGGTATTAACATTGCGCCGGACCATCTGCACAGCGAAAGCTTTAAAGCAGATATcc                                          >  1:3212364/1‑68 (MQ=255)
     gTCAAATTTGGTATTAACATTGCGCCGGACCATCTGCACAGCGAAAGCTTTAAAGCAGATATCCAGaaa                                     >  1:3369020/1‑69 (MQ=255)
     gTCAAATTTGGTATTAACATTGCGCCGGACCATCTGCACAGCGAAAGCTTTAAAGCAGATATCCAGaaa                                     >  1:3193555/1‑69 (MQ=255)
     gTCAAATTTGGTATTAACATTGCGCCGGACCATCTGCACAGCGAAAGCTTTAAAGCAGATATCCAGaaa                                     >  1:1796584/1‑69 (MQ=255)
       cAAATTTGGCATTAACATTGCGCCAGCCCACTTGCACAGCGAAAGCTTTAAAGAGATATCCAGAAACTg                                  <  1:1474096/69‑1 (MQ=255)
        aaaTTTGGCATTAACATTGCGCCAGCCCACTTGCACAGCGAAAGCTTTAAAGCGGATATCCAGAAACTg                                  <  1:1891949/69‑1 (MQ=255)
        aaaTTTGGCATTAACATTGCGCCAGCCCACTTGCACAGCGAAAGCTTTAAAGCGGATATCCAGAAACTg                                  <  1:2755300/69‑1 (MQ=255)
        aaaTTTGGCATTAACATTGCGCCAGCCCACTTGCACAGCGAAAGCTTTAAAGCGGATATCCAGAAACTg                                  <  1:3287216/69‑1 (MQ=255)
        aaaTTTGGCATTAACATTGCGCCAGCCCACTTGCACAGCGAAAGCTTTAAAGCGGATATCCAGAAACTg                                  <  1:1875353/69‑1 (MQ=255)
        aaaTTTGGCATTAACATTGCGCCAGCCCACTTGCACAGCGAAAGCTTTAAAGCGGATATCCAGAAACTg                                  <  1:430281/69‑1 (MQ=255)
         aaTTTGGCATTAACATTGCGCCAGCCCACTTGCACAGCGAAAGCTTTAAAGCGGATATCCAGAAACTg                                  <  1:1924519/68‑1 (MQ=255)
         aaTTTGGCATTAACATTGCGCCAGCCCACTTGCACAGCGAAAGCTTTAAAGCGGATATCCAGAAACTg                                  <  1:2808129/68‑1 (MQ=255)
         aaTTTGGCATTAACATTGCGCCAGCCCACTTGCACAGCGAAAGCTTTAAAGCGGATATCCAGAAACTg                                  <  1:13799/68‑1 (MQ=255)
            ttGGCATTAACATTGCGCCAGCCCACTTGCACAGCGAAAGCTTTAAAGCGGATATCCAGAAACTg                                  <  1:841413/65‑1 (MQ=255)
                   tAACATTGCGCCAGCCCACTTGCACAGCGAAAGCTTTAAAGCGGATATCCAGAAACTGCTCACTTccc                        >  1:2273725/1‑68 (MQ=255)
                   tAACATTGCGCCAGCCCACTTGCACAGCGAAAGCTTTAAAGCGGATATCCAGAAACTGCTCACTTCCCt                       >  1:2410533/1‑69 (MQ=255)
                   tAACATTGCGCCAGCCCACTTGCACAGCGAAAGCTTTAAAGCGGATATCCAGAAACTGCTCACTTCCCt                       >  1:3348648/1‑69 (MQ=255)
                   tAACATTGCGCCAGCCCACTTGCACAGCGAAAGCTTTAAAGCGGATATCCAGAAACTGCTCACTTCCCt                       >  1:813402/1‑69 (MQ=255)
                   tAACATTGCGCCAGCCCACTTGCACAGCGAAAGCTTTAAAGCGGATATCCAGAAACTGCTCACTTCCCt                       >  1:86670/1‑69 (MQ=255)
                   tAACATTGCGCCAGCCCACTTGCACAGCGAAAGCTTTAAAGCGGATATCCAGAAACTGCTCACTTCCCt                       >  1:2221960/1‑69 (MQ=255)
                   tAACATTGCGCCAGCCCACTTGCACAGCGAAAGCTTTAAAGCGGATATCCAGAAACTGCTCACTTCCCt                       >  1:2196782/1‑69 (MQ=255)
                   tAACATTGCGCCAGCCCACTTGCACAGCGAAAGCTTTAAAGCGGATATCCAGAAACTGCTCACTTCCCt                       >  1:1399758/1‑69 (MQ=255)
                   tAACATTGCGCCAGCCCACTTGCACAGCGAAAGCATTAAAGCGGaaat                                            >  1:2885354/1‑48 (MQ=255)
                    aaCATTGCGCCGGACCATCTGCACAGCGAGAGCTTTAAAGCAGATATCCAGa                                       >  1:1358038/1‑52 (MQ=255)
                    aaCATTGCGCCGGACCATCTGCACAGCGAAAGCTTTAAAGCAGATATCCAGa                                       >  1:2478745/1‑52 (MQ=255)
                    aaCATTGCGCCGGACCATCTGCACAGCGAAAGCTTTAAAGCAGATATCCAGa                                       >  1:2643961/1‑52 (MQ=255)
                    aaCATTGCGCCGGACCATCTGCACAGCGAAAGCTTTAAAGCAGATATCCAGa                                       >  1:3236588/1‑52 (MQ=255)
                    aaCATTGCGCCGGACCATCTGCACAGCGAAAGCTTTAAAGCAGATATCCAGAAACTGCTCACTTCCCt                       >  1:1358902/1‑68 (MQ=255)
                    aaCATTGCGCCGGACCATCTGCACAGCGAAAGCTTTAAAGCAGATATCCAGAAACTGCTCACTTCCCTg                      >  1:1069973/1‑69 (MQ=255)
                    aaCATTGCGCCGGACCATCTGCACAGCGAAAGCTTTAAAGCAGATATCCAGAAACTGCTCACTTCCCTg                      >  1:2210302/1‑69 (MQ=255)
                    aaCATTGCGCCGGACCATCTGCACAGCGAAAGCTTTAAAGCAGATATCCAGAAACTGCTCACTTCCCTg                      >  1:626870/1‑69 (MQ=255)
                    aaCATTGCGCCGGACCATCTGCACAGCGAAAGCTTTAAAGCAGATATCCAGAAACTGCTCACTTCCCTg                      >  1:1411727/1‑69 (MQ=255)
                    aaCATTGCGCCGGACCATCTGCACAGCGAAAGCTATAAAGCAGATATCCAGa                                       >  1:3247264/1‑52 (MQ=255)
                    aaCATTGCGCCGGACCATCTGAACAGCGAAAGCTTTAAAGCAGATATcc                                          >  1:2636501/1‑49 (MQ=255)
                          gcgcCGGACCATCTGCACAGCGAAAGCTTTAAAGCAGATATCCAGACACTGCTCACTTCCCTGCCCGc                 >  1:1085938/1‑68 (MQ=255)
                          gcgcCGGACCATCTGCACAGCGAAAGCTTTAAAGCAGATATCCAGAAACTGCTCACTTCCCTGCCCGca                >  1:1931908/1‑69 (MQ=255)
                          gcgcCGGACCATCTGCACAGCGAAAGCTTTAAAGCAGATATCCAGAAACTGCTCACTTCCCTGCCCGca                >  1:1987304/1‑69 (MQ=255)
                          gcgcCGGACCATCTGCACAGCGAAAGCTTTAAAGCAGATATCCAGAAACTGCTCACTTCCCTGCCCGca                >  1:2156103/1‑69 (MQ=255)
                          gcgcCGGACCATCTGCACAGCGAAAGCTTTAAAGCAGATATCCAGAAACTGCTCACTTCCCTGCCCGc                 >  1:2940371/1‑68 (MQ=255)
                                       tGCACAGCGAAAGCTTTAAAGCAGATATCCAGAAACTGCTCACTTCCCTGCCCGCACACCATTTCCAGa   >  1:1014638/1‑69 (MQ=255)
                                        gCACAGCGAAAGCTTTAAAGCGGATATCCAGAAACTg                                  <  1:1652148/37‑1 (MQ=255)
                                           caGCGAAAGCTTTAAAGCAGATATCCAGAAACTGCTCACTTCCCTGCCCGCACACCATTTCCAGAt  <  1:3008258/66‑1 (MQ=255)
                                                             |                                               
TAGGCGTCAAATTTGGTATTAACATTGCGCCGGACCATCTGCACAGCGAAAGCTTTAAAGCAGATATCCAGAAACTGCTCACTTCCCTGCCCGCACACCATTTCCAGAT  >  minE/1403200‑1403308

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