Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,897,694 T→C 7.7% E127E (GAA→GAG ygcU ← predicted FAD containing dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,897,6940TC 7.7% 131.1 / 3.3 52E127E (GAA→GAGygcUpredicted FAD containing dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (18/30);  new base C (2/2);  total (20/32)
Fisher's exact test for biased strand distribution p-value = 6.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

CGCCAGCGGCTTTGACTGCGGAGAATGCCCCGTGGTGTAACCTTTTTCACGCAACGCGTTTTCCAGCACTTCCAGCGGAACACCACATTGCGCCGTCGCCTGCATATTCTCAATATCAATATTAATGATTTGATTCAT  >  W3110S.gb/2897625‑2897762
                                                                     |                                                                    
cgcCAGCGGCTTTGACTGCGGAGAATG‑‑CCGTGGTGTAACCTTTTTCACGCAACGCGTTTTCCAGCACTTcc                                                                   <  1:2313170/71‑1 (MQ=255)
 gcCAGCGGCTTTGACTGCGGAGAATGCCCCGTGGTGTAACCTTTTTCACGCAACGCGTTTTCCAGCACTTc                                                                    >  1:1281668/1‑71 (MQ=255)
  ccAGCGGCTTTGACTGCGGAGAATGCCCCGTGGTGTAACCTTTTTCACGCAACGCGTTTTCCAGCACttac                                                                   <  1:2960820/71‑3 (MQ=255)
  ccAGCGGCTTTGACTGCGGAGAATGCCCCGTGGTGTAACCTTTTTCACGCAACGCGTTTTCCAGCACTTc                                                                    >  1:1382601/1‑70 (MQ=255)
    aGCGGCTTTGACTGCGGAGAATGCCCCGTGGTGTAACCTTTTTCACGCAACGCGTTTTCCAGCACTTCCAg                                                                 <  1:1698583/71‑1 (MQ=255)
     gCGGCTTTGACTGCGGAGAATGCCCCGTGGTGTAACCTTTTTCACGCAACGCGTTTTCCAGCACTTCCAga                                                                <  1:148446/71‑2 (MQ=255)
     gCGGCTTTGACTGCGGAGAATGCCCCGTGGTGTAACCTTTTTCACGCAACGCGTTTTCCAGCACTTCCAg                                                                 <  1:866969/70‑1 (MQ=255)
       ggCTTTGACTGCGGAGAATGCCCCGTGGTGTAACCTTTTTCACGCAACGCGTTTTCCAGCACTTCCAGCgg                                                              <  1:707822/71‑1 (MQ=255)
       ggCTTTGACTGC‑‑AGAATGCCCCGTGGTGTAACCTTTTTCACGCAACGCGTTTTCCAGCACTTCCAGCGGa                                                             >  1:1623171/1‑70 (MQ=255)
            tGACTGCGGAGAATGCCCCGTGGTGTAACCTTTTTCACGCAACGCGTTTTCCAGCACTTCCAg                                                                 <  1:1052378/63‑1 (MQ=255)
            tGACTGCGGAGAATGCCCCGTGGTGTAACCTTTTTCACGCAACGCGTTTTCCAGCACTTCCAGCGGAaca                                                          <  1:819952/70‑1 (MQ=255)
             gACTGCGGAGAATGCCCCGTGGTGTAACCTTTTTCACGCAACGCGTTTTCCAGCACTTCCAGCGGAacac                                                         <  1:804918/70‑1 (MQ=255)
              aCTGCGGAGAATGCCCCGTGGTGTAACCTTTTTCACGCAACGCGTTTTCCAGCACTTCCAGCg                                                               >  1:1289295/1‑63 (MQ=255)
                tGCGGAGAATGCCCCGTGGTGTAACCTTTTTCACGCAACGCGTTTTCCAGCACTTCCAGCGGaa                                                            <  1:2213778/64‑1 (MQ=255)
                  cGGAGAATGCCCCGTGGTGTAACCTTTTTCACGCAACGCGTTTTCCAGCACTTCCAGCGGAACACcaca                                                     <  1:1810594/69‑1 (MQ=255)
                  cGGAGAATGCCCCGTGGTGTAACCTTTTTCACGCAACGCGTTTTCCAGCACTTCCAGCGGAACACCACAtt                                                   <  1:3121492/71‑1 (MQ=255)
                            cccGTGGTGTAACCTTTTTCACGCAACGCGTTTTCCAGCACCTCCAGCGGAACGCCACATTGCGCCGTGGc                                         <  1:260626/71‑1 (MQ=255)
                            cccGTGGTGTAACCTTTTTCACGCAACGCGTTTTCCAGCACCTCCAGCGGAACGCCACATTGCGCCGTGGc                                         <  1:1161142/71‑1 (MQ=255)
                               gtggtgTAACCTTTTTCACGCAACGCGTTTTCCAGCACCTCCAGCGGAACGCCACATTGCGCCGTGGCCTg                                      >  1:222469/1‑71 (MQ=255)
                               gtggtgTAACCTTTTTCACGCAACGCGTTTTCCAGCACCTCCAGCGGAACGCCACATTGCGCCGTGGCCTg                                      >  1:29706/1‑71 (MQ=255)
                                tggtgTAACCTTTTTCACGCAACGCGTTTTCCAGCACTTCCAGCGGAACACCACATTGCGCCg                                             <  1:4416041/63‑1 (MQ=255)
                                tggtgTAACCTTTTTCACGCAACGCGTTTTCCAGCACTTCCAGCGGAACACCACATTGCGCCGTCGCCTg                                      >  1:1728696/1‑70 (MQ=255)
                                tggtgTAACCTTTTTCACGCAACGCGTTTTCCAGCACTTCCAGCGGAACACCACATTGCGCCGTCGCCTg                                      >  1:4347363/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:1122412/1‑63 (MQ=255)
                                            tttCACGCAACGCGTTTTCCAGCACTTCCAGCGGAACACCACATTGCGCCGTCGCCTGCatat                                 >  1:1530306/1‑63 (MQ=255)
                                              tCACGCAACGCGTTTTCCAGCACTTCCAGCGGAACACCACATTGCGCCGTCGCCTGCATATTCtcaatat                        <  1:898369/70‑1 (MQ=255)
                                               cACGCAACGCGTTTTCCAGCACTTCCAGCGGAACACCACATTGCGCCGTCGCCTGCATATTCtcaatatc                       <  1:2701145/70‑1 (MQ=255)
                                               cACGCAACGCGTTTTCCAGCACTTCCAGCGGAACACCACATTGCGCCGTCGCCTGCATATTCtcaatatc                       >  1:2001845/1‑70 (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)
                                                                     |                                                                    
CGCCAGCGGCTTTGACTGCGGAGAATGCCCCGTGGTGTAACCTTTTTCACGCAACGCGTTTTCCAGCACTTCCAGCGGAACACCACATTGCGCCGTCGCCTGCATATTCTCAATATCAATATTAATGATTTGATTCAT  >  W3110S.gb/2897625‑2897762

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: