Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,403,233 A→C 42.9% D359A (GAC→GCC)  rtn → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,403,2330AC42.9% 1.7 / 49.5 49D359A (GAC→GCC) rtnconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (26/2);  new base C (9/12);  total (35/14)
Fisher's exact test for biased strand distribution p-value = 2.33e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.63e-01

GCGATGCCGCAGAATTAGAAAAAGTGCTGCCGGTAGGCGTCAAATTTGGTATTAACATTGCGCCGGACCATCTGCACAGCGAAAGCTTTAAAGCAGATATCCAGAAACTGCTCACTTCCCTGCCCGCA  >  minE/1403167‑1403294
                                                                  |                                                             
gcgATGCCGCAGAATTAGAAAAAGTGCTGCCGGTAGGCGTCAAATTTGGTATTAACATTGCGCCGGAcc                                                             <  1:3073422/69‑1 (MQ=255)
 cgATGCCGCAGATTTAGAAAAAGTGCTGCCGGTAGGCGTCAAATTTGGTATTAACATTGCGCCGGAcc                                                             <  1:3043012/68‑1 (MQ=255)
            aaTTAGAAAAAGTACTGCCGGTAGGCGTCAAATTTGGCATTAACATTGCGCCAGCCCACTTGCACAGCg                                                 <  1:3341182/69‑1 (MQ=255)
            aaTTAGAAAAAGTACTGCCGGTAGGCGTCAAATTTGGCATTAACATTGCGCCAGCCCACTTGCACAGCg                                                 <  1:1377287/69‑1 (MQ=255)
                                 tAGGCGTCAAATTTGGTATTAACATTGCGCCGGACCATCTGCACAGCGAAAGCTTTAAAGCAGata                               >  1:2893130/1‑66 (MQ=255)
                                 tAGGCGTCAAATTTGGTATTAACATTGCGCCGGACCATCTGCACAGCGAAAGCTTTAAAGCAGata                               >  1:2114244/1‑66 (MQ=255)
                                 tAGGCGTCAAATTTGGTATTAACATTGCGCCGGACCATCTGCACAGCGAAAGCTTTAAAGCAGata                               >  1:584142/1‑66 (MQ=255)
                                  aGGCGTCAAATTTGGTATTAAGATTGCGCCGGACCATCTGCACAGCGAAAGCTTTAAAGCAGATATcc                            >  1:2294002/1‑68 (MQ=255)
                                  aGGCGTCAAATTTGGTATTAACATTGCGCCGGACCATCTGCACAGCGAAAGCTTTAAAGCAGATATcc                            >  1:2297790/1‑68 (MQ=255)
                                  aGGCGTCAAATTTGGTATTAACATTGCGCCGGACCATCTGCACAGCGAAAGCTTTAAAGCAGATATcc                            >  1:3212364/1‑68 (MQ=255)
                                  aGGCGTCAAATTTGGTATTAACATTGCGCCGGACCATCTGCACAGCGAAAGCTTTAAAGCAGATATcc                            >  1:2782946/1‑68 (MQ=255)
                                      gTCAAATTTGGTATTAACATTGCGCCGGACCATCTGCACAGCGAAAGCTTTAAAGCAGATATCCAGaaa                       >  1:3193555/1‑69 (MQ=255)
                                      gTCAAATTTGGTATTAACATTGCGCCGGACCATCTGCACAGCGAAAGCTTTAAAGCAGATATCCAGaaa                       >  1:3369020/1‑69 (MQ=255)
                                      gTCAAATTTGGTATTAACATTGCGCCGGACCATCTGCACAGCGAAAGCTTTAAAGCAGATATCCAGaaa                       >  1:1796584/1‑69 (MQ=255)
                                        cAAATTTGGCATTAACATTGCGCCAGCCCACTTGCACAGCGAAAGCTTTAAAGAGATATCCAGAAACTg                    <  1:1474096/69‑1 (MQ=255)
                                         aaaTTTGGCATTAACATTGCGCCAGCCCACTTGCACAGCGAAAGCTTTAAAGCGGATATCCAGAAACTg                    <  1:430281/69‑1 (MQ=255)
                                         aaaTTTGGCATTAACATTGCGCCAGCCCACTTGCACAGCGAAAGCTTTAAAGCGGATATCCAGAAACTg                    <  1:1891949/69‑1 (MQ=255)
                                         aaaTTTGGCATTAACATTGCGCCAGCCCACTTGCACAGCGAAAGCTTTAAAGCGGATATCCAGAAACTg                    <  1:1875353/69‑1 (MQ=255)
                                         aaaTTTGGCATTAACATTGCGCCAGCCCACTTGCACAGCGAAAGCTTTAAAGCGGATATCCAGAAACTg                    <  1:2755300/69‑1 (MQ=255)
                                         aaaTTTGGCATTAACATTGCGCCAGCCCACTTGCACAGCGAAAGCTTTAAAGCGGATATCCAGAAACTg                    <  1:3287216/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:3348648/1‑69 (MQ=255)
                                                    tAACATTGCGCCAGCCCACTTGCACAGCGAAAGCTTTAAAGCGGATATCCAGAAACTGCTCACTTCCCt         >  1:2410533/1‑69 (MQ=255)
                                                    tAACATTGCGCCAGCCCACTTGCACAGCGAAAGCTTTAAAGCGGATATCCAGAAACTGCTCACTTCCCt         >  1:2221960/1‑69 (MQ=255)
                                                    tAACATTGCGCCAGCCCACTTGCACAGCGAAAGCTTTAAAGCGGATATCCAGAAACTGCTCACTTCCCt         >  1:2196782/1‑69 (MQ=255)
                                                    tAACATTGCGCCAGCCCACTTGCACAGCGAAAGCTTTAAAGCGGATATCCAGAAACTGCTCACTTCCCt         >  1:813402/1‑69 (MQ=255)
                                                    tAACATTGCGCCAGCCCACTTGCACAGCGAAAGCTTTAAAGCGGATATCCAGAAACTGCTCACTTCCCt         >  1:1399758/1‑69 (MQ=255)
                                                    tAACATTGCGCCAGCCCACTTGCACAGCGAAAGCTTTAAAGCGGATATCCAGAAACTGCTCACTTCCCt         >  1:86670/1‑69 (MQ=255)
                                                    tAACATTGCGCCAGCCCACTTGCACAGCGAAAGCATTAAAGCGGaaat                              >  1:2885354/1‑48 (MQ=255)
                                                     aaCATTGCGCCGGACCATCTGCACAGCGAGAGCTTTAAAGCAGATATCCAGa                         >  1:1358038/1‑52 (MQ=255)
                                                     aaCATTGCGCCGGACCATCTGCACAGCGAAAGCTTTAAAGCAGATATCCAGa                         >  1:3236588/1‑52 (MQ=255)
                                                     aaCATTGCGCCGGACCATCTGCACAGCGAAAGCTTTAAAGCAGATATCCAGa                         >  1:2478745/1‑52 (MQ=255)
                                                     aaCATTGCGCCGGACCATCTGCACAGCGAAAGCTTTAAAGCAGATATCCAGa                         >  1:2643961/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:2156103/1‑69 (MQ=255)
                                                           gcgcCGGACCATCTGCACAGCGAAAGCTTTAAAGCAGATATCCAGAAACTGCTCACTTCCCTGCCCGca  >  1:1987304/1‑69 (MQ=255)
                                                           gcgcCGGACCATCTGCACAGCGAAAGCTTTAAAGCAGATATCCAGAAACTGCTCACTTCCCTGCCCGca  >  1:1931908/1‑69 (MQ=255)
                                                           gcgcCGGACCATCTGCACAGCGAAAGCTTTAAAGCAGATATCCAGAAACTGCTCACTTCCCTGCCCGc   >  1:2940371/1‑68 (MQ=255)
                                                                  |                                                             
GCGATGCCGCAGAATTAGAAAAAGTGCTGCCGGTAGGCGTCAAATTTGGTATTAACATTGCGCCGGACCATCTGCACAGCGAAAGCTTTAAAGCAGATATCCAGAAACTGCTCACTTCCCTGCCCGCA  >  minE/1403167‑1403294

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: