Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,688,001 A→G 40.0% S440G (AGT→GGT)  yifB → predicted bifunctional protein, enzyme and transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,688,0010AG40.0% 37.5 / 70.2 60S440G (AGT→GGT) yifBpredicted bifunctional protein, enzyme and transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (17/19);  new base G (16/8);  total (33/27)
Fisher's exact test for biased strand distribution p-value = 1.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAACAACGCGTAATGGCCGCCAGAGAGCGCCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGCGATGTGGTTGGAAGGAACACTGAT  >  W3110S.gb/3687932‑3688068
                                                                     |                                                                   
aaaCAACGCGTAATGGCCGCCAGAGAGCGCCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATAg                                                                    >  1:2760543/1‑71 (MQ=255)
    aacGCGTAATGGCCGCCAGAGAGCGCCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATAGTCCg                                                                <  1:2521794/71‑1 (MQ=255)
      cgcgTAATGGCCGCCAGAGAGCGCCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATa                                                                     >  1:482845/1‑64 (MQ=255)
        cgTAATGGCCGCCAGAGAGCGCCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATAGTCCGGaaa                                                            <  1:135852/71‑1 (MQ=255)
          tAATGGCCGCCAGAGAGCGCCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATAGTCCGGaaa                                                            <  1:4096580/69‑1 (MQ=255)
            aTGGCCGCCAGAGAGCGCCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATAGTCCGGAAATag                                                         >  1:2766827/1‑69 (MQ=255)
             tGGCCGCCAGAGAGCGCCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATAGTCCGGAAATACg                                                        <  1:4730690/70‑1 (MQ=255)
             tGGCCGCCAGAGAGCGCCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGc                                                       >  1:818002/1‑71 (MQ=255)
              ggCCGCCAGAGAGCGCCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGc                                                       >  1:117982/1‑70 (MQ=255)
               gccgccAGAGAGCGCCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATGGTCCGGAAATACGcc                                                      <  1:262680/70‑1 (MQ=255)
               gccgccAGAGAGCGCCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATGGTCCGGAAATACGcc                                                      <  1:1323151/70‑1 (MQ=255)
               gccgccAGAGAGCGCCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATGGTCCGGAAATACGcc                                                      <  1:2744131/70‑1 (MQ=255)
                  gccAGAGAGCGCCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATGGTCCGGAAATACGCCAAtt                                                  <  1:4446910/71‑1 (MQ=255)
                   ccAGAGAGCGCCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATAGTCCGGa                                                              <  1:333683/58‑1 (MQ=255)
                     agagagCGCCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAAtt                                                  <  1:152825/68‑1 (MQ=255)
                        gagCGCCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATGGTCCGGAAATACGCCAATTCTGCaa                                            <  1:4570235/71‑1 (MQ=255)
                        gagCGCCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATGGTCCGGAAATACGCCAATTCTGCa                                             >  1:2230429/1‑70 (MQ=255)
                             ccAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTg                                       <  1:1701311/71‑1 (MQ=255)
                              cAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTg                                       <  1:697553/70‑1 (MQ=255)
                                aTTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTc                                                 >  1:4061511/1‑58 (MQ=255)
                                aTTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTga                                      >  1:3878103/1‑69 (MQ=255)
                                 tttAAGCGGCAGAATAAACTGAATGCCTGGCTGGATGGTCCGGAAATACGCCAATTCTGCAAACTTga                                      >  1:4698452/1‑68 (MQ=255)
                                  ttAAGCGGCAGAATAAACTGAATGCCTGGCTGGATGGTCCGGAAATACGCCAATTCTGCAAACTTGagag                                   >  1:2983790/1‑70 (MQ=255)
                                   tAAGCGGCAGAATAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGCCTATTCTGCAAACTTGAGAGc                                  >  1:3950502/1‑70 (MQ=255)
                                     aGCGGCAGAATAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGaa                               <  1:3231009/71‑1 (MQ=255)
                                       cGGCAGAATAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCg                                 >  1:996407/1‑67 (MQ=255)
                                            gAATAAACTGAATGCCTGGCTGGATGGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGCGa                        >  1:3493928/1‑71 (MQ=255)
                                            gAATAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAAc                                          <  1:1372497/53‑1 (MQ=255)
                                            gAATAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTg                                       <  1:3232423/56‑1 (MQ=255)
                                            gAATAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGCGa                        >  1:1148145/1‑71 (MQ=255)
                                              aTAAACTGAATGCCTGGCTGGATGGTCCGGAAATACGCCAATTCTGCAAACtt                                        <  1:2192248/53‑1 (MQ=255)
                                              aTAAACTGAATGCCTGGCTGGATGGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCg                                 >  1:3080967/1‑60 (MQ=255)
                                              aTAAACTGAATGCCTGGCTGGATGGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCg                                 >  1:4502958/1‑60 (MQ=255)
                                              aTAAACTGAATGCCTGGCTGGATGGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCg                                 >  1:3881898/1‑60 (MQ=255)
                                              aTAAACTGAATGCCTGGCTGGATAGTCCGGAAATAc                                                         >  1:3111110/1‑36 (MQ=255)
                                               tAAACTGAATGCCTGGCTGGATGGTCCGGAAATACGCAAATTCTGCAAACTTGAGAGCGAAGATGCGAtgt                     >  1:2003694/1‑71 (MQ=255)
                                               tAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAAtt                                                  >  1:3063023/1‑42 (MQ=255)
                                               tAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGCg                         <  1:2445221/67‑1 (MQ=255)
                                               tAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGCGAtgt                     >  1:2233729/1‑71 (MQ=255)
                                                aaaCTGAATGCCTGGCTGGATGGTCCGGAAATACGc                                                       >  1:850825/1‑36 (MQ=255)
                                                aaaCTGAATGCCTGGCTGGATGGTCCGGAAATACGc                                                       >  1:1530504/1‑36 (MQ=255)
                                                  aCTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTga                                      >  1:1050961/1‑51 (MQ=255)
                                                     gAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGaa                               <  1:89006/55‑1 (MQ=255)
                                                      aaTGCCTGGCTGGATGGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGCGATGTGGTTgg               <  1:3603015/70‑1 (MQ=255)
                                                      aaTGCCTGGCTGGATGGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGCGATGTGGTTgg               <  1:4687335/70‑1 (MQ=255)
                                                      aaTGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATg                           >  1:3362564/1‑58 (MQ=255)
                                                      aaTGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATg                           >  1:3234059/1‑58 (MQ=255)
                                                       aTGCCTGGCTGGATGGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGCGATGTGGTTgg               >  1:4703415/1‑69 (MQ=255)
                                                       aTGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGa                             >  1:1839419/1‑55 (MQ=255)
                                                        tGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGCGATGTGGTTggaag            <  1:2547858/71‑1 (MQ=255)
                                                         gCCTGGCTGGATGGTCCGGAAATACGCCAATTCTGGAAACTTGAGAGCGaa                               >  1:2744348/1‑51 (MQ=39)
                                                         gCCTGGCTGGATGGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGaa                               >  1:1765114/1‑51 (MQ=255)
                                                          cccggctggATGGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGCGATGTGGTTggaagg           >  1:417345/4‑70 (MQ=255)
                                                          ccTGGCTGGATGGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGCGATGTGGt                  >  1:1608127/1‑63 (MQ=255)
                                                          ccTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGCGATGTGGTTggaagg           <  1:3190305/70‑1 (MQ=255)
                                                           ctggctggATGGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGa                             >  1:1892389/1‑51 (MQ=255)
                                                           ctggctggATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGaa                               <  1:822821/49‑1 (MQ=255)
                                                            tggctggATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGCGATGTGGTTggaaggaa         <  1:1584824/70‑1 (MQ=255)
                                                                  gATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGCGATGTGGTTGGAAGGAACACTg    <  1:1888636/69‑1 (MQ=255)
                                                                  gATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGCGATGTGGTTGGAAGGAACACTg    <  1:2138186/69‑1 (MQ=255)
                                                                  gATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGCGATGTGGTTGGAAGGAACACTGAt  >  1:4730766/1‑71 (MQ=255)
                                                                     |                                                                   
AAACAACGCGTAATGGCCGCCAGAGAGCGCCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGCGATGTGGTTGGAAGGAACACTGAT  >  W3110S.gb/3687932‑3688068

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: