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

TGCTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGTAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGT  >  W3110S.gb/3824220‑3824343
                                                            |                                                               
tGCTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTa                                                       >  1:1765379/1‑71 (MQ=255)
tGCTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTa                                                       >  1:3011148/1‑71 (MQ=255)
  cTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTAc                                                      <  1:685833/70‑1 (MQ=255)
     cAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGa                                                   <  1:2452107/70‑1 (MQ=255)
     cAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGa                                                   <  1:3499151/70‑1 (MQ=255)
     cAGGCTGATGGTGGCACCCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGa                                                   >  1:3402402/1‑70 (MQ=255)
      aGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAg                                                 >  1:1185656/1‑71 (MQ=255)
      aGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAg                                                 >  1:1552006/1‑71 (MQ=255)
       ggCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGc                                                <  1:785079/71‑1 (MQ=255)
          tGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCtg                                              >  1:2898407/1‑70 (MQ=255)
            aTGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCtggtg                                           >  1:347732/1‑71 (MQ=255)
             tggtggCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCtggtg                                           <  1:2396549/70‑1 (MQ=255)
               gtggCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAg                                                    <  1:1343257/59‑1 (MQ=255)
               gtggCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGaa                                        <  1:3816593/71‑1 (MQ=255)
               gtggCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGaa                                        <  1:804424/71‑1 (MQ=255)
                   cACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCtggt                                            <  1:682079/63‑1 (MQ=255)
                     cgcgTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTAc                                                      <  1:1064647/51‑1 (MQ=255)
                     cgcgTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTAc                                                      <  1:2942408/51‑1 (MQ=255)
                     cgcgTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGc                                  <  1:149842/71‑1 (MQ=255)
                      gcgTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTAc                                                      <  1:2741423/50‑1 (MQ=255)
                      gcgTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGc                                  >  1:319993/1‑70 (MQ=255)
                      gcgTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGc                                  >  1:1035214/1‑70 (MQ=255)
                       cgTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGa                                                              >  1:1782944/1‑41 (MQ=255)
                       cgTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGa                                                              >  1:2940762/1‑41 (MQ=255)
                        gTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAg                               >  1:3529328/1‑71 (MQ=255)
                          aCCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCt                             >  1:31388/1‑71 (MQ=255)
                           ccGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCt                             >  1:260223/1‑70 (MQ=255)
                            cgcgTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGa                           <  1:1679606/71‑1 (MQ=255)
                             gcgTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGaa                          <  1:200650/71‑1 (MQ=255)
                             gcgTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGaa                          <  1:3543527/71‑1 (MQ=255)
                              cgTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGaa                                        <  1:1018445/56‑1 (MQ=255)
                              cgTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGaa                          >  1:882730/1‑70 (MQ=255)
                               gTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGaaaa                        <  1:1626646/71‑1 (MQ=255)
                                  gATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGaaaa                        >  1:1043797/1‑68 (MQ=255)
                                   aTTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTg                            <  1:3025680/63‑1 (MQ=255)
                                   aTTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTg                            <  1:462902/63‑1 (MQ=255)
                                    ttACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCaac                   >  1:355790/1‑70 (MQ=255)
                                     tACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCaa                                <  1:639065/57‑1 (MQ=255)
                                     tACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGaaa                         >  1:3078708/1‑64 (MQ=255)
                                     tACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCaacc                  >  1:3086974/1‑69 (MQ=255)
                                      aCGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCtggtgg                                          >  1:1259883/1‑46 (MQ=255)
                                       cGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAAccc                 <  1:2118285/70‑1 (MQ=255)
                                       cGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAAccc                 <  1:740737/70‑1 (MQ=255)
                                       cGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCg                <  1:725980/71‑1 (MQ=255)
                                         ggTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGAt              <  1:395804/71‑1 (MQ=255)
                                           tGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGAt              <  1:784379/69‑1 (MQ=255)
                                           tGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa            <  1:1435681/71‑1 (MQ=255)
                                             ccTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGaa                                        <  1:1645113/41‑1 (MQ=255)
                                             ccTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAAccc                 >  1:1995111/1‑64 (MQ=255)
                                             ccTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGaa           <  1:109752/70‑1 (MQ=255)
                                             ccTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa            >  1:3587834/1‑69 (MQ=255)
                                              cTGCGTGGCGCTGGCAGATCCGCTACAGAAGCtggt                                            >  1:1309915/1‑36 (MQ=25)
                                                 cGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGaaaa                        >  1:123999/1‑53 (MQ=255)
                                                 cGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAggg       <  1:3551429/70‑1 (MQ=255)
                                                   tGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCa                                 <  1:1026938/42‑1 (MQ=255)
                                                     gcgcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATgg   <  1:487478/70‑1 (MQ=255)
                                                     gcgcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGt  <  1:215592/71‑1 (MQ=255)
                                                      cgcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATgg   >  1:189376/1‑69 (MQ=255)
                                                      cgcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATgg   >  1:2831145/1‑69 (MQ=255)
                                                      cgcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGt  <  1:509620/70‑1 (MQ=255)
                                                       gcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGa                           <  1:2469428/44‑1 (MQ=255)
                                                            |                                                               
TGCTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGTAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGT  >  W3110S.gb/3824220‑3824343

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: