Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,824,533:1 +G 100% pseudogene (666/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,5331.G97.7% 133.2 / ‑4.7 42pseudogene (666/687 nt)rphECK3633:JW3618:b3643; ribonuclease PH
Reads supporting (aligned to +/- strand):  ref base . (1/0);  major base G (14/26);  minor base T (1/0);  total (16/26)
Fisher's exact test for biased strand distribution p-value = 3.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-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.
Rejected as polymorphism: Polymorphic base substitution creates a homopolymer stretch.

AAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGAGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGCGACT  >  W3110S.gb/3824463‑3824597
                                                                       |                                                                 
aaGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGa                                                                    <  1:2856689/71‑1 (MQ=255)
    gCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGa                                                                    >  1:869407/1‑67 (MQ=255)
     cAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAggggg                                                              <  1:4417846/71‑2 (MQ=255)
     cAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAgggg                                                               >  1:4168668/1‑70 (MQ=255)
      aGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAggggga                                                             <  1:3947732/71‑3 (MQ=255)
      aGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAggggga                                                             <  1:3838589/71‑3 (MQ=255)
          gACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGG‑‑CTGGCCCGAGGGGgaatcga                                                        <  1:2021028/70‑1 (MQ=255)
                   aGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTg                                                >  1:2467755/1‑71 (MQ=255)
                      cGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTa                                              >  1:1603847/1‑70 (MQ=255)
                        aGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCg                                           <  1:3029150/71‑1 (MQ=255)
                              tCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGa                                                                    >  1:3974503/1‑41 (MQ=38)
                               cACCCATGAAGAGCTAGTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAg                                     >  1:129467/1‑70 (MQ=255)
                                 cccATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACTCAGAAg                                  <  1:616061/71‑1 (MQ=255)
                                    aTGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcg                                                         <  1:1848610/45‑1 (MQ=255)
                                    aTGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcg                                                         <  1:3536623/45‑1 (MQ=255)
                                    aTGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcg                                                         <  1:3221245/45‑1 (MQ=255)
                                    aTGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGCaatcgaa                                                       >  1:3466098/1‑47 (MQ=255)
                                       aaGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAggcggc                             <  1:2969131/70‑1 (MQ=255)
                                       aaGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGgcggcg                            <  1:68340/71‑1 (MQ=255)
                                        agagCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCAtt                                                 <  1:952001/49‑1 (MQ=255)
                                         gagCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCt                          <  1:3006448/71‑1 (MQ=255)
                                         gagCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCt                          <  1:2953942/71‑1 (MQ=255)
                                            cTACTCTCCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTgg                        >  1:2357963/1‑70 (MQ=255)
                                            cTACTCACCTTGTTGGCGCTGGCCCGATGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTgg                        >  1:2830441/1‑70 (MQ=255)
                                            cTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCa                                      <  1:2204419/56‑1 (MQ=255)
                                               cTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCaaa                    >  1:357430/1‑71 (MQ=255)
                                               cTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCaa                     <  1:2515145/70‑1 (MQ=255)
                                                tCACCTTGTTGGCGCTGGCCCGAGGGGgaatcgaatc                                                     <  1:4597311/37‑1 (MQ=21)
                                                  accTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGa                                          >  1:1957595/3‑46 (MQ=255)
                                                   ccTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTg                                                <  1:4732197/38‑1 (MQ=25)
                                                   ccTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGa                >  1:3288559/2‑71 (MQ=255)
                                                     ttgttgGCGCTGGCCCGAGGGGGAATCGAATCCAt                                                  <  1:4597953/35‑1 (MQ=25)
                                                     ttgttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACt                  <  1:3414716/67‑1 (MQ=255)
                                                      tgttgGCGCTGGCCCGAGGGGGAATCGAATCCATTg                                                >  1:3753846/1‑36 (MQ=37)
                                                      tgttgGCGCTGGCCCGAGGGGGAATCGAATCCATTg                                                <  1:143866/36‑1 (MQ=37)
                                                      tgttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAg                                     <  1:391787/47‑1 (MQ=255)
                                                       gttgGCGCTGGCCCGAGGGGGAATCGAATCCATTg                                                >  1:2490355/1‑35 (MQ=25)
                                                       gttgGCGCTGGCCCGAGGGGGAATCGAATCCATTg                                                >  1:392927/1‑35 (MQ=25)
                                                       gttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCaaa                    >  1:3266515/1‑63 (MQ=255)
                                                        ttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGAttttt           <  1:1687347/71‑1 (MQ=255)
                                                          ggCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGAtttt            >  1:1921205/1‑68 (MQ=255)
                                                          ggCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTa          <  1:1147306/70‑1 (MQ=255)
                                                           gcgctGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTaa         <  1:1493546/67‑1 (MQ=255)
                                                             gctGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGc      <  1:1195669/70‑1 (MQ=255)
                                                              ctGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGc      >  1:873362/1‑70 (MQ=255)
                                                                 gCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGa                >  1:774983/1‑57 (MQ=255)
                                                                 gCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGCGa    <  1:553185/69‑1 (MQ=255)
                                                                 gCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGCGACt  <  1:3801411/71‑1 (MQ=255)
                                                                 gCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGCGACt  <  1:3000336/71‑1 (MQ=255)
                                                                       |                                                                 
AAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGAGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGCGACT  >  W3110S.gb/3824463‑3824597

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