Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,031,910 A→G 29.2% N64D (AAC→GAC)  idi → isopentenyl diphosphate isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,031,9100AG29.2% 30.3 / 15.4 24N64D (AAC→GAC) idiisopentenyl diphosphate isomerase
Reads supporting (aligned to +/- strand):  ref base A (9/8);  new base G (3/4);  total (12/12)
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

ATGCCAAAGGACAATTATTAGTTACCCGCCGCGCACTGAGCAAAAAAGCATGGCCTGGCGTGTGGACTAACTCGGTTTGTGGGCACCCACAACTGGGAGAAAGCAACGAAGACGCAGTGATCCGCCGTTGCCGT  >  W3110S.gb/3031842‑3031975
                                                                    |                                                                 
aTGCCAAAGGACAATTATTAGTTACCCGCCGCGCACTGAGCAAAAAAGCATGGCCTGGCGTGTGGACTaa                                                                  <  1:2400584/70‑1 (MQ=255)
   ccAAAGGACAATTATTAGTTACCCGCCGCGCACTGAGCAAAAAAGCATGGCCTGGCGTGTGGACTAACTCg                                                              >  1:237206/1‑71 (MQ=255)
        ggACAATTATTAGTTACCCGCCGCGCACTGAGCAAAAAAGCATGGCCTGGCGTGTGGACTAACTCGGTTtg                                                         <  1:2212211/71‑1 (MQ=255)
         gACAATTATTAGTTACCCGCCGCGCACTGAGCAAAAAAGCATGGCCTGGCGTGTGGACTAACTCGGTTtgt                                                        >  1:1403042/1‑71 (MQ=255)
          aCAATTATTAGTTACCCGCCGCGCACTGAGCAAAAAAGCATGGCCTGGCGTGTGGACTGACTCGGTTtgtg                                                       <  1:1961411/71‑1 (MQ=255)
          aCAATTATTAGTTACCCGCCGCGCACTGAGCAAAAAAGCATGGCCTGGCGTGTGGACTAACTCGGTTtgtg                                                       <  1:2221959/71‑1 (MQ=255)
           cAATTATTAGTTACCCGCCGCGCACTGAGCAAAAAAGCATGGCCTGGCGTGTGGACTGACTCGGTTTGTgg                                                      >  1:294090/1‑71 (MQ=255)
               tattaGTTACCCGCCGCGCACTGAGCAAAAAAGCATGGCCTGGCGTGTGGACTGACTCGGTTTGTgg                                                      >  1:2834152/1‑67 (MQ=255)
               tattaGTTACCCGCCGCGCACTGAGCAAAAAAGCATGGCCTGGCGTGTGGACTGACTCGGTTTGTGGGCAc                                                  >  1:1918195/1‑71 (MQ=255)
                attaGTTACCCGCCGCGCACTGAGCAAAAAAGCATGGCCTGGCGTGTGGACTAACTCGGTTTGTGGGCAcc                                                 >  1:602100/1‑71 (MQ=255)
                attaGTTACCCGCCGCGCACTGAGCAAAAAAGCATGGCCTGGCGTGTGGACTAACTCGGTTTGTGGGCAcc                                                 >  1:331761/1‑71 (MQ=255)
                attaGTTACCCGCCGCGCACTGAGCAAAAAAGCATGGCCTGGCGTGTGGACTAACTCGGTTTGTGGGCAcc                                                 >  1:1601255/1‑71 (MQ=255)
                    gTTACCCGCCGCGCACTGAGCAAAAAAGCATGGCCTGGCGTGTGGACTAACTCGGTTTGTGGGCACCcac                                              <  1:2255412/70‑1 (MQ=255)
                       aCCCGCCGCGCACTGAGCAAAAAAGCATGGCCTGGCGTGTGGACTAACTCGGTTTGTGGGCACCCACAACt                                          >  1:3116351/1‑71 (MQ=255)
                                 cACTGAGCAAAAAAGCATGGCCTGGCGTGTGGACTGACTCGGTTTGTGGGCACCCACAACTGGGAGAAAGc                                <  1:1888042/71‑1 (MQ=255)
                                   cTGAGCAAAAAAGCATGGCCTGGCGTGTGGACTAACTCGGTTtg                                                         <  1:3218722/44‑1 (MQ=255)
                                   cTGAGCAAAAAAGCATGGCCTGGCGTGTGGACTAACTCGGTTtg                                                         <  1:295870/44‑1 (MQ=255)
                                      aGCAAAAAAGCATGGCCTGGCGTGTGGACTAACTCGGTTTGTGGGCACCCACAACTGGGAGAAAGCAACg                            <  1:364693/70‑1 (MQ=255)
                                           aaaaGCATGGCCTGGCGTGTGGACTGACTCGGTTTGTGGGCACCCACAACTGGGAGAAAGCAACGAAGAc                       <  1:1337067/70‑1 (MQ=255)
                                            aaaGCATGGCCTGGCGTGTGGACTAACTCGGTTTGTGGGCACCCACAACTGGGAGAAAGCAACGAAGAc                       >  1:2101447/1‑69 (MQ=255)
                                            aaaGCATGGCCTGGCGTGTGGACTAACTCGGTTTGTGGGCACCCACAACTGGGAGAAAGCAACGAAGACg                      >  1:1129978/1‑70 (MQ=255)
                                                cATGGCCTGGCGTGTGGACTAACTCGGTTTGTGGGCACCCACAACTGGGAGAAAGCAACGAAGACGCAGTg                 <  1:2061219/71‑1 (MQ=255)
                                                 aTGGCCTGGCGTGTGGACTGACTCGGTTTGTGGGCACCCACAACTGGGAGAAAGCAACGAAGACGCAGTGa                <  1:2091746/71‑1 (MQ=255)
                                                    gCCTGGCGTGTGGACTAACTCGGTTTGTGGGCACCCACAACTGGGAGAAAGCAACGAAGACGCAGTGATcc             <  1:3526708/71‑1 (MQ=255)
                                                        ggCGTGTGGACTAACTCGGTTTGTGGGCACCCACAACTGGGAGAAAGCAACGAAGACGCAGTGATccgccg         >  1:772279/1‑71 (MQ=255)
                                                                gactgactCGGTTTGTGGGCACCCACAACTGGGAGAAAGCAACGAAGACGCAGTGATCCgccgttgccgt  >  1:270088/1‑70 (MQ=255)
                                                                gactgactCGGTTTGTGGGCACCCACAACTGGGAGAAAGCAACGAAGACGCAGTGATCCgccgttgccgt  >  1:3261844/1‑70 (MQ=255)
                                                                gactgactCGGTTTGTGGGCACCCACAACTGGGAGAAAGCAACGAAGACGCAGTGATCCgccgttgccgt  >  1:1900874/1‑70 (MQ=255)
                                                                    |                                                                 
ATGCCAAAGGACAATTATTAGTTACCCGCCGCGCACTGAGCAAAAAAGCATGGCCTGGCGTGTGGACTAACTCGGTTTGTGGGCACCCACAACTGGGAGAAAGCAACGAAGACGCAGTGATCCGCCGTTGCCGT  >  W3110S.gb/3031842‑3031975

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: