Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 115,012 C→T 25.0% G330D (GGC→GAC)  yacH ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE115,0120CT25.0% 189.5 / 55.5 151G330D (GGC→GAC) yacHhypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (30/83);  new base T (14/24);  total (44/107)
Fisher's exact test for biased strand distribution p-value = 3.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.18e-01

Only 100 of 194 total aligned reads displayed.

ATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATG  >  minE/114945‑115079
                                                                   |                                                                   
aTGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTgccagcc                                                                    <  1:2208708/69‑1 (MQ=255)
aTGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTgccagcc                                                                    <  1:2237007/69‑1 (MQ=255)
aTGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTgccagcc                                                                    <  1:559016/69‑1 (MQ=255)
aTGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTgccagcc                                                                    <  1:1170359/69‑1 (MQ=255)
aTGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTgccagcc                                                                    <  1:507164/69‑1 (MQ=255)
aTGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTgccagcc                                                                    <  1:1436772/69‑1 (MQ=255)
aTGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTgccagcc                                                                    <  1:3073983/69‑1 (MQ=255)
aTGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTgccagcc                                                                    <  1:1386390/69‑1 (MQ=255)
aTGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTgccagc                                                                     <  1:1552474/68‑1 (MQ=255)
   cTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCAtt                                                                 >  1:2345764/1‑69 (MQ=255)
   cTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCAtt                                                                 >  1:271746/1‑69 (MQ=255)
   cTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCAtt                                                                 >  1:3182313/1‑69 (MQ=255)
       ccGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTgccagcc                                                                    <  1:3420053/62‑1 (MQ=255)
             aTACGGTTAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGAt                                                      <  1:485836/69‑1 (MQ=255)
               aCGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGAt                                                      <  1:1708893/68‑1 (MQ=255)
                    tGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCTTGTCATTACCGTCACGa                                                       >  1:2384072/1‑62 (MQ=255)
                    tGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGTCATTACCGTCACGATGACCGc                                                >  1:3057605/1‑69 (MQ=255)
                    tGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGTCATTACCGTCACGATGACCGc                                                >  1:642526/1‑69 (MQ=255)
                    tGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGTCATTACCGTCACGATGACCGc                                                >  1:2491520/1‑69 (MQ=255)
                    tGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGTCATTACCGTCACGATGACCGc                                                >  1:702393/1‑69 (MQ=255)
                    tGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGTCATTACCGTCACGATGACCGc                                                >  1:2185309/1‑69 (MQ=255)
                    tGAAATTATTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGTCATTACCGTCACGATGACCGc                                                >  1:3384323/1‑69 (MQ=255)
                     gAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACg                                                        >  1:843937/1‑60 (MQ=255)
                        aTTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATc                                            >  1:90331/1‑69 (MQ=255)
                        aTTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATc                                            >  1:930872/1‑69 (MQ=255)
                        aTTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATc                                            >  1:1987882/1‑69 (MQ=255)
                        aTTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATc                                            >  1:972467/1‑69 (MQ=255)
                        aTTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATc                                            >  1:1289656/1‑69 (MQ=255)
                        aTTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATc                                            >  1:1502408/1‑69 (MQ=255)
                        aTTGTTGACGTCGATATTGATATTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATc                                            >  1:2086222/1‑69 (MQ=255)
                         ttgttgACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCg                                           >  1:2554561/1‑69 (MQ=255)
                           gttgACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCgtg                                         <  1:3163503/69‑1 (MQ=255)
                           gttgACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCgtg                                         <  1:1586519/69‑1 (MQ=255)
                           gttgACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCgtg                                         <  1:1047978/69‑1 (MQ=255)
                           gttgACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCgtg                                         <  1:1009083/69‑1 (MQ=255)
                           gttgACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCgtg                                         <  1:1890019/69‑1 (MQ=255)
                                    gATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATa                                <  1:2496286/69‑1 (MQ=255)
                                    gATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATa                                <  1:1104342/69‑1 (MQ=255)
                                    gATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATa                                <  1:2067587/69‑1 (MQ=255)
                                    gATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATa                                <  1:597361/69‑1 (MQ=255)
                                    gATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATa                                <  1:53794/69‑1 (MQ=255)
                                    gATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATa                                <  1:2692444/69‑1 (MQ=255)
                                          gATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGAt                                                      <  1:3336679/41‑1 (MQ=255)
                                              ttgtCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAatcatc                          >  1:81848/1‑65 (MQ=255)
                                                gtCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGt                     <  1:2538106/68‑1 (MQ=255)
                                                gtCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGt                     <  1:3288298/68‑1 (MQ=255)
                                                gtCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGt                     <  1:1682473/68‑1 (MQ=255)
                                                gtCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGt                     <  1:1519685/68‑1 (MQ=255)
                                                gtCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGt                     <  1:186291/68‑1 (MQ=255)
                                                gtCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGt                     <  1:1790984/68‑1 (MQ=255)
                                                gtCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTc                    >  1:2888045/1‑69 (MQ=255)
                                                 tCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGt                     <  1:1228897/67‑1 (MQ=255)
                                                   gCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATa                                <  1:1030700/54‑1 (MQ=255)
                                                     cGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCatgat               <  1:436357/69‑1 (MQ=255)
                                                     cGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCatgat               <  1:1123283/69‑1 (MQ=255)
                                                     cGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCatgat               <  1:784421/69‑1 (MQ=255)
                                                     cGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCatgat               <  1:1005936/69‑1 (MQ=255)
                                                     cGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCatgat               <  1:3165147/69‑1 (MQ=255)
                                                     cGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCatgat               <  1:1671110/69‑1 (MQ=255)
                                                     cGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCatgat               <  1:2142976/69‑1 (MQ=255)
                                                      gTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCatgat               <  1:672516/68‑1 (MQ=255)
                                                      gTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCatgat               <  1:1359865/68‑1 (MQ=255)
                                                        tgtgTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCatgatgat            <  1:2876296/69‑1 (MQ=255)
                                                        tgtgTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCatgatgat            <  1:1313446/69‑1 (MQ=255)
                                                        tgtgTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCatgatgat            <  1:2057360/69‑1 (MQ=255)
                                                        tgtgTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCatgatgat            <  1:1349372/69‑1 (MQ=255)
                                                        tgtgTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCatgatgat            <  1:2765861/69‑1 (MQ=255)
                                                        tgtgTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCatgatgat            <  1:2038075/69‑1 (MQ=255)
                                                        tgtgTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCatgatgat            <  1:1651270/69‑1 (MQ=255)
                                                        tgtgTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCatgatgat            <  1:848125/69‑1 (MQ=255)
                                                        tgtgTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCatgatgat            <  1:29062/69‑1 (MQ=255)
                                                         gtgtTGCCAGCCATTAGCGTCACGATGACCGCCATCgtg                                         <  1:1968010/39‑1 (MQ=255)
                                                           gtTGCCAGTCATTACCGTCACGATGACCGCCATCg                                           <  1:768832/35‑1 (MQ=255)
                                                           gtTGCCAGTCATTACCGTCACGATGACCGCCATCg                                           <  1:1951714/35‑1 (MQ=255)
                                                           gtTGCCAGTCATTACCGTCACGATGACCGCCATCg                                           <  1:1068140/35‑1 (MQ=255)
                                                           gtTGCCAGTCATTACCGTCACGATGACCGCCATCg                                           <  1:657771/35‑1 (MQ=255)
                                                           gtTGCCAGTCATTACCGTCACGATGACCGCCATCg                                           <  1:96718/35‑1 (MQ=255)
                                                           gtTGCCAGTCATTACCGTCACGATGACCGCCATCg                                           <  1:2843838/35‑1 (MQ=255)
                                                           gtTGCCAGTCATTACCGTCACGATGACCGCCATCg                                           <  1:1740043/35‑1 (MQ=255)
                                                           gtTGCCAGTCATTACCGTCACGATGACCGCCATCg                                           <  1:2240483/35‑1 (MQ=255)
                                                           gtTGCCAGTCATTACCGTCACGATGACCGCCATCg                                           <  1:383299/35‑1 (MQ=255)
                                                           gtTGCCAGTCATTACCGTCACGATGACCGCCATCg                                           <  1:2360354/35‑1 (MQ=255)
                                                           gtTGCCAGTCATTACCGTCACGATGACCGCCATCg                                           <  1:2711932/35‑1 (MQ=255)
                                                           gtTGCCAGTCATTACCGTCACGATGACCGCCATCg                                           <  1:1406633/35‑1 (MQ=255)
                                                           gtTGCCAGTCATTACCGTCACGATGACCGCCATCg                                           <  1:1473022/35‑1 (MQ=255)
                                                           gtTGCCAGTCATTACCGTCACGATGACCGCCATCg                                           <  1:3157229/35‑1 (MQ=255)
                                                              gccagtcaTTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATCGTCatgatgatgatg        >  1:3206339/1‑67 (MQ=255)
                                                              gccagtcaTTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATCGTCATGATGATGATGGt      >  1:2810971/1‑69 (MQ=255)
                                                              gccagtcaTTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATCGTCATGATGATGATGGt      >  1:2902949/1‑69 (MQ=255)
                                                              gccagtcaTTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATCGTCATGATGATGATGGt      >  1:29947/1‑69 (MQ=255)
                                                              gccagtcaTTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATCGTCATGATGATGATGGt      >  1:1954175/1‑69 (MQ=255)
                                                              gccagtcaTTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATCGTCATGATGATGATGGt      >  1:868478/1‑69 (MQ=255)
                                                              gccagtcaTTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATCGTCATGATGATGATGGt      >  1:3072812/1‑69 (MQ=255)
                                                              gccagtcaTTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATCGTCATGATGATGATGGt      >  1:3282131/1‑69 (MQ=255)
                                                              gccagtcaTTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATCGTCATGATGATGATGGt      >  1:1485409/1‑69 (MQ=255)
                                                              gccagtcaTTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATCGTCATGATGATGATGGt      >  1:2133680/1‑69 (MQ=255)
                                                              gccagtcaTTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATCGTCATGATGATGATGGt      >  1:1630653/1‑69 (MQ=255)
                                                              gccagtcaTTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATCGTCATGATGATGATGGt      >  1:1545452/1‑69 (MQ=255)
                                                              gccagtcaTTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATCGTCATGATGATGATGGt      >  1:1242278/1‑69 (MQ=255)
                                                              gccagtcaTTACCCTCACGATGACCGCCATCGTGGTGATGATAATCATCATCGTCATGATGATGATGGt      >  1:3100732/1‑69 (MQ=255)
                                                                   |                                                                   
ATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATG  >  minE/114945‑115079

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: