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

GTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAGGGGACGG  >  W3110S.gb/3824342‑3824477
                                                                    |                                                                   
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)
   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)
        aGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAgaga                                                           >  1:1288725/1‑71 (MQ=255)
          tttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                                                               <  1:1046803/65‑1 (MQ=255)
          tttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                                                               <  1:660060/65‑1 (MQ=255)
            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:3677195/1‑48 (MQ=37)
                             cGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCaga                                                             >  1:3773098/1‑48 (MQ=37)
                               gtgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGc                                                                <  1:2791894/41‑1 (MQ=25)
                                tgAAGCAATTTGCGATCTGGAATACTTTGAAGACTCAGCCg                                                                 <  1:2120230/40‑1 (MQ=21)
                                       atttGCGATCTGGAATACGTTGAAGACTCATCCGCa                                                               <  1:2218332/35‑1 (MQ=37)
                                       atttGCGATCTGGAATACGTTGAAGACTCAGCCGCa                                                               <  1:1695353/35‑1 (MQ=38)
                                         ttGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGa                                                  >  1:1206707/1‑47 (MQ=255)
                                         ttGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGgc                          >  1:2940418/1‑71 (MQ=255)
                                          tGCGATCTGGAATACGTTGAAGACTCATCCGCAGAGACCGa                                                       >  1:3835102/1‑41 (MQ=255)
                                          tGCGATCTGGAATACGTTGAAGACTCAGCCGCaga                                                             >  1:1215656/1‑35 (MQ=255)
                                          tGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa                                                       >  1:1055575/1‑41 (MQ=255)
                                          tGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa                                                       >  1:1604701/1‑41 (MQ=255)
                                          tGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGAcc                                    >  1:1308395/1‑60 (MQ=255)
                                          tGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGAcc                                    >  1:4156377/1‑60 (MQ=255)
                                          tGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGgcg                         >  1:993577/1‑71 (MQ=255)
                                          tGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGgcg                         >  1:81460/1‑71 (MQ=255)
                                          tGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGgc                          <  1:2618532/70‑1 (MQ=255)
                                              aTCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGa                               >  1:1804375/1‑61 (MQ=255)
                                                cTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGcatcat                   >  1:2587670/1‑71 (MQ=255)
                                                cTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGcatcat                   >  1:2273313/1‑71 (MQ=255)
                                                 tGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtgatg                                       <  1:2911980/50‑1 (MQ=255)
                                                 tGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                                           <  1:3758498/46‑1 (MQ=255)
                                                 tGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCAtt                  >  1:4281149/1‑71 (MQ=255)
                                                 tGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCAtt                  >  1:2608204/1‑71 (MQ=255)
                                                   gAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGcatca                    <  1:3108948/67‑1 (MQ=255)
                                                   gAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTg                 <  1:3994253/70‑1 (MQ=255)
                                                   gAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTg                 >  1:1096093/1‑70 (MQ=255)
                                                   gAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTg                 <  1:859356/70‑1 (MQ=255)
                                                   gAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTg                 <  1:3421862/70‑1 (MQ=255)
                                                   gAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGAGCATCATTg                 <  1:570732/70‑1 (MQ=255)
                                                       aCGTTGAAGACTCAGCCGCAGAGACCGACATGAACg                                               <  1:1443354/36‑1 (MQ=255)
                                                        cGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGcatc                     >  1:2516912/1‑61 (MQ=255)
                                                        cGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTg            >  1:2750721/1‑70 (MQ=255)
                                                        cGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGc           <  1:3391473/71‑1 (MQ=255)
                                                        cGTAGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGa                               >  1:2424677/1‑51 (MQ=255)
                                                         gTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGc           <  1:123634/70‑1 (MQ=255)
                                                          ttGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGa                               <  1:4287188/49‑1 (MQ=255)
                                                          ttGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAg         <  1:1559376/71‑1 (MQ=255)
                                                           tGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtga                                         >  1:2709826/1‑38 (MQ=255)
                                                            gAAGACTCAGCCGCAGAGACCGACATGAACGTAGtga                                         >  1:2562774/1‑37 (MQ=255)
                                                                aCTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGGATCATTGAAGTGCAGGGGACg   <  1:2355888/71‑1 (MQ=255)
                                                                 ctcaGCCGCAGAGACCGTCATGAACGTAGTGATGACCg                                   >  1:3235794/1‑38 (MQ=38)
                                                                 ctcaGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAg         >  1:658493/1‑64 (MQ=255)
                                                                 ctcaGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAg         >  1:1903645/1‑64 (MQ=255)
                                                                 ctcaGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAGGGGACgg  <  1:3007404/71‑1 (MQ=255)
                                                                    |                                                                   
GTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAGGGGACGG  >  W3110S.gb/3824342‑3824477

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: