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

GAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACGTAGTGATGACCG  >  W3110S.gb/3824317‑3824444
                                                               ||                                                               
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)
     ccAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTg                                                      >  1:3457091/1‑71 (MQ=255)
     ccAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTg                                                      >  1:732897/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:4018923/71‑1 (MQ=255)
                 aaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                                          <  1:4401845/71‑1 (MQ=255)
                 aaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                                          <  1:4555976/71‑1 (MQ=255)
                  aGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                                            >  1:1971432/1‑68 (MQ=255)
                  aGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                                            >  1:199667/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:4541497/69‑1 (MQ=255)
                        ggTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGActc                                     <  1:4439705/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: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:2171570/1‑43 (MQ=19)
                                    ttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGaa                                                   >  1:2782281/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:3605688/4‑71 (MQ=255)
                                           ggtaTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGaa                >  1:2567568/4‑71 (MQ=255)
                                            gtaTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGc                               <  1:2546241/53‑1 (MQ=255)
                                            gtaTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                              >  1:1657567/3‑56 (MQ=255)
                                            gtaTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                              >  1:451038/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)
                                                 gtaaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt          >  1:2802439/4‑71 (MQ=255)
                                                 gtaaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt          >  1:3588257/4‑71 (MQ=255)
                                                 gtaaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt          >  1:755263/4‑71 (MQ=255)
                                                  taaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGa                                        >  1:3369696/3‑40 (MQ=25)
                                                   aaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtga        >  1:1728009/2‑71 (MQ=255)
                                                    aaCGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                                          >  1:2851310/1‑36 (MQ=25)
                                                    aaCGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                                          >  1:1337480/1‑36 (MQ=25)
                                                    aaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtga        <  1:4152279/70‑1 (MQ=255)
                                                     aCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtgat       >  1:613142/1‑70 (MQ=255)
                                                     aCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtgat       >  1:3651296/1‑70 (MQ=255)
                                                     aCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtgat       >  1:976764/1‑70 (MQ=255)
                                                       ggtgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGAc    >  1:28994/4‑71 (MQ=255)
                                                         tgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCg  <  1:4042876/70‑1 (MQ=255)
                                                         tgAAGCAATTTGCGAGCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGAcc   <  1:1028428/69‑1 (MQ=255)
                                                               ||                                                               
GAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACGTAGTGATGACCG  >  W3110S.gb/3824317‑3824444

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: