Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,824,380 2 bp→AA 100% pseudogene (513‑514/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,3800GA100.0% 111.0 / NA 41pseudogene (513/687 nt)rphECK3633:JW3618:b3643; ribonuclease PH
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (20/21);  total (20/21)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
*W3110S.gb3,824,3810GA100.0% 112.3 / NA 41pseudogene (514/687 nt)rphECK3633:JW3618:b3643; ribonuclease PH
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (20/21);  total (20/21)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACGTAGT  >  W3110S.gb/3824318‑3824436
                                                              ||                                                       
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:2015998/69‑1 (MQ=255)
       acccGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGaa                                           <  1:886792/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: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)
                                                              ||                                                       
AAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACGTAGT  >  W3110S.gb/3824318‑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: