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.

AAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGAGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGCGACT  >  W3110S.gb/3824463‑3824597
                                                                        |                                                                  
aaGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGa                                                                     <  1:842064/71‑1 (MQ=255)
aaGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGa                                                                     <  1:3172517/71‑1 (MQ=255)
 aGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAg                                                                   <  1:1550367/71‑1 (MQ=255)
 aGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAg                                                                   <  1:877079/71‑1 (MQ=255)
     cAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAgggg                                                              >  1:3408409/1‑70 (MQ=255)
     cAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAgggg                                                              >  1:2333825/1‑70 (MQ=255)
     cAGGGGACGGCAGAAGGCGAACCGTTCACCCACGAAGAGTTACTCACCTTGTTGGCGCTGGCCCGAg                                                                   >  1:2178739/1‑67 (MQ=255)
     cAGGGGACGGCAGAAGGCGAACCGTTCACCCACGAAGAGTTACTCACCTTGTTGGCGCTGGCCCGAg                                                                   >  1:1758876/1‑67 (MQ=255)
      aGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAggggga                                                            <  1:510270/71‑3 (MQ=255)
           aCGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcga                                                        <  1:3463240/71‑1 (MQ=255)
           aCGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcga                                                        <  1:2082046/71‑1 (MQ=255)
           aCGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcga                                                        <  1:3158464/71‑1 (MQ=255)
              gCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcgaatc                                                     <  1:2005220/71‑1 (MQ=255)
               cAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATcc                                                    <  1:3321471/71‑1 (MQ=255)
                   aGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTg                                                >  1:2754826/1‑71 (MQ=255)
                       gAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAg                                             <  1:2643313/70‑1 (MQ=255)
                        aGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCg                                           <  1:3456541/71‑1 (MQ=255)
                        aGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCg                                           <  1:3634911/71‑1 (MQ=255)
                         gCCGTTCACCCATGAAGAGCTACTCACATTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGa                                          <  1:1607274/71‑1 (MQ=255)
                          ccGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGAc                                         <  1:388967/71‑1 (MQ=255)
                             ttCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcgaatc                                                     <  1:1171780/56‑1 (MQ=255)
                             ttCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAg                                                                   <  1:662987/43‑1 (MQ=38)
                               cACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGa                                          >  1:1817844/1‑65 (MQ=255)
                               cACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAg                                     <  1:1281680/70‑1 (MQ=255)
                               cACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAg                                     <  1:398828/70‑1 (MQ=255)
                                aCCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAg                                     <  1:1725043/69‑1 (MQ=255)
                                aCCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGaa                                   <  1:1908664/71‑1 (MQ=255)
                                aCCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGaa                                   <  1:2670935/71‑1 (MQ=255)
                                aCCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGa                                    <  1:1146178/70‑1 (MQ=255)
                                aCCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGCGGAATCGAATCCATTGTAGCGACGCAGa                                    <  1:3100012/70‑1 (MQ=255)
                                  ccATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGaa                                   >  1:66574/1‑69 (MQ=255)
                                   cATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcgaa                                                       >  1:3761213/1‑48 (MQ=255)
                                    aTGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAgggc                               >  1:872260/1‑69 (MQ=255)
                                    aTGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAggcg                               >  1:2043925/1‑71 (MQ=255)
                                     tGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAggggga                                                            <  1:1406776/40‑3 (MQ=37)
                                       aaGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGaa                                   >  1:896045/1‑64 (MQ=255)
                                       aaGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGaa                                   >  1:1317158/1‑64 (MQ=255)
                                       aaGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGaa                                   >  1:3399611/1‑64 (MQ=255)
                                       aaGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGaa                                   >  1:3490247/1‑64 (MQ=255)
                                       aaGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGaa                                   >  1:3548122/1‑64 (MQ=255)
                                       aaGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGgcggcg                            >  1:2730550/1‑71 (MQ=255)
                                            cTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTgg                        >  1:3361413/1‑70 (MQ=255)
                                               cTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCa                                                   <  1:682923/40‑1 (MQ=25)
                                               cTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTg                                                <  1:68162/43‑1 (MQ=25)
                                               cTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTg                                                <  1:1821237/43‑1 (MQ=25)
                                               cTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTg                                                <  1:78031/43‑1 (MQ=25)
                                                tCACCTTGTTGGCGCTGGCCCGAAGGGGAATCGAATCCATTGt                                               <  1:1408013/43‑1 (MQ=25)
                                                  accTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGgcgc                           >  1:1826232/3‑61 (MQ=255)
                                                  accTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTgg                        >  1:989526/3‑64 (MQ=255)
                                                  accTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTgg                        >  1:2646424/3‑64 (MQ=255)
                                                  accTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACt                  >  1:2518020/3‑70 (MQ=255)
                                                  accTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACt                  >  1:1153426/3‑70 (MQ=255)
                                                  accTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACt                  >  1:3221753/3‑70 (MQ=255)
                                                    cTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGc                                            <  1:3347901/42‑1 (MQ=37)
                                                    cTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGc                                            <  1:7739/42‑1 (MQ=37)
                                                     ttgttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGAtt              >  1:3632924/1‑71 (MQ=255)
                                                     ttgttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGAtt              >  1:717324/1‑71 (MQ=255)
                                                     ttgttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGAtt              >  1:1641703/1‑71 (MQ=255)
                                                        ttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGc                                            <  1:84170/38‑1 (MQ=37)
                                                          ggCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGAtt              >  1:1526648/1‑66 (MQ=255)
                                                          ggCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTa          <  1:2938920/70‑1 (MQ=255)
                                                          ggCGCTGGCCCGAGTGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTa          <  1:3075069/70‑1 (MQ=255)
                                                           gcgctGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTa          <  1:569/66‑1 (MQ=255)
                                                                 gCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGCGACt  <  1:450894/71‑1 (MQ=255)
                                                                 gCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGCGACt  <  1:2070924/71‑1 (MQ=255)
                                                                 gCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGCGACt  <  1:1908807/71‑1 (MQ=255)
                                                                 gCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGCGACt  <  1:1677434/71‑1 (MQ=255)
                                                                  cccGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGCGACt  <  1:34546/70‑1 (MQ=255)
                                                                        |                                                                  
AAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGAGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGCGACT  >  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: