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

TGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGA  >  W3110S.gb/3824301‑3824429
                                                                   |                                                             
tggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAAc                                                            <  1:4527159/71‑1 (MQ=255)
tggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAAc                                                            <  1:3213000/71‑1 (MQ=255)
    aaaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGTTATTGTAAACggtg                                                        <  1:2983537/71‑3 (MQ=255)
    aaaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggtg                                                        <  1:3241473/71‑3 (MQ=255)
                 aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCg                                            <  1:5439552/70‑1 (MQ=255)
                 aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCg                                            <  1:1199641/70‑1 (MQ=255)
                 aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCg                                            <  1:3109191/70‑1 (MQ=255)
                 aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCg                                            <  1:1668572/70‑1 (MQ=255)
                 aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCg                                            <  1:2756239/70‑1 (MQ=255)
                 aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCg                                            <  1:1852075/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:3929690/66‑1 (MQ=255)
                    aCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGc                                             <  1:3892492/66‑1 (MQ=255)
                      cAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTgg                                      <  1:3035787/71‑1 (MQ=255)
                           cGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATa                                  >  1:1388672/1‑70 (MQ=255)
                           cGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATa                                  >  1:2595266/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:1140904/1‑71 (MQ=255)
                             aTGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGt                               >  1:1363834/1‑71 (MQ=255)
                             aTGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGt                               >  1:486208/1‑71 (MQ=255)
                             aTGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGt                               >  1:3348652/1‑71 (MQ=255)
                             aTGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGt                               >  1:3144794/1‑71 (MQ=255)
                               gAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGtt                              >  1:1932099/1‑70 (MQ=255)
                                aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                             <  1:4608954/70‑1 (MQ=255)
                                aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                             <  1:4698148/70‑1 (MQ=255)
                                aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                             >  1:3761830/1‑70 (MQ=255)
                                aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                             <  1:4979676/70‑1 (MQ=255)
                                aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                             <  1:5279873/70‑1 (MQ=255)
                                aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                             >  1:936637/1‑70 (MQ=255)
                                aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                             <  1:2544988/70‑1 (MQ=255)
                                aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                             <  1:2598845/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:672901/71‑1 (MQ=255)
                                             ccGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa               <  1:2209275/71‑1 (MQ=255)
                                             ccGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa               <  1:3947362/71‑1 (MQ=255)
                                             ccGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa               <  1:2209887/71‑1 (MQ=255)
                                             ccGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa               <  1:2309774/71‑1 (MQ=255)
                                             ccGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa               <  1:2338975/71‑1 (MQ=255)
                                             ccGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa               <  1:5195137/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:605332/1‑54 (MQ=255)
                                                  gTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                           >  1:3794095/1‑54 (MQ=255)
                                                  gTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                           >  1:5389457/1‑54 (MQ=255)
                                                  gTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                           >  1:4603232/1‑54 (MQ=255)
                                                   tttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGAcc         >  1:1762656/1‑71 (MQ=255)
                                                   tttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGAcc         >  1:2656689/1‑71 (MQ=255)
                                                   tttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGAcc         <  1:4836370/71‑1 (MQ=255)
                                                   tttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGAcc         >  1:3203247/1‑71 (MQ=255)
                                                   tttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGAcc         >  1:3339334/1‑71 (MQ=255)
                                                   tttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGAcc         <  1:1098119/71‑1 (MQ=255)
                                                   tttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGAc          >  1:1580514/1‑70 (MQ=255)
                                                    ttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGAcc         >  1:3756879/1‑70 (MQ=255)
                                                    ttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGAcc         >  1:2296878/1‑70 (MQ=255)
                                                    ttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGAcc         >  1:198331/1‑70 (MQ=255)
                                                    ttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGAcc         >  1:4370683/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:4818327/60‑1 (MQ=255)
                                                               ttGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCg        <  1:1012350/60‑1 (MQ=255)
                                                               ttGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCg        <  1:5491904/60‑1 (MQ=255)
                                                               ttGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCg        <  1:5587661/60‑1 (MQ=255)
                                                               ttGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGAAGACCGACATGa  <  1:69279/64‑1 (MQ=255)
                                                                   |                                                             
TGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGA  >  W3110S.gb/3824301‑3824429

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: