Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 811,947 G→T 23.3% G382V (GGC→GTC) 
A2S (GCA→TCA) 
bioF →
bioC →
8‑amino‑7‑oxononanoate synthase
predicted methltransferase, enzyme of biotin synthesis

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb811,9470GT23.3% 163.9 / 24.2 60G382V (GGC→GTC) 
A2S (GCA→TCA) 
bioF
bioC
8‑amino‑7‑oxononanoate synthase
predicted methltransferase, enzyme of biotin synthesis
Reads supporting (aligned to +/- strand):  ref base G (21/25);  new base T (4/10);  total (25/35)
Fisher's exact test for biased strand distribution p-value = 3.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

GACTGCGCTTAACGCTAACCGCTGCGCATGAAATGCAGGATATCGACCGTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATG  >  W3110S.gb/811878‑812013
                                                                     |                                                                  
gACTGCGCTTAACGCTAACCGCTGCGCATGAAATGCAGGATATCGACCGTCTGCTGGAGGTGCTGCATgg                                                                    >  1:2520823/1‑70 (MQ=255)
gACTGCGCTTAACGCTAACCGCTGCGCATGAAATGCAGGATATCGACCGTCTGCTGGAGGTGCTGCATGGc                                                                   <  1:3664353/71‑1 (MQ=255)
    gcgcTTAACGCTAACCGCTGCGCATGAAATGCAGGATATCGACCGTCTGCTGGAGGTGATGAATGGCAACg                                                               >  1:1614002/1‑71 (MQ=255)
     cgcTTAACGCTAACCGCTGCGCATGAAATGCAGGATATCGACCGTCTGCTGGAGGTGCTGCATGGCAACgg                                                              >  1:19956/1‑71 (MQ=255)
     cgcTTAACGCTAACCGCTGCGCATGAAATGCAGGATATCGACCGTCTGCTGGAGGTGCTGCATGGCAACgg                                                              >  1:3945641/1‑71 (MQ=255)
       cTTAACGCTAACCGCTGCGCATGAAATGCAGGATATCGACCGTCTGCTGGAGGTGCTGCATGGCAACGGt                                                             >  1:4199842/1‑70 (MQ=255)
       cTTAACGCTAACCGCTGCGCATGAAATGCAGGATATCGACCGTCTGCTGGAGGTGCTGCATGGCAACGGt                                                             >  1:3512302/1‑70 (MQ=255)
          aaCGCTAACCGCTGCGCATGAAATGCAGGATATCGACCGTCTGCTGGAGGTGCTGCATGGCAACGGTtaa                                                          >  1:2474120/1‑70 (MQ=255)
           aCGCTAACCGCTGCGCATGAAATGCAGGATATCGACCGTCTGCTGGAGGTGCTGCATGGCAACGGTtaata                                                        <  1:2939528/71‑1 (MQ=255)
           aCGCTAACCGCTGCGCATGAAATGCAGGATATCGACCGTCTGCTGGAGGTGCTGCATGGCAACGGTtaat                                                         <  1:1853794/70‑1 (MQ=255)
             gCTAACCGCTGCGCATGAAATGCAGGATATCGACCGTCTGCTGGAGGTGCTGCATGGCAACGGTTAATaaa                                                      <  1:3068751/71‑1 (MQ=255)
             gCTAACCGCTGCGCATGAAATGCAGGATATCGACCGTCTGCTGGAGGTGCTGCATGGCAACGGTTAATaaa                                                      <  1:4532462/71‑1 (MQ=255)
                aaCCGCTGCGCATGAAATGCAGGATATCGACCGTCTGCTGGAGGTGCTGCATGGCaa                                                                 >  1:967097/1‑57 (MQ=255)
                aaCCGCTGCGCATGAAATGCAGGATATCGACCGTCTGCTGGAGGTGCTGCATGGCaa                                                                 >  1:4738048/1‑57 (MQ=255)
                 aCCGCTGCGCATGAAATGCAGGATATCGACCGTCTGCTGGAGGTGCTGCATGTCaa                                                                 >  1:1911088/1‑56 (MQ=255)
                 aCCGCTGCGCATGAAATGCAGGATATCGACCGTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACa                                                    <  1:3457867/69‑1 (MQ=255)
                  ccGCTGCTCATGAAATGCAGGATATCGACCGTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGc                                                 <  1:2547142/71‑1 (MQ=255)
                   cGCTGCGCATGAAATGCAGGATATCGACCGTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGc                                                 >  1:2027567/1‑70 (MQ=255)
                        cgcATGAAATGCAGGATATCGACCGTCTGCTGGAGGTGCTGCATGGCAACGGTTAGTAAACAAGCCATTg                                            <  1:3587196/70‑1 (MQ=255)
                        cgcATGAAATGCAGGATATCGACCGTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTg                                            <  1:4366305/70‑1 (MQ=255)
                              aaaTGCAGGATATCGACCGTCTGCTGGAGGTGCTGCATGGCAACGGTTAATaaa                                                      >  1:4776199/1‑54 (MQ=255)
                                aTGCAGGATATCGACCGTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCAt                                   <  1:2301235/71‑1 (MQ=255)
                                aTGCAGGATATCGACCGTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCAt                                   <  1:1526200/71‑1 (MQ=255)
                                aTGCAGGATATCGACCGTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCAt                                   <  1:4039421/71‑1 (MQ=255)
                                   cAGGATATCGACCGTCTGCTGGAGGTGCTGCATGTCAACGGTTAATAAACAAGCCATTg                                            >  1:1029522/1‑59 (MQ=255)
                                   cAGGATATCGACCGTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTg                                >  1:4229150/1‑71 (MQ=255)
                                   cAGGATATCGACCGTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTg                                >  1:4534334/1‑71 (MQ=255)
                                    aGGATATCGACCGTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTgg                               <  1:667598/71‑1 (MQ=255)
                                     ggATATCGACCGTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGcc                                                <  1:4746170/53‑1 (MQ=255)
                                      gATATCGACCGTCTGCTGGAGGTGCTGCATGTCAACGGTtaataa                                                       <  1:3709286/45‑1 (MQ=255)
                                      gATATCGACCGTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTg                                <  1:293124/68‑1 (MQ=255)
                                       atatCGACCGTCTGCTGGAGGTGCTGCATGTCAACGGTtaataa                                                       <  1:1855611/44‑1 (MQ=255)
                                       atatCGACCGTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCAt                                              <  1:1206193/53‑1 (MQ=255)
                                        tatCGACCGTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTg                                <  1:1370903/66‑1 (MQ=255)
                                         atCGACCGTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCggg                          >  1:4548264/1‑71 (MQ=255)
                                          tCGACCGTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGc                         <  1:1615037/71‑1 (MQ=255)
                                            gACCGTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAAc                                                     >  1:1273236/1‑41 (MQ=255)
                                            gACCGTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCg                                       >  1:1107461/1‑55 (MQ=255)
                                            gACCGTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAg                       >  1:1368789/1‑71 (MQ=255)
                                               cGTCTGCTGGAGGTGCTGCATGTCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCg                    <  1:5017398/71‑1 (MQ=255)
                                               cGTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCg                    <  1:2456786/71‑1 (MQ=255)
                                               cGTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCg                    <  1:2846066/71‑1 (MQ=255)
                                                 tCTGCTGGAGGTGCTGCATGTCAACGGTTAATAAACAAGcc                                                <  1:5019753/41‑1 (MQ=255)
                                                 tCTGCTGGAGGTGCTGCATGTCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCg                            <  1:3102682/61‑1 (MQ=255)
                                                 tCTGCTGGAGGTGCTGCATGTCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCg                            <  1:4412860/61‑1 (MQ=255)
                                                   tgctgGAGGTGCTGCATGTCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGcaca                >  1:1986738/1‑71 (MQ=255)
                                                   tgctgGAGGTGCTGCATGTCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGcaca                >  1:1535427/1‑71 (MQ=255)
                                                      tgGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTa             >  1:2996330/1‑71 (MQ=255)
                                                       ggAGGTGCTGCATGGCAACGGTTAATAAACAAGCCAt                                              >  1:3956921/1‑37 (MQ=255)
                                                        gAGGTGCTGCATGGCAACGGTTAATAAACAAGCCa                                               <  1:1611792/35‑1 (MQ=255)
                                                        gAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTg                                            >  1:4071199/1‑38 (MQ=255)
                                                           gTGCTGCATGTCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTg                                <  1:3202133/47‑1 (MQ=255)
                                                           gTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAg         >  1:355294/1‑70 (MQ=255)
                                                            tgctgcATGTCAACGGTTAATAACCAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCa       <  1:2774922/71‑1 (MQ=255)
                                                            tgctgcATGTCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCa       <  1:3572144/71‑1 (MQ=255)
                                                            tgctgcATGTCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCa       <  1:1687792/71‑1 (MQ=255)
                                                              ctgcATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCaa      <  1:569095/70‑1 (MQ=255)
                                                                 cATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACAt   <  1:2623209/70‑1 (MQ=255)
                                                                  aTGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATg  >  1:854289/1‑70 (MQ=255)
                                                                   tGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTg                                <  1:3625577/39‑1 (MQ=255)
                                                                   tGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTg                                <  1:1875549/39‑1 (MQ=255)
                                                                    ggCAACGGTTAATAAACAAGCCATTGCAGTGGCAt                                   <  1:1355623/35‑1 (MQ=255)
                                                                     gCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAg                       >  1:780468/1‑46 (MQ=255)
                                                                     gCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATg  >  1:4747430/1‑67 (MQ=255)
                                                                     |                                                                  
GACTGCGCTTAACGCTAACCGCTGCGCATGAAATGCAGGATATCGACCGTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATG  >  W3110S.gb/811878‑812013

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