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

GTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACG  >  W3110S.gb/3824300‑3824432
                                                                    |                                                                
gtggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTaaa                                                                >  1:1640399/1‑71 (MQ=255)
 tggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTaaa                                                                >  1:2743344/1‑70 (MQ=255)
 tggAAAACGGCAAGCCGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAAc                                                               <  1:2246619/71‑1 (MQ=255)
  ggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAAc                                                               <  1:1785498/70‑1 (MQ=255)
            aaGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACg                                                              >  1:835850/1‑61 (MQ=255)
              gCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggtg                                                           <  1:713490/62‑3 (MQ=255)
                 gAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTaaa                                                                <  1:233126/54‑1 (MQ=255)
                 gAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGc                                                >  1:2322201/1‑70 (MQ=255)
                  aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGca                                                      >  1:1680602/1‑62 (MQ=255)
                  aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCg                                               <  1:2978456/70‑1 (MQ=255)
                  aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCg                                               <  1:2938087/70‑1 (MQ=255)
                  aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGa                                              >  1:952803/1‑71 (MQ=255)
                   aaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAAc                                                               >  1:1595356/1‑53 (MQ=255)
                   aaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGa                                              <  1:2198851/70‑1 (MQ=255)
                    aaccaaccCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATc                                            >  1:1453574/1‑71 (MQ=255)
                     aCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAg                                                        >  1:1311104/1‑58 (MQ=255)
                           ccGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATa                                     <  1:2311671/71‑1 (MQ=255)
                              aTGAAAGGGATGGTAGCCTCAGTTTCTGTCGGTATTGTAAACGGTGAAGcaa                                                     >  1:3100767/1‑50 (MQ=255)
                              aTGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTgg                                         <  1:1344436/64‑1 (MQ=255)
                                gAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACg                                                              >  1:1485738/1‑41 (MQ=38)
                                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:1637724/70‑1 (MQ=255)
                                             gCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCg                    >  1:3338537/1‑70 (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:769137/1‑71 (MQ=255)
                                                      tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa          >  1:297211/1‑71 (MQ=255)
                                                      tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa          >  1:2055310/1‑71 (MQ=255)
                                                      tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa          >  1:3337848/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:1795261/70‑1 (MQ=255)
                                                              taTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg  <  1:1160668/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)
                                                                    |                                                                
GTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACG  >  W3110S.gb/3824300‑3824432

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: