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.

ATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGA‑GGGG‑AATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTA  >  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)
         cAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑ggggg                                                      <  1:4417846/71‑2 (MQ=255)
         cAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑gggg                                                       >  1:4168668/1‑70 (MQ=255)
          aGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑ggggga                                                     <  1:3838589/71‑3 (MQ=255)
          aGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑ggggga                                                     <  1:3947732/71‑3 (MQ=255)
            gggACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGccc                                                              <  1:1172937/61‑1 (MQ=255)
            gggACGGCAAAAGACGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCg                                                             <  1:2700233/62‑1 (MQ=255)
              gACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGG‑‑CTGGCCCGAGGGGg‑aatcga                                                <  1:2021028/70‑1 (MQ=255)
                       aGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGG‑AATCGAATCCATTg                                        >  1:2467755/1‑71 (MQ=255)
                        ggCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCg                                                             <  1:4379975/50‑1 (MQ=255)
                          cGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGG‑AATCGAATCCATTGTa                                      >  1:1603847/1‑70 (MQ=255)
                            aGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGG‑AATCGAATCCATTGTAGCg                                   <  1:3029150/71‑1 (MQ=255)
                                gTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGccc                                                              <  1:3403448/41‑1 (MQ=38)
                                  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:3221245/45‑1 (MQ=255)
                                        aTGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGg‑aatcg                                                 <  1:3536623/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)
                                                                 |                                                                   
ATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGA‑GGGG‑AATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTA  >  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: ‑