Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 809,043 T→C 13.0% I213V (ATT→GTT)  bioA ← 7,8‑diaminopelargonic acid synthase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb809,0430TC13.0% 116.7 / 10.1 46I213V (ATT→GTT) bioA7,8‑diaminopelargonic acid synthase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base T (24/16);  new base C (2/4);  total (26/20)
Fisher's exact test for biased strand distribution p-value = 3.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.50e-01

TTTTGCGGATTCGTTTTAACCATTCCGGATGGTACATGCGCATCCCGCCTGCGCCCTGGACAATCGGCTCAATGATCACCGCCGCGATTTCATGACGATGCGCCGCCATCAGGCGGGCAAAGCCCACCA  >  W3110S.gb/808980‑809108
                                                               |                                                                 
ttttGCGGATTCGTTTTAACCATTCCGGATGGTACATGCGCATCCCGCCTGCGCCCTGGACAATCGGCTc                                                             >  1:3452447/1‑70 (MQ=255)
ttttGCGGATTCGTTTTAACCATTCCGGATGGTACATGCGCATCCCGCCTGCGCCCTGGACAATCGGCTc                                                             >  1:191550/1‑70 (MQ=255)
ttttGCGGATTCGTTTTAACCATTCCGGATGGTACATGCGCATCCCGCCTGCGCCCTGGACAATCGGCTc                                                             >  1:4004673/1‑70 (MQ=255)
ttttGCGGATTCGTTTTAACCATTCCGGATGGTACATGCGCATCCCGCCTGCGCCCTGGACAATCGGCTCa                                                            <  1:1510808/71‑1 (MQ=255)
ttttGCGGATTCGTTTTAACCATTCCGGATGGTACATGCGCATCCCGCCTGCGCCCTGGACAACCGGCTc                                                             >  1:3606619/1‑70 (MQ=255)
     cGGATTCGTTTTAACCATTCCGGATGGTACATGCGCATCCCGCCTGCGCCCTGGACAATCGGCTCAATGAt                                                       <  1:1042654/71‑1 (MQ=255)
     cGGATTCGTTTTAACCATTCCGGATGGTACATGCGCATCCCGCCTGCGCCCTGGACAATCGGCTCAATGAt                                                       <  1:384808/71‑1 (MQ=255)
       gATTCGTTTTAACCATTCCGGATGGTACATGCGCATCCCGCCTGCGCCCTGGACAATCGGCTCAATGATCa                                                     >  1:1322561/1‑71 (MQ=255)
       gATTCGTTTTAACCATTCCGGATGGTACATGCGCATCCCGCCTGCGCCCTGGACAATCGGCTCAATGATCa                                                     >  1:3045451/1‑71 (MQ=255)
          tCGTTTTAACCATTCCGGATGGTACATGCGCATCCCGCCTGCGCCCTGGACAATCGGCTCAATGATCAccg                                                  >  1:3854197/1‑71 (MQ=255)
          tCGTTTTAACCATTCCGGATGGTACATGCGCATCCCGCCTGCGCCCTGGACAATCGGCTCAATGATCAcc                                                   >  1:3760793/1‑70 (MQ=255)
          tCGTTTTAACCATTCCGGATGGTACATGCGCATCCCGCCTGCGCCCTGGACAACCGGCTCAATGATCAccg                                                  >  1:2335891/1‑71 (MQ=255)
               ttAACCATTCCGGATGGTACATGCGCATCCCGCCTGCGCCCTGGACAATCGGCTCAATGATCACcgccgc                                              <  1:1055558/70‑1 (MQ=255)
               ttAACCATTCCGGATGGTACATGCGCATCCCGCCTGCGCCCTGGACAATCGGCTCAATGATCACcgccgc                                              >  1:2412359/1‑70 (MQ=255)
                tAACCATTCCGGATGGTACATGCGCATCCCGCCTGCGCCCTGGACAATCGGCTCAATGATCACCGCCGCGa                                            >  1:1969089/1‑71 (MQ=255)
                   ccATTCCGGATTGTACATGCGCATCCCGCCTGCGCCCTGGACAACCGGCTCAATGATCACCGCCGCGAttt                                         <  1:2045376/71‑1 (MQ=255)
                   ccATTCCGGATGGTACATGCGCATCCCGCCTGCGCCCTGGACAATCGGCTCAATGATCACCGCCTCGAttt                                         <  1:528184/71‑1 (MQ=255)
                   ccATTCCGGATGGTACATGCGCATCCCGCCTGCGCCCTGGACAACCGGCTCAATGATCACCGCCGCGAttt                                         <  1:1008590/71‑1 (MQ=255)
                        ccGGATGGTACATGCGCATCCCGCCTGCGCCCTGGACAATCGGCTCAATGATCACCGCCGCGATTTCATg                                     >  1:623425/1‑70 (MQ=255)
                        ccGGATGGTACATGCGCATCCCGCCTGCGCCCTGGACAATCGGCTCAATGATCACCGCCGCGATTTCATGa                                    >  1:2600787/1‑71 (MQ=255)
                        ccGGATGGTACATGCGCATCCCGCCTGCGCCCTGGACAATCGGCTCAATGATCACCGCCGCGATTTCATGa                                    >  1:3001068/1‑71 (MQ=255)
                        ccGGATGGTACATGCGCATCCCGCCTGCGCCCTGGACAATCGGCTCAATGATCACCGCCGCGATTTCATGa                                    >  1:2854215/1‑71 (MQ=255)
                         cGGATGGTACATGCGCATCCCGCCTGCGCCCTGGACAATCGGCTCAATGATCACCGCCGCGATTTCATGAc                                   <  1:1320577/71‑1 (MQ=255)
                         cGGATGGTACATGCGCATCCCGCCTGCGCCCTGGACAATCGGCTCAATGATCACCGCCGCGATTTCATGAc                                   <  1:2664342/71‑1 (MQ=255)
                         cGGATGGTACATGCGCATCCCGCCTGCGCCCTGGACAATCGGCTCAATGATCACCGCCGCGATTTCATGAc                                   <  1:3973240/71‑1 (MQ=255)
                         cGGATGGTACATGCGCATCCCGCCTGCGCCCTGGACAATCGGCTCAATGATCACCGCCGCGATTTCATGAc                                   <  1:1601465/71‑1 (MQ=255)
                          ggATGGTACATGCGCATCCCGCCTGCGCCCTGGACAATCGGCTCAATGATCACCGCCGCGATTTCATGAc                                   <  1:193385/70‑1 (MQ=255)
                                  cATGCGCATCCCGCCTGCGCCCTGGACAATCGGCTCAATGATCACCGCCGCGATTTCATGAc                                   <  1:1843379/62‑1 (MQ=255)
                                     gcgcATCCCGCCTGCGCCCTGGACAATCGGCTCAATGATCACCGCCGCGATTTCATGACGATGCgccgcc                        >  1:2266210/1‑70 (MQ=255)
                                      cgcATCCCGCCTGCGCCCTGGACAATCGGCTCAATGATCACCGCCGCGATTTCATGACGa                                 >  1:1953269/1‑60 (MQ=255)
                                            ccGCCTGCGCCCTGGACAATCGGCTCAATGATCAc                                                    >  1:3036071/1‑35 (MQ=255)
                                            ccGCCTGCGCCCTGGACAATCGGCTCAATGATCAc                                                    >  1:612748/1‑35 (MQ=255)
                                            ccGCCTGCGCCCTGGACAATCGGCTCAATGATCACCGCCGCGATTTCa                                       >  1:3261802/1‑48 (MQ=255)
                                                cTGCGCCCTGGACAACCGGCTCAATGATCACCGCCTCGATTTCATGACGATGCGCCGCCATCAGGCGGGCa            <  1:2320059/71‑1 (MQ=255)
                                                cTGCGCCCTGGACAACCGGCTCAATGATCACCGCCGCGATTTCATGACGATGCGCCGCCATCAGGCGGCCa            <  1:1464848/71‑1 (MQ=255)
                                                 tGCGCCCTGGACAATCGGCTCAATGATCACCGCCGCGATTTCATTACGATGCGCCGCCATCAGGCGGGCa            >  1:3929371/1‑70 (MQ=255)
                                                  gcgcCCTGGACAATCGGCTCAATGATCACCGCCGCGa                                            >  1:1256247/1‑37 (MQ=255)
                                                  gcgcCCTGGACAATCGGCTCAATGATCACCGCCGCGAttt                                         <  1:1708254/40‑1 (MQ=255)
                                                   cgcCCTGGACAATCGGCTCAATGATCACCGCCGCGATTTCATGACGATGCGCCGCCATCAGGCGGGCAAAg         >  1:2839776/1‑71 (MQ=255)
                                                          gACAATCGGCTCAATGATCACCGCCGCGATTTCATGACGATGCGCCGCCATCAGGCGGGCAAAGCccacca  >  1:2948519/1‑71 (MQ=255)
                                                          gACAATCGGCTCAATGATCACCGCCGCGATTTCATGACGATGCGCCGCCATCAGGCGGGCAAAGCccacca  >  1:474406/1‑71 (MQ=255)
                                                          gACAATCGGCTCAATGATCACCGCCGCGATTTCATGACGATGCGCCGCCATCAGGCCGGCAAAGCccac    <  1:232698/69‑1 (MQ=255)
                                                          gACAATCGGCTCAATGATCACCGCCGCGATTTCATGACGATGCGCCGCCATAAGGCGGGCAAAGCccacca  >  1:24942/1‑71 (MQ=255)
                                                              aTCGGCTCAATGATCACCGCCGCGATTTCATGACGATgc                              <  1:2143463/39‑1 (MQ=255)
                                                              aTCGGCTCAATGATCACCGCCGCGATTTCATGACGATgc                              <  1:762648/39‑1 (MQ=255)
                                                              aTCGGCTCAATGATCACCGCCGCGATTTCATGAAGATGCGCCGCCATCa                    <  1:444095/49‑1 (MQ=255)
                                                               |                                                                 
TTTTGCGGATTCGTTTTAACCATTCCGGATGGTACATGCGCATCCCGCCTGCGCCCTGGACAATCGGCTCAATGATCACCGCCGCGATTTCATGACGATGCGCCGCCATCAGGCGGGCAAAGCCCACCA  >  W3110S.gb/808980‑809108

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: