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% 148.3 / NA 51pseudogene (513/687 nt)rphECK3633:JW3618:b3643; ribonuclease PH
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (23/28);  total (23/28)
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% 147.4 / NA 51pseudogene (514/687 nt)rphECK3633:JW3618:b3643; ribonuclease PH
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (23/28);  total (23/28)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACGTAGTGAT  >  W3110S.gb/3824317‑3824439
                                                               ||                                                          
gAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGc                                                       >  1:2322201/1‑70 (MQ=255)
 aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCg                                                      <  1:2978456/70‑1 (MQ=255)
 aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCg                                                      <  1:2938087/70‑1 (MQ=255)
 aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGa                                                     >  1:952803/1‑71 (MQ=255)
  aaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGa                                                     <  1:2198851/70‑1 (MQ=255)
   aaccaaccCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATc                                                   >  1:1453574/1‑71 (MQ=255)
          ccGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATa                                            <  1:2311671/71‑1 (MQ=255)
             aTGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTgg                                                <  1:1344436/64‑1 (MQ=255)
               gAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGtt                                        >  1:124476/1‑70 (MQ=255)
                aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                                       <  1:3311843/70‑1 (MQ=255)
                aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                                       <  1:3238901/70‑1 (MQ=255)
                 aaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                                       >  1:2753446/1‑69 (MQ=255)
                 aaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                                     <  1:2148535/71‑1 (MQ=255)
                  aGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAg                                    <  1:359915/71‑1 (MQ=255)
                   gggATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGAGATCTg                                                 >  1:80718/1‑57 (MQ=255)
                    ggATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGAc                                  <  1:2153997/71‑1 (MQ=255)
                            gCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCg                           >  1:3338537/1‑70 (MQ=255)
                            gCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCg                           <  1:1637724/70‑1 (MQ=255)
                             ccGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCg                           <  1:3126838/69‑1 (MQ=255)
                             ccGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                         <  1:268550/71‑1 (MQ=255)
                                   tttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTgg                                                <  1:324375/42‑1 (MQ=19)
                                    ttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGAcc                   >  1:3429196/1‑70 (MQ=255)
                                     tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGaa                                              >  1:2728424/1‑42 (MQ=19)
                                     tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa                 >  1:2055310/1‑71 (MQ=255)
                                     tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa                 >  1:769137/1‑71 (MQ=255)
                                     tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa                 >  1:3337848/1‑71 (MQ=255)
                                     tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa                 >  1:297211/1‑71 (MQ=255)
                                      cTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                                       <  1:3000520/48‑1 (MQ=25)
                                      cTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                         >  1:2608084/1‑62 (MQ=255)
                                         tCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGc                          <  1:3194150/58‑1 (MQ=255)
                                             taTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                         <  1:625125/54‑1 (MQ=255)
                                             taTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg         <  1:2951511/70‑1 (MQ=255)
                                             taTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg         <  1:1160668/70‑1 (MQ=255)
                                             taTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg         <  1:1795261/70‑1 (MQ=255)
                                               ttGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGa                                   <  1:894895/43‑1 (MQ=25)
                                               ttGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                         <  1:2365629/53‑1 (MQ=255)
                                               ttGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg         >  1:3573473/1‑69 (MQ=255)
                                                 gtaaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt     >  1:2707807/4‑71 (MQ=255)
                                                 gtaaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt     >  1:2701837/4‑71 (MQ=255)
                                                 gtaaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt     >  1:2421184/4‑71 (MQ=255)
                                                 gtaaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt     >  1:232639/4‑71 (MQ=255)
                                                  taaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg         <  1:2844892/64‑1 (MQ=255)
                                                  taaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtg    <  1:1720416/69‑1 (MQ=255)
                                                  taaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtg    <  1:2218113/69‑1 (MQ=255)
                                                  taaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt     <  1:3650140/68‑1 (MQ=255)
                                                    aaCGGTGAAGCAATTTGCGATCTGGAATACGTTGTAGa                                   >  1:847218/1‑38 (MQ=21)
                                                    aaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCaga                       >  1:1244153/1‑50 (MQ=255)
                                                    aaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                         <  1:1083781/48‑1 (MQ=37)
                                                    aaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATg             >  1:160377/1‑60 (MQ=255)
                                                    aaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtgat  >  1:3351123/1‑71 (MQ=255)
                                                     aCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg         <  1:463985/63‑1 (MQ=255)
                                                               ||                                                          
GAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACGTAGTGAT  >  W3110S.gb/3824317‑3824439

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: