Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,824,362 A→T 100% pseudogene (495/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,3620AT100.0% 96.8 / NA 35pseudogene (495/687 nt)rphECK3633:JW3618:b3643; ribonuclease PH
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (13/22);  total (13/22)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GCTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAG  >  W3110S.gb/3824296‑3824419
                                                                  |                                                         
gCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTAt                                                         <  1:2070653/69‑1 (MQ=255)
  tggtggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTAtt                                                        >  1:2257668/1‑68 (MQ=255)
      ggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtgta                                                     <  1:1430028/67‑2 (MQ=255)
      ggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtgta                                                     <  1:2385950/67‑2 (MQ=255)
      ggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTaaa                                                   >  1:2711419/1‑69 (MQ=255)
                       aaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGc                                   >  1:2665610/1‑68 (MQ=255)
                       aaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCg                                  >  1:1076050/1‑69 (MQ=255)
                       aaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCg                                  >  1:3133106/1‑69 (MQ=255)
                             acccGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTgg                            <  1:1211448/67‑1 (MQ=255)
                             acccGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTgg                            <  1:3128159/67‑1 (MQ=255)
                                cGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAAt                         >  1:385139/1‑69 (MQ=255)
                                    gAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCg                                  <  1:2746670/56‑1 (MQ=255)
                                      aaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCg                                  <  1:2244614/54‑1 (MQ=255)
                                      aaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCg                                  <  1:2433668/54‑1 (MQ=255)
                                       aGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGt                     >  1:3158052/1‑66 (MQ=255)
                                       aGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                   <  1:2272051/68‑1 (MQ=255)
                                        gggATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                 <  1:75732/69‑1 (MQ=255)
                                        gggATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                 >  1:64447/1‑69 (MQ=255)
                                        gggATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                 <  1:1900428/69‑1 (MQ=255)
                                        gggATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                 >  1:401819/1‑69 (MQ=255)
                                        gggATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                 >  1:3468036/1‑69 (MQ=255)
                                                    gCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa     <  1:525453/69‑1 (MQ=255)
                                                         ttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCagat  <  1:885086/67‑2 (MQ=255)
                                                         ttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCagag  <  1:3194085/67‑1 (MQ=255)
                                                         ttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCagag  <  1:2743509/67‑1 (MQ=255)
                                                         ttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCagag  <  1:2285045/67‑1 (MQ=255)
                                                         ttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCagag  <  1:770426/67‑1 (MQ=255)
                                                          tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTgg                            <  1:1338804/40‑1 (MQ=16)
                                                          tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTgg                            <  1:1742436/40‑1 (MQ=16)
                                                          tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTgg                            <  1:1840713/40‑1 (MQ=16)
                                                          tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGcc        >  1:1230306/1‑60 (MQ=255)
                                                          tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGcc        >  1:55162/1‑60 (MQ=255)
                                                          tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGcc        >  1:1914151/1‑60 (MQ=255)
                                                          tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGc         <  1:2779168/59‑1 (MQ=255)
                                                          tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGc         <  1:2080452/59‑1 (MQ=255)
                                                                  |                                                         
GCTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAG  >  W3110S.gb/3824296‑3824419

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