Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,876,422 A→G 41.1% G334G (GGT→GGC eno ← enolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,876,4220AG41.1% 24.2 / 52.3 73G334G (GGT→GGCenoenolase
Reads supporting (aligned to +/- strand):  ref base A (17/26);  new base G (15/15);  total (32/41)
Fisher's exact test for biased strand distribution p-value = 4.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGATTGCAGCCAGAGTTTCGGTCAGAGAACCGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACAGGTCGTCACCAACCAG  >  minE/1876354‑1876484
                                                                    |                                                              
ttGATTGCAGCCAGAGTTTCGGTCAGAGAACCGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATa                                                                >  1:1133409/1‑69 (MQ=255)
ttGATTGCAGCCAGAGTTTCGGTCAGAGAACCGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATa                                                                >  1:1133747/1‑69 (MQ=255)
ttGATTGCAGCCAGAGTTTCGGTCAGAGAACCGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATa                                                                >  1:781148/1‑69 (MQ=255)
ttGATTGCAGCCAGAGTTTCGGTCAGAGAACCGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATa                                                                >  1:35409/1‑69 (MQ=255)
ttGATTGCAGCCAGAGTTTCGGTCAGAGAACCGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATa                                                                >  1:1847757/1‑69 (MQ=255)
ttGATTGCAGCCAGAGTTTCGGTCAGAGAACCGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATa                                                                >  1:2595503/1‑69 (MQ=255)
ttGATTGCAGCCAGAGTTTCGGTCAGAGAACCGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATa                                                                >  1:3046917/1‑69 (MQ=255)
ttGATTGCAGCCAGAGTTTCGGTCAGAGAACCGATCTGGTTGAATTAGATCAGGATGGAGTTAGCGATa                                                                >  1:159495/1‑69 (MQ=255)
         gCCAGAGTTTCGGTCAGAGAACCTATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATACCTTTTTc                                                        <  1:1464272/68‑1 (MQ=255)
         gCCAGAGTTTCGGTCAGAGAACCGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATACCTTTTTc                                                        <  1:2858731/68‑1 (MQ=255)
         gCCAGAGTTTCGGTCAGAGAACCGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATACCTTTTTc                                                        <  1:2976125/68‑1 (MQ=255)
         gCCAGAGTTTCGGTCAGAGAACCGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATACCTTTTTc                                                        <  1:2063430/68‑1 (MQ=255)
         gCCAGAGTTTCGGTCAGAGAACCGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATACCTTTTTc                                                        <  1:3051903/68‑1 (MQ=255)
         gCCAGAGTTTCGGTCAGAGAACCGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATACCTTTTTc                                                        <  1:2778924/68‑1 (MQ=255)
         gCCAGAGTTTCGGTCAGAGAACCGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATACCTTTTTc                                                        <  1:2144334/68‑1 (MQ=255)
         gCCAGAGTTTCGGTCAGAGAACCGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATACCTTTTTc                                                        <  1:2370504/68‑1 (MQ=255)
         gCCAGAGTTTCGGTCAGAGAACCGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATACCTTTTTc                                                        <  1:1761117/68‑1 (MQ=255)
         gCCAGAGTTTCGGTCAGAGAACCGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATACCTTTTTc                                                        <  1:2264451/68‑1 (MQ=255)
         gCCAGAGTTTCGGTCAGAGAACCGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATACCTTTTTc                                                        <  1:3367726/68‑1 (MQ=255)
         gCCAGAGTTTCGGTCAGAGAACCGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATACCTTTTTc                                                        <  1:345821/68‑1 (MQ=255)
         gCCAGAGTTTCGGTCAGAGAACCGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATACCTTTTTc                                                        <  1:383187/68‑1 (MQ=255)
         gCCAGAGTTTCGGTCAGAGAACCGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATACCTTTTTc                                                        <  1:2396592/68‑1 (MQ=255)
         gCCAGAGTTTCGGTCAGAGAACCGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATACCTTTTTc                                                        <  1:1493667/68‑1 (MQ=255)
         gCCAGAGTTTCGGTCAGAGAACCGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATACCTTTTTc                                                        <  1:2470151/68‑1 (MQ=255)
         gCCAGAGTTTCGGTCAGAGAACCGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATACCTTTTTc                                                        <  1:1253278/68‑1 (MQ=255)
         gCCAGAGTTTCGGTCAGAGAACCGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATACCTTTTTc                                                        <  1:650924/68‑1 (MQ=255)
         gCCAGAGTTTCGGTCAGAGAACCGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATACCTTTTTc                                                        <  1:1315255/68‑1 (MQ=255)
              agTTTCGGTCAGAGAACCGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATACCTTTTTCGATAcc                                                  <  1:1445294/69‑1 (MQ=255)
              agTTTCGGTCAGAGAACCGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATACCTTTTTCGATAcc                                                  <  1:2893355/69‑1 (MQ=255)
                        agagAACCGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATACCTTTTTc                                                        <  1:1662391/53‑1 (MQ=255)
                                           aTTTGATCAGGATGGAGTTAGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACg                      <  1:2114592/68‑1 (MQ=255)
                                           aTTTGATCAGGATGGAGTTAGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACg                      <  1:2033288/68‑1 (MQ=255)
                                           aTTTGATCAGGATGGAGTTAGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACg                      <  1:1970374/68‑1 (MQ=255)
                                            tttGATCAGGATGGAGTTAGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACg                      <  1:671784/67‑1 (MQ=255)
                                              tGATCAGGATGGAGTTAGCGATGCCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTAc                       >  1:2226962/1‑64 (MQ=255)
                                              tGATCAGGATGGAGTTAGCGATGCCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAAc                   >  1:2858264/1‑68 (MQ=255)
                                              tGATCAGGATGGAGTTAGCGATGCCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAAc                   >  1:2653551/1‑68 (MQ=255)
                                              tGATCAGGATGGAGTTAGCGATGCCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACa                  >  1:819559/1‑69 (MQ=255)
                                              tGATCAGGATGGAGTTAGCGATGCCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACa                  >  1:2494292/1‑69 (MQ=255)
                                              tGATCAGGATGGAGTTAGCGATGCCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACa                  >  1:2711504/1‑69 (MQ=255)
                                              tGATCAGGATGGAGTTAGCGATGCCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACa                  >  1:3299307/1‑69 (MQ=255)
                                              tGATCAGGATGGAGTTAGCGATGCCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACa                  >  1:386242/1‑69 (MQ=255)
                                              tGATCAGGATGGAGTTAGCGATGCCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACa                  >  1:425656/1‑69 (MQ=255)
                                              tGATCAGGATGGAGTTAGCGATGCCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACa                  >  1:653834/1‑69 (MQ=255)
                                              tGATCAGGATGGAGTTAGCGATGCCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACa                  >  1:224368/1‑69 (MQ=255)
                                              tGATCAGGATGGAGTTAGCGATGCCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACa                  >  1:1325567/1‑69 (MQ=255)
                                              tGATCAGGATGGAGTTAGCGATGCCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACa                  >  1:151448/1‑69 (MQ=255)
                                              tGATCAGGATGGAGTTAGCGATGCCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACa                  >  1:1871569/1‑69 (MQ=255)
                                              tGATCAGGATGGAGTTAGCGATGCCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACa                  >  1:1966956/1‑69 (MQ=255)
                                                 tCAGGATGGAGTTAGCGATGCCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACAgg                <  1:2960419/68‑1 (MQ=255)
                                                 tCAGGATGGAGTTAGCGATGCCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACAgg                <  1:1565517/68‑1 (MQ=255)
                                                 tCAGGATGGAGTTAGCGATGCCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACAgg                <  1:1852878/68‑1 (MQ=255)
                                                 tCAGGATGGAGTTAGCGATGCCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACAgg                <  1:3256479/68‑1 (MQ=255)
                                                 tCAGGATGGAGTTAGCGATGCCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACAgg                <  1:3255996/68‑1 (MQ=255)
                                                 tCAGGATGGAGTTAGCGATGCCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACAgg                <  1:1524891/68‑1 (MQ=255)
                                                 tCAGGATGGAGTTAGCGATGCCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACAgg                <  1:3150808/68‑1 (MQ=255)
                                                 tCAGGATGGAGTTAGCGATGCCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACAgg                <  1:2936710/68‑1 (MQ=255)
                                                 tCAGGATGGAGTTAGCGATGCCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACAgg                <  1:1446038/68‑1 (MQ=255)
                                                 tCAGGATGGAGTTAGCGATGCCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACAgg                <  1:473068/68‑1 (MQ=255)
                                                 tCAGGATGGAGTTAGCGATGCCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACAgg                <  1:2853160/68‑1 (MQ=255)
                                                 tCAGGATGGAGTTAGCGATGCCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACAgg                <  1:1321661/68‑1 (MQ=255)
                                                 tCAGGATGGAGTTAGCGATGCCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACAgg                <  1:217687/68‑1 (MQ=255)
                                                 tCAGGATGGAGTTAGCGATGCCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACAgg                <  1:712265/68‑1 (MQ=255)
                                                  cAGGATGGAGTTAGCGATGCCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACAgg                <  1:1907698/67‑1 (MQ=255)
                                                              aGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGt                          >  1:2555928/1‑45 (MQ=255)
                                                              aGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGt                          >  1:720576/1‑45 (MQ=255)
                                                              aGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGt                          >  1:370081/1‑45 (MQ=255)
                                                              aGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGt                          >  1:124057/1‑45 (MQ=255)
                                                              aGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGt                          >  1:1549297/1‑45 (MQ=255)
                                                              aGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGt                          >  1:1572817/1‑45 (MQ=255)
                                                              aGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGt                          >  1:3159068/1‑45 (MQ=255)
                                                              aGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGt                          >  1:3097041/1‑45 (MQ=255)
                                                              aGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGt                          >  1:1988072/1‑45 (MQ=255)
                                                              aGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGt                          >  1:2001478/1‑45 (MQ=255)
                                                              aGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGt                          >  1:2064883/1‑45 (MQ=255)
                                                              aGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGt                          >  1:2474173/1‑45 (MQ=255)
                                                              aGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACAGGTCGTCACCAACCAg  >  1:2141385/1‑69 (MQ=255)
                                                              aGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACAGGTCGTCACCAACCAg  >  1:1273558/1‑69 (MQ=255)
                                                              aGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACAGGTCGTCACCAACCAg  >  1:2600063/1‑69 (MQ=255)
                                                              aGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACAGGTCGTCACCAACCAg  >  1:1107602/1‑69 (MQ=255)
                                                              aGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACAGCTCGTCACCAACCAg  >  1:238009/1‑69 (MQ=255)
                                                                    |                                                              
TTGATTGCAGCCAGAGTTTCGGTCAGAGAACCGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACAGGTCGTCACCAACCAG  >  minE/1876354‑1876484

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