Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,824,380 2 bp→AA 100% pseudogene (513‑514/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,3800GA100.0% 193.6 / NA 69pseudogene (513/687 nt)rphECK3633:JW3618:b3643; ribonuclease PH
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (37/32);  total (37/32)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
*W3110S.gb3,824,3810GA100.0% 195.7 / NA 69pseudogene (514/687 nt)rphECK3633:JW3618:b3643; ribonuclease PH
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (37/32);  total (37/32)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACGTAGTGATG  >  W3110S.gb/3824317‑3824440
                                                               ||                                                           
gAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGc                                                        >  1:232767/1‑70 (MQ=255)
gAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGc                                                        >  1:3694425/1‑70 (MQ=255)
 aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGa                                                      >  1:2329814/1‑71 (MQ=255)
 aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGa                                                      <  1:1813908/71‑1 (MQ=255)
 aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGa                                                      <  1:2884486/71‑1 (MQ=255)
     ccAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTg                                                  <  1:202790/71‑1 (MQ=255)
      cAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTgg                                                 <  1:1846119/71‑1 (MQ=255)
      cAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTgg                                                 <  1:2392574/71‑1 (MQ=255)
      cAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTgg                                                 <  1:3809293/71‑1 (MQ=255)
      cAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTg                                                  >  1:1456726/1‑70 (MQ=255)
         cccGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGaa                                               >  1:2823998/2‑70 (MQ=255)
           cGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGc                                                        >  1:2842749/1‑59 (MQ=255)
            gATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACg                                           <  1:4138096/71‑1 (MQ=255)
               gAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATa                                             >  1:3574559/1‑66 (MQ=255)
               gAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATa                                             >  1:4186287/1‑66 (MQ=255)
               gAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGtt                                         >  1:195445/1‑70 (MQ=255)
                aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                                        >  1:4378864/1‑70 (MQ=255)
                aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                                        <  1:2759955/70‑1 (MQ=255)
                aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                                        <  1:3877002/70‑1 (MQ=255)
                aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                                        <  1:2387139/70‑1 (MQ=255)
                aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                                        <  1:1080610/70‑1 (MQ=255)
                 aaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                                      <  1:3541650/71‑1 (MQ=255)
                      aTGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGActc                                 <  1:4418413/71‑1 (MQ=255)
                        ggTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGAct                                  >  1:1583316/1‑68 (MQ=255)
                         gTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGc                              >  1:3423536/1‑71 (MQ=255)
                         gTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGc                              >  1:1650945/1‑71 (MQ=255)
                          tAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGc                              >  1:1458063/1‑70 (MQ=255)
                           aGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCg                            >  1:2958070/1‑71 (MQ=255)
                            gCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTgg                                                 <  1:4310905/49‑1 (MQ=25)
                            gCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCg                            >  1:3188825/1‑70 (MQ=255)
                            gCCGCAGTTTCTGTCGGTATTGTAAACCGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCg                            <  1:3261125/70‑1 (MQ=255)
                              cGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCg                            >  1:2442978/1‑68 (MQ=255)
                               gCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                          <  1:3875356/69‑1 (MQ=255)
                                cAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                          >  1:2427079/1‑68 (MQ=255)
                                 aGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTg                                                         >  1:49513/1‑36 (MQ=16)
                                 aGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGc                                                        <  1:487192/37‑1 (MQ=16)
                                 aGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGc                                                        <  1:1146740/37‑1 (MQ=16)
                                 aGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAgaga                      >  1:1288725/1‑71 (MQ=255)
                                   tttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTgg                                                 <  1:1899496/42‑1 (MQ=19)
                                   tttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATa                                             <  1:3957182/46‑1 (MQ=25)
                                   tttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATa                                             <  1:695985/46‑1 (MQ=25)
                                   tttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                          <  1:1046803/65‑1 (MQ=255)
                                   tttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                          <  1:660060/65‑1 (MQ=255)
                                    ttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGActcag                               >  1:474601/1‑57 (MQ=255)
                                    ttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGActcag                               >  1:452878/1‑57 (MQ=255)
                                     tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGt                                          >  1:225996/1‑47 (MQ=25)
                                     tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCaga                        >  1:269608/1‑65 (MQ=255)
                                      ctgtctgtATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa                  <  1:1316784/70‑1 (MQ=255)
                                      cTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa                  <  1:3379863/70‑1 (MQ=255)
                                           ggtaTTGTAAACGGTAAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGa             <  1:456305/67‑1 (MQ=255)
                                             taTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg          <  1:3738145/70‑1 (MQ=255)
                                               ttGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGt         <  1:4272502/70‑1 (MQ=255)
                                                 gtaaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGATTCAGCCGCAGAGACCGACATGAACGTAGt      >  1:319440/4‑71 (MQ=255)
                                                 gtaaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTa        >  1:3088690/4‑69 (MQ=255)
                                                 gtaaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTa        >  1:3093857/4‑69 (MQ=255)
                                                 gtaaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAg       >  1:1960744/4‑70 (MQ=255)
                                                 gtaaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAg       <  1:735082/67‑1 (MQ=255)
                                                 gtaaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt      >  1:2800096/4‑71 (MQ=255)
                                                  taaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa                  >  1:3023341/3‑58 (MQ=255)
                                                  taaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa                  >  1:3018234/3‑58 (MQ=255)
                                                   aaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtga    >  1:1856261/2‑71 (MQ=255)
                                                    aaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTa        >  1:4188759/1‑66 (MQ=255)
                                                    aaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtgat   >  1:364552/1‑71 (MQ=255)
                                                     aCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtgatg  <  1:2866903/71‑1 (MQ=255)
                                                     aCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtga    >  1:7743/1‑69 (MQ=255)
                                                      cGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCaga                        >  1:3773098/1‑48 (MQ=37)
                                                      cGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCaga                        >  1:3677195/1‑48 (MQ=37)
                                                        gtgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGc                           <  1:2791894/41‑1 (MQ=25)
                                                         tgAAGCAATTTGCGATCTGGAATACTTTGAAGACTCAGCCg                            <  1:2120230/40‑1 (MQ=21)
                                                               ||                                                           
GAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACGTAGTGATG  >  W3110S.gb/3824317‑3824440

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