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

AAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAAC  >  W3110S.gb/3824304‑3824431
                                                                |                                                               
aaaaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggt                                                           <  1:855431/71‑2 (MQ=255)
 aaaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggtg                                                          <  1:1762149/71‑3 (MQ=255)
 aaaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggt                                                           >  1:866780/1‑69 (MQ=255)
 aaaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggt                                                           >  1:1592648/1‑69 (MQ=255)
     ggCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGc                                                      <  1:700526/71‑1 (MQ=255)
      gCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGca                                                     <  1:1241079/71‑2 (MQ=255)
       cAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGca                                                     >  1:278509/1‑69 (MQ=255)
        aaGCTGAAAACCAACCCGATTAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGcaat                                                   <  1:2031320/71‑4 (MQ=255)
               aaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggtg                                                          >  1:315460/1‑55 (MQ=255)
                aaccaaccCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATc                                           <  1:777151/71‑1 (MQ=255)
                              aaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                             <  1:77281/71‑1 (MQ=255)
                                 ggATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGAc                          >  1:2010852/1‑71 (MQ=255)
                                  gATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGaa                                      >  1:637724/1‑58 (MQ=255)
                                     ggTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGActcag                      >  1:1333260/1‑69 (MQ=255)
                                          ccGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                 <  1:2065419/71‑1 (MQ=255)
                                           cGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGActc                        <  1:1407855/63‑1 (MQ=255)
                                           cGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCag                <  1:1950099/71‑1 (MQ=255)
                                           cGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                 <  1:1604829/70‑1 (MQ=255)
                                              aGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGc                  >  1:1585674/1‑66 (MQ=255)
                                              aGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAgaga             >  1:1213613/1‑71 (MQ=255)
                                                      tCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGActca                       <  1:1927840/53‑2 (MQ=255)
                                                         gtaTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCaga               >  1:1870513/3‑58 (MQ=255)
                                                         gtaTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAAc  <  1:1325359/69‑1 (MQ=255)
                                                                |                                                               
AAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAAC  >  W3110S.gb/3824304‑3824431

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