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.

AAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGA‑GGGG‑AATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGCGACT  >  W3110S.gb/3824463‑3824597
                                                                       |                                                                 
aaGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGa                                                                    <  1:2856689/71‑1 (MQ=255)
    gCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGa                                                                    >  1:869407/1‑67 (MQ=255)
     cAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑ggggg                                                              <  1:4417846/71‑2 (MQ=255)
     cAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑gggg                                                               >  1:4168668/1‑70 (MQ=255)
      aGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑ggggga                                                             <  1:3947732/71‑3 (MQ=255)
      aGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑ggggga                                                             <  1:3838589/71‑3 (MQ=255)
          gACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGG‑‑CTGGCCCGAGGGGg‑aatcga                                                        <  1:2021028/70‑1 (MQ=255)
                   aGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGG‑AATCGAATCCATTg                                                >  1:2467755/1‑71 (MQ=255)
                      cGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGG‑AATCGAATCCATTGTa                                              >  1:1603847/1‑70 (MQ=255)
                        aGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGG‑AATCGAATCCATTGTAGCg                                           <  1:3029150/71‑1 (MQ=255)
                              tCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGa                                                                    >  1:3974503/1‑41 (MQ=38)
                               cACCCATGAAGAGCTAGTCACCTTGTTGGCGCTGGCCCGAGGGGG‑AATCGAATCCATTGTAGCGACGCAg                                     >  1:129467/1‑70 (MQ=255)
                                 cccATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGG‑AATCGAATCCATTGTAGCGACTCAGAAg                                  <  1:616061/71‑1 (MQ=255)
                                    aTGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGg‑aatcg                                                         <  1:1848610/45‑1 (MQ=255)
                                    aTGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGg‑aatcg                                                         <  1:3536623/45‑1 (MQ=255)
                                    aTGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGg‑aatcg                                                         <  1:3221245/45‑1 (MQ=255)
                                    aTGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑GGGGCaatcgaa                                                       >  1:3466098/1‑47 (MQ=255)
                                       aaGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGG‑AATCGAATCCATTGTAGCGACGCAGAAggcggc                             <  1:2969131/70‑1 (MQ=255)
                                       aaGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGG‑AATCGAATCCATTGTAGCGACGCAGAAGgcggcg                            <  1:68340/71‑1 (MQ=255)
                                        agagCTACTCACCTTGTTGGCGCTGGCCCGAGGGGG‑AATCGAATCCAtt                                                 <  1:952001/49‑1 (MQ=255)
                                         gagCTACTCACCTTGTTGGCGCTGGCCCGAGGGGG‑AATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCt                          <  1:3006448/71‑1 (MQ=255)
                                         gagCTACTCACCTTGTTGGCGCTGGCCCGAGGGGG‑AATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCt                          <  1:2953942/71‑1 (MQ=255)
                                            cTACTCTCCTTGTTGGCGCTGGCCCGAGGGGG‑AATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTgg                        >  1:2357963/1‑70 (MQ=255)
                                            cTACTCACCTTGTTGGCGCTGGCCCGATGGGG‑AATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTgg                        >  1:2830441/1‑70 (MQ=255)
                                            cTACTCACCTTGTTGGCGCTGGCCCGAGGGGG‑AATCGAATCCATTGTAGCGACGCa                                      <  1:2204419/56‑1 (MQ=255)
                                               cTCACCTTGTTGGCGCTGGCCCGAGGGGG‑AATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCaaa                    >  1:357430/1‑71 (MQ=255)
                                               cTCACCTTGTTGGCGCTGGCCCGAGGGGG‑AATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCaa                     <  1:2515145/70‑1 (MQ=255)
                                                tCACCTTGTTGGCGCTGGCCCGAGGGGg‑aatcgaatc                                                     <  1:4597311/37‑1 (MQ=21)
                                                  accTTGTTGGCGCTGGCCCGAGGGGG‑AATCGAATCCATTGTAGCGa                                          >  1:1957595/3‑46 (MQ=255)
                                                   ccTTGTTGGCGCTGGCCCGAGGGGG‑AATCGAATCCATTg                                                <  1:4732197/38‑1 (MQ=25)
                                                   ccTTGTTGGCGCTGGCCCGAGGGGG‑AATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGa                >  1:3288559/2‑71 (MQ=255)
                                                     ttgttgGCGCTGGCCCGAGGGGG‑AATCGAATCCAt                                                  <  1:4597953/35‑1 (MQ=25)
                                                     ttgttgGCGCTGGCCCGAGGGGG‑AATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACt                  <  1:3414716/67‑1 (MQ=255)
                                                      tgttgGCGCTGGCCCGAGGGGG‑AATCGAATCCATTg                                                >  1:3753846/1‑36 (MQ=37)
                                                      tgttgGCGCTGGCCCGAGGGGG‑AATCGAATCCATTg                                                <  1:143866/36‑1 (MQ=37)
                                                      tgttgGCGCTGGCCCGAGGGGG‑AATCGAATCCATTGTAGCGACGCAg                                     <  1:391787/47‑1 (MQ=255)
                                                       gttgGCGCTGGCCCGAGGGGG‑AATCGAATCCATTg                                                >  1:2490355/1‑35 (MQ=25)
                                                       gttgGCGCTGGCCCGAGGGGG‑AATCGAATCCATTg                                                >  1:392927/1‑35 (MQ=25)
                                                       gttgGCGCTGGCCCGAGGGGG‑AATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCaaa                    >  1:3266515/1‑63 (MQ=255)
                                                        ttgGCGCTGGCCCGAGGGGG‑AATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGAttttt           <  1:1687347/71‑1 (MQ=255)
                                                          ggCGCTGGCCCGAGGGGG‑AATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGAtttt            >  1:1921205/1‑68 (MQ=255)
                                                          ggCGCTGGCCCGAGGGGG‑AATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTa          <  1:1147306/70‑1 (MQ=255)
                                                           gcgctGGCCCGAGGGGG‑AATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTaa         <  1:1493546/67‑1 (MQ=255)
                                                             gctGGCCCGAGGGGG‑AATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGc      <  1:1195669/70‑1 (MQ=255)
                                                              ctGGCCCGAGGGGG‑AATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGc      >  1:873362/1‑70 (MQ=255)
                                                                 gCCCGAGGGGG‑AATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGa                >  1:774983/1‑57 (MQ=255)
                                                                 gCCCGAGGGGG‑AATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGCGa    <  1:553185/69‑1 (MQ=255)
                                                                 gCCCGAGGGGG‑AATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGCGACt  <  1:3801411/71‑1 (MQ=255)
                                                                 gCCCGAGGGGG‑AATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGCGACt  <  1:3000336/71‑1 (MQ=255)
                                                                       |                                                                 
AAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGA‑GGGG‑AATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGCGACT  >  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: ‑