Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,009,916 A→T 43.5% pseudogene (111/411 nt) yghX ← ECK2993:JW5926:b2999; predicted hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,009,9160AT43.5% 21.6 / 64.0 62pseudogene (111/411 nt)yghXECK2993:JW5926:b2999; predicted hydrolase
Reads supporting (aligned to +/- strand):  ref base A (16/19);  new base T (19/8);  total (35/27)
Fisher's exact test for biased strand distribution p-value = 7.19e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.66e-01

CCACATCGGCAGTGGGTGCCTGACGACCATAAAACGGCACCGCGCAGGCCAGTTCCGGATACGCGACAGCCGCCGCGTTCGATACGCCACCGCCATAGCAAAATCCGGTAATACCCACTTTGCCTGTCGCTTGCG  >  minE/2009849‑2009983
                                                                   |                                                                   
ccACATCGGCAGTGGGTGCCTGACGACCATAAAACGGCACCGCGCAGGCCAGTTCCGGATACGCGACa                                                                     >  1:1160555/1‑68 (MQ=255)
ccACATCGGCAGTGGGTGCCTGACGACCATAAAACGGCACCGCGCAGGCCAGTTCCGGATACGCGACAg                                                                    >  1:1036146/1‑69 (MQ=255)
ccACATCGGCAGTGGGTGCCTGACGACCATAAAACGGCACCGCGCAGGCCAGTTCCGGATACGCGACAg                                                                    >  1:1574477/1‑69 (MQ=255)
ccACATCGGCAGTGGGTGCCTGACGACCATAAAACGGCACCGCGCAGGCCAGTTCCGGATACGCGACAg                                                                    >  1:1015727/1‑69 (MQ=255)
ccACATCGGCAGTGGGTGCCTGACGACCATAAAACGGCACCGCGCAGGCCAGTTCCGGATACGCGACAg                                                                    >  1:3390948/1‑69 (MQ=255)
     tCGGCAGTGGGTGCCTGACGACCATAAAACGGCACCGCGCAGGCCAGTTCCGGATACGCGACAgccgcc                                                               >  1:1364156/1‑69 (MQ=255)
     tCGGCAGTGGGTGCCTGACGACCATAAAACGGCACCGCGCAGGCCAGTTCCGGATACGCGACAgccgcc                                                               >  1:2114490/1‑69 (MQ=255)
     tCGGCAGTGGGTGCCTGACGACCATAAAACGGCACCGCGCAGGCCAGTTCCGGATACGCGACAgccgcc                                                               >  1:378099/1‑69 (MQ=255)
     tCGGCAGTGGGTGCCTGACGACCATAAAACGGCACCGCGCAGGCCAGTTCCGGATACGCGACAgccgcc                                                               >  1:895654/1‑69 (MQ=255)
             gggTGCCTGACGACCATAAAACGGCACCGCGCAGGCCAGTTCCGGATACGCGACAGCCGCCGCGTTc                                                         <  1:953306/67‑1 (MQ=255)
             gggTGCCTGACGACCATAAAACGGCACCGCGCAGGCCAGTTCCGGATACGCGACAGCCGCCGCGTTc                                                         <  1:1227533/67‑1 (MQ=255)
             gggTGCCTGACGACCATAAAACGGCACCGCGCAGGCCAGTTCCGGATACGCGACAGCCGCCGCGTTc                                                         <  1:2702948/67‑1 (MQ=255)
             gggTGCCTGACGACCATAAAACGGCACCGCGCAGGCCAGTTCCGGATACGCGACAGCCGCCGCGTTc                                                         <  1:476323/67‑1 (MQ=255)
             gggTGCCTGACGACCATAAAACGGCACCGCGCAGGCCAGTTCCGGATACGCGACAGCCGCCGCGTTCGa                                                       <  1:2669978/69‑1 (MQ=255)
             gggTGCCTGACGACCATAAAACGGCACCGCGCAGGCCAGTTCCGGATACGCGACAGCCGCCGCGTTCGa                                                       <  1:3315177/69‑1 (MQ=255)
             gggTGCCTGACGACCATAAAACGGCACCGCGCAGGCCAGTTCCGGATACGCGACAGCCGCCGCGTTCGa                                                       <  1:1497989/69‑1 (MQ=255)
             gggTGCCTGACGACCATAAAACGGCACCGCGCAGGCCAGTTCCGGATACGCGACAGCCGCCGCGTTCGa                                                       <  1:1947348/69‑1 (MQ=255)
             gggTGCCTGACGACCATAAAACGGCACCGCGCAGGCCAGTTCCGGATACGCGACAGCCGCCGCGTTCGa                                                       <  1:2252645/69‑1 (MQ=255)
              ggTGCCTGACGACCATAAAACGGCACCGCGCAGGCCAGTTCCGGATACGCGACAGCCGCCGCGTTCGAt                                                      >  1:3380699/1‑69 (MQ=255)
              ggTGCCTGACGACCATAAAACGGCACCGCGCAGGCCAGTTCCGGATACGCGACAGCCGCCGCGTTCGAt                                                      >  1:1585629/1‑69 (MQ=255)
                      acgacCTTAAAACGGCACCGCGCAGGCCAGTTCCGGATACGCGACAGCCGCCGCGTTc                                                         <  1:2340077/58‑1 (MQ=255)
                          ccATAAAACGGCACCGCGCAGGCCAGTTCCGGATACGCGACTGCTGCCGCGTTCGATACGCCCCCGCCa                                          <  1:2579099/69‑1 (MQ=255)
                          ccATAAAACGGCACCGCGCAGGCCAGTTCCGGATACGCGACTGCTGCCGCGTTCGATACGCCCCCGCCa                                          <  1:2912647/69‑1 (MQ=255)
                          ccATAAAACGGCACCGCGCAGGCCAGTTCCGGATACGCGACTGCTGCCGCGTTCGATACGCCCCCGCCa                                          <  1:2949205/69‑1 (MQ=255)
                          ccATAAAACGGCACCGCGCAGGCCAGTTCCGGATACGCGACTGCTGCCGCGTTCGATACGCCCCCGCCa                                          <  1:1780203/69‑1 (MQ=255)
                          ccATAAAACGGCACCGCGCAGGCCAGTTCCGGATACGCGACTGCTGCCGCGTTCGATACGCCCCCGCCa                                          <  1:1381936/69‑1 (MQ=255)
                          ccATAAAACGGCACCGCGCAGGCCAGTTCCGGATACGCGACTGCTGCCGCGTTCGATACGCCCCCGCCa                                          <  1:401526/69‑1 (MQ=255)
                           cATAAAACGGCACCGCGCAGGCCAGTTCCGGATACGCGACTGCTGCCGCGTTCGATACGCCCCCGCCa                                          <  1:683938/68‑1 (MQ=255)
                               aaaCGGCACCGCGCAGGCCAGTTCCGGATACGCGACAGCCGCCGCGTTCGATACGCCACCGCCATAGc                                      >  1:3383086/1‑68 (MQ=255)
                               aaaCGGCACCGCGCAGGCCAGTTCCGGATACGCGACAGCCGCCGCGTTCGATACGCCACCGCCATAGc                                      >  1:1023712/1‑68 (MQ=255)
                               aaaCGGCACCGCGCAGGCCAGTTCCGGATACGCGACAGCCGCCGCGTTCGATACGCCACCGCCATAGc                                      >  1:942607/1‑68 (MQ=255)
                                 aCGGCACCGCGCAGGCCAGTTCCGGATACGCGACAGCCGCCGCGTTCGATACGCCACCGCCATAGCaa                                    <  1:1437565/68‑1 (MQ=255)
                                        cgcgCAGGCCAGTTCCGGATACGCGACTGCTGCCGCGTTCGATACGCCCCCGCCATAGGAAAATcc                               >  1:2250304/1‑66 (MQ=255)
                                        cgcgCAGGCCAGTTCCGGATACGCGACTGCTGCCGCGTTCGATACGCCCCCGCCATAGCAAAATcc                               >  1:3010789/1‑66 (MQ=255)
                                        cgcgCAGGCCAGTTCCGGATACGCGACTGCTGCCGCGTTCGATACGCCCCCGCCATAGCAAAATcc                               >  1:405478/1‑66 (MQ=255)
                                        cgcgCAGGCCAGTTCCGGATACGCGACTGCTGCCGCGTTCGATACGCCCCCGCCATAGCAAAATcc                               >  1:2589361/1‑66 (MQ=255)
                                        cgcgCAGGCCAGTTCCGGATACGCGACTGCTGCCGCGTTCGATACGCCCCCGCCATAGCAAAATcc                               >  1:206417/1‑66 (MQ=255)
                                        cgcgCAGGCCAGTTCCGGATACGCGACTGCTGCCGCGTTCGATACGCCCCCGCCATAGCAAAATcc                               >  1:2156744/1‑66 (MQ=255)
                                        cgcgCAGGCCAGTTCCGGATACGCGACTGCTGCCGCGTTCGATACGCCCCCGCCATAGCAAAATcc                               >  1:2272008/1‑66 (MQ=255)
                                        cgcgCAGGCCAGTTCCGGATACGCGACTGCTGCCGCGTTCGATACGCCCCCGCCATAGCAAAATcc                               >  1:911269/1‑66 (MQ=255)
                                        cgcgCAGGCCAGTTCCGGATACGCGACAGCCGCCGCGTTCGATACGCCACCGCCATAGCAAAATCCGGt                            <  1:2935195/69‑1 (MQ=255)
                                        cgcgCAGGCCAGTTCCGGATACGCGACAGCCGCCGCGTTCGATACGCCACCGCCATAGCAAAATCCGGt                            <  1:2097641/69‑1 (MQ=255)
                                        cgcgCAGGCCAGTTCCGGATACGCGACAGCCGCCGCGTTCGATACGCCACCGCCATAGCAAAATCCGGt                            <  1:3335104/69‑1 (MQ=255)
                                        cgcgCAGGCCAGTTCCGGATACGCGACAGCCGCCGCGTTCGATACGCCACCGCCATAGCAAAATCCGGt                            <  1:1049247/69‑1 (MQ=255)
                                         gcgcAGGCCAGTTCCGGATACGCGACAGCCGCCGCGTTCGATACGCCACCGCCATAGCAAAATCCGGt                            >  1:369885/1‑68 (MQ=255)
                                              ggCCAGTTCCGGATACGCGACAGCCGCCGCGTTCGa                                                       <  1:485490/36‑1 (MQ=255)
                                                    ttCCGGATACGCGACTGCTGCCGCGTTCGATACGCCCCCGCCa                                          <  1:451315/43‑1 (MQ=255)
                                                      ccGGATACGCGACAGCCGCCGCGTTCGATACGCCACCGCCATAGCAAAATCCGGTAATACCCACTTTGc              >  1:1257310/1‑69 (MQ=255)
                                                      ccGGATACGCGACAGCCGCCGCGTTCGATACGCCACCGCCATAGCAAAATCCGGTAATACCCACTTTGc              >  1:1111753/1‑69 (MQ=255)
                                                              gcgACTGCTGCCGCGTTCGATACGCCCCCGCCATATCAAAATcc                               >  1:1326441/1‑44 (MQ=25)
                                                              gcgACTGCTGCCGCGTTCGATACGCCCCCGCCATAGCAAAAt                                 >  1:9762/1‑42 (MQ=37)
                                                              gcgACTGCTGCCGCGTTCGATACGCCCCCGCCATAGCAAAATcc                               >  1:298104/1‑44 (MQ=38)
                                                              gcgACTGCTGCCGCGTTCGATACGCCCCCGCCATAGCAAAATcc                               >  1:536393/1‑44 (MQ=38)
                                                              gcgACTGCTGCCGCGTTCGATACGCCCCCGCCATAGCAAAATcc                               >  1:2839332/1‑44 (MQ=38)
                                                              gcgACTGCTGCCGCGTTCGATACGCCCCCGCCATAGCAAAATcc                               >  1:781017/1‑44 (MQ=38)
                                                              gcgACTGCTGCCGCGTTCGATACGCCCCCGCCATAGCAAAATcc                               >  1:2620553/1‑44 (MQ=38)
                                                              gcgACTGCTGCCGCGTTCGATACGCCCCCGCCATAGCAAAATcc                               >  1:2386801/1‑44 (MQ=38)
                                                              gcgACTGCTGCCGCGTTCGATACGCCCCCGCCATAGCAAAATcc                               >  1:2257189/1‑44 (MQ=38)
                                                              gcgACTGCTGCCGCGTTCGATACGCCCCCGCCATAGCAAAATcc                               >  1:1166805/1‑44 (MQ=38)
                                                              gcgACTGCTGCCGCGTTCGATACGCCCCCGCCATAGCAAAATcc                               >  1:1079669/1‑44 (MQ=38)
                                                                  cAGCCGCCGCGTTCGATACGCCACCGCCATAGCAAAATCCGGTAATACCCACTTTGCCTGTCGCTTGCg  <  1:354797/69‑1 (MQ=255)
                                                                  cAGCCGCCGCGTTCGATACGCCACCGCCATAGCAAAATCCGGTAATACCCACTTTGCCTGTCGCTTGCg  <  1:2878955/69‑1 (MQ=255)
                                                                  cAGCCGCCGCGTTCGATACGCCACCGCCATAGCAAAATCCGGTAATACCCACTTTGCCTGTCGCTTGCg  <  1:653060/69‑1 (MQ=255)
                                                                   |                                                                   
CCACATCGGCAGTGGGTGCCTGACGACCATAAAACGGCACCGCGCAGGCCAGTTCCGGATACGCGACAGCCGCCGCGTTCGATACGCCACCGCCATAGCAAAATCCGGTAATACCCACTTTGCCTGTCGCTTGCG  >  minE/2009849‑2009983

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: