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

GCTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACC  >  W3110S.gb/3824296‑3824422
                                                                  |                                                            
gCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtg                                                          <  1:1088996/71‑1 (MQ=255)
 cTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtgt                                                         <  1:805309/71‑1 (MQ=255)
    gtggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtgta                                                        <  1:3960378/69‑2 (MQ=255)
    gtggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtgta                                                        <  1:430018/69‑2 (MQ=255)
    gtggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtgta                                                        <  1:1413347/69‑2 (MQ=255)
     tggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtgta                                                        <  1:3588770/68‑2 (MQ=255)
      ggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACg                                                    >  1:4462847/1‑71 (MQ=255)
        aaaaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggt                                                  <  1:3798818/71‑2 (MQ=255)
         aaaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggtg                                                 <  1:4347890/71‑3 (MQ=255)
         aaaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggtg                                                 <  1:3342841/71‑3 (MQ=255)
           aCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggtg                                                 >  1:1020035/1‑67 (MQ=255)
             ggCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAAc                                                     <  1:2945560/63‑1 (MQ=255)
              gCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGca                                            <  1:2034617/71‑2 (MQ=255)
              gCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGca                                            <  1:1235893/71‑2 (MQ=255)
              gCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGca                                            <  1:1370020/71‑2 (MQ=255)
                aaGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGcaat                                          <  1:4112534/71‑4 (MQ=255)
                   cTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtgta                                                        <  1:1418172/54‑2 (MQ=255)
                    tGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtg                                                          <  1:2521594/51‑1 (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:3457091/1‑71 (MQ=255)
                          ccAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTg                                >  1:732897/1‑71 (MQ=255)
                           cAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTg                                       <  1:2319686/63‑1 (MQ=255)
                           cAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTg                                       <  1:4711422/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:2285206/1‑70 (MQ=255)
                                     aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                      <  1:1055704/70‑1 (MQ=255)
                                      aaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                    <  1:4018923/71‑1 (MQ=255)
                                      aaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                    <  1:4401845/71‑1 (MQ=255)
                                      aaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                    <  1:4555976/71‑1 (MQ=255)
                                      aaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                    <  1:2084026/71‑1 (MQ=255)
                                       aGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                      >  1:199667/1‑68 (MQ=255)
                                       aGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                      >  1:1971432/1‑68 (MQ=255)
                                        gggATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGa                  >  1:4418962/1‑71 (MQ=255)
                                        gggATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGa                  <  1:1812582/71‑1 (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:2973066/64‑1 (MQ=255)
                                              gTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAg                   <  1:4142666/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:4014277/71‑1 (MQ=255)
                                                    gCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCaga      <  1:2382290/71‑1 (MQ=255)
                                                     cAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCagag     <  1:811571/71‑1 (MQ=255)
                                                     cAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa        >  1:3347983/1‑68 (MQ=255)
                                                     cAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa        >  1:1796857/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:4521166/55‑1 (MQ=255)
                                                           cTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGActc               <  1:3631139/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)
                                                                  |                                                            
GCTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACC  >  W3110S.gb/3824296‑3824422

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: