Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,940,955 T→C 10.5% D318G (GAT→GGT)  lpxM ← myristoyl‑acyl carrier protein (ACP)‑dependent acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,940,9550TC10.5% 97.0 / 4.0 38D318G (GAT→GGT) lpxMmyristoyl‑acyl carrier protein (ACP)‑dependent acyltransferase
Reads supporting (aligned to +/- strand):  ref base T (17/17);  new base C (2/2);  total (19/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTGAACTTATCATCAGGCGAAGGCCTCTCCTCGCGAGAGGCTTTTTTATTTGATGGGATAAAGATCTTTGCGCTTATACGGCTGGATTTCGCCCGGTTTGCGAGTTTTCAGCAATTTTAGTATCCAGGTGTA  >  W3110S.gb/1940890‑1941022
                                                                 |                                                                   
cTTGAACTTATCATCAGGCGAAGGCCTCTCCTCGCGAGAGGCTTTTTTATTTGATGGGATAAAGATCTTTg                                                                <  1:2975723/71‑1 (MQ=255)
  tGAACTTATCATCAGGCGAAGGCCTCTCCTCGCGAGAGGCTTTTTTATTTGATGGGATAAAGATCTTTgc                                                               >  1:1772525/1‑70 (MQ=255)
   gAACTTATCATCAGGCGAAGGCCTCTCCTCGCGAGAGGCTTTTTTATTTGATGGGATAAAGATCTTTgcgc                                                             >  1:4367610/1‑71 (MQ=255)
    aaCTTATCATCAGGCGAAGGCCTCTCCTCGCGAGAGGCTTTTTTATTTGATGGGATAAAGATCTTTGCGCt                                                            <  1:2342423/71‑1 (MQ=255)
      cTTATCATCAGGCGAAGGCCTCTCCTCGCGAGAGGCTTTTTTATTTGATGGGATAAAGATCTTTGCGCTta                                                          <  1:4350450/71‑1 (MQ=255)
      cTTATCATCAGGCGAAGGCCTCTCCTCGCGAGAGGCTTTTTTATTTGATGGGATAAAGATCTTTGCGCTta                                                          <  1:2416615/71‑1 (MQ=255)
           catcaGGCGAAGGCCTCTCCTCGCGAGAGGCTTTTTTATTTGATGGGATAAAGATCTTTGCGCTTATACgg                                                     >  1:3738599/1‑71 (MQ=255)
                ggCGAAGGCCTCTCCTCGCGAGAGGCTTTTTTATTTGATGGGATAAAGATCTTTGCGCTTATACGGCTGGa                                                <  1:3568309/71‑1 (MQ=255)
                         ctctCCTCGCGAGAGGCTTTTTTATTTGATGGGATAAAGATCTTTGCGCTTATACGGCTGGATTTCGCCCg                                       >  1:2600575/1‑71 (MQ=255)
                         ctctCCTCGCGAGAGGCTTTTTTATTTGATGGGATAAAGACCTTTGCGCTta                                                          >  1:2627393/1‑52 (MQ=255)
                          tctcCTCGCGAGAGGCTTTTTTATTTGATGGGATAAAGATCTTTGCGCTTATACGGCTGGATTTCGcc                                         <  1:2719635/68‑1 (MQ=255)
                            tcctcGCGAGAGGCTTTTTTATTTGATGGGATAAAGATCt                                                                   <  1:1216016/40‑1 (MQ=255)
                            tcctcGCGAGAGGCTTTTTTATTTGATGGGATAAAGATCt                                                                   >  1:1348827/1‑40 (MQ=255)
                            tcctcGCGAGAGGCTTTTTTATTTGATGGGATAAAGATCTTTGCGCTtata                                                        >  1:4115050/1‑51 (MQ=255)
                            tcctcGCGAGAGGCTTTTTTATTTGATGGGATAAAGATCTTTGCGCTTATACGGCTGGATTTCGCCCGGta                                    >  1:1644980/1‑70 (MQ=255)
                              ctcGCGAGAGGCTTTTTTATTTGATGGGATAAAGATCTTTGCGCTTATACGGCTGGATTTCGCCCGGTTTg                                  >  1:2645904/1‑71 (MQ=255)
                              ctcGCGAGAGGCTTTTTTATTTGATGGGATAAAGACCTTTGCGCTTATACGGCTGGAt                                               <  1:909087/58‑1 (MQ=255)
                              ctcGCGAGAGGCTTTTTTATTTGATGGGATAAAGACCTTTGCGCTTATACGGCTGGAt                                               <  1:2770012/58‑1 (MQ=255)
                                 gcgAGAGGCTTTTTTATTTGATGGGATAAAGATCTTTGCGCt                                                            >  1:1043868/1‑42 (MQ=255)
                                 gcgAGAGGCTTTTTTATTTGATGGGATAAAGATCTTTGCGCTTATACGGCTGGa                                                >  1:4007530/1‑54 (MQ=255)
                                 gcgAGAGGCTTTTTTATTTGATGGGATAAAGATCTTTGCGCTTATACGGCTGGATTTCGCCCGGTTTGc                                 >  1:1362073/1‑69 (MQ=255)
                                    agagGCTTTTTTATTTGATGGGATAAAGATCTTTGCGCTTATACGGCTGGATTTCGCCCGGTTTGCGAGtt                            <  1:784187/71‑1 (MQ=255)
                                     gagGCTTTTTTATTTGATGGGATAAAGACCTTTGCGCTTATACGGCTGGATTTCGCCCGGTTTGCGAGttt                           >  1:2021491/1‑71 (MQ=255)
                                      agGCTTTTTTATTTGATGGGATAAAGATCTTTGCGCTTATACGGCTg                                                  <  1:1391524/47‑1 (MQ=255)
                                          ttttttATTTGATGGGATAAAGATCTTTGCGCTTATACGGCTGGATTTCGCCCGGTTTGCGAGTTTTCAg                       >  1:3806414/1‑70 (MQ=255)
                                           tttttATTTGATGGGATAAAGATCTTTGCGCTTATACGGCTGGATTTCGCCCGGTTTGCGAGtt                            <  1:1871973/64‑1 (MQ=255)
                                            ttttATTTGATGGGATAAAGATCTTTGCGCTTATACGGCTGGATTTCGCCCGGTTTGCGAGtt                            <  1:99859/63‑1 (MQ=255)
                                            ttttATTTGATGGGATAAAGATCTTTGCGCTTATACGGCTGGATTTCGCCCGGTTTGCGAGTTTTCAGCaa                    >  1:2222801/1‑71 (MQ=255)
                                                aTTTGATGGGATAAAGATCTTTGCGCTTATACGGCTGGATTTCGCCCGGTTTGCGAGTTTTCAGCAAttt                 <  1:3176011/70‑1 (MQ=255)
                                                aTTTGATGGGATAAAGATCTTTGCGCTTATACGGCTGGATTTCGCCCGGTTTGCGAGTTTTCAGCAAttt                 <  1:731009/70‑1 (MQ=255)
                                                    gATGGGATAAAGATCTTTGCGCTTATACGGCTGGATTTCGCCCGGTTTGCGAGTTTTCAGCAATTTTAGTa            >  1:2301253/1‑71 (MQ=255)
                                                    gATGGGATAAAGATCTTTGCGCTTATACGGCTGGATTTCGCCCGGTTTGCGAGTATTCAGCAATTTTAGt             >  1:3702145/1‑70 (MQ=255)
                                                       gggATAAAGATCTTTGCGCTTATACGGCTGGATTTCGCCCGGTTTGc                                 <  1:1947336/47‑1 (MQ=255)
                                                        ggATAAAGATCTTTGCGCTTATAAGGCTGGATTTCGCCCGGTTTGCGAGTTTTCAGCAATTTTAGTATct         <  1:3072711/70‑2 (MQ=255)
                                                          aTAAAGATCTTTGCGCTTATACGGCTGGATTTCGCCCTGTTTGCGAg                              <  1:4340374/47‑1 (MQ=255)
                                                             aaGATCTTTGCGCTTATACGGCTGGATTTCGCCCGGTTTGCGAGTTTTCAGCAATTTTAGTATCCAg       >  1:2942390/1‑67 (MQ=255)
                                                              aGATCTTTGCGCTTATACGGCTGGATTTCGCCCGGTTTGCGAGTTTTCAGCAATTTTAGTATcc         <  1:3353166/64‑1 (MQ=255)
                                                              aGATCTTTGCGCTTATACGGCTGGATTTCGCCCGGTTTGCGAGTTTTCAGCAATTTTAGTATCCAGGTGTa  >  1:2333613/1‑71 (MQ=255)
                                                                 tCTTTGCGCTTATACGGCTGGATTTCGCCCGGTTTGCGAGTTTTCAGCAAt                   <  1:1320936/51‑1 (MQ=255)
                                                                 |                                                                   
CTTGAACTTATCATCAGGCGAAGGCCTCTCCTCGCGAGAGGCTTTTTTATTTGATGGGATAAAGATCTTTGCGCTTATACGGCTGGATTTCGCCCGGTTTGCGAGTTTTCAGCAATTTTAGTATCCAGGTGTA  >  W3110S.gb/1940890‑1941022

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: