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

AAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACGTAGTGATGA  >  W3110S.gb/3824318‑3824441
                                                        |                                                                   
aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCg                                                        <  1:2756239/70‑1 (MQ=255)
aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCg                                                        <  1:5439552/70‑1 (MQ=255)
aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCg                                                        <  1:3109191/70‑1 (MQ=255)
aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCg                                                        <  1:1199641/70‑1 (MQ=255)
aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCg                                                        <  1:1852075/70‑1 (MQ=255)
aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCg                                                        <  1:1668572/70‑1 (MQ=255)
 aaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCg                                                        <  1:4127087/69‑1 (MQ=255)
 aaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGGCGGTATTGTAAACGGTGAAGCAATTTGCg                                                        <  1:195709/69‑1 (MQ=255)
   aCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGc                                                         <  1:1014694/66‑1 (MQ=255)
   aCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGc                                                         <  1:5028288/66‑1 (MQ=255)
   aCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGc                                                         <  1:3892492/66‑1 (MQ=255)
   aCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGc                                                         <  1:3929690/66‑1 (MQ=255)
     cAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTgg                                                  <  1:3035787/71‑1 (MQ=255)
          cGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATa                                              >  1:2595266/1‑70 (MQ=255)
          cGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATa                                              >  1:1388672/1‑70 (MQ=255)
          agATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATa                                              >  1:2788512/2‑70 (MQ=255)
            aTGATAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGt                                           >  1:4059505/1‑71 (MQ=255)
            aTGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACg                                            >  1:5228621/1‑70 (MQ=255)
            aTGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGt                                           >  1:1363834/1‑71 (MQ=255)
            aTGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGt                                           >  1:3348652/1‑71 (MQ=255)
            aTGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGt                                           >  1:3144794/1‑71 (MQ=255)
            aTGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGt                                           >  1:1140904/1‑71 (MQ=255)
            aTGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGt                                           >  1:486208/1‑71 (MQ=255)
              gAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGtt                                          >  1:1932099/1‑70 (MQ=255)
               aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                                         >  1:936637/1‑70 (MQ=255)
               aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                                         <  1:4608954/70‑1 (MQ=255)
               aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                                         >  1:3761830/1‑70 (MQ=255)
               aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                                         <  1:4698148/70‑1 (MQ=255)
               aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                                         <  1:4979676/70‑1 (MQ=255)
               aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                                         <  1:5279873/70‑1 (MQ=255)
               aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                                         <  1:2598845/70‑1 (MQ=255)
               aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                                         <  1:2544988/70‑1 (MQ=255)
                  gggATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTgg                                                  <  1:1454512/58‑1 (MQ=255)
                    gATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGAct                                   <  1:4493781/71‑1 (MQ=255)
                    gATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGAct                                   <  1:4868930/71‑1 (MQ=255)
                    gATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGAct                                   <  1:1786200/71‑1 (MQ=255)
                     aTGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTTAAGCAATTTGCGATCTGGAATACGTTGAAGActc                                  <  1:1438562/71‑1 (MQ=255)
                     aTGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGActc                                  <  1:1450717/71‑1 (MQ=255)
                      tGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGActca                                 <  1:777687/71‑2 (MQ=255)
                      tGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGActca                                 <  1:3788150/71‑2 (MQ=255)
                            ccGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                           <  1:5195137/71‑1 (MQ=255)
                            ccGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                           <  1:3947362/71‑1 (MQ=255)
                            ccGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                           <  1:2209275/71‑1 (MQ=255)
                            ccGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                           <  1:2209887/71‑1 (MQ=255)
                            ccGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                           <  1:672901/71‑1 (MQ=255)
                            ccGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                           <  1:2309774/71‑1 (MQ=255)
                            ccGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                           <  1:2338975/71‑1 (MQ=255)
                             cGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                           <  1:1402501/70‑1 (MQ=255)
                               cAGTTTCTGTGGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                                         <  1:5050411/54‑1 (MQ=255)
                               cAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                                         <  1:1719573/54‑1 (MQ=255)
                                 gTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                                       >  1:4603232/1‑54 (MQ=255)
                                 gTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                                       >  1:5389457/1‑54 (MQ=255)
                                 gTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                                       >  1:605332/1‑54 (MQ=255)
                                 gTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                                       >  1:3794095/1‑54 (MQ=255)
                                  tttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGAcc                     <  1:4836370/71‑1 (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)
                                  tttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGAcc                     >  1:3339334/1‑71 (MQ=255)
                                  tttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGAc                      >  1:1580514/1‑70 (MQ=255)
                                   ttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGAcc                     >  1:2296878/1‑70 (MQ=255)
                                   ttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGAcc                     >  1:4370683/1‑70 (MQ=255)
                                   ttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGAcc                     >  1:198331/1‑70 (MQ=255)
                                   ttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGAcc                     >  1:3756879/1‑70 (MQ=255)
                                    tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGaa                                                >  1:1690742/1‑42 (MQ=19)
                                     cTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa                   <  1:2274009/70‑1 (MQ=255)
                                      tGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa                   <  1:3501349/69‑1 (MQ=255)
                                              ttGTAAACTGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                           <  1:3701528/53‑1 (MQ=255)
                                              ttGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCg                    <  1:5587661/60‑1 (MQ=255)
                                              ttGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCg                    <  1:5491904/60‑1 (MQ=255)
                                              ttGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCg                    <  1:4818327/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)
                                                        |                                                                   
AAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACGTAGTGATGA  >  W3110S.gb/3824318‑3824441

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: