Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,197,638 G→A 67.6% L313L (CTG→CTA icd → isocitrate dehydrogenase, specific for NADP+

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,197,6380GA67.6% 26.1 / 19.5 26L313L (CTG→CTAicdisocitrate dehydrogenase, specific for NADP+
Reads supporting (aligned to +/- strand):  ref base G (2/6);  new base A (8/9);  total (11/15)
Fisher's exact test for biased strand distribution p-value = 4.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.60e-01

CGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAAT  >  W3110S.gb/1197577‑1197693
                                                             |                                                       
cGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTAGCAGCGCa                                                 <  1:88543/70‑1 (MQ=255)
          tgatgTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTAGCAGCGCAGGTTGGcggta                                      >  1:665354/1‑71 (MQ=255)
            atgTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTAGCAGCGCAGGTTGGcggtatc                                    >  1:1676332/1‑71 (MQ=255)
               ttATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTAGCAGCGCAGGTTGGcggtatcggt                                 <  1:2130338/71‑1 (MQ=255)
               ttATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTAGCAGCGCAGGTTGGcggtatcgg                                  <  1:1925538/70‑1 (MQ=255)
                    gCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTAGCAGCGCAGGTTGGCGGTATCGGTATCGc                            <  1:454139/71‑1 (MQ=255)
                    gCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTAGCAGCGCAGGTTGGCGGTATCGGTATCGc                            <  1:1774324/71‑1 (MQ=255)
                       tGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTTGCAg                                                     >  1:173478/1‑43 (MQ=255)
                         tATGAACCTGAACGGTGACTACATTTCTGCCGCCCTAGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCt                        >  1:1233953/1‑70 (MQ=255)
                           tGAACCTGAACGGTGACTACATTTCTGACGCCCTAGCAGCGCa                                                 >  1:2252544/1‑43 (MQ=255)
                           tGAACCTGAACGGTGACTACATTTCTGACGCCCTAGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGt                     >  1:710892/1‑71 (MQ=255)
                             aaCCTGAACGGTGACTACATTTCTGACGCCCTAGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTg                    >  1:948058/1‑70 (MQ=255)
                              aCCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGc                            <  1:2606949/61‑1 (MQ=255)
                              aCCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGCGGTATCGGTATCGc                            <  1:2399764/60‑1 (MQ=255)
                                      ggTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGcggtat                                     >  1:118670/1‑44 (MQ=255)
                                         gACTACATTTCTGACGCCCTAGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCgg          >  1:1095185/1‑68 (MQ=255)
                                           cTACATTTCTGACGCCCTAGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACg     <  1:1544962/71‑1 (MQ=255)
                                           cTACATTTCTGACGCCCTAGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACg     <  1:1204785/71‑1 (MQ=255)
                                              cATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACg     <  1:1957039/68‑1 (MQ=255)
                                              cATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAAt  >  1:2221478/1‑71 (MQ=255)
                                              cATTTCTGACGCCCTAGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGaa   <  1:37179/70‑1 (MQ=255)
                                              cATTTCTGACGCCCTAGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGa    <  1:926149/69‑1 (MQ=255)
                                               aTTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACg     <  1:1775364/67‑1 (MQ=255)
                                                 ttCTGACGCCCTAGCAGCGCAGGTTGGcggtatcggtat                               >  1:1535459/1‑39 (MQ=255)
                                                      aCGCCCTGGCAGCGCAGGTTGGCGGTATCGTTATCGCCCCTGGTGCaaa                <  1:1395673/49‑1 (MQ=255)
                                                      aCGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCaaa                <  1:1216675/49‑1 (MQ=255)
                                                             |                                                       
CGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAAT  >  W3110S.gb/1197577‑1197693

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