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

ACGGGTGCCTGCGTGGCGCTGGTAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGG  >  W3110S.gb/3824258‑3824393
                                                                    |                                                                   
aCGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAAccc                                                                   >  1:4316506/1‑71 (MQ=255)
aCGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAAccc                                                                   >  1:1457664/1‑71 (MQ=255)
aCGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAAccc                                                                   >  1:4475720/1‑71 (MQ=255)
     tGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa                                                              <  1:3633971/71‑1 (MQ=255)
     tGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa                                                              <  1:1322077/71‑1 (MQ=255)
     tGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa                                                              <  1:5544971/71‑1 (MQ=255)
     tGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa                                                              <  1:5390511/71‑1 (MQ=255)
      gCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa                                                              >  1:4423688/1‑70 (MQ=255)
      gCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa                                                              <  1:2946350/70‑1 (MQ=255)
            gTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGAt                                                       <  1:3780255/71‑1 (MQ=255)
                       aGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAg                                            >  1:3928452/1‑71 (MQ=255)
                       aGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAg                                            >  1:231778/1‑71 (MQ=255)
                       aGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAg                                            >  1:3972682/1‑71 (MQ=255)
                                 cAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCACCCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTc                                    <  1:2203000/69‑1 (MQ=255)
                                 cAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTc                                    <  1:5312511/69‑1 (MQ=255)
                                 cAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTc                                    <  1:797574/69‑1 (MQ=255)
                                 cAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTc                                    <  1:4144762/69‑1 (MQ=255)
                                   gAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAg                                                           <  1:1188227/44‑1 (MQ=255)
                                     aGCTGGTTGAAAACTGCAAGCTGAAAACCAACCCGATGAAAg                                                           <  1:4466118/42‑1 (MQ=37)
                                       cTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa                                                              <  1:4086378/37‑1 (MQ=255)
                                           tggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAAc                        <  1:4527159/71‑1 (MQ=255)
                                           tggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAAc                        <  1:3213000/71‑1 (MQ=255)
                                             gAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGcc                                                >  1:4566214/1‑45 (MQ=255)
                                             gAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGcc                                                >  1:5383577/1‑45 (MQ=255)
                                             gAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGcc                                                >  1:944690/1‑45 (MQ=255)
                                               aaaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGTTATTGTAAACggtg                    <  1:2983537/71‑3 (MQ=255)
                                               aaaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggtg                    <  1:3241473/71‑3 (MQ=255)
                                                         cTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTc                                        <  1:859790/41‑1 (MQ=255)
                                                         cTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTc                                        <  1:1762030/41‑1 (MQ=255)
                                                         cTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTc                                        <  1:1942744/41‑1 (MQ=255)
                                                         cTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTc                                        <  1:2460073/41‑1 (MQ=255)
                                                         cTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTc                                        <  1:435331/41‑1 (MQ=255)
                                                          tGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTc                                        <  1:4557147/40‑1 (MQ=255)
                                                            aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGtt                                          <  1:5130720/36‑1 (MQ=255)
                                                            aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGtt                                          <  1:163107/36‑1 (MQ=255)
                                                            aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGtt                                          <  1:2465436/36‑1 (MQ=255)
                                                            aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGtt                                          <  1:4265235/36‑1 (MQ=255)
                                                            aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGtt                                          <  1:1567793/36‑1 (MQ=255)
                                                            aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGtt                                          <  1:1563507/36‑1 (MQ=255)
                                                            aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCt                                       <  1:5383897/39‑1 (MQ=255)
                                                            aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCt                                       <  1:3077526/39‑1 (MQ=255)
                                                            aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCg        <  1:5439552/70‑1 (MQ=255)
                                                            aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCg        <  1:1668572/70‑1 (MQ=255)
                                                            aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCg        <  1:1199641/70‑1 (MQ=255)
                                                            aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCg        <  1:3109191/70‑1 (MQ=255)
                                                            aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCg        <  1:1852075/70‑1 (MQ=255)
                                                            aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCg        <  1:2756239/70‑1 (MQ=255)
                                                             aaaCCAACCCGATGAAAGGGATGGTAGCCGCAGtt                                          <  1:4176591/35‑1 (MQ=255)
                                                             aaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCg        <  1:4127087/69‑1 (MQ=255)
                                                             aaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGGCGGTATTGTAAACGGTGAAGCAATTTGCg        <  1:195709/69‑1 (MQ=255)
                                                               aCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGc         <  1:5028288/66‑1 (MQ=255)
                                                               aCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGc         <  1:3929690/66‑1 (MQ=255)
                                                               aCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGc         <  1:3892492/66‑1 (MQ=255)
                                                               aCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGc         <  1:1014694/66‑1 (MQ=255)
                                                                 cAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTgg  <  1:3035787/71‑1 (MQ=255)
                                                                    |                                                                   
ACGGGTGCCTGCGTGGCGCTGGTAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGG  >  W3110S.gb/3824258‑3824393

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: