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

CGGGTGCCTGCGTGGCGCTGGTAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTG  >  W3110S.gb/3824259‑3824392
                                                                   |                                                                  
cGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAAccc                                                                  <  1:1811421/70‑1 (MQ=255)
cGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAAccc                                                                  <  1:3254817/70‑1 (MQ=255)
    tGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa                                                             <  1:3921587/71‑1 (MQ=255)
    tGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa                                                             <  1:2313830/71‑1 (MQ=255)
         gCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAgg                                                         >  1:2893815/1‑70 (MQ=255)
          cGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAggg                                                        <  1:399752/70‑1 (MQ=255)
          cGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGa                                                       >  1:2343706/1‑71 (MQ=255)
              gcgcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATgg                                                    >  1:3842272/1‑70 (MQ=255)
                 cTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAgg                                                         <  1:3047220/62‑1 (MQ=255)
                 cTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAgg                                                         <  1:543844/62‑1 (MQ=255)
                 cTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAgg                                                         <  1:765549/62‑1 (MQ=255)
                      aGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGaaa                                                           >  1:3736737/1‑55 (MQ=255)
                         tGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGttt                                        >  1:5120582/1‑71 (MQ=255)
                         tGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGttt                                        >  1:3416512/1‑71 (MQ=255)
                         tGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGt                                          >  1:95136/1‑69 (MQ=255)
                          gcgcTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTc                                       >  1:2155605/1‑71 (MQ=255)
                           cgcTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCt                                      <  1:2184062/71‑1 (MQ=255)
                            gcTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa                                                             <  1:733709/47‑1 (MQ=255)
                            gcTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa                                                             <  1:1389982/47‑1 (MQ=255)
                            gcTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTg                                     <  1:4194897/71‑1 (MQ=255)
                             cTACAGAAGCTGGTTGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGAt                                                      <  1:77512/53‑1 (MQ=255)
                             cTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGAt                                                      <  1:898708/53‑1 (MQ=255)
                             cTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGAt                                                      <  1:860702/53‑1 (MQ=255)
                             cTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGAt                                                      <  1:3798643/53‑1 (MQ=255)
                             cTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGAGGAAAGGGATGGTAGCCGCAGTTTCTg                                     <  1:2566869/70‑1 (MQ=255)
                               aCAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATg                                                              <  1:4815302/43‑1 (MQ=255)
                                cAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTc                                   <  1:2988874/69‑1 (MQ=255)
                                cAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTc                                   <  1:955911/69‑1 (MQ=255)
                                    aGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGaaa                                                           >  1:1523607/1‑41 (MQ=255)
                                    aGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTATCCGCAGTTTCTGTCGGTAt                              >  1:3931491/1‑70 (MQ=255)
                                    aGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTAtt                             >  1:3415096/1‑71 (MQ=255)
                                         gtggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGa                                                       >  1:3225332/1‑40 (MQ=255)
                                         gtggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtgta                          <  1:4032490/69‑2 (MQ=255)
                                          tggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAAc                       <  1:321152/71‑1 (MQ=255)
                                           ggAAAACGGCAAGCTGAAAACCAACCAGATGAAAGGGATGGTAGCCGCAGtt                                         >  1:127877/1‑52 (MQ=255)
                                            gAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGa                                                       >  1:2463116/1‑37 (MQ=255)
                                            gAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGa                                                       >  1:553514/1‑37 (MQ=255)
                                             aaaaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggt                    <  1:92736/71‑2 (MQ=255)
                                              aaaCGGCAAGCTGAAAACCAACCCGATGAAAGGGa                                                       >  1:4004231/1‑35 (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)
                                                       gCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGt                                          <  1:2733854/39‑1 (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:4152608/71‑1 (MQ=255)
                                                               ccAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTg  <  1:2390956/71‑1 (MQ=255)
                                                               ccAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTg  <  1:5107184/71‑1 (MQ=255)
                                                               ccAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTg  <  1:4912872/71‑1 (MQ=255)
                                                                   |                                                                  
CGGGTGCCTGCGTGGCGCTGGTAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTG  >  W3110S.gb/3824259‑3824392

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: