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

GGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACGT  >  W3110S.gb/3824302‑3824433
                                                                  |                                                                 
ggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACg                                                               >  1:4462847/1‑71 (MQ=255)
  aaaaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggt                                                             <  1:3798818/71‑2 (MQ=255)
   aaaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggtg                                                            <  1:3342841/71‑3 (MQ=255)
   aaaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggtg                                                            <  1:4347890/71‑3 (MQ=255)
     aCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggtg                                                            >  1:1020035/1‑67 (MQ=255)
       ggCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAAc                                                                <  1:2945560/63‑1 (MQ=255)
        gCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGca                                                       <  1:1370020/71‑2 (MQ=255)
        gCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGca                                                       <  1:2034617/71‑2 (MQ=255)
        gCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGca                                                       <  1:1235893/71‑2 (MQ=255)
          aaGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGcaat                                                     <  1:4112534/71‑4 (MQ=255)
               gAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGc                                                 >  1:3574604/1‑70 (MQ=255)
               gAAAACCAACCCGATGAAAGGGATGGTAGACGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGc                                                 >  1:3988098/1‑70 (MQ=255)
                aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGc                                                 <  1:3142755/69‑1 (MQ=255)
                 aaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTCTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGAt                                              >  1:4097017/1‑71 (MQ=255)
                    ccAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAAc                                                                <  1:3958616/50‑1 (MQ=255)
                    ccAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTg                                           >  1:732897/1‑71 (MQ=255)
                    ccAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTg                                           >  1:3457091/1‑71 (MQ=255)
                     cAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTg                                                  <  1:4711422/63‑1 (MQ=255)
                     cAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTg                                                  <  1:2319686/63‑1 (MQ=255)
                     cAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTg                                           >  1:3041745/1‑70 (MQ=255)
                       acccGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGc                                                 >  1:1195222/3‑62 (MQ=255)
                       acccGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCg                                                >  1:1835890/3‑63 (MQ=255)
                         ccGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATa                                      <  1:3488775/71‑1 (MQ=255)
                             tGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGtt                                  >  1:715/1‑71 (MQ=255)
                               aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                                 <  1:1055704/70‑1 (MQ=255)
                               aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                                 >  1:2285206/1‑70 (MQ=255)
                                aaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                               <  1:2084026/71‑1 (MQ=255)
                                aaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                               <  1:4555976/71‑1 (MQ=255)
                                aaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                               <  1:4018923/71‑1 (MQ=255)
                                aaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                               <  1:4401845/71‑1 (MQ=255)
                                 aGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                                 >  1:199667/1‑68 (MQ=255)
                                 aGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                                 >  1:1971432/1‑68 (MQ=255)
                                  gggATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGa                             <  1:1812582/71‑1 (MQ=255)
                                  gggATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGa                             >  1:4418962/1‑71 (MQ=255)
                                   ggATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGa                             <  1:2600756/70‑1 (MQ=255)
                                     aTGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                               <  1:2046483/66‑1 (MQ=255)
                                     aTGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGAct                           >  1:1413903/1‑70 (MQ=255)
                                     aTGGTAGCCGCAGGTTCTGTCTGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                               <  1:350990/66‑1 (MQ=255)
                                       ggTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGActc                          <  1:4439705/69‑1 (MQ=255)
                                       ggTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGActc                          <  1:4541497/69‑1 (MQ=255)
                                        gTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAg                              <  1:4142666/64‑1 (MQ=255)
                                        gTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAg                              <  1:2973066/64‑1 (MQ=255)
                                         tAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGcc                      >  1:2644508/1‑71 (MQ=255)
                                          aGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGa                                >  1:344769/1‑60 (MQ=255)
                                             cGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGtt                                  >  1:515557/1‑55 (MQ=255)
                                             cGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCg                     >  1:2098201/1‑68 (MQ=255)
                                              gCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCaga                 <  1:2382290/71‑1 (MQ=255)
                                              gCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCaga                 <  1:4014277/71‑1 (MQ=255)
                                               cAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCagag                <  1:811571/71‑1 (MQ=255)
                                               cAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                   >  1:1796857/1‑68 (MQ=255)
                                               cAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                   >  1:3347983/1‑68 (MQ=255)
                                                 gTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGAc              >  1:1708692/1‑71 (MQ=255)
                                                  tttCTGTCGGTATTGTAAACGGTGAAGCAACTTGCGATCTGGAATACGTTGAAGAc                            <  1:4157827/56‑1 (MQ=255)
                                                   ttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGaa                                        >  1:2782281/1‑43 (MQ=19)
                                                   ttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGaa                                        >  1:2171570/1‑43 (MQ=19)
                                                   ttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGAcc             >  1:3281685/1‑70 (MQ=255)
                                                   ttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGAcc             >  1:4556967/1‑70 (MQ=255)
                                                    tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                               >  1:2310230/1‑51 (MQ=255)
                                                    tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                               >  1:4549621/1‑51 (MQ=255)
                                                     cTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGActc                          <  1:3631139/55‑1 (MQ=255)
                                                     cTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGActc                          <  1:4521166/55‑1 (MQ=255)
                                                     cTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGActc                          <  1:203206/55‑1 (MQ=255)
                                                      tGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTgg                                          <  1:586474/38‑1 (MQ=16)
                                                      tGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                   <  1:749720/61‑1 (MQ=255)
                                                          ggtaTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGaa     >  1:2567568/4‑71 (MQ=255)
                                                          ggtaTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGaa     >  1:3605688/4‑71 (MQ=255)
                                                           gtaTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGc                    <  1:2546241/53‑1 (MQ=255)
                                                           gtaTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                   >  1:451038/3‑56 (MQ=255)
                                                           gtaTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                   >  1:1657567/3‑56 (MQ=255)
                                                             aTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACAt        >  1:299284/1‑65 (MQ=255)
                                                              ttGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCg                     <  1:2915673/51‑1 (MQ=255)
                                                              ttGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                   <  1:742734/53‑1 (MQ=255)
                                                              ttGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGt  <  1:3199156/70‑1 (MQ=255)
                                                                  |                                                                 
GGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACGT  >  W3110S.gb/3824302‑3824433

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: