Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,721,287 T→C 5.8% intergenic (+18/‑96) yfiQ → / → pssA fused predicted acyl‑CoA synthetase NAD(P)‑binding subunit and ATP‑binding subunit/phosphatidylserine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,721,2870TC 5.8% 187.8 / 2.3 69intergenic (+18/‑96)yfiQ/pssAfused predicted acyl‑CoA synthetase NAD(P)‑binding subunit and ATP‑binding subunit/phosphatidylserine synthase
Reads supporting (aligned to +/- strand):  ref base T (29/36);  new base C (2/2);  total (31/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

GAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGTTATTATTGCCCGCTTATGTCGTCTGCATTGCACAGAGGACCCTTC  >  W3110S.gb/2721219‑2721356
                                                                    |                                                                     
gagGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGtt                                                                      >  1:455392/1‑70 (MQ=255)
         gTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTa                                                             <  1:2990692/70‑1 (MQ=255)
          ttGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAAt                                                           <  1:3853395/71‑1 (MQ=255)
          ttGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAAt                                                           <  1:3667687/71‑1 (MQ=255)
             ggCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGAGCACTTAATCgg                                                        <  1:3358814/71‑1 (MQ=255)
             ggCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCgg                                                        <  1:3076846/71‑1 (MQ=255)
             ggCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCgg                                                        <  1:3415828/71‑1 (MQ=255)
             ggCTTACGCTAAATCTTGCCCAGCGCGAGGAACCATGAGTAAGGTACTGGAAATGTTGACCACTTAATctg                                                        <  1:2668910/71‑3 (MQ=255)
                ttACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGt                                                                       >  1:197875/1‑53 (MQ=255)
                  aCGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGa                                                                    >  1:3689072/1‑54 (MQ=255)
                  aCGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGAc                                                     <  1:1611314/69‑1 (MQ=255)
                   cGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGtt                                                                      <  1:2712750/51‑1 (MQ=255)
                      tAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAAt                                                           >  1:1453682/1‑59 (MQ=255)
                      tAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAAt                                                           >  1:273391/1‑59 (MQ=255)
                      tAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGgtgt                                               >  1:1665012/1‑71 (MQ=255)
                      tAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGgtgt                                               >  1:1909947/1‑71 (MQ=255)
                      tAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGgtgt                                               >  1:1989498/1‑71 (MQ=255)
                        aaTCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGa                                                                    >  1:2078896/1‑48 (MQ=255)
                        aaTCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGtta                                             >  1:3098005/1‑71 (MQ=255)
                          tCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCgg                                                        <  1:1080620/58‑1 (MQ=255)
                            ttGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGa                                                                    >  1:2971607/1‑44 (MQ=255)
                            ttGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGa                                                      >  1:1807335/1‑58 (MQ=255)
                            ttGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGTtattat                                         >  1:2897352/1‑71 (MQ=255)
                            ttGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGTtattat                                         >  1:2254195/1‑71 (MQ=255)
                              gCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATAGGGACTGgtgt                                               >  1:1635915/1‑63 (MQ=255)
                                ccAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGTTATTATTGcc                                     >  1:3480970/1‑71 (MQ=255)
                                  aGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGTTATTATTGccc                                    >  1:2407345/1‑70 (MQ=255)
                                   gcgcGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCgg                                                        <  1:565380/49‑1 (MQ=255)
                                    cgcgAGGATTCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTTGTGTTATTATTGCCCGCt                                 <  1:2985406/71‑1 (MQ=255)
                                    cgcgAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGTTATTATTGCCCGCt                                 <  1:2799822/71‑1 (MQ=255)
                                    cgcgAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGTTATTATTGCCCGCt                                 <  1:444471/71‑1 (MQ=255)
                                    cgcgAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGTTATTATTGCCCGCt                                 <  1:1157554/71‑1 (MQ=255)
                                    cgcgAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGTTATTATTGCCCGCt                                 <  1:2626628/71‑1 (MQ=255)
                                    cgcgAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGTTATTATTGCCCGCt                                 <  1:1703098/71‑1 (MQ=255)
                                    cgcgAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGTTATTATTGCCCGCt                                 <  1:2419149/71‑1 (MQ=255)
                                    cgcgAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGTTATTATTGCCCGCt                                 <  1:1976297/71‑1 (MQ=255)
                                     gcgAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCTGGACTGGTGTTATTATTGCCCGCt                                 <  1:1673723/70‑1 (MQ=255)
                                     gcgAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGTTATTATTGCCCGCt                                 <  1:3742848/70‑1 (MQ=255)
                                       gAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGTTATTATTGCCCGCTTAt                              >  1:128340/1‑71 (MQ=255)
                                       gAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGTTATTATTGCCCGCTTAt                              >  1:311914/1‑71 (MQ=255)
                                        aGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGTTATTATTGCCCGCTTATg                             <  1:2394341/71‑1 (MQ=255)
                                         ggAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCgg                                                        <  1:3247752/43‑1 (MQ=255)
                                         ggAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGTTATTATTGCCCGCTTATgt                            <  1:3535695/71‑1 (MQ=255)
                                          gAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGTTATTATTGCCCGCTTATgtc                           <  1:837501/71‑1 (MQ=255)
                                          gAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGTTATTATTGCCCGCTTATgtc                           <  1:264170/71‑1 (MQ=255)
                                           aaTCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGttatt                                           <  1:2139416/54‑1 (MQ=255)
                                           aaTCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGttatt                                           <  1:3262464/54‑1 (MQ=255)
                                            aTCATGAGTAAGGTAATGGAAATGTTGACCACTTAATCGGGACTGGTGTTATTATTGCCCGCTTATgtcgt                         >  1:1505910/1‑71 (MQ=255)
                                             tCATGAGTAAGGTACTGGAAATGTTGACCACTTAATc                                                          <  1:2430991/37‑1 (MQ=255)
                                             tCATGAGTAAGGTACTGGAAATGCTGACCACTTAATCgg                                                        <  1:547175/39‑1 (MQ=255)
                                             tCATGAGTAAGGTACTGGAAATGCTGACCACTTAATCgg                                                        <  1:3360830/39‑1 (MQ=255)
                                               aTGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGTTATTATTGc                                      >  1:2052786/1‑55 (MQ=255)
                                               aTGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGTTATTATTGc                                      >  1:1630450/1‑55 (MQ=255)
                                               aTGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGTTATTATTGc                                      >  1:1259020/1‑55 (MQ=255)
                                               aTGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGTTATTATTGCCCGCTTa                               >  1:1867411/1‑62 (MQ=255)
                                               aTGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGTTATTATTGCCCGCTTATGTCGTCTg                      <  1:3481030/71‑1 (MQ=255)
                                               aTGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGTTATTATTGCCCGCTTATGTCGTCTg                      <  1:950603/71‑1 (MQ=255)
                                               aTGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGTGTTATTATTGCCCGCTTATGTCGTCTg                      >  1:792262/1‑70 (MQ=255)
                                               aTGAGTAAGGTACTGGAAATGCTGACCACTTAATCGGGa                                                      >  1:3349788/1‑39 (MQ=255)
                                               aTGAGTAAGGTACTGGAAATGCTGACCACTTAATCGGGa                                                      >  1:2394615/1‑39 (MQ=255)
                                                 gAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGTTATTATTGccc                                    >  1:2125042/1‑55 (MQ=255)
                                                 gAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGTTATTATTGccc                                    >  1:2125078/1‑55 (MQ=255)
                                                 gAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGTTATTATTGCCCGCTTATGTCGTCTGCa                    >  1:2377758/1‑71 (MQ=255)
                                                  aGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGTTATTATTGCCCGCTTATgtcgtc                        >  1:1355219/1‑66 (MQ=255)
                                                   gTAAGGTACTGGATATGTTGACCACTTAATCGGGACt                                                    >  1:746565/1‑37 (MQ=255)
                                                   gTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGttatta                                          <  1:3426016/47‑1 (MQ=255)
                                                   gTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGttatta                                          <  1:2658374/47‑1 (MQ=255)
                                                        gTACTGGAAATGTTGACCACTTAATCGGGACTGGTGTTATTATTGCCCGCTTATGTCGTCTGCATTGcaca             <  1:2497894/71‑1 (MQ=255)
                                                         tACTGGAAATGTTGACCACTTAATCGGGACTGGTGTTATTATTGCCCGCTTATGTCGTCTGCATTGCACag            <  1:2265629/71‑1 (MQ=255)
                                                           cTGGAAATGTTGACCACTTAATCGGGACTGGTGTTATTATTGCCCGCTTATGTCGTCTGCATTGc                >  1:2264265/1‑65 (MQ=255)
                                                            tGGAAATGTTGACCACTTAATCGGGACTGGTGTTATTATTGCCCGCt                                 <  1:1477345/47‑1 (MQ=255)
                                                            tGGAAATGTTGACCACTTAATCGGGACTGGTGTTATTATTGCCCGCt                                 <  1:2147266/47‑1 (MQ=255)
                                                                    ttGACCACTTAATCGGGACTGGTGTTATTATTGCCCGCTTATGTCGTCTGCATTTCACAGAGGACCCTTc  >  1:2517244/1‑70 (MQ=255)
                                                                    ttGACCACTTAATCGGGACTGGTGTTATTATTGCCCGCTTATGTCGTCTGCATTGCACAGAGGACCCTTc  >  1:3393825/1‑70 (MQ=255)
                                                                    |                                                                     
GAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGTTATTATTGCCCGCTTATGTCGTCTGCATTGCACAGAGGACCCTTC  >  W3110S.gb/2721219‑2721356

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