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

CGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACGTAGT  >  W3110S.gb/3824308‑3824436
                                                                  |                                                              
cGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAg                                                            >  1:2116206/1‑71 (MQ=255)
 ggCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGc                                                           <  1:1979421/71‑1 (MQ=255)
  gCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGca                                                          <  1:1287921/71‑2 (MQ=255)
   cAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGcaa                                                         <  1:2634669/71‑3 (MQ=255)
   cAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGcaa                                                         >  1:1312015/1‑69 (MQ=255)
   cAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGcaa                                                         <  1:875142/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:1957157/1‑71 (MQ=255)
            aaccaaccCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATc                                                >  1:1383021/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:2361219/70‑1 (MQ=255)
                         aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                                    <  1:2583787/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:830817/1‑40 (MQ=16)
                                          aGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGAt                                                 >  1:272810/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)
                                                        ttGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                                  <  1:1540805/41‑1 (MQ=25)
                                                         tgtaAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAg   >  1:1626945/1‑71 (MQ=255)
                                                          gtaaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACAt           >  1:80629/4‑62 (MQ=255)
                                                          gtaaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt  >  1:518804/4‑71 (MQ=255)
                                                          gtaaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt  >  1:1643029/4‑71 (MQ=255)
                                                          gtaaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt  >  1:1669191/4‑71 (MQ=255)
                                                                  |                                                              
CGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACGTAGT  >  W3110S.gb/3824308‑3824436

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: