Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,181,288 T→G 76.7% D64E (GAT→GAG cobB → deacetylase of acetyl‑CoA synthetase, NAD‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,181,2880TG76.7% 78.9 / 9.0 30D64E (GAT→GAGcobBdeacetylase of acetyl‑CoA synthetase, NAD‑dependent
Reads supporting (aligned to +/- strand):  ref base T (2/5);  new base G (14/9);  total (16/14)
Fisher's exact test for biased strand distribution p-value = 2.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

TACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCG  >  W3110S.gb/1181224‑1181346
                                                                |                                                          
tACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGcc                                                        <  1:538743/69‑1 (MQ=255)
  cTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCtg                                                      <  1:3172119/69‑1 (MQ=255)
          gACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCCCAGAGGGCCTGTGGgaag                                               >  1:2647713/1‑68 (MQ=255)
                 gCAGGAATTTCTGCGGAATCAGGTTTTCGTACCTTTCGCGCCGCAGAGGGCCtg                                                      <  1:1985158/54‑1 (MQ=255)
                 gCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCtg                                                      <  1:1085228/54‑1 (MQ=255)
                   aGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGAGGGCCTGTGGGAAGAACATCggg                                      >  1:2618187/1‑68 (MQ=255)
                      aaTTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGt                                     >  1:3333473/1‑66 (MQ=255)
                       aTTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGAGGGCCTGTGGgaag                                               >  1:1575591/1‑55 (MQ=255)
                       aTTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGAGGGCCTGTGGGAAGAACATCGGGt                                     >  1:2847561/1‑65 (MQ=255)
                       aTTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGAGGGCCTGTGGGAAGAACATCGGGt                                     >  1:1007099/1‑65 (MQ=255)
                       aTTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGAGGGCCTGTGGGAAGAACATCGGGt                                     >  1:2396940/1‑65 (MQ=255)
                       aTTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGAGGGCCTGTGGGAAGAACATCGGGt                                     >  1:211127/1‑65 (MQ=255)
                       aTTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGAGGGCCTGTGGGAAGAACATCGGGt                                     >  1:1189024/1‑65 (MQ=255)
                       aTTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGAGGGCCTGAGGGAAGAACATCGGGt                                     >  1:1638331/1‑65 (MQ=255)
                         ttCTGCGTAATCAGGTATTCGTACCTTTCGCGCCGCAGAGGGCCTGTGGGAAGAACATCGGGTTGAAGa                               <  1:1038752/69‑1 (MQ=255)
                         ttCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCtg                                                      <  1:302323/46‑1 (MQ=255)
                          tCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGAt                              >  1:3611994/1‑69 (MQ=255)
                          tCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGAGGGCCTGTGGgaag                                               <  1:2223795/52‑1 (MQ=255)
                          tCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGAGGGCCTGTGGGAAGAACATCGGGTTGAAGAt                              >  1:342301/1‑69 (MQ=255)
                                gAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGc                         <  1:3482964/68‑1 (MQ=255)
                                gAATCAGGTATTCGTACCTTTCGCGCCGCAGAGGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGc                         <  1:1306553/68‑1 (MQ=255)
                                     aGGTATTCGTACCTTTCGCGCCGCAGAGGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGc                         >  1:3474960/1‑63 (MQ=255)
                                         aTTCGTACCTTTCGCGCCGCAGAGGGCCTGTGGGAAGAACATCGGGTTGaa                                 <  1:1993493/51‑1 (MQ=255)
                                         aTTCGTACCTTTCGCGCCGCAGAGGGCCTGTGGGAAGAACATCGGGTTGaa                                 <  1:762813/51‑1 (MQ=255)
                                            cGTACCTTTCGCGCCGCAGAGGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGt            >  1:3370088/1‑69 (MQ=255)
                                            cGTACCTTTCGCGCCGCAGAGGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGt            >  1:3501482/1‑69 (MQ=255)
                                                ccTTTCGCGCCGCAGAGGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCg        <  1:3685827/69‑1 (MQ=255)
                                                    tCGCGCCGCAGAGGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATc     >  1:2847585/1‑68 (MQ=255)
                                                       cgccgcAGAGGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATcgcg  <  1:2938981/68‑1 (MQ=255)
                                                       cgccgcAGAGGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATcgcg  <  1:855557/68‑1 (MQ=255)
                                                                |                                                          
TACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCG  >  W3110S.gb/1181224‑1181346

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: