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.G94.7% 178.0 / 1.5 57pseudogene (666/687 nt)rphECK3633:JW3618:b3643; ribonuclease PH
Reads supporting (aligned to +/- strand):  ref base . (0/3);  new base G (23/31);  total (23/34)
Fisher's exact test for biased strand distribution p-value = 2.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.83e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

AAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCG‑A‑G‑G‑GGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGCGACT  >  W3110S.gb/3824463‑3824597
                                                                        |                                                                  
aaGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCG‑a                                                                     <  1:842064/71‑1 (MQ=255)
aaGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCG‑a                                                                     <  1:3172517/71‑1 (MQ=255)
 aGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCG‑A‑g                                                                   <  1:1550367/71‑1 (MQ=255)
 aGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCG‑A‑g                                                                   <  1:877079/71‑1 (MQ=255)
     cAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCG‑A‑g‑g‑gg                                                              >  1:3408409/1‑70 (MQ=255)
     cAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCG‑A‑g‑g‑gg                                                              >  1:2333825/1‑70 (MQ=255)
     cAGGGGACGGCAGAAGGCGAACCGTTCACCCACGAAGAGTTACTCACCTTGTTGGCGCTGGCCCG‑A‑g                                                                   >  1:2178739/1‑67 (MQ=255)
     cAGGGGACGGCAGAAGGCGAACCGTTCACCCACGAAGAGTTACTCACCTTGTTGGCGCTGGCCCG‑A‑g                                                                   >  1:1758876/1‑67 (MQ=255)
      aGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCG‑A‑g‑g‑ggga                                                            <  1:510270/71‑3 (MQ=255)
           aCGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCG‑AGG‑G‑Ggaatcga                                                        <  1:3463240/71‑1 (MQ=255)
           aCGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCG‑AGG‑G‑Ggaatcga                                                        <  1:2082046/71‑1 (MQ=255)
           aCGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCG‑AGG‑G‑Ggaatcga                                                        <  1:3158464/71‑1 (MQ=255)
              gCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCG‑AGG‑G‑Ggaatcgaatc                                                     <  1:2005220/71‑1 (MQ=255)
               cAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCG‑AGG‑G‑GGAATCGAATcc                                                    <  1:3321471/71‑1 (MQ=255)
                   aGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCG‑AGG‑G‑GGAATCGAATCCATTg                                                >  1:2754826/1‑71 (MQ=255)
                       gAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCG‑AGG‑G‑GGAATCGAATCCATTGTAg                                             <  1:2643313/70‑1 (MQ=255)
                        aGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCG‑AGG‑G‑GGAATCGAATCCATTGTAGCg                                           <  1:3456541/71‑1 (MQ=255)
                        aGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCG‑AGG‑G‑GGAATCGAATCCATTGTAGCg                                           <  1:3634911/71‑1 (MQ=255)
                         gCCGTTCACCCATGAAGAGCTACTCACATTGTTGGCGCTGGCCCG‑AGG‑G‑GGAATCGAATCCATTGTAGCGa                                          <  1:1607274/71‑1 (MQ=255)
                          ccGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCG‑AGG‑G‑GGAATCGAATCCATTGTAGCGAc                                         <  1:388967/71‑1 (MQ=255)
                             ttCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCG‑AGG‑G‑Ggaatcgaatc                                                     <  1:1171780/56‑1 (MQ=255)
                             ttCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCG‑A‑g                                                                   <  1:662987/43‑1 (MQ=38)
                               cACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCG‑AGG‑G‑GGAATCGAATCCATTGTAGCGa                                          >  1:1817844/1‑65 (MQ=255)
                               cACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCG‑AGG‑G‑GGAATCGAATCCATTGTAGCGACGCAg                                     <  1:1281680/70‑1 (MQ=255)
                               cACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCG‑AGG‑G‑GGAATCGAATCCATTGTAGCGACGCAg                                     <  1:398828/70‑1 (MQ=255)
                                aCCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCG‑AGG‑G‑GGAATCGAATCCATTGTAGCGACGCAg                                     <  1:1725043/69‑1 (MQ=255)
                                aCCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCG‑AGG‑G‑GGAATCGAATCCATTGTAGCGACGCAGaa                                   <  1:1908664/71‑1 (MQ=255)
                                aCCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCG‑AGG‑G‑GGAATCGAATCCATTGTAGCGACGCAGaa                                   <  1:2670935/71‑1 (MQ=255)
                                aCCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCG‑AGG‑G‑GGAATCGAATCCATTGTAGCGACGCAGa                                    <  1:1146178/70‑1 (MQ=255)
                                aCCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCG‑A‑G‑GCGGAATCGAATCCATTGTAGCGACGCAGa                                    <  1:3100012/70‑1 (MQ=255)
                                  ccATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCG‑AGG‑G‑GGAATCGAATCCATTGTAGCGACGCAGaa                                   >  1:66574/1‑69 (MQ=255)
                                   cATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCG‑AGG‑G‑Ggaatcgaa                                                       >  1:3761213/1‑48 (MQ=255)
                                    aTGAAGAGCTACTCACCTTGTTGGCGCTGGCCCG‑AGG‑G‑GGAATCGAATCCATTGTAGCGACGCAGAAgggc                               >  1:872260/1‑69 (MQ=255)
                                    aTGAAGAGCTACTCACCTTGTTGGCGCTGGCCCG‑AGG‑G‑GGAATCGAATCCATTGTAGCGACGCAGAAggcg                               >  1:2043925/1‑71 (MQ=255)
                                     tGAAGAGCTACTCACCTTGTTGGCGCTGGCCCG‑A‑g‑g‑ggga                                                            <  1:1406776/40‑3 (MQ=37)
                                       aaGAGCTACTCACCTTGTTGGCGCTGGCCCG‑AGG‑G‑GGAATCGAATCCATTGTAGCGACGCAGaa                                   >  1:896045/1‑64 (MQ=255)
                                       aaGAGCTACTCACCTTGTTGGCGCTGGCCCG‑AGG‑G‑GGAATCGAATCCATTGTAGCGACGCAGaa                                   >  1:1317158/1‑64 (MQ=255)
                                       aaGAGCTACTCACCTTGTTGGCGCTGGCCCG‑AGG‑G‑GGAATCGAATCCATTGTAGCGACGCAGaa                                   >  1:3399611/1‑64 (MQ=255)
                                       aaGAGCTACTCACCTTGTTGGCGCTGGCCCG‑AGG‑G‑GGAATCGAATCCATTGTAGCGACGCAGaa                                   >  1:3490247/1‑64 (MQ=255)
                                       aaGAGCTACTCACCTTGTTGGCGCTGGCCCG‑AGG‑G‑GGAATCGAATCCATTGTAGCGACGCAGaa                                   >  1:3548122/1‑64 (MQ=255)
                                       aaGAGCTACTCACCTTGTTGGCGCTGGCCCG‑AGG‑G‑GGAATCGAATCCATTGTAGCGACGCAGAAGgcggcg                            >  1:2730550/1‑71 (MQ=255)
                                            cTACTCACCTTGTTGGCGCTGGCCCG‑AGG‑G‑GGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTgg                        >  1:3361413/1‑70 (MQ=255)
                                               cTCACCTTGTTGGCGCTGGCCCG‑AGG‑G‑GGAATCGAATCCa                                                   <  1:682923/40‑1 (MQ=25)
                                               cTCACCTTGTTGGCGCTGGCCCG‑AGG‑G‑GGAATCGAATCCATTg                                                <  1:68162/43‑1 (MQ=25)
                                               cTCACCTTGTTGGCGCTGGCCCG‑AGG‑G‑GGAATCGAATCCATTg                                                <  1:1821237/43‑1 (MQ=25)
                                               cTCACCTTGTTGGCGCTGGCCCG‑AGG‑G‑GGAATCGAATCCATTg                                                <  1:78031/43‑1 (MQ=25)
                                                tCACCTTGTTGGCGCTGGCCCGAA‑G‑G‑GGAATCGAATCCATTGt                                               <  1:1408013/43‑1 (MQ=25)
                                                  accTTGTTGGCGCTGGCCCG‑AGG‑G‑GGAATCGAATCCATTGTAGCGACGCAGAAGGCGgcgc                           >  1:1826232/3‑61 (MQ=255)
                                                  accTTGTTGGCGCTGGCCCG‑AGG‑G‑GGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTgg                        >  1:989526/3‑64 (MQ=255)
                                                  accTTGTTGGCGCTGGCCCG‑AGG‑G‑GGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTgg                        >  1:2646424/3‑64 (MQ=255)
                                                  accTTGTTGGCGCTGGCCCG‑AGG‑G‑GGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACt                  >  1:2518020/3‑70 (MQ=255)
                                                  accTTGTTGGCGCTGGCCCG‑AGG‑G‑GGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACt                  >  1:1153426/3‑70 (MQ=255)
                                                  accTTGTTGGCGCTGGCCCG‑AGG‑G‑GGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACt                  >  1:3221753/3‑70 (MQ=255)
                                                    cTTGTTGGCGCTGGCCCG‑AGG‑G‑GGAATCGAATCCATTGTAGc                                            <  1:3347901/42‑1 (MQ=37)
                                                    cTTGTTGGCGCTGGCCCG‑AGG‑G‑GGAATCGAATCCATTGTAGc                                            <  1:7739/42‑1 (MQ=37)
                                                     ttgttgGCGCTGGCCCG‑AGG‑G‑GGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGAtt              >  1:3632924/1‑71 (MQ=255)
                                                     ttgttgGCGCTGGCCCG‑AGG‑G‑GGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGAtt              >  1:717324/1‑71 (MQ=255)
                                                     ttgttgGCGCTGGCCCG‑AGG‑G‑GGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGAtt              >  1:1641703/1‑71 (MQ=255)
                                                        ttgGCGCTGGCCCG‑AGG‑G‑GGAATCGAATCCATTGTAGc                                            <  1:84170/38‑1 (MQ=37)
                                                          ggCGCTGGCCCG‑AGG‑G‑GGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGAtt              >  1:1526648/1‑66 (MQ=255)
                                                          ggCGCTGGCCCG‑AGG‑G‑GGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTa          <  1:2938920/70‑1 (MQ=255)
                                                          ggCGCTGGCCCG‑A‑GTG‑GGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTa          <  1:3075069/70‑1 (MQ=255)
                                                           gcgctGGCCCG‑AGG‑G‑GGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTa          <  1:569/66‑1 (MQ=255)
                                                                 gCCCG‑AGG‑G‑GGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGCGACt  <  1:450894/71‑1 (MQ=255)
                                                                 gCCCG‑AGG‑G‑GGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGCGACt  <  1:2070924/71‑1 (MQ=255)
                                                                 gCCCG‑AGG‑G‑GGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGCGACt  <  1:1908807/71‑1 (MQ=255)
                                                                 gCCCG‑AGG‑G‑GGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGCGACt  <  1:1677434/71‑1 (MQ=255)
                                                                  cccG‑AGG‑G‑GGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGCGACt  <  1:34546/70‑1 (MQ=255)
                                                                        |                                                                  
AAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCG‑A‑G‑G‑GGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGCGACT  >  W3110S.gb/3824463‑3824597

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