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% 118.9 / NA 43pseudogene (495/687 nt)rphECK3633:JW3618:b3643; ribonuclease PH
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (18/25);  total (18/25)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACAT  >  W3110S.gb/3824302‑3824427
                                                            |                                                                 
ggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACg                                                         >  1:2076975/1‑71 (MQ=255)
      cGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAg                                                   >  1:2116206/1‑71 (MQ=255)
       ggCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtg                                                               <  1:2609078/58‑1 (MQ=255)
       ggCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGc                                                  <  1:1979421/71‑1 (MQ=255)
        gCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGca                                                 <  1:1287921/71‑2 (MQ=255)
         cAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGcaa                                                <  1:875142/71‑3 (MQ=255)
         cAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGcaa                                                >  1:1312015/1‑69 (MQ=255)
         cAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGcaa                                                <  1:2634669/71‑3 (MQ=255)
                aaaaCCAACCCGATGAAAGGGATTGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGa                                         <  1:384074/71‑1 (MQ=255)
                aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGc                                           <  1:2029839/69‑1 (MQ=255)
                aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGc                                           <  1:1657110/69‑1 (MQ=255)
                 aaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGAt                                        <  1:2339737/71‑1 (MQ=255)
                  aaccaaccCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATc                                       >  1:1383021/1‑71 (MQ=255)
                  aaccaaccCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATc                                       >  1:1957157/1‑71 (MQ=255)
                       acccGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGaa                                  <  1:886792/69‑1 (MQ=255)
                       acccGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGaa                                  <  1:2015998/69‑1 (MQ=255)
                            aTGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATa                                <  1:1381082/68‑1 (MQ=255)
                              gAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGtt                            >  1:1498527/1‑70 (MQ=255)
                              gAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGtt                            >  1:207934/1‑70 (MQ=255)
                               aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGa                                         <  1:593433/56‑1 (MQ=255)
                               aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                           <  1:2583787/70‑1 (MQ=255)
                               aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                           <  1:2361219/70‑1 (MQ=255)
                                aaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                         <  1:105870/71‑1 (MQ=255)
                                 aGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGa                                   >  1:106911/1‑60 (MQ=255)
                                 aGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAg                        <  1:103834/71‑1 (MQ=255)
                                 aGGGATGGTAGCCGAAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTgg                                    >  1:1679366/1‑59 (MQ=255)
                                      tGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTgg                                    <  1:1018677/54‑1 (MQ=255)
                                        gTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCt                                      <  1:1445567/50‑1 (MQ=255)
                                        gTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGc                 >  1:967319/1‑71 (MQ=255)
                                         tAGCCGCAGTTTCTGTCTGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGtt                            <  1:404111/59‑1 (MQ=255)
                                            ccGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGc                                           <  1:1006490/41‑1 (MQ=19)
                                            ccGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                           <  1:1698249/57‑1 (MQ=255)
                                                aGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGAt                                        >  1:272810/1‑40 (MQ=16)
                                                aGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGAt                                        >  1:830817/1‑40 (MQ=16)
                                                 gTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGAc        >  1:1306627/1‑71 (MQ=255)
                                                  tttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCg                                          <  1:2173982/36‑1 (MQ=12)
                                                  tttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGAt                                        >  1:2582116/1‑38 (MQ=16)
                                                  tttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCt                                      <  1:85347/40‑1 (MQ=16)
                                                  tttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATa                                <  1:1019973/46‑1 (MQ=25)
                                                    tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                         >  1:1772663/1‑51 (MQ=255)
                                                    tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa     >  1:782049/1‑71 (MQ=255)
                                                    tCTGTCGGTATTGTAAACGGTGAAGAAATTTGCGATCTGGAATACGTTg                           >  1:1831823/1‑49 (MQ=19)
                                                        tCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACAt  >  1:1376429/1‑70 (MQ=255)
                                                            |                                                                 
GGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACAT  >  W3110S.gb/3824302‑3824427

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: