Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,992,923 C→T 25.6% pseudogene (227/375 nt) pbl → ECK2852:JW5457:b2854; predicted peptidoglycan‑binding enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,992,9230CT25.6% 157.8 / 59.1 90pseudogene (227/375 nt)pblECK2852:JW5457:b2854; predicted peptidoglycan‑binding enzyme
Reads supporting (aligned to +/- strand):  ref base C (19/47);  new base T (9/14);  total (29/61)
Fisher's exact test for biased strand distribution p-value = 4.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.22e-01

AAACGTGGAATAATAAAAGATGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGG  >  W3110S.gb/2992856‑2992992
                                                                   |                                                                     
aaaCGTGGAATAATAAAAGATGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCga                                                                    <  1:2639683/71‑1 (MQ=255)
     tgtaataatAAAAGATGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGAttt                                                               <  1:3481345/68‑1 (MQ=255)
     tGGAATAATAAAAGATGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGAttt                                                               >  1:859336/1‑71 (MQ=255)
      ggAATAATAAAAGATGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGAttt                                                               >  1:3903994/1‑70 (MQ=255)
       gAATAATAAAAGATGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAATCGAGATTTta                                                             <  1:2628538/71‑1 (MQ=255)
            ataaAAGATGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAATCGAGATTTTATATaa                                                        >  1:4018937/1‑71 (MQ=255)
            ataaAAGATGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATaa                                                        >  1:633640/1‑71 (MQ=255)
              aaaaGATGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAAc                                                       >  1:3376514/1‑70 (MQ=255)
              aaaaGATGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAAc                                                       >  1:3666247/1‑70 (MQ=255)
              aaaaGATGAACGAGATTTGTTAGATAACCCTTGGCTGAAGATTTAAATAGGAATCgag                                                                   <  1:1257198/58‑1 (MQ=255)
                aaGATGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCAt                                                    <  1:3212331/71‑1 (MQ=255)
                  gATGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCAttt                                                  <  1:1488711/71‑1 (MQ=255)
                    tGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAATCGAGATTTTATATAACCAttttt                                                >  1:569876/1‑71 (MQ=255)
                    tGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCga                                                                    <  1:1776392/51‑1 (MQ=255)
                    tGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCga                                                                    <  1:3456/51‑1 (MQ=255)
                    tGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCAttttt                                                <  1:892741/71‑1 (MQ=255)
                    tGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCAttttt                                                <  1:1400740/71‑1 (MQ=255)
                    tGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCAttttt                                                >  1:1986533/1‑71 (MQ=255)
                         gagaTTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTTCACGt                                           >  1:3737757/1‑71 (MQ=255)
                            aTTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCga                                                                    <  1:1713099/43‑1 (MQ=255)
                            aTTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCga                                                                    <  1:2775406/43‑1 (MQ=255)
                            aTTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTTCAc                                             <  1:908790/66‑1 (MQ=255)
                             tttGTTAGATAACCCTTGTCTGAACATTAAAATAGGAATCGAGATTTTATATAACCATTTTTCACGTtgtg                                       <  1:1284907/71‑1 (MQ=255)
                             tttGTTAGATAACCCTTGTCTGAACATTAAAATAGGAATCGAGATTTTATATAACCATTTTTCACGTtgtg                                       <  1:2024537/71‑1 (MQ=255)
                             tttGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTTCACGTtgtg                                       <  1:2554336/71‑1 (MQ=255)
                              ttGTTAGATAACCCTTGTCTGAACATTAAAATAGGAATCGAGATTTTATATAACCATTTTTCACGTtgtg                                       >  1:3924766/1‑70 (MQ=255)
                              ttGTTAGATAACCCTTGTCTGAACATTAAAATAGGAATCGAGATTTTATATAACCATTTTTCACGTtgtg                                       >  1:1472134/1‑70 (MQ=255)
                              ttGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTTCACGt                                           >  1:1303366/1‑66 (MQ=255)
                               tGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTTCACGTtgtg                                       >  1:1204128/1‑69 (MQ=255)
                                     aTAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTTCACGTtgtg                                       <  1:3118256/63‑1 (MQ=255)
                                     aTAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTTCACGTtgtg                                       <  1:1982032/63‑1 (MQ=255)
                                     aTAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTTCACGTtgtg                                       <  1:341307/63‑1 (MQ=255)
                                       aaCCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTTCACGTTGTGGAGTAACGTg                             >  1:852311/1‑71 (MQ=255)
                                         cccTTGTCTGAACATTAAAATAGGAATCGAGATTTTATATAACCATTTTTCACGTTGTGGAGTAACGTGGc                           <  1:2283887/71‑1 (MQ=255)
                                         cccTTGTCTGAACATTAAAATAGGAATCGAGATTTTATATAACCATTTTTCACGTTGTGGAGTAACGTGGc                           <  1:2593497/71‑1 (MQ=255)
                                         cccTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTTCACGTTGTGGAGTAACGTGGc                           <  1:3394910/71‑1 (MQ=255)
                                         cccTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTTCACGTTGTGGAGTAACGTGGc                           <  1:1163851/71‑1 (MQ=255)
                                         cccTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTTCACGTTGTGGAGTAACGTGGc                           <  1:3842544/71‑1 (MQ=255)
                                         cccTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTTCACGTTGTGGAGTAACGTGGc                           <  1:3877050/71‑1 (MQ=255)
                                         cccTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTTCACGTTGTGGAGTAACGTGGc                           <  1:1952267/71‑1 (MQ=255)
                                         cccTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTTCACGTTGTGGAGTAACGTGGc                           <  1:1942402/71‑1 (MQ=255)
                                         cccTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTTCACGTTGTGGAGTAACGTGGc                           <  1:14527/71‑1 (MQ=255)
                                          ccTTGTCTGAACATTAAAATAGGAATCGAGATTTTATATAACCATTTTTCACGTTGTGGAGTAACGTGGc                           <  1:1438263/70‑1 (MQ=255)
                                          ccTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTTCACGTtgtg                                       <  1:3337436/58‑1 (MQ=255)
                                          ccTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTTCACGTTGTGGAGTaa                                 >  1:3492039/1‑64 (MQ=255)
                                          ccTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTTCACGTTGTGGAGTAACGTGGc                           <  1:857702/70‑1 (MQ=255)
                                           cTTGTCTGAACATTAAAATAGGAATCGAGATTTTATATAACCATTTTTCACGTTGTGGAGTAACGTGGCaa                         <  1:1490041/71‑1 (MQ=255)
                                           cTTGTCTGAACATTAAAATAGGAACCGAGATTTta                                                             <  1:840247/35‑1 (MQ=255)
                                           cTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTTCACGTTGTGGa                                     >  1:2292441/1‑59 (MQ=255)
                                           cTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTTCACGTTGTGGAGTAACGTGGCaa                         <  1:2680660/71‑1 (MQ=255)
                                           cTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTTCACGTTGTGGAGTAACGTGGCaa                         <  1:3381499/71‑1 (MQ=255)
                                              gTCTGAACATTAAAATAGGAATCGAGATTTTATATAACCATTTTTCACGTTGTGGAGTAACGTGGCAATg                       <  1:964266/70‑1 (MQ=255)
                                              gTCTGAACATTAAAATAGGAATCGAGATTTTATATAACCATTTTTCACGTTGTGGAGTAACGTGGCAATGc                      >  1:3926043/1‑71 (MQ=255)
                                              gTCTGAACATTAAAATAGGAATCGAGATTTTATATAACCATTTTTCACGTTGTGGAGTAACGTGGCAATGc                      >  1:2222113/1‑71 (MQ=255)
                                              gTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTTCACGTTGTGGAGTAACGTGGCAAt                        <  1:735309/69‑1 (MQ=255)
                                              gTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTTCACGTTGTGGAGTAACGTGGCAATGc                      >  1:2192288/1‑71 (MQ=255)
                                              gTCTGAACATTAAAATAGGAAACGAGATTTTATATAACCATTTTTCACGTTGTGGAGTAACGTGGCAATGc                      >  1:2395228/1‑71 (MQ=255)
                                               tCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTTCACGTtgtg                                       <  1:2219602/53‑1 (MQ=255)
                                                cTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTTCACGTtgtg                                       <  1:3230032/52‑1 (MQ=255)
                                                cTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTTCACGTtgtg                                       <  1:1364036/52‑1 (MQ=255)
                                                   aaCATTAAAATAGGAACCGAGATTTTATATAACCATTTTTCACGTTGTGGAGTAACGTGGCAATGCCTggg                 >  1:3777978/1‑71 (MQ=255)
                                                   aaCATTAAAATAGGAACCGAGATTTTATATAACCATTTTTCACGTTGTGGAGTAACGTGGCAATGCCTggg                 >  1:2880151/1‑71 (MQ=255)
                                                      aTTAAAATAGGAATCGAGATTTTATATAACCATTTTTCACGTTGTGGAGTAACGTGGc                           <  1:818692/58‑1 (MQ=255)
                                                       ttAAAATAGGAATCGAGATTTTATATAACCATTTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACt             <  1:689885/71‑1 (MQ=255)
                                                       ttAAAATAGGAACCGAGATTTTATATAACCATTTTTCACGt                                           <  1:1294067/41‑1 (MQ=255)
                                                       ttAAAATAGGAACCGAGATTTTATATAACCATTTTTCACGt                                           <  1:227747/41‑1 (MQ=255)
                                                       ttAAAATAGGAACCGAGATTTTATATAACCATTTTTCACGt                                           <  1:469902/41‑1 (MQ=255)
                                                        tAAAATAGGAACCGAGATTTTATATAACCATTTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACt             >  1:1839036/1‑70 (MQ=255)
                                                         aaaaTAGGAATCGAGATTTTATATAACCATTTTTCACGTTGTGGa                                     >  1:419373/1‑45 (MQ=255)
                                                         aaaaTAGGAATCGAGATTTTATATAACCATTTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTta           >  1:1709963/1‑71 (MQ=255)
                                                         aaaaTAGGAATCGAGATTTTATATAACCATTTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTta           >  1:489150/1‑71 (MQ=255)
                                                         aaaaTAGGAACCGAGATTTTATATAACCATTTTTCACGTTGTGGAGTAACGTGGCAATGcc                     >  1:2481209/1‑61 (MQ=255)
                                                         aaaaTAGGAACCGAGATTTTATATAACCATTTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTta           >  1:1144833/1‑71 (MQ=255)
                                                         aaaaTAGGAACCGAGATTTTATATAACCATTTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTta           >  1:2966949/1‑71 (MQ=255)
                                                          aaaTAGGAATCGAGATTTTATATAACCATTTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTtat          <  1:3665599/71‑1 (MQ=255)
                                                          aaaTAGGAACCGAGATTTTATATAACCATTTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTtat          <  1:2125309/71‑1 (MQ=255)
                                                          aaaTAGGAACCGAGATTTTATATAACCATTTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTtat          <  1:2365276/71‑1 (MQ=255)
                                                          aaaTAGGAACCGAGATTTTATATAACCATTTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTtat          <  1:3777402/71‑1 (MQ=255)
                                                          aaaTAGGAACCGAGATTTTATATAACCATTTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTtat          <  1:2418794/71‑1 (MQ=255)
                                                          aaaTAGGAACCGAGATTTTATATAACCATTTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTtat          <  1:1483504/71‑1 (MQ=255)
                                                            aTAGGAACCGAGATTTTATATAACCATTTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATaa        <  1:3277813/71‑1 (MQ=255)
                                                              aGGAATCGAGATTTTATATAACCATTTTTCACGTTGTGGAGTAACGTGGc                           <  1:4058096/50‑1 (MQ=255)
                                                              aGGAATCGAGATTTTATATAACCATTTTTCACGTTGTGGAGTAACGTGGc                           <  1:2472952/50‑1 (MQ=255)
                                                                 aaCCGAGATTTTATATAACCATTTTTCACGTTGTGGAGTAACGTGGCAATGcc                     <  1:3147748/53‑1 (MQ=255)
                                                                 aaCCGAGATTTTATATAACCATTTTTCACGTTGTGGAGTAACGTGGCAATGcc                     <  1:3507509/53‑1 (MQ=255)
                                                                 aaCCGAGATTTTATATAACCATTTTTCACGTTGTGGAGTAACGTGGCAATGcc                     <  1:1204134/53‑1 (MQ=255)
                                                                 aaCCGAGATTTTATATAACCATTTTTCACGTTGTGGAGTAACGTGGCAATGcc                     <  1:597607/53‑1 (MQ=255)
                                                                 aaCCGAGATTTTATATAACCATTTTTCACGTTGTGGAGTAACGTGGCAATGcc                     <  1:3099476/53‑1 (MQ=255)
                                                                 aaCCGAGATTTTATATAACCATTTTTCACGTTGTGGAGTAACGTGGCAATGcc                     <  1:4039584/53‑1 (MQ=255)
                                                                  aCCGAGATTTTATATAACCATTTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCggg  <  1:280215/71‑1 (MQ=255)
                                                                   |                                                                     
AAACGTGGAATAATAAAAGATGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGG  >  W3110S.gb/2992856‑2992992

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: