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% 172.8 / NA 58pseudogene (495/687 nt)rphECK3633:JW3618:b3643; ribonuclease PH
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (28/30);  total (28/30)
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
                                                                   |                                                                  
aGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTATCCGCAGTTTCTGTCGGTAt                                                                  >  1:3931491/1‑70 (MQ=255)
aGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTAtt                                                                 >  1:3415096/1‑71 (MQ=255)
     gtggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtgta                                                              <  1:4032490/69‑2 (MQ=255)
      tggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAAc                                                           <  1:321152/71‑1 (MQ=255)
         aaaaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggt                                                        <  1:92736/71‑2 (MQ=255)
          aaaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggtg                                                       <  1:4492691/71‑3 (MQ=255)
           aaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggt                                                        <  1:1106145/69‑2 (MQ=255)
           aaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGa                                                      <  1:3796428/71‑1 (MQ=255)
             cGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGa                                                      >  1:3934646/1‑69 (MQ=255)
              ggCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAg                                                    >  1:2153373/1‑70 (MQ=255)
               gCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGca                                                  <  1:4512878/71‑2 (MQ=255)
                cAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGca                                                  <  1:2789782/70‑2 (MQ=255)
                    cTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTg                                             >  1:2994338/1‑71 (MQ=255)
                       aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGa                                          >  1:2824550/1‑71 (MQ=255)
                         aaccaaccCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATc                                        >  1:562265/1‑71 (MQ=255)
                          aCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCg                                           <  1:2761375/67‑1 (MQ=255)
                          aCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCt                                       >  1:3451634/1‑71 (MQ=255)
                           ccAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggt                                                        <  1:24226/53‑2 (MQ=255)
                           ccAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCt                                       <  1:578141/70‑1 (MQ=255)
                           ccAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTg                                      <  1:5107184/71‑1 (MQ=255)
                           ccAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTg                                      <  1:4912872/71‑1 (MQ=255)
                           ccAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTg                                      <  1:2390956/71‑1 (MQ=255)
                           ccAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTg                                      <  1:4152608/71‑1 (MQ=255)
                              acccGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGaa                                   <  1:1333910/69‑1 (MQ=255)
                                ccGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACTGTGAAGCAATTTGCGATCTGGAATa                                 >  1:4762662/1‑71 (MQ=255)
                                ccGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATa                                 >  1:173370/1‑71 (MQ=255)
                                 cGATGAAAGGGATGGTAGCCGCAGTTTCCGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATAc                                <  1:4271093/71‑1 (MQ=255)
                                  gATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACg                               >  1:90940/1‑71 (MQ=255)
                                   aTGAAAGTGATGGTAGCCGCAGTTTCTGTCGGTATTGTaaa                                                            >  1:5033687/1‑41 (MQ=37)
                                   aTGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTaaa                                                            >  1:415179/1‑41 (MQ=38)
                                   aTGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGca                                                  <  1:1937946/51‑2 (MQ=255)
                                     gAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATa                                 >  1:4400866/1‑66 (MQ=255)
                                     gAAAGGGATGGTAGACGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATa                                 >  1:104257/1‑66 (MQ=255)
                                      aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                            <  1:241735/70‑1 (MQ=255)
                                       aaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                            >  1:2291937/1‑69 (MQ=255)
                                       aaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                          <  1:2220482/71‑1 (MQ=255)
                                       aaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                          <  1:634891/71‑1 (MQ=255)
                                        aGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGca                                                  <  1:4760620/46‑2 (MQ=255)
                                         gggATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGa                        <  1:2845717/71‑1 (MQ=255)
                                          ggATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATAc                                <  1:3977246/62‑1 (MQ=255)
                                            aTGGTAGCCGCAGTTTCTGTCGGTATTGTAAACgg                                                         >  1:2658063/1‑35 (MQ=37)
                                               gTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGtt                             >  1:4703860/1‑60 (MQ=255)
                                               gTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                            <  1:4977454/61‑1 (MQ=255)
                                               gTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                            <  1:3081076/61‑1 (MQ=255)
                                               gTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGActcag                   <  1:509473/70‑3 (MQ=255)
                                                 aGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCg                >  1:2543418/1‑71 (MQ=255)
                                                  gCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGa                                          >  1:3760291/1‑44 (MQ=21)
                                                  gCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGa                                          >  1:4224767/1‑44 (MQ=21)
                                                  gCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGa                                          >  1:3796225/1‑44 (MQ=21)
                                                          ttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGtt                             <  1:3987026/49‑1 (MQ=25)
                                                          ttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGtt                             <  1:2549385/49‑1 (MQ=25)
                                                           tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCg       >  1:3131981/1‑70 (MQ=255)
                                                            cTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa      >  1:4881398/1‑70 (MQ=255)
                                                            cTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa      >  1:4875778/1‑70 (MQ=255)
                                                            cTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa      >  1:4284490/1‑70 (MQ=255)
                                                            cTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa      <  1:2047520/70‑1 (MQ=255)
                                                                cGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATg  >  1:3211571/1‑70 (MQ=255)
                                                                cGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATg  >  1:3517274/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: