Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,824,524 T→G 100% pseudogene (657/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,5240TG98.0% 165.4 / ‑3.9 50pseudogene (657/687 nt)rphECK3633:JW3618:b3643; ribonuclease PH
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base G (21/28);  minor base . (0/1);  total (21/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-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.

ATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGAGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTA  >  W3110S.gb/3824459‑3824589
                                                                 |                                                                   
atTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGc                                                                >  1:2730131/1‑71 (MQ=255)
 ttGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGcc                                                               >  1:3247575/1‑71 (MQ=255)
 ttGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGcc                                                               >  1:460214/1‑71 (MQ=255)
   gAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGccc                                                              >  1:820586/1‑70 (MQ=255)
    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:3838589/71‑3 (MQ=255)
          aGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAggggga                                                     <  1:3947732/71‑3 (MQ=255)
            gggACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGccc                                                              <  1:1172937/61‑1 (MQ=255)
            gggACGGCAAAAGACGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCg                                                             <  1:2700233/62‑1 (MQ=255)
              gACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGG‑‑CTGGCCCGAGGGGgaatcga                                                <  1:2021028/70‑1 (MQ=255)
                       aGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTg                                        >  1:2467755/1‑71 (MQ=255)
                        ggCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCg                                                             <  1:4379975/50‑1 (MQ=255)
                          cGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTa                                      >  1:1603847/1‑70 (MQ=255)
                            aGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCg                                   <  1:3029150/71‑1 (MQ=255)
                                gTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGccc                                                              <  1:3403448/41‑1 (MQ=38)
                                  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:3221245/45‑1 (MQ=255)
                                        aTGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcg                                                 <  1:3536623/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)
                                                                 |                                                                   
ATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGAGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTA  >  W3110S.gb/3824459‑3824589

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: