Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 163,666 C→T 25.9% E77K (GAA→AAA)  yaeI ← predicted phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE163,6660CT25.9% 45.5 / 14.5 27E77K (GAA→AAA) yaeIpredicted phosphatase
Reads supporting (aligned to +/- strand):  ref base C (5/15);  new base T (2/5);  total (7/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01

GCCGAAAAATTCAGCGACATATCAAATAATACGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACGAGAGTAATGGAGATC  >  minE/163596‑163729
                                                                      |                                                               
gCCGAAAAATTCAGCGACATATCAAATAATACGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTc                                                                 <  1:995066/71‑1 (MQ=255)
gCCGAAAAATTCAGCGACATATCAAATAATACGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTc                                                                 <  1:321299/71‑1 (MQ=255)
 ccGAAAAATTCAGCGACATATCAAATAATACGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTCt                                                                <  1:596721/71‑1 (MQ=255)
 ccGAAAAATTCAGCGACATATCAAATAATACGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTCt                                                                <  1:1006481/71‑1 (MQ=255)
      aaaTTCAGCGACATATCAAATAATACGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTc                                                                 >  1:764694/1‑65 (MQ=255)
      aaaTTCAGCGACATATCAAATAATACGTAATCGCCGCCCAGTAATATCAAATCGGGATTTTGTTCTATGcc                                                           <  1:1726775/71‑1 (MQ=255)
                     tCAAATAATACGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTTTATGCCAAGAGCAATc                                                 <  1:342461/66‑1 (MQ=255)
                     tCAAATAATACGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTCTATGCCAAGAGCAATCGCGt                                             >  1:1314379/1‑70 (MQ=255)
                          taataCGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTTTATGCCAAGAGCAATCGCGTCAGaaa                                       <  1:249385/71‑1 (MQ=255)
                            ataCGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAAt                                      <  1:968233/70‑1 (MQ=255)
                                    tCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTa                              >  1:283303/1‑70 (MQ=255)
                                            cAGTAATATCAAATCGGGCTTTTGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACaa                     <  1:531981/71‑1 (MQ=255)
                                            cAGTAATATCAAATCGGGCTTTTGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACaa                     <  1:43173/71‑1 (MQ=255)
                                            cAGTAATATCAAATCGGGCTTTTGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACaa                     <  1:1504962/71‑1 (MQ=255)
                                             aGTAATATCAAATCGGGCTTTTGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTaa                             <  1:963241/62‑1 (MQ=255)
                                                    tCAAATCGGGCTTTTGTTTTATGCCAATAGCAATCGCGTCAGAAATCa                                    <  1:1174133/48‑1 (MQ=255)
                                                    tCAAATCGGGCTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTa                              <  1:2173239/54‑1 (MQ=255)
                                                    tCAAATCGGGCTTTTGTTTTATGCCAAGAGCAATCGCGTAAGGAATCAGGCTTAAAGGaa                        >  1:2160551/1‑60 (MQ=255)
                                                    tCAAATCGGGCTTTTGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAgg                                  <  1:229342/50‑1 (MQ=255)
                                                    tCAAATCGGGCTTTTGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAgg                                  <  1:1514360/50‑1 (MQ=255)
                                                      aaaTCGGGCTTTTGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACGAGAGt            <  1:345412/70‑1 (MQ=255)
                                                      aaaTCGGGCTTTTGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACGAGAGt            <  1:1615543/70‑1 (MQ=255)
                                                      aaaTCGGGCTTTTGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACGAGAGt            <  1:97447/70‑1 (MQ=255)
                                                      aaaTCGGGCTTTTGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACGAGAGTa           <  1:1037005/71‑1 (MQ=255)
                                                          cGGGCTTTTGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACGAGAGTAATg        >  1:1780560/1‑70 (MQ=255)
                                                           gggCTTTTGTTCTATGCCAAGAGCAATCGCGTCAGaaa                                       <  1:65231/38‑1 (MQ=255)
                                                                tttGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACgaga              >  1:2014627/1‑58 (MQ=255)
                                                                tttGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACgaga              >  1:1889174/1‑58 (MQ=255)
                                                                 ttGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACGAGAGTAATGga      >  1:1922242/1‑65 (MQ=255)
                                                                 ttGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACGAGAGTAATGGAGATc  <  1:1485207/69‑1 (MQ=255)
                                                                      |                                                               
GCCGAAAAATTCAGCGACATATCAAATAATACGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACGAGAGTAATGGAGATC  >  minE/163596‑163729

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: