Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,688,001 A→G 41.5% 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,0010AG41.5% 17.3 / 50.2 41S440G (AGT→GGT) yifBpredicted bifunctional protein, enzyme and transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (14/10);  new base G (12/5);  total (26/15)
Fisher's exact test for biased strand distribution p-value = 5.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCGTAATGGCCGCCAGAGAGCGCCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGCGATGTGGTTGGAAGGAACACTGA  >  W3110S.gb/3687938‑3688067
                                                               |                                                                  
cgcgTAATGGCCGCCAGAGAGCGCCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATAGTCCGGa                                                             <  1:1992782/71‑1 (MQ=255)
 gcgTAATGGCCGCCAGAGAGCGCCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATGGTCCGGa                                                             >  1:2017462/1‑70 (MQ=255)
      aTGGCCGCCAGAGAGCGCCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATAGTCCGGAAATACg                                                       >  1:151872/1‑71 (MQ=255)
           cgccAGAGAGCGCCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTTGATAGTCCGGaa                                                            <  1:793746/61‑1 (MQ=255)
            gccAGAGAGCGCCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATGGTCCGGAAATa                                                         >  1:634443/1‑63 (MQ=255)
            gccAGAGAGCGCCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAAtt                                                 >  1:886602/1‑71 (MQ=255)
            gccAGAGAGCGCCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAAtt                                                 >  1:2337252/1‑71 (MQ=255)
               agagagCGCCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATGGTCCGGaaa                                                           >  1:2392178/1‑58 (MQ=255)
               agagagCGCCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCt                                               >  1:1563629/1‑70 (MQ=255)
                  gagCGCCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCa                                            >  1:318130/1‑70 (MQ=255)
                   agCGCCAATTTAAGCGGCAGAATATACTGAATGCCTGGCTGGATGGTCCGGAAATACGCCAATTCTGCaaa                                          >  1:537147/1‑71 (MQ=255)
                   agCGCCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATGGTCCGGAAATACGCCAATTCTGCaaa                                          <  1:2378935/71‑1 (MQ=255)
                   agCGCCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATAGTCCGGaaa                                                           <  1:626241/54‑1 (MQ=255)
                   agCGCCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATAGTCCGGAAATa                                                         >  1:39484/1‑56 (MQ=255)
                   agCGCCAATTTAAGAGGCAGAATAAACTGAATGCCTGGCTGGATa                                                                    >  1:1757223/1‑45 (MQ=255)
                     cgcCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATGGTCCGGAAATACGCCAATTCTGCAAACt                                        >  1:1184604/1‑71 (MQ=255)
                      gcCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATGGTCCGGAAATACGCCa                                                    <  1:1839329/58‑1 (MQ=255)
                      gcCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACtt                                       <  1:714413/71‑1 (MQ=255)
                          aTTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATGGTCCGGAAATACGCCAATTCTGCAAACTTga                                     >  1:2028609/1‑69 (MQ=255)
                           tttAAGCGGCAGAATAAACTGAATGCCTGGCTTGATAGTCCGGaaa                                                           <  1:969530/46‑1 (MQ=255)
                           tttAAGCGGCAGAATAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTGagag                                  <  1:1426593/71‑1 (MQ=255)
                                    cAGAATAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGc                         <  1:610805/71‑1 (MQ=255)
                                     aGAATAAACTGAATGCCTGGCTGGATGGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGc                         >  1:1251926/1‑70 (MQ=255)
                                      gAATAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGCGa                       >  1:61891/1‑71 (MQ=255)
                                         tAAACTGAATGCCTGGCTGGATGGTCCGGAAATACGc                                                      >  1:1004063/1‑37 (MQ=255)
                                         tAAACTGAATGCCTGGCTGGATGGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGa                               <  1:2250597/60‑1 (MQ=255)
                                         tAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGCCaa                                                   >  1:975085/1‑40 (MQ=255)
                                         tAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAAc                                         <  1:1326692/50‑1 (MQ=255)
                                         tAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGc                         <  1:1270960/66‑1 (MQ=255)
                                         tAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGCGAtgt                    >  1:2205039/1‑71 (MQ=255)
                                           aaCTGAATGCCTGGCTGGATGGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGCGATGTgg                  >  1:473079/1‑71 (MQ=255)
                                           aaCTGAATGCCTGGCTGGATGGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGCGATGTgg                  >  1:823619/1‑71 (MQ=255)
                                           aaCTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGCGATGTgg                  >  1:735865/1‑71 (MQ=255)
                                              tGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGCGATGTGGTTg               >  1:2372714/1‑71 (MQ=255)
                                                   gCCTGGCTGGATGGTCCGGAAATACGCCAATTCTGCAAACtt                                       <  1:1093302/42‑1 (MQ=255)
                                                    ccTGGCTGGATGGTCCGGAAATACTCCAATTCTGCAAACTTGAGAGCGAAGATGCGATGTGGTTggaagg          <  1:2044397/70‑1 (MQ=255)
                                                    ccTGGCTGGATAGTCCGGAAATACGCAAATTCTGCAAACTTGAGAGCGAAGATGAGATGTGGTTggaagg          >  1:2160374/1‑70 (MQ=255)
                                                     ctggctggATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGCGATGTGGTTggaaggaa        <  1:2258963/71‑1 (MQ=255)
                                                        gctggATAGTCCGGAAATACGCCAATTCTGCAAACTTgaga                                   >  1:1819273/1‑41 (MQ=255)
                                                          tggATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGCGATGTGGTTGGAAGGAACACt    >  1:950003/1‑70 (MQ=255)
                                                           ggATGGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGCGATGTGGTTGGAAGGAACACTGa  >  1:1579643/1‑71 (MQ=255)
                                                            gATGGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGa                            >  1:637233/1‑44 (MQ=255)
                                                               |                                                                  
CGCGTAATGGCCGCCAGAGAGCGCCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGCGATGTGGTTGGAAGGAACACTGA  >  W3110S.gb/3687938‑3688067

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