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

CTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACGTAGTGATGACCG  >  W3110S.gb/3824315‑3824444
                                                                 ||                                                               
cTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTg                                                             <  1:2035815/71‑1 (MQ=255)
 tGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGc                                                            >  1:329634/1‑71 (MQ=255)
 tGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGc                                                            >  1:2328291/1‑71 (MQ=255)
  gAAAACCAACCCGATGAAAGGGATGGTAGCCGAAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGc                                                            >  1:3915910/1‑70 (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:4411400/1‑71 (MQ=255)
       ccAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTg                                                      <  1:4296311/71‑1 (MQ=255)
        cAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCt                                                       >  1:2810242/1‑69 (MQ=255)
          acccGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCg                                                           <  1:1873219/61‑1 (MQ=255)
             cGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTgg                                                     <  1:2621737/66‑1 (MQ=255)
                tGAAAGGGATGGTAGCCGCAGTTTCTGTCTGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGtt                                             <  1:4020840/71‑1 (MQ=255)
                tGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGc                                                            <  1:1537975/56‑1 (MQ=255)
                tGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCg                                                           <  1:391973/57‑1 (MQ=255)
                 gAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGAt                                                         >  1:2731633/1‑58 (MQ=255)
                  aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                                            <  1:1280933/70‑1 (MQ=255)
                  aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                                            <  1:797022/70‑1 (MQ=255)
                   aaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                                          <  1:113090/71‑1 (MQ=255)
                    aGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                                            >  1:3318424/1‑68 (MQ=255)
                     gggATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGa                                        >  1:3361526/1‑71 (MQ=255)
                      ggATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGAc                                       >  1:1352855/1‑71 (MQ=255)
                      ggATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGAc                                       <  1:4664406/71‑1 (MQ=255)
                      ggATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGAc                                       >  1:550385/1‑71 (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:2726925/1‑53 (MQ=255)
                           gTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGa                                                    >  1:265191/1‑53 (MQ=255)
                           gTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGActcaga                                  >  1:3309461/1‑68 (MQ=255)
                           gTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGc                                  >  1:1692238/1‑71 (MQ=255)
                           gTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGc                                  >  1:4437080/1‑71 (MQ=255)
                             aGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCg                                >  1:2969710/1‑71 (MQ=255)
                             aGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCg                                >  1:4756674/1‑71 (MQ=255)
                             aGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCg                                >  1:969216/1‑71 (MQ=255)
                             aGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCg                                >  1:1403140/1‑71 (MQ=255)
                               ccGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                              <  1:2343926/71‑1 (MQ=255)
                                 gCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                              >  1:1084617/1‑69 (MQ=255)
                                  cAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGAt                                                         >  1:4783037/1‑41 (MQ=19)
                                  cAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGAt                                                         >  1:4676068/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:1164690/51‑1 (MQ=255)
                                     tttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                                            <  1:508682/51‑1 (MQ=255)
                                     tttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                                            <  1:2752766/51‑1 (MQ=255)
                                     tttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGc                                  <  1:88627/61‑1 (MQ=255)
                                     tttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                              <  1:2808130/65‑1 (MQ=255)
                                      ttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAgaga                          >  1:4257332/1‑68 (MQ=255)
                                      ttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAgaga                          >  1:4251805/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:31162/1‑51 (MQ=255)
                                       tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                                          >  1:705350/1‑51 (MQ=255)
                                       tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                                          >  1:1282908/1‑51 (MQ=255)
                                       tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGa                                        >  1:822230/1‑53 (MQ=255)
                                       tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                              >  1:61881/1‑63 (MQ=255)
                                       tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                              >  1:3595046/1‑63 (MQ=255)
                                       tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCg                       <  1:3095358/70‑1 (MQ=255)
                                       tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCg                       <  1:3101014/70‑1 (MQ=255)
                                       tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCg                       <  1:3410738/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:2158607/71‑1 (MQ=255)
                                           tCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATg                  <  1:2876150/71‑1 (MQ=255)
                                            cGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATg                  >  1:4333780/1‑70 (MQ=255)
                                             ggtaTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGa                 <  1:4503438/67‑1 (MQ=255)
                                                 ttGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg              >  1:3350290/1‑69 (MQ=255)
                                                 ttGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg              >  1:2619281/1‑69 (MQ=255)
                                                 ttGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg              >  1:2963666/1‑69 (MQ=255)
                                                   gtaaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt          >  1:4469711/4‑71 (MQ=255)
                                                   gtaaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt          >  1:1365419/4‑71 (MQ=255)
                                                    taaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtg         <  1:1075244/69‑1 (MQ=255)
                                                    taaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt          <  1:1027049/68‑1 (MQ=255)
                                                    taaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt          <  1:679653/68‑1 (MQ=255)
                                                    taaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt          <  1:1887434/68‑1 (MQ=255)
                                                      aaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtgat       <  1:3792754/71‑1 (MQ=255)
                                                         ggtgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt          <  1:4449869/62‑1 (MQ=255)
                                                           tgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGc                               <  1:1629147/41‑1 (MQ=37)
                                                           tgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCaga                            >  1:3259276/2‑45 (MQ=37)
                                                            gAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCg  >  1:1341514/1‑70 (MQ=255)
                                                                 ||                                                               
CTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACGTAGTGATGACCG  >  W3110S.gb/3824315‑3824444

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: