Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 978,430 A→G 26.5% E561E (GAA→GAG mukB → fused chromosome partitioning proteins

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb978,4300AG26.5% 53.6 / 21.7 34E561E (GAA→GAGmukBfused chromosome partitioning proteins
Reads supporting (aligned to +/- strand):  ref base A (12/13);  new base G (4/5);  total (16/18)
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

GTCAGGGCAAGAATTTTGATATCGACGAACTGGAAGCCCTGCATCAGGAACTGGAAGCACGCATTGCCTCTCTTTCCGATAGCGTGTCTAACGCCCGTGAAGAGCGCATGG  >  W3110S.gb/978381‑978491
                                                 |                                                             
gTCAGGGCAAGAATTTTGATATCGACGAACTGGAAGCCCTGCATCAGGAACTGGAAGCACGCATTGCctct                                          >  1:1092460/1‑71 (MQ=255)
gTCAGGGCAAGAATTTTGATATCGACGAACTGGAAGCCCTGCATCAGGAACTGGAAGCACGCATTGCctct                                          >  1:1644578/1‑71 (MQ=255)
gTCAGGGCAAGAATTTTGATATCGACGAACTGGAAGCCCTGCATCAGGAACTGGAAGCACGCATTGCctc                                           >  1:329053/1‑70 (MQ=255)
                  atatCGACGAACTGGAAGCCCTGCATCAGGAACTGGAAGCACGCATTGCCTCTCTTTCCGATAGCGTGTCt                        <  1:519522/71‑1 (MQ=255)
                  atatCGACGAACTGGAAGCCCTGCATCAGGAACTGGAAGCACGCATTGCCTCTCTTTCCGATAGCGTGTCt                        <  1:858278/71‑1 (MQ=255)
                  atatCGACGAACTGGAAGCCCTGCATCAGGAACTGGAAGCACGCATTGCCTCTCTTTCCGATAGCGTGTCt                        <  1:242406/71‑1 (MQ=255)
                  atatCGACGAACTGGAAGCCCTGCATCAGGAACTGGAAGCACGCATTGCCTCTCTTTCCGATAGCGTGTCt                        <  1:259917/71‑1 (MQ=255)
                  atatCGACGAACTGGAAGCCCTGCATCAGGAACTGGAAGCACGCATTGCCTCTCTTTCCGATAGCGTGTCt                        <  1:314359/71‑1 (MQ=255)
                   tatCGACGAACTGGAAGCCCTGCATCAGGAACTGGAAGCACGCATTGCCTCTCTTTCCGATAGCGTGTCt                        <  1:50942/70‑1 (MQ=255)
                        acgaACTGGAAGCCCTGCATCAGGAACTGGAAGCACGCATTGCCTCTCTTTCCGATAGCGGGTCt                        <  1:1307504/65‑1 (MQ=255)
                                   gCCCTGCATCAGGAACTGGAAGCACGCATTGCCTCTCTTTCCGATAGCGTGTCTAACGCCCGTGAAGAgcg       <  1:143709/71‑1 (MQ=255)
                                   gCCCTGCATCAGGAACTGGAAGCACGCATTGCCTCTCTTTCCGATAGCGTGTCTAACGCCCGTGAAGAgcg       <  1:649601/71‑1 (MQ=255)
                                   gCCCTGCATCAGGAACTGGAAGCACGCATTGCCTCTCTTTCCGATAGCGTGTCTAACGCCCGTGAAGAgcg       <  1:838034/71‑1 (MQ=255)
                                   gCCCTGCATCAGGAACTGGAAGCACGCATTGCCTCTCTTTCCGATAGCGTGTCTAACGCCCGTGAAGAgcg       <  1:517053/71‑1 (MQ=255)
                                   gCCCTGCATCAGGAACTGGAAGCACGCATTGCCTCTCTTTCCGATAGCGTGTCTAACGCCCGTGAAGAgcg       <  1:932633/71‑1 (MQ=255)
                                   gCCCTGCATCAGGAACTGGAAGCACGCATTGCCTCTCTTTCCGATAGCGTGTCTAACGCCCGTGAAGAgcg       <  1:1014965/71‑1 (MQ=255)
                                      cTGCATCAGGAACTGGAAGCACGCATTGCCTCTCTTTCCGATAGCGTGTCTAACGCCCGTGAAGAGCGCAt    >  1:838896/1‑71 (MQ=255)
                                         cATCAGGAGCTGGAAGCACGCATTGCCTCTCTTTCCGATAGCGTGTCTAACGCCCGTGAAGAGCGCATgg  <  1:1307752/70‑1 (MQ=255)
                                         cATCAGGAGCTGGAAGCACGCATTGCCTCTCTTTCCGATAGCGTGTCTAACGCCCGTGAAGAGCGCATgg  <  1:739023/70‑1 (MQ=255)
                                         cATCAGGAGCTGGAAGCACGCATTGCCTCTCTTTCCGATAGCGTGTCTAACGCCCGTGAAGAGCGCATgg  <  1:987139/70‑1 (MQ=255)
                                         cATCAGGAGCTGGAAGCACGCATTGCCTCTCTTTCCGATAGCGTGTCTAACGCCCGTGAAGAGCGCATgg  <  1:1532135/70‑1 (MQ=255)
                                          aTCAGGAGCTGGAAGCACGCATTGCCTCTCTTTCCGATAGCGTGTCTAACGCCCGTGAAGAGCGCATgt  <  1:426534/69‑2 (MQ=255)
                                             aGGAGCTGGAAGCACGCATTGCCTCTCTTTCCGATAGCGTGTCTAACGCCCGTGAAGAGCGCa     >  1:196792/1‑63 (MQ=255)
                                             aGGAGCTGGAAGCACGCATTGCCTCTCTTTCCGATAGCGTGTCTAACGCCCGTGAAGAGCGCa     >  1:1681708/1‑63 (MQ=255)
                                             aGGAGCTGGAAGCACGCATTGCCTCTCTTTCCGATAGCGTGTCTAACGCCCGTGAAGAGCGCa     >  1:144365/1‑63 (MQ=255)
                                             aGGAGCTGGAAGCACGCATTGCCTCTCTTTCCGATAGCGTGTCTAACGCCCGTGAAGAGCGCa     >  1:752886/1‑63 (MQ=255)
                                             aGGAACTGGAAGCACGCATTGCCTCTCTTTCCGATAGCgtg                           >  1:1402878/1‑41 (MQ=255)
                                             aGGAACTGGAAGCACGCATTGCCTCTCTTTCCGATAGCgtg                           >  1:541891/1‑41 (MQ=255)
                                             aGGAACTGGAAGCACGCATTGCCTCTCTTTCCGATAGCgtg                           >  1:457310/1‑41 (MQ=255)
                                             aGGAACTGGAAGCACGCATTGCCTCTCTTTCCGATAGCgtg                           >  1:227754/1‑41 (MQ=255)
                                             aGGAACTGGAAGCACGCATTGCCTCTCTTTCCGATAGCgtg                           >  1:1853221/1‑41 (MQ=255)
                                             aGGAACTGGAAGCACGCATTGCCTCTCTTTCCGATAGCgtg                           >  1:1546712/1‑41 (MQ=255)
                                             aGGAACTGGAAGCACGCATTGCCTCTCTTTCCGATAGCg                             >  1:801269/1‑39 (MQ=255)
                                             aGGAACTGGAAGCACGCATTGCCTCTCTTTCCGATAGAgtg                           >  1:999756/1‑41 (MQ=255)
                                                 |                                                             
GTCAGGGCAAGAATTTTGATATCGACGAACTGGAAGCCCTGCATCAGGAACTGGAAGCACGCATTGCCTCTCTTTCCGATAGCGTGTCTAACGCCCGTGAAGAGCGCATGG  >  W3110S.gb/978381‑978491

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: