Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,181,288 T→G 100% 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,2880TG78.4% 97.9 / 12.9 37D64E (GAT→GAGcobBdeacetylase of acetyl‑CoA synthetase, NAD‑dependent
Reads supporting (aligned to +/- strand):  ref base T (7/1);  new base G (21/8);  total (28/9)
Fisher's exact test for biased strand distribution p-value = 6.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.74e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTG  >  W3110S.gb/1181226‑1181352
                                                              |                                                                
cTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCtgtg                                                          >  1:1387818/1‑71 (MQ=255)
cTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCtgtg                                                          >  1:3454985/1‑71 (MQ=255)
cTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCtgt                                                           >  1:3458420/1‑70 (MQ=255)
cTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGAGGGCCtgtg                                                          >  1:411800/1‑71 (MQ=255)
cTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGAGGGCCtgtg                                                          >  1:2782875/1‑71 (MQ=255)
   gTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGAGGGCCTGTggg                                                        >  1:2001623/1‑70 (MQ=255)
   gTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGAGGGCCTGTggg                                                        >  1:2421410/1‑70 (MQ=255)
                       ttCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGAt                                    <  1:1920874/70‑1 (MQ=255)
                                       aTTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAgg                   >  1:125046/1‑71 (MQ=255)
                                       aTTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAgg                   >  1:39556/1‑71 (MQ=255)
                                       aTTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAgg                   >  1:3204945/1‑71 (MQ=255)
                                       aTTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAgg                   >  1:2621154/1‑71 (MQ=255)
                                        ttCGTACCTTTCGCGCCGCAGAGGGCCTGTGGGAAGAACATCGGGTTGaa                                       <  1:691508/50‑1 (MQ=255)
                                        ttCGTACCTTTCGCGCCGCAGAGGGCCTGTGGGAAGAACATCGGGTTGaa                                       <  1:218869/50‑1 (MQ=255)
                                          cGTACCTTTCGCGCCGCAGAGGGCCTGTGGGAAGAACATc                                               >  1:2992448/1‑40 (MQ=255)
                                          cGTACCTTTCGCGCCGCAGAGGGCCTGTGGGAAGAACATCGGGTTGAAGa                                     >  1:2800371/1‑50 (MQ=255)
                                          cGTACCTTTCGCGCCGCAGAGGGCCTGTGGGAAGAACATCGGGTTGAAGATGTgg                                >  1:2945513/1‑55 (MQ=255)
                                          cGTACCTTTCGCGCCGCAGAGGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGttt                >  1:2430854/1‑71 (MQ=255)
                                          cGTACCTTTCGCGCCGCAGAGGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGttt                >  1:2448203/1‑71 (MQ=255)
                                          cGTACCTTTCGCGCCGCAGAGGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGttt                >  1:3209332/1‑71 (MQ=255)
                                          cGTACCTTTCGCGCCGCAGAGGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGtt                 >  1:220122/1‑70 (MQ=255)
                                                tttCGCGCCGCAGAGGGCCTGTGGGAAGAACATCGGGTTGAAGa                                     >  1:2355138/1‑44 (MQ=255)
                                                tttCGCGCCGCAGAGGGCCTGTGGGAAGAACATCGGGTTGAAGa                                     >  1:2556735/1‑44 (MQ=255)
                                                tttCGCGCCGCAGAGGGCCTGTGGGAAGAACATCGGGTTGAAGa                                     >  1:2197432/1‑44 (MQ=255)
                                                 ttCGCGCCGCAGAGGGCCTGTGGGAAGAACATCGGGTTGAAGAtg                                   >  1:2387452/1‑45 (MQ=255)
                                                 ttCGCGCCGCAGAGGGCCTGTGGGAAGAACATCGGGTTGAAGAtg                                   >  1:1485646/1‑45 (MQ=255)
                                                 ttCGCGCCGCAGAGGGCCTGTGGGAAGAACATCGGGTTGAAGAtg                                   >  1:1674290/1‑45 (MQ=255)
                                                 ttCGCGCCGCAGAGGGCCTGTGGGAAGAACATCGGGTTGAAGAtg                                   >  1:2900859/1‑45 (MQ=255)
                                                 ttCGCGCCGCAGAGGGCCTGTGGGAAGAACATCGGGTTGAAGATGTgg                                >  1:622608/1‑48 (MQ=255)
                                                 ttCGCGCCGCAGAGGGCCTGTGGGAAGAACATCGGGTTGAAGATGTgg                                >  1:1217654/1‑48 (MQ=255)
                                                  tCGCGCCGCAGAGGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATcgc         >  1:342199/1‑70 (MQ=255)
                                                   cgcgCCGCAGAGGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATcgcg        <  1:1971484/70‑1 (MQ=255)
                                                   cgcgCCGCAGAGGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATcgcg        <  1:1403581/70‑1 (MQ=255)
                                                    gcgcCGCAGAGGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGa       <  1:309512/70‑1 (MQ=255)
                                                    gcgcCGCAGAGGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGa       <  1:1834449/70‑1 (MQ=255)
                                                       ccgcAGAGGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGAAGGTTTCGATCGCGATCCTg  <  1:2866506/71‑1 (MQ=255)
                                                        cgcAGAGGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTg  <  1:3437475/71‑1 (MQ=255)
                                                              |                                                                
CTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTG  >  W3110S.gb/1181226‑1181352

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