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

AGCTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATG  >  W3110S.gb/3824295‑3824428
                                                                   |                                                                  
aGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTAtt                                                                 >  1:3680275/1‑71 (MQ=255)
      tggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtgta                                                              <  1:2144706/68‑2 (MQ=255)
        gAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACgg                                                         >  1:3952797/1‑71 (MQ=255)
          aaaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggtg                                                       <  1:4236526/71‑3 (MQ=255)
           aaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGa                                                      <  1:1135845/71‑1 (MQ=255)
              ggCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGaa                                                     >  1:756269/1‑69 (MQ=255)
              ggCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGc                                                   <  1:79219/71‑1 (MQ=255)
              ggCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGc                                                   <  1:2292187/71‑1 (MQ=255)
                cAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtgta                                                              <  1:2525647/58‑2 (MQ=255)
                 aaGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGcaat                                                >  1:4226816/1‑68 (MQ=255)
                    cTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggtg                                                       <  1:1375978/61‑3 (MQ=255)
                    cTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAg                                                    <  1:494591/64‑1 (MQ=255)
                    cTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTg                                             <  1:2035815/71‑1 (MQ=255)
                     tGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGc                                            >  1:2328291/1‑71 (MQ=255)
                     tGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGc                                            >  1:329634/1‑71 (MQ=255)
                      gAAAACCAACCCGATGAAAGGGATGGTAGCCGAAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGc                                            >  1:3915910/1‑70 (MQ=255)
                       aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggtg                                                       <  1:1856311/58‑3 (MQ=255)
                       aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCg                                           <  1:383498/70‑1 (MQ=255)
                         aaccaaccCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGAt                                         <  1:620597/70‑1 (MQ=255)
                         aaccaaccCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATc                                        <  1:1702977/71‑1 (MQ=255)
                         aaccaaccCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATc                                        >  1:4108682/1‑71 (MQ=255)
                          aCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCt                                       <  1:319233/71‑1 (MQ=255)
                          aCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCt                                       >  1:4502242/1‑71 (MQ=255)
                           ccAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTg                                      <  1:4296311/71‑1 (MQ=255)
                           ccAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTg                                      >  1:4411400/1‑71 (MQ=255)
                            cAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCt                                       >  1:2810242/1‑69 (MQ=255)
                              acccGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCg                                           <  1:1873219/61‑1 (MQ=255)
                                 cGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTgg                                     <  1:2621737/66‑1 (MQ=255)
                                   aTGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtgta                                                              <  1:3740567/39‑2 (MQ=39)
                                    tGAAAGGGATGGTAGCCGCAGTTTCTGTCTGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGtt                             <  1:4020840/71‑1 (MQ=255)
                                    tGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGc                                            <  1:1537975/56‑1 (MQ=255)
                                    tGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCg                                           <  1:391973/57‑1 (MQ=255)
                                     gAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTaaa                                                            >  1:4309337/1‑39 (MQ=38)
                                     gAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTaaa                                                            >  1:1276784/1‑39 (MQ=38)
                                     gAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGAt                                         >  1:2731633/1‑58 (MQ=255)
                                     gAAAGGGATGGTAGCCGCAGTTTATGTCGGTATTGTaaa                                                            >  1:1209041/1‑39 (MQ=37)
                                      aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                            <  1:797022/70‑1 (MQ=255)
                                      aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                            <  1:1280933/70‑1 (MQ=255)
                                       aaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGaa                                                     >  1:4291110/1‑44 (MQ=38)
                                       aaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                          <  1:113090/71‑1 (MQ=255)
                                        aGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                            >  1:3318424/1‑68 (MQ=255)
                                         gggATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGa                        >  1:3361526/1‑71 (MQ=255)
                                          ggATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACgg                                                         <  1:1544957/37‑1 (MQ=38)
                                          ggATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACgg                                                         <  1:624535/37‑1 (MQ=38)
                                          ggATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGAc                       >  1:550385/1‑71 (MQ=255)
                                          ggATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGAc                       >  1:1352855/1‑71 (MQ=255)
                                          ggATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGAc                       <  1:4664406/71‑1 (MQ=255)
                                             tGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCg                                           <  1:1613280/48‑1 (MQ=25)
                                             tGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCg                                           <  1:3832069/48‑1 (MQ=25)
                                             tGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                            <  1:990900/63‑1 (MQ=255)
                                               gTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGa                                    >  1:265191/1‑53 (MQ=255)
                                               gTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGa                                    >  1:2726925/1‑53 (MQ=255)
                                               gTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGActcaga                  >  1:3309461/1‑68 (MQ=255)
                                               gTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGc                  >  1:4437080/1‑71 (MQ=255)
                                               gTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGc                  >  1:1692238/1‑71 (MQ=255)
                                                 aGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCg                >  1:969216/1‑71 (MQ=255)
                                                 aGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCg                >  1:4756674/1‑71 (MQ=255)
                                                 aGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCg                >  1:1403140/1‑71 (MQ=255)
                                                 aGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCg                >  1:2969710/1‑71 (MQ=255)
                                                   ccGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa              <  1:2343926/71‑1 (MQ=255)
                                                     gCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa              >  1:1084617/1‑69 (MQ=255)
                                                      cAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGAt                                         >  1:4676068/1‑41 (MQ=19)
                                                      cAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGAt                                         >  1:4783037/1‑41 (MQ=19)
                                                      cAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGAt                                         >  1:138282/1‑41 (MQ=19)
                                                        gTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTgg                                     >  1:4299779/1‑43 (MQ=19)
                                                         tttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTgg                                     <  1:351855/42‑1 (MQ=19)
                                                         tttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                            <  1:2752766/51‑1 (MQ=255)
                                                         tttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                            <  1:508682/51‑1 (MQ=255)
                                                         tttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                            <  1:1164690/51‑1 (MQ=255)
                                                         tttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGc                  <  1:88627/61‑1 (MQ=255)
                                                         tttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa              <  1:2808130/65‑1 (MQ=255)
                                                          ttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAgaga          >  1:4251805/1‑68 (MQ=255)
                                                          ttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAgaga          >  1:4257332/1‑68 (MQ=255)
                                                          ttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGAcc        >  1:4250726/1‑70 (MQ=255)
                                                           tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGaa                                   >  1:3722172/1‑42 (MQ=19)
                                                           tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTTAAGACTCAGCCGCAGAGACCGa      >  1:4375180/1‑71 (MQ=255)
                                                           tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                          >  1:2096770/1‑51 (MQ=255)
                                                           tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                          >  1:705350/1‑51 (MQ=255)
                                                           tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                          >  1:31162/1‑51 (MQ=255)
                                                           tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                          >  1:1282908/1‑51 (MQ=255)
                                                           tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGa                        >  1:822230/1‑53 (MQ=255)
                                                           tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa              >  1:3595046/1‑63 (MQ=255)
                                                           tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa              >  1:61881/1‑63 (MQ=255)
                                                           tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCg       <  1:3410738/70‑1 (MQ=255)
                                                           tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCg       <  1:3095358/70‑1 (MQ=255)
                                                           tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCg       <  1:3101014/70‑1 (MQ=255)
                                                           tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa      >  1:2294059/1‑71 (MQ=255)
                                                           tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa      >  1:3160188/1‑71 (MQ=255)
                                                             tGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACa    >  1:3662788/1‑71 (MQ=255)
                                                               tCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACAt   <  1:3682119/70‑1 (MQ=255)
                                                               tCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATg  <  1:2876150/71‑1 (MQ=255)
                                                               tCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATg  <  1:2158607/71‑1 (MQ=255)
                                                                cGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATg  >  1:4333780/1‑70 (MQ=255)
                                                                   |                                                                  
AGCTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATG  >  W3110S.gb/3824295‑3824428

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: