Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 647,863 A→G 9.1% 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,8630AG 9.1% 116.0 / 4.7 44M311T (ATG→ACG) citFcitrate lyase, citrate‑ACP transferase (alpha) subunit
Reads supporting (aligned to +/- strand):  ref base A (17/23);  new base G (2/2);  total (19/25)
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.82e-01

GCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAATATCGCGGCTACGCATTTTGTCTTCCAGGA  >  W3110S.gb/647794‑647932
                                                                     |                                                                     
gcagcaTGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGTTCAACCAt                                                                      <  1:3784602/71‑1 (MQ=255)
gcagcaTGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCAt                                                                      >  1:2285424/1‑71 (MQ=255)
  agcaTGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCAtcg                                                                    >  1:3305938/1‑71 (MQ=255)
  agcaTGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCAtcg                                                                    >  1:1666418/1‑71 (MQ=255)
        gCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCgg                                                              <  1:604675/71‑1 (MQ=255)
         cTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGt                                                             <  1:2164206/71‑1 (MQ=255)
          tGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGt                                                             >  1:2039768/1‑70 (MQ=255)
            tCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTaa                                                           >  1:720967/1‑70 (MQ=255)
            tCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTaa                                                           >  1:909119/1‑70 (MQ=255)
            tCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTaa                                                           <  1:1287654/70‑1 (MQ=255)
              aaaGCTCTGCAGATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTa                                                            >  1:1665193/1‑67 (MQ=255)
                  ctctGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCGTCGTCGCGGTAATACCGcc                                                    <  1:608713/71‑1 (MQ=255)
                     tGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAg                                                 >  1:2992163/1‑71 (MQ=255)
                      gCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGcc                                                    <  1:2868206/67‑1 (MQ=255)
                       cacaTCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCaa                                                  <  1:652245/68‑1 (MQ=255)
                       cacaTCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCaa                                                  <  1:734307/68‑1 (MQ=255)
                        acaTCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGc                                              <  1:3648255/71‑1 (MQ=255)
                           tCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATAcc                                                       <  1:1146306/59‑1 (MQ=255)
                            ccAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCAtc                                                                     >  1:2599201/1‑44 (MQ=255)
                              agcagTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAggg                                               >  1:4079433/1‑64 (MQ=255)
                               gcagTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCAtcgt                                                                   <  1:117714/43‑1 (MQ=255)
                               gcagTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCg                                       >  1:1204566/1‑71 (MQ=255)
                                cagTTTGCGGATCAGACTTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCgg                                      <  1:52426/71‑1 (MQ=255)
                                cagTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCAtcgt                                                                   <  1:934285/42‑1 (MQ=255)
                                cagTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCgg                                                              >  1:2099627/1‑47 (MQ=255)
                                     tGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGcgcg                                  >  1:3185145/1‑70 (MQ=255)
                                     tGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGcgcg                                  >  1:311483/1‑70 (MQ=255)
                                     tGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGcgcg                                  <  1:391719/70‑1 (MQ=255)
                                         gATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGaa                                           <  1:2084509/57‑1 (MQ=255)
                                          aTCAGACCTTTTTCGTGCAGGTCAACCGTCGTCGCGGt                                                             <  1:1230511/38‑1 (MQ=255)
                                          aTCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAATATc                            <  1:3694905/71‑1 (MQ=255)
                                           tCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAATATcg                           <  1:2823518/71‑1 (MQ=255)
                                            cAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAAt                                                          <  1:3226880/39‑1 (MQ=255)
                                                 cTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAATATCGCGGCTa                     <  1:1323471/71‑1 (MQ=255)
                                                   ttttCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGaa                                >  1:1599675/1‑58 (MQ=255)
                                                    tttCGTGCAGGTCAACCGTCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAATATCGCGGCTACg                   >  1:3318870/1‑70 (MQ=255)
                                                        gTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAATATCGCGGCTACGCAtt               <  1:642915/70‑1 (MQ=255)
                                                           cAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAATATCGCGGCTACGCATTTTg            <  1:1961871/70‑1 (MQ=255)
                                                               tCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAATATCGCGGCTACGCATTTTGTCTTc       <  1:379573/71‑1 (MQ=255)
                                                                cAACCGTCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAATATCGCGGCTACGCATTTTGTCTTc       >  1:1449713/1‑70 (MQ=255)
                                                                 aaCCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGTCGCGAATATCGCGGCTACGCATTTTGTCTTc       >  1:3626529/1‑69 (MQ=255)
                                                                 aaCCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAATATCGCGGCTACGCATTTTGTCTTcc      <  1:1550680/70‑1 (MQ=255)
                                                                 aaCCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAATATCGCGGCTACGCATTTTGTCTTc       >  1:272064/1‑69 (MQ=255)
                                                                    cATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAATATCGCGGCTACGCATTTTGTCTTCCAGGa  <  1:3889959/71‑1 (MQ=255)
                                                                     |                                                                     
GCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAATATCGCGGCTACGCATTTTGTCTTCCAGGA  >  W3110S.gb/647794‑647932

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: