Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 647,863 A→G 25.5% M311T (ATG→ACG)  citF ← citrate lyase, citrate‑ACP transferase (alpha) subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb647,8630AG25.5% 79.1 / 29.3 47M311T (ATG→ACG) citFcitrate lyase, citrate‑ACP transferase (alpha) subunit
Reads supporting (aligned to +/- strand):  ref base A (18/17);  new base G (2/10);  total (20/27)
Fisher's exact test for biased strand distribution p-value = 4.67e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAATATCGCGGCTACGCATTTTGTCTTCCAGGA  >  W3110S.gb/647795‑647932
                                                                    |                                                                     
cagcaTGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCAtc                                                                     <  1:3590723/71‑1 (MQ=255)
cagcaTGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCAtc                                                                     >  1:2320523/1‑71 (MQ=255)
  gcaTGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCAtcgt                                                                   >  1:1655783/1‑71 (MQ=255)
    aTGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCGtcgtcg                                                                 >  1:4183566/1‑71 (MQ=255)
    aTGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCAtcgtcg                                                                 >  1:1178802/1‑71 (MQ=255)
     tGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCAtcgtc                                                                  >  1:3842208/1‑69 (MQ=255)
     tGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTcgc                                                                >  1:2739157/1‑71 (MQ=255)
       gCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTcgcg                                                               <  1:1094397/70‑1 (MQ=255)
       gCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTcgcg                                                               <  1:3822582/70‑1 (MQ=255)
       gCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTcgcg                                                               <  1:4311493/70‑1 (MQ=255)
       gCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCgg                                                              >  1:3839951/1‑71 (MQ=255)
           tCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTaa                                                           <  1:4456902/70‑1 (MQ=255)
                 ctctGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGcc                                                    <  1:2180689/71‑1 (MQ=255)
                    tGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCa                                                   <  1:1430540/69‑1 (MQ=255)
                              gcagTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCaa                                                  >  1:2451247/1‑60 (MQ=255)
                              gcagTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCaa                                                  >  1:1400574/1‑60 (MQ=255)
                                  tttGCGGATCAGACCTTTTTCGTGCAGGTCAACCGTCGTCGCGGt                                                             <  1:1604637/45‑1 (MQ=255)
                                  tttGCGGATCAGACCTTTTTCGTGCAGGTCAACCGTCGTCGCGGt                                                             <  1:1121511/45‑1 (MQ=255)
                                    tGCGGATCAGACCTTTTTCGTGCAGGTCAACCGTCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGcgcg                                  <  1:4762047/70‑1 (MQ=255)
                                    tGCGGATCAGACCTTTTTCGTGCAGGTCAACCGTCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGcgcg                                  <  1:3972511/70‑1 (MQ=255)
                                    tGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTaa                                                           >  1:4772470/1‑45 (MQ=255)
                                    tGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAAt                                                          >  1:1489057/1‑46 (MQ=255)
                                    tGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGcgcg                                  >  1:2189854/1‑70 (MQ=255)
                                    tGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGcgcg                                  >  1:1960387/1‑70 (MQ=255)
                                    tGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGcgcg                                  >  1:2345680/1‑70 (MQ=255)
                                    tGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGcgcg                                  >  1:1718966/1‑70 (MQ=255)
                                     gCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGcgcg                                  <  1:2092024/69‑1 (MQ=255)
                                      cGGATCAGACCTTTTTCGTGCAGGTCAACCGTCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGaa                                <  1:4984017/70‑1 (MQ=255)
                                        gATCAGACCTTTTTCGTGCAGGTCAACCGTCGTCGCGGTa                                                            <  1:4743626/40‑1 (MQ=255)
                                        gATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAatat                             <  1:3625094/71‑1 (MQ=255)
                                        gATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAatat                             <  1:2641309/71‑1 (MQ=255)
                                           cAGACCTTTTTCGTGCAGGTCAACCGTCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAATATcgc                          >  1:3831733/1‑71 (MQ=255)
                                            aGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGt                                                             <  1:1410519/35‑1 (MQ=255)
                                             gACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAgg                                                >  1:3099939/1‑47 (MQ=255)
                                             gACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAgg                                                >  1:3435107/1‑47 (MQ=255)
                                             gACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAgg                                                >  1:4718718/1‑47 (MQ=255)
                                                  ttttCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGcgcg                                  <  1:3335140/56‑1 (MQ=255)
                                                  ttttCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAATATCGCGGCt                      <  1:3790173/68‑1 (MQ=255)
                                                        tGCAGGTCAACCGTCGTCGCGGTAATACCGCCAAGGGCGAAGTCGgcgct                                  <  1:3348175/50‑2 (MQ=255)
                                                         gCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAata                              >  1:1575180/1‑53 (MQ=255)
                                                          cAGGTCAACCGTCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAata                              <  1:3882266/52‑1 (MQ=255)
                                                            ggTCAACCGTCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAATATCGCGGCTACGCATTTTGTCt         <  1:3938551/71‑1 (MQ=255)
                                                             gTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAATATc                            <  1:1331977/51‑1 (MQ=255)
                                                              tCAACCGTCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAATATCGCGGCTACGCATTTTGTCTTc       <  1:2730510/71‑1 (MQ=255)
                                                                   cATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAATATCGCGGCTACGCATTTTGTCTTCCAGGa  <  1:2383577/71‑1 (MQ=255)
                                                                   cATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAATATCGCGGCTACGCATTTTGTCTTCCAGGa  <  1:2550288/71‑1 (MQ=255)
                                                                   cATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAATATCGCGGCTACGCATTTTGTCTTCCAGGa  <  1:3813398/71‑1 (MQ=255)
                                                                    aTCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAATATCGCGGCTAc                    >  1:4528504/1‑52 (MQ=255)
                                                                    aTCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAATATCGCGGCTAc                    >  1:3773932/1‑52 (MQ=255)
                                                                    |                                                                     
CAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAATATCGCGGCTACGCATTTTGTCTTCCAGGA  >  W3110S.gb/647795‑647932

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