Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,894,340 C→T 73.7% W32* (TGG→TGA ygcQ ← predicted flavoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,894,3400CT73.7% 56.2 / 24.7 38W32* (TGG→TGAygcQpredicted flavoprotein
Reads supporting (aligned to +/- strand):  ref base C (8/2);  new base T (16/12);  total (24/14)
Fisher's exact test for biased strand distribution p-value = 2.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.40e-01

GCTGCCACTGTTCCACCAACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTC  >  W3110S.gb/2894276‑2894403
                                                                |                                                               
gCTGCCACTGTTCCACCAACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGTCAGTgc                                                           <  1:2669828/71‑1 (MQ=255)
 cTGCCACTGTTCCACCAACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGTCAGTgcg                                                          <  1:1275395/71‑1 (MQ=255)
    ccACTGTTCCACCAACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGTCAGTGCGcc                                                        <  1:1572057/70‑1 (MQ=255)
            ccaccaACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGTCAGTGCGCCAGGGt                                                   >  1:783059/1‑67 (MQ=255)
                   cGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGTCAGTGCGCCAGGGTACAACCGCTGa                                        <  1:851064/71‑1 (MQ=255)
                     cATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGTCAGTGCGCCAGGGTACAACCGCTGaaa                                      >  1:2123543/1‑71 (MQ=255)
                       tCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGTCAGTGCGCCAGGGTACAACCGCTGAAATc                                    <  1:1974045/71‑1 (MQ=255)
                         cAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGTCAGTGCGCCAGGGTACAACCGCTGAAATCCt                                  >  1:712472/1‑71 (MQ=255)
                         cAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGTCAGTGCGCCAGGGTACAACCGCTGAAATCCt                                  >  1:2602436/1‑71 (MQ=255)
                         cAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGTCAGTGCGCCAGGGTACAACCGCTGAAATCCt                                  >  1:2599247/1‑71 (MQ=255)
                         cAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAg                                                              >  1:563881/1‑43 (MQ=255)
                         cAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACaa                                               >  1:2781829/1‑58 (MQ=255)
                         cAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACaa                                               >  1:1768344/1‑58 (MQ=255)
                           gCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGc                                           >  1:2672667/1‑60 (MQ=255)
                            cACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGTCAGTGCGCCAGGGTACAACCGCTGAAATCCTgcg                               >  1:395/1‑71 (MQ=255)
                            cACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTgcg                               >  1:1363587/1‑71 (MQ=255)
                              ccTGTTCCGCCACTACAGGCTGAGGTTCAATCTGTCa                                                               >  1:153145/1‑37 (MQ=255)
                              ccTGTTCCGCCACTACAGGCTGAGGTTCAATCTGTCa                                                               >  1:628936/1‑37 (MQ=255)
                              ccTGTTCCGCCACTACAGGCTGAGGTTCAATCTGTCa                                                               >  1:696850/1‑37 (MQ=255)
                               cTGTTCCGCCACTACAGGCTGAGGTTCAATCTGTCAGTGCGCCAGGGTACAACCGCTGaa                                       <  1:466385/60‑1 (MQ=255)
                                tGTTCCGCCACTACAGGCTGAGGTTCAATCTGTCAGTGCGCCAGGGTACAACCGCTGAAATCCTGcgccgc                           >  1:2264604/1‑71 (MQ=255)
                                    ccGCCACTACAGGCTGAGGTTCAATCTGTCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCgccagc                       >  1:2296695/1‑71 (MQ=255)
                                      gCCACTACAGGCTGAGGTTCAATCTGTCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCgccagcca                     <  1:1647052/71‑1 (MQ=255)
                                       ccACTACAGGCTGAGGTTCAATCTGTCAGTGCGCCAGGGTACGACCGCTGAAATCCTGCGCCGccagccag                    <  1:1816354/71‑1 (MQ=255)
                                       ccACTACAGGCTGAGGTTCAATCTGTCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGccagccag                    <  1:2775445/71‑1 (MQ=255)
                                            aCAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATc                                    >  1:2202147/1‑50 (MQ=255)
                                             cAGGCTGAGGTTCAATCTGTCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGc               >  1:1096917/1‑70 (MQ=255)
                                             cAGGCTGAGGTTCAATCTGTCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGc               >  1:1392198/1‑70 (MQ=255)
                                             cAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGc               >  1:2597115/1‑70 (MQ=255)
                                             cAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCt              <  1:1607316/71‑1 (MQ=255)
                                                gCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTAt            <  1:444861/70‑1 (MQ=255)
                                                    aGGTTCAATCTGTCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCgg       <  1:1787372/71‑1 (MQ=255)
                                                      gTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTa             >  1:803415/1‑63 (MQ=255)
                                                       ttCAATCTGTCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCt              <  1:246450/61‑1 (MQ=255)
                                                        tCAATCTGTCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCAtt   >  1:582452/1‑71 (MQ=255)
                                                        tCAATCTGTCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCAtt   >  1:777875/1‑71 (MQ=255)
                                                         cAATCTGTCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATtc  >  1:1514126/1‑71 (MQ=255)
                                                           aTCTGTCAGTGCGCCAGGGTACGACCGCTGAAATCCTGCGCCGccagccag                    <  1:1819651/51‑1 (MQ=39)
                                                                |                                                               
GCTGCCACTGTTCCACCAACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTC  >  W3110S.gb/2894276‑2894403

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: