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

CCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAGGGGACG  >  W3110S.gb/3824346‑3824476
                                                                |                                                                  
ccGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                                                              <  1:2309774/71‑1 (MQ=255)
ccGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                                                              <  1:2338975/71‑1 (MQ=255)
ccGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                                                              <  1:672901/71‑1 (MQ=255)
ccGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                                                              <  1:2209887/71‑1 (MQ=255)
ccGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                                                              <  1:2209275/71‑1 (MQ=255)
ccGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                                                              <  1:5195137/71‑1 (MQ=255)
ccGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                                                              <  1:3947362/71‑1 (MQ=255)
 cGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                                                              <  1:1402501/70‑1 (MQ=255)
      tttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGAcc                                                        <  1:4836370/71‑1 (MQ=255)
      tttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGAcc                                                        >  1:3339334/1‑71 (MQ=255)
      tttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGAcc                                                        >  1:1762656/1‑71 (MQ=255)
      tttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGAcc                                                        >  1:3203247/1‑71 (MQ=255)
      tttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGAcc                                                        >  1:2656689/1‑71 (MQ=255)
      tttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGAcc                                                        <  1:1098119/71‑1 (MQ=255)
      tttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGAc                                                         >  1:1580514/1‑70 (MQ=255)
       ttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGAcc                                                        >  1:2296878/1‑70 (MQ=255)
       ttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGAcc                                                        >  1:198331/1‑70 (MQ=255)
       ttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGAcc                                                        >  1:3756879/1‑70 (MQ=255)
       ttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGAcc                                                        >  1:4370683/1‑70 (MQ=255)
         cTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa                                                      <  1:2274009/70‑1 (MQ=255)
          tGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa                                                      <  1:3501349/69‑1 (MQ=255)
                  ttGTAAACTGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                                                              <  1:3701528/53‑1 (MQ=255)
                  ttGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCg                                                       <  1:4818327/60‑1 (MQ=255)
                  ttGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCg                                                       <  1:5491904/60‑1 (MQ=255)
                  ttGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCg                                                       <  1:5587661/60‑1 (MQ=255)
                  ttGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCg                                                       <  1:1012350/60‑1 (MQ=255)
                  ttGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGAAGACCGACATGa                                                 <  1:69279/64‑1 (MQ=255)
                         cGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtgatga                                     <  1:930740/71‑1 (MQ=255)
                            tgAAGCAATTTGCGATCTGGAATACGTTGAAGTCTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCg                                  <  1:4579464/70‑1 (MQ=255)
                             gAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGTCGCAGAGACCGACATGAACGTAGtgat                                       <  1:2387505/65‑1 (MQ=255)
                             gAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtgat                                       <  1:4385124/65‑1 (MQ=255)
                             gAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtgat                                       <  1:5271530/65‑1 (MQ=255)
                             gAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtgat                                       <  1:3855815/65‑1 (MQ=255)
                             gAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtgat                                       <  1:1984699/65‑1 (MQ=255)
                             gAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtgat                                       <  1:336934/65‑1 (MQ=255)
                                    tttGCGATCTGGAATACGTTGAAGACTCAGCCGCaga                                                            <  1:4132763/37‑1 (MQ=255)
                                    tttGCGATCTGGAATACGTTGAAGACTCAGCCGCaga                                                            <  1:5413073/37‑1 (MQ=255)
                                    tttGCGATCTGGAATACGTTGAAGACTCAGCCGCaga                                                            <  1:2784929/37‑1 (MQ=255)
                                    tttGCGATCTGGAATACGTTGAAGACTCAGCCGCaga                                                            <  1:1834478/37‑1 (MQ=255)
                                    tttGCGATCTGGAATACGTTGAAGACTCAGCCGCaga                                                            <  1:1069731/37‑1 (MQ=255)
                                    tttGCGATCTGGAATACGTTGAAGACTCAGCCGCaga                                                            <  1:4666033/37‑1 (MQ=255)
                                      tGCGATCTGGAATACGTTGAAGACTCAGCCGCAgaga                                                          >  1:581370/1‑37 (MQ=255)
                                      tGCGATCTGGAATACGTTGAAGACTCAGCCGCAgaga                                                          >  1:660204/1‑37 (MQ=255)
                                      tGCGATCTGGAATACGTTGAAGACTCAGCCGCAgaga                                                          >  1:4687222/1‑37 (MQ=255)
                                      tGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGgcg                        >  1:2834898/1‑71 (MQ=255)
                                      tGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGgc                         >  1:4947766/1‑70 (MQ=255)
                                        cGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGgcgc                       <  1:1285637/70‑1 (MQ=255)
                                        cGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGgcgc                       <  1:1604287/70‑1 (MQ=255)
                                        cGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGgcgc                       <  1:5450524/70‑1 (MQ=255)
                                        cGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGgcgc                       <  1:5522814/70‑1 (MQ=255)
                                           tCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGt                                             >  1:2891884/1‑45 (MQ=255)
                                           tCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGcatca                   >  1:5286127/1‑71 (MQ=255)
                                           tCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGcatca                   >  1:1957113/1‑71 (MQ=255)
                                           tCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGcatca                   >  1:27681/1‑71 (MQ=255)
                                           tCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGcatca                   >  1:3992839/1‑71 (MQ=255)
                                           tCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGcatca                   >  1:2466437/1‑71 (MQ=255)
                                               gAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTg                <  1:1329709/70‑1 (MQ=255)
                                               gAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTg                <  1:4955333/70‑1 (MQ=255)
                                               gAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTg                <  1:3079835/70‑1 (MQ=255)
                                               gAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTg                <  1:5365841/70‑1 (MQ=255)
                                               gAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTg                <  1:2286519/70‑1 (MQ=255)
                                                    cGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtgatga                                     >  1:5111900/1‑44 (MQ=255)
                                                    cGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGt            >  1:3030550/1‑69 (MQ=255)
                                                    cGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTg           <  1:5032830/70‑1 (MQ=255)
                                                    cGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTg           >  1:5055302/1‑70 (MQ=255)
                                                    cGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTg           >  1:1004961/1‑70 (MQ=255)
                                                    cGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTg           >  1:1040208/1‑70 (MQ=255)
                                                    cGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTg           >  1:785811/1‑70 (MQ=255)
                                                    cGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTg           >  1:2289368/1‑70 (MQ=255)
                                                    cGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTg           >  1:2936597/1‑70 (MQ=255)
                                                    cGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTg           <  1:2277894/70‑1 (MQ=255)
                                                    cGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTg           >  1:2409671/1‑70 (MQ=255)
                                                    cGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTg           <  1:2138712/70‑1 (MQ=255)
                                                    cGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTg           <  1:276786/70‑1 (MQ=255)
                                                    cGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTg           >  1:3010145/1‑70 (MQ=255)
                                                    cGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTg           <  1:3225065/70‑1 (MQ=255)
                                                    cGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTg           >  1:4116911/1‑70 (MQ=255)
                                                     gTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCTTCATTGAAGTGCa         <  1:2134734/71‑1 (MQ=255)
                                                       tGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAgg       >  1:1778764/1‑71 (MQ=255)
                                                       tGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAgg       >  1:1441299/1‑71 (MQ=255)
                                                       tGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAgg       >  1:2738999/1‑71 (MQ=255)
                                                       tGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAgg       >  1:421926/1‑71 (MQ=255)
                                                       tGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAgg       >  1:289900/1‑71 (MQ=255)
                                                        gAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAggg      >  1:1278118/1‑71 (MQ=255)
                                                         aaGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAggg      <  1:4221147/70‑1 (MQ=255)
                                                         aaGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAggg      <  1:4051100/70‑1 (MQ=255)
                                                         aaGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAggg      <  1:4220853/70‑1 (MQ=255)
                                                         aaGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAggg      <  1:2413383/70‑1 (MQ=255)
                                                         aaGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAggg      <  1:5436506/70‑1 (MQ=255)
                                                         aaGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAggg      <  1:2908094/70‑1 (MQ=255)
                                                         aaGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAggg      <  1:4922806/70‑1 (MQ=255)
                                                           gACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGa                              >  1:1021353/1‑44 (MQ=255)
                                                           gACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGa                              >  1:3959740/1‑44 (MQ=255)
                                                           gACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTg           <  1:1464480/63‑1 (MQ=255)
                                                            aCTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTg           <  1:2815278/62‑1 (MQ=255)
                                                            aCTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAGGGGACg  <  1:57015/71‑1 (MQ=255)
                                                             ctcaGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAGGGGACg  <  1:3622862/70‑1 (MQ=255)
                                                                |                                                                  
CCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAGGGGACG  >  W3110S.gb/3824346‑3824476

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: