Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 810,113 G→A 27.5% M116I (ATG→ATA bioB → biotin synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb810,1130GA27.5% 66.9 / 26.4 40M116I (ATG→ATAbioBbiotin synthase
Reads supporting (aligned to +/- strand):  ref base G (17/12);  new base A (4/7);  total (21/19)
Fisher's exact test for biased strand distribution p-value = 2.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGTATGGGCGCGGCGTGGAAGAATCCCCACGAACGCGATATGCCGTACCTGGAACAAATGGTGCAGGGGGTAAAAGCGATGGGGCTGGAGGCGTGTATGACGCTGGGCACGTTGAGTGAATCT  >  W3110S.gb/810053‑810176
                                                            |                                                               
cTGTATGGGCGCGGCGTGGAAGAATCCCCACGAACGCGATATGCCGTACCTGGAACAAATGGTGCAggggg                                                       >  1:303089/1‑71 (MQ=255)
     tGGGCGCGGCGTGGAAGAATCCCCAGGAACGCGATATGCCGTACCTGGAACAAATGGTGCAGGGGGTaa                                                    >  1:1084091/1‑69 (MQ=255)
     tGGGCGCGGCGTGGAAGAATCCCCACGAACGCGATATGCCGTACCTGGAACAAATGGTGCAGGGGGTaa                                                    >  1:3089189/1‑69 (MQ=255)
     tGGGCGCGGCGTGGAAGAATCCCCACGAACGCGATATGCCGTACCTGGAACAAATGGTGCAGGGGGTaa                                                    >  1:1517825/1‑69 (MQ=255)
     tGGGCGCGGCGTGGAAGAATCCCCACGAACGCGATATGCCGTACCTGGAACAAATGGTGCAGGGGGTaa                                                    >  1:2303764/1‑69 (MQ=255)
       ggCGCGGCGTGGAAGAATCCCCACGAACGCGATATGCCGTACCTGGAACAAATGGTGCAGGGGGTAAAAGc                                                >  1:4642177/1‑71 (MQ=255)
           cggcgTGGAAGAATCCCCACGAACGCGATATGCCGTACCTGGAACAAATAGTg                                                              <  1:829111/53‑1 (MQ=255)
                tGGAAGAATCCCCACGAACGCGATATGCCGTACCTGGAACAAATGGTGCAGGGGGTa                                                     <  1:2469414/57‑1 (MQ=255)
                tGGAAGAATCCCCACGAACGCGATATGCCGTACCTGGAACAAATGGTGCAGGGGGTa                                                     <  1:2365619/57‑1 (MQ=255)
                  gaagaaTCCCCACGAACGCGATATGCCGTACCTGGAACAAATGGTGCAGGGGGTaaa                                                   >  1:2602753/1‑57 (MQ=255)
                  gaagaaTCCCCACGAACGCGATATGCCGTACCTGGAACAAATGGTGCAGGGGGTAAAAGCGATGGGGCTg                                      <  1:1914460/70‑1 (MQ=255)
                  gaagaaTCCCCACGAACGCGATATGCCGTACCTGGAACAAATGGTGCAGGGGGTAAAAGCGATGGGGCTg                                      <  1:2731461/70‑1 (MQ=255)
                  gaagaaTCCCCACGAACGCGATATGCCGTACCTGGAACAAATGGTGCAGGGGGTAAAAGCGATGGGGCTg                                      >  1:805628/1‑70 (MQ=255)
                  gaagaaTCCCCACGAACGCGATATGCCGTACCTGGAACAAATAGTGCAGGGGGTAAAAGCGATGGGGCTg                                      >  1:1400568/1‑70 (MQ=255)
                  gaagaaTCCCCACGAACGCGATATGCCGTACCTGGAACAAATAGTGCAGGGGGTAAAAGCGATGGGGCTg                                      >  1:3917291/1‑70 (MQ=255)
                    agaaTCCCCACGAACGCGATATGCCGTACCTGGAACAAATGGTGCAGGGGGTAAAAGCGATGGGGCTGGAg                                   >  1:2047253/1‑71 (MQ=255)
                      aaTCCCCACGAACGCGATATGCCGTACCTGGAACAAATGGTGCAGGGGGt                                                      <  1:4766288/50‑1 (MQ=255)
                           ccACGAACGCGATATGCCGTACCTGGAACAAATAGTGCAGGGGGTAAAAGc                                                <  1:4109011/51‑1 (MQ=255)
                            cACGAACGCGATATGCCGTACCTGGAACAAATGGTGCAGGGGGTAAAAGc                                                <  1:773767/50‑1 (MQ=255)
                              cGAACGCGATATGCCGTACCTGGAACAAATGGTGCAGGGGGTAAAAGCg                                               <  1:4184007/49‑1 (MQ=255)
                              cGAACGCGATATGCCGTACCTGGAACAAATGGTGCAGGGGGTAAAAGCg                                               <  1:349755/49‑1 (MQ=255)
                                    cgATATGCCGTACCTGGAACAAATGGTGCAGGGGGTAAAAGCGATGGGGCTgg                                     >  1:2173063/1‑53 (MQ=255)
                                    cgATATGCCGTACCTGGAACAAATGGTGCAGGGGGTAAAAGCGATGGGGCTgg                                     >  1:1230199/1‑53 (MQ=255)
                                    cgATATGCCGTACCTGGAACAAATGGTGCAGGGGGTAAAAGCGATGGGGCTgg                                     >  1:3426617/1‑53 (MQ=255)
                                    cgATATGCCGTACCTGGAACAAATGGTGCAGGGGGTAAAAGCGATGGGGCTgg                                     >  1:2289836/1‑53 (MQ=255)
                                        atGCCGTACCTGGAACAAATGGTGCAGGGGGTAAAAGCGATGGGGCTGGAGGCGTGTATGACGCTGGGCAc               >  1:528187/1‑71 (MQ=255)
                                        atGCCGTACCTGGAACAAATGGTGCAGGGGGTAAAAGCGATGGGGCTGGAGGCGTGTATGACGCTGGGCAc               >  1:4588281/1‑71 (MQ=255)
                                        atGCCGTACCTGGAACAAATAGTGCAGGGGGTAAAAGCGATGGGGCTGGAGGCGTGTATGACGCTGGGCa                >  1:4461979/1‑70 (MQ=255)
                                         tGCCGTACCTGGAACAAATGGTGCAGGGGGTAAAAGCGATGGGGCTGGAGGc                                 <  1:4159720/52‑1 (MQ=255)
                                         tGCCGTACCTGGAACAAATGGTGCAGGGGGTAAAAGCGATGGGGCTGGAGGCGTGTATGACGCTg                    >  1:903731/1‑65 (MQ=255)
                                         tGCCGTACCTGGAACAAATAGTGCAGGGGGTAAAAGc                                                <  1:4047265/37‑1 (MQ=255)
                                         tGCCGTACCTGGAACAAATAGTGCAGGGGGTAAAAGc                                                <  1:2408150/37‑1 (MQ=255)
                                         tGCCGTACCTGGAACAAATAGTGCAGGGGGTAAAAGc                                                <  1:1455639/37‑1 (MQ=255)
                                            cGTACCTGGAACAAATAGTGCAGGGGGTAAAAGCGATGGGGCTGGAGGCGTGTATGACGCTGGGCACGtt            >  1:4338993/1‑70 (MQ=255)
                                                ccTGGAACAAATGGTGCAGGGGGTAAAAGCGATGGGGCTGGAGGCGTGTATGACGCTGGGCACGTTGAGt        >  1:589857/1‑70 (MQ=255)
                                                  tGGAACAAATGGTGCAGGGGGTAAAAGCGATGGGGCTGGAGGCgtgt                             <  1:1219645/47‑1 (MQ=255)
                                                    gAACAAATAGTGCAGGGGGTAAAAGCGATGGGGCTGGAGGCGTGTATGACGCTGGGCACGTTGAGTGAAt    <  1:2934490/70‑1 (MQ=255)
                                                      aCAAATGGTGCAGCGGGTAAAAGCGATGGGGCTGGAGGCGTGTATGACGCTGGGCACGTTGAGTGAAtct  <  1:846769/70‑1 (MQ=255)
                                                       cAAATGGTGCAGGGGGTAAAAGCGATGGGGCTGGAGGCgtg                              <  1:4466738/41‑1 (MQ=255)
                                                        aaaTAGTGCAGGGGGTAAAAGCGATGGGGCTGGAGGCGTGTAt                           <  1:1327052/43‑1 (MQ=255)
                                                            |                                                               
CTGTATGGGCGCGGCGTGGAAGAATCCCCACGAACGCGATATGCCGTACCTGGAACAAATGGTGCAGGGGGTAAAAGCGATGGGGCTGGAGGCGTGTATGACGCTGGGCACGTTGAGTGAATCT  >  W3110S.gb/810053‑810176

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: