Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,824,280 T→C 100% pseudogene (413/687 nt) rph → ECK3633:JW3618:b3643; ribonuclease PH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,824,2800TC98.0% 161.1 / ‑4.8 48pseudogene (413/687 nt)rphECK3633:JW3618:b3643; ribonuclease PH
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (26/21);  minor base A (0/1);  total (26/22)
Fisher's exact test for biased strand distribution p-value = 4.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CGACGTGCTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGTAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCC  >  W3110S.gb/3824215‑3824347
                                                                 |                                                                   
cGACGTGCTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGa                                                                  <  1:4442520/69‑1 (MQ=255)
cGACGTGCTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATg                                                                <  1:2510459/71‑1 (MQ=255)
    gTGCTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCt                                                            <  1:3216489/71‑1 (MQ=255)
    gTGCTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCt                                                            <  1:749325/71‑1 (MQ=255)
    gTGCTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCt                                                            <  1:3401425/71‑1 (MQ=255)
     tGCTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTa                                                           >  1:2728365/1‑71 (MQ=255)
       cTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTa                                                           >  1:1478379/1‑69 (MQ=255)
          cAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAg                                                        >  1:2404812/1‑69 (MQ=255)
          cAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAg                                                        >  1:4620922/1‑69 (MQ=255)
             gCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCt                                                   <  1:1597469/71‑1 (MQ=255)
                 aTGGTGGCACGCGTACCGCGTCGATTAAGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCtggtg                                               >  1:3407632/1‑71 (MQ=255)
                   ggtggCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCtggtgg                                              >  1:1276608/1‑70 (MQ=255)
                    gtggCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGaa                                            <  1:969610/71‑1 (MQ=255)
                       gCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAAc                                         >  1:879685/1‑71 (MQ=255)
                            cgTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATg                                                                >  1:3724019/1‑43 (MQ=255)
                            cgTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATg                                                                >  1:1203100/1‑43 (MQ=255)
                                ccGCGTCGATTACGGTTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGc                                                    >  1:3508036/1‑51 (MQ=38)
                                ccGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCa                                     <  1:4284502/66‑1 (MQ=255)
                                ccGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCt                                 >  1:3503667/1‑70 (MQ=255)
                                ccGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCt                                 >  1:3385424/1‑70 (MQ=255)
                                 cgcgTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTg                                >  1:327451/1‑70 (MQ=255)
                                  gcgTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGaa                              <  1:130872/71‑1 (MQ=255)
                                  gcgTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGaa                              <  1:1006153/71‑1 (MQ=255)
                                    gTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGaaaa                            >  1:1060268/1‑71 (MQ=255)
                                    gTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGaa                              <  1:3282792/69‑1 (MQ=255)
                                     tCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAAc                           <  1:158723/71‑1 (MQ=255)
                                     tCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAAc                           <  1:3737806/71‑1 (MQ=255)
                                     tCGATTACGGGTGCCTGCGTGGCGCTGGAAGATgc                                                               <  1:3961640/35‑1 (MQ=255)
                                       gATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGaaa                                           >  1:375083/1‑53 (MQ=255)
                                       gATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAAcc                          >  1:280055/1‑70 (MQ=255)
                                          tACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAg                                                        >  1:221045/1‑37 (MQ=255)
                                          tACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCaacc                      >  1:67198/1‑69 (MQ=255)
                                                tGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa                <  1:1189798/71‑1 (MQ=255)
                                                tGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa                <  1:3355073/71‑1 (MQ=255)
                                                tGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa                <  1:2228531/71‑1 (MQ=255)
                                                    tGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGa                               <  1:3560108/52‑1 (MQ=255)
                                                     gCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGaa                                            <  1:177250/38‑1 (MQ=255)
                                                      cGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAAc                                         >  1:2756446/1‑40 (MQ=255)
                                                      cGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAAc                           >  1:519569/1‑54 (MQ=255)
                                                      cGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGa                   >  1:1106436/1‑62 (MQ=255)
                                                      cGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGa          >  1:1422445/1‑71 (MQ=255)
                                                      cGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGa          >  1:238631/1‑71 (MQ=255)
                                                       gTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGaaaa                                          >  1:1073106/1‑38 (MQ=255)
                                                        tGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGAt         <  1:1365718/70‑1 (MQ=255)
                                                          gcgcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGa                   >  1:2023847/1‑58 (MQ=255)
                                                           cgcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGt      <  1:3602334/70‑1 (MQ=255)
                                                             cTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTg                                <  1:2322653/42‑1 (MQ=255)
                                                              tGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGcc  >  1:2702474/1‑71 (MQ=255)
                                                                 |                                                                   
CGACGTGCTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGTAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCC  >  W3110S.gb/3824215‑3824347

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: