Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,897,706 A→G 9.7% G123G (GGT→GGC ygcU ← predicted FAD containing dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,897,7060AG 9.7% 99.5 / 3.9 41G123G (GGT→GGCygcUpredicted FAD containing dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (19/18);  new base G (2/2);  total (21/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01

TGACTGCGGAGAATGCCCCGTGGTGTAACCTTTTTCACGCAACGCGTTTTCCAGCACTTCCAGCGGAACACCACATTGCGCCGTCGCCTGCATATTCTCAATATCAATATTAATGATTTGATTCATGGCGGAGC  >  W3110S.gb/2897637‑2897770
                                                                     |                                                                
tGACTGCGGAGAATGCCCCGTGGTGTAACCTTTTTCACGCAACGCGTTTTCCAGCACTTCCAGCGGAaca                                                                  <  1:819952/70‑1 (MQ=255)
 gACTGCGGAGAATGCCCCGTGGTGTAACCTTTTTCACGCAACGCGTTTTCCAGCACTTCCAGCGGAacac                                                                 <  1:804918/70‑1 (MQ=255)
      cGGAGAATGCCCCGTGGTGTAACCTTTTTCACGCAACGCGTTTTCCAGCACTTCCAGCGGAACACcaca                                                             <  1:1810594/69‑1 (MQ=255)
      cGGAGAATGCCCCGTGGTGTAACCTTTTTCACGCAACGCGTTTTCCAGCACTTCCAGCGGAACACCACAtt                                                           <  1:3121492/71‑1 (MQ=255)
                cccGTGGTGTAACCTTTTTCACGCAACGCGTTTTCCAGCACCTCCAGCGGAACGCCACATTGCGCCGTGGc                                                 <  1:1161142/71‑1 (MQ=255)
                cccGTGGTGTAACCTTTTTCACGCAACGCGTTTTCCAGCACCTCCAGCGGAACGCCACATTGCGCCGTGGc                                                 <  1:260626/71‑1 (MQ=255)
                   gtggtgTAACCTTTTTCACGCAACGCGTTTTCCAGCACCTCCAGCGGAACGCCACATTGCGCCGTGGCCTg                                              >  1:29706/1‑71 (MQ=255)
                   gtggtgTAACCTTTTTCACGCAACGCGTTTTCCAGCACCTCCAGCGGAACGCCACATTGCGCCGTGGCCTg                                              >  1:222469/1‑71 (MQ=255)
                    tggtgTAACCTTTTTCACGCAACGCGTTTTCCAGCACTTCCAGCGGAACACCACATTGCGCCg                                                     <  1:4416041/63‑1 (MQ=255)
                    tggtgTAACCTTTTTCACGCAACGCGTTTTCCAGCACTTCCAGCGGAACACCACATTGCGCCGTCGCCTg                                              >  1:4347363/1‑70 (MQ=255)
                    tggtgTAACCTTTTTCACGCAACGCGTTTTCCAGCACTTCCAGCGGAACACCACATTGCGCCGTCGCCTg                                              >  1:1728696/1‑70 (MQ=255)
                      gtgtAACCTTTTTCACGCAACGCGTTTTCCAGCACTTCCAGCGGAACACCACATTGCGCCGTCGCCTGCa                                            <  1:2776690/70‑1 (MQ=255)
                         tAACCTTTTTCACGCAACGCGTTTTCCAGCACTTCCAGCGGAAcacca                                                               <  1:2008696/48‑1 (MQ=255)
                         tAACCTTTTTCACGCAACGCGTTTTCCAGCACTTCCAGCGGAAcacca                                                               <  1:1655881/48‑1 (MQ=255)
                         tAACCTTTTTCACGCAACGCGTTTTCCAGCACTTCCAGCGGAAcacca                                                               <  1:1610150/48‑1 (MQ=255)
                         tAACCTTTTTCACGCAACGCGTTTTCCAGCACTTCCAGCGGAACACCACATTGCGCCGTCGcc                                                >  1:3125274/1‑63 (MQ=255)
                           aCCTTTTTCACGCAACGCGTTTTCCAGCACTTCCAGCGGAacac                                                                 <  1:2119816/44‑1 (MQ=255)
                           aCCTTTTTCACGCAACGCGTTTTCCAGCACTTCCAGCGGAACACcaca                                                             >  1:556919/1‑48 (MQ=255)
                           aCCTTTTTCACGCAACGCGTTTTCCAGCACTTCCAGCGGAACACCACATTGCGCCGTCGCCTGCATATtct                                      >  1:2234318/1‑71 (MQ=255)
                           aCCTTTTTCACGCAACGCGTTTTCCAGCACTTCCAGCGGAACACCACATTGCGCCGTCGCCTGCATATtct                                      >  1:3024536/1‑71 (MQ=255)
                            ccTTTTTCACGCAACGCGTTTTCCAGCACTTCCAGCGGAACACCACATTGCGCCGTCGCCTGCATATtctc                                     <  1:3429951/71‑1 (MQ=255)
                             cTTTTTCACGCAACGCGTTTTCCAGCACTTCCAGCGGAACACCACATTGCGCCGTCGCCTg                                              <  1:2240662/61‑1 (MQ=255)
                             cTTTTTCACGCAACGCGTTTTCCAGCACTTCCAGCGGAACACCACATTGCGCCGTCGCCTGCATATtctc                                     >  1:2380698/1‑70 (MQ=255)
                              tttttCACGCAACGCGTTTTCCAGCACTTCCAGCGGAacac                                                                 >  1:827767/1‑41 (MQ=255)
                               ttttCACGCAACGCGTTTTCCAGCACTTCCAGCGGAACACCACATTGCGCCGTCGCCTGCa                                            >  1:487777/1‑61 (MQ=255)
                                tttCACGCAACGCGTTTTCCAGCACTTCCAGCGGAACACCACATTGCGCCGTCGCCTGCatat                                         >  1:1530306/1‑63 (MQ=255)
                                tttCACGCAACGCGTTTTCCAGCACTTCCAGCGGAACACCACATTGCGCCGTCGCCTGCatat                                         >  1:1122412/1‑63 (MQ=255)
                                  tCACGCAACGCGTTTTCCAGCACTTCCAGCGGAACACCACATTGCGCCGTCGCCTGCATATTCtcaatat                                <  1:898369/70‑1 (MQ=255)
                                   cACGCAACGCGTTTTCCAGCACTTCCAGCGGAACACCACATTGCGCCGTCGCCTGCATATTCtcaatatc                               >  1:2001845/1‑70 (MQ=255)
                                   cACGCAACGCGTTTTCCAGCACTTCCAGCGGAACACCACATTGCGCCGTCGCCTGCATATTCtcaatatc                               <  1:2701145/70‑1 (MQ=255)
                                       cAACGCGTTTTCCAGCACTTCCAGCGGAACACCACATTGCGCCGTCGCCTGCATATTCTcaatatcaatat                          <  1:4277807/71‑1 (MQ=255)
                                         aCGCGTTTTCCAGCACTTCCAGCGGAACACCACATTGCGCCGTCGCCTGCATATtctc                                     <  1:811008/58‑1 (MQ=255)
                                         aCGCGTTTTCCAGCACTTCAAGCGGAACACCACATTGCGCCGTCGCCTGCATATtctc                                     <  1:4077669/58‑1 (MQ=255)
                                          cgcgTTTTCCAGCACTTCCAGCGGAACACCACATTGCGCCGTCGCCTGCATATTCTCAATATCAATATTaa                       >  1:2601193/1‑71 (MQ=255)
                                            cgTTTTCCAGCACTTCCAGCGGAACACCACATTGCGCCGTCGCCTGCATATTCTCAATATCAATATTa                        <  1:3264145/68‑1 (MQ=255)
                                              ttttCCAGCACTTCCAGCGGAACACCACATTGCGCCGTCGCCTGCATATTCtcaatatcaat                            <  1:2322407/62‑1 (MQ=255)
                                                ttCCAGCACTTCCAGCGGAACACCACATTGCGCCGTCGCCTGCATATTCTCAATATCAATATTAAtgattt                 >  1:4326314/1‑71 (MQ=255)
                                                 tCCAGCACTTCCAGCGGAACACCACATTGCGCCGTCGCCTGCATATTCTCAATATCAATATTAAtgattt                 <  1:2878113/70‑1 (MQ=255)
                                                 tCCAGCAATTCCAGCGGAACACCACATTGCGCCGTCGCCTGCATATTCTCAATATCAATATTAAtgattt                 <  1:1234986/70‑1 (MQ=255)
                                                       aCTTCCAGCGGAACACCACATTGCGCCGTCGCCTGCATATTCTCAATATCAATATTAATGATTTGATTCAt          <  1:1797788/71‑1 (MQ=255)
                                                        cTTCCAGCGGAACACCACATTGCGCCGTCGCCTGCATATTCTCAATATCAATATTAAtgatttgatt             <  1:1582520/67‑1 (MQ=255)
                                                         ttCCAGCGGAACACCACATTGCGCCGTCTCCTGCATATTCTCAATATCAATATTAAtgatttg                <  1:2148540/63‑1 (MQ=255)
                                                          tCCAGCGGAACACCACATTGCGCCGTCGCCTGCata                                          >  1:3616281/1‑36 (MQ=255)
                                                          tCCAGCGGAACACCACATTGCGCCGTCGCCTGCata                                          >  1:3740625/1‑36 (MQ=255)
                                                           ccAGCGGAACACCACATTGCGCCGTCGCCTGCATATTCTCAATATCAATATTAATGATTTGATTCATGGCg      >  1:1565575/1‑71 (MQ=255)
                                                              gCGGAACACCACATTGCGCCGTCGCCTGCATATTCTCAATATCAATATTAATGATTTGATTCATGGCGGAg   >  1:3739125/1‑71 (MQ=255)
                                                               cGGGACACCACATTGCGCCGTCGCCTGCATATTCTCAATATCAATATTAATGATTTGATTCATGGCGGAGc  >  1:3360079/1‑71 (MQ=255)
                                                               cGGAACACCACATTGCGCCGTCGCCTGCATATTCTCAATATCAATATTAATGATTTGATTCATGGCGGAGc  >  1:1268521/1‑71 (MQ=255)
                                                                   acacCACATTGCGCCGTCGCCTGCATATTCtcaatatcaata                           >  1:3245571/1‑42 (MQ=255)
                                                                     |                                                                
TGACTGCGGAGAATGCCCCGTGGTGTAACCTTTTTCACGCAACGCGTTTTCCAGCACTTCCAGCGGAACACCACATTGCGCCGTCGCCTGCATATTCTCAATATCAATATTAATGATTTGATTCATGGCGGAGC  >  W3110S.gb/2897637‑2897770

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: