Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,339,708 G→A 100% A722T (GCC→ACC)  acnA → aconitate hydratase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,339,7080GA97.1% 95.9 / ‑3.8 34A722T (GCC→ACC) acnAaconitate hydratase 1
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (14/19);  total (14/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CGAAAAGACTTTAACTCCTACGGTTCGCGGCGTGGTAACCATGAAGTGATGATGCGCGGCACCTTCGCCAATATTCGCATCCGTAATGAAATGGTGCCTGGCGTTGAAGGGGGGATGACGCGGCATTTACCTGA  >  W3110S.gb/1339642‑1339775
                                                                  |                                                                   
cGAAAAGACTTTAACTCCTACGGTTCGCGGCGTGGTAACCATGAAGTGATGATGCGCGGCACCTTCACCaa                                                                 <  1:2710489/71‑1 (MQ=255)
cGAAAAGACTTTAACTCCTACGGTTCGCGGCGTGGTAACCATGAAGTGATGATGCGCGGCACCTTCACCaa                                                                 <  1:1618044/71‑1 (MQ=255)
     aGACTTTAACTCCTACGGTTCGCGGCGTGGTAACCATGAAGTGATGATGCGCGGCACCTTCACCAATATTc                                                            >  1:3216805/1‑71 (MQ=255)
     aGACTTTAACTCCTACGGTTCGCGGCGTGGTAACCATGAAGTGATGATGCGCGGCACCTTCACCAATATTc                                                            >  1:3580480/1‑71 (MQ=255)
      gACTTTAACTCCTACGGTTCGCGGCGTGGTAACCATGAAGTGATGATGCGCGGCACCTTCACCAATATTCg                                                           <  1:1137597/71‑1 (MQ=255)
      gACTTTAACTCCTACGGTTCGCGGCGTGGTAACCATGAAGTGATGATGCGCGGCACCTTCACCAATATTCg                                                           <  1:2637956/71‑1 (MQ=255)
      gACTTTAACTCCTACGGTTCGCGGCGTGGTAACCATGAAGTGATGATGCGCGGCACCTTCACCAATATTCg                                                           <  1:2336718/71‑1 (MQ=255)
       aCTTTAACTCCTACGGTTCGCGGCGTGGTAACCATGAAGTGATGATGCGCGGCACCTTCACCAATATTCg                                                           >  1:1581356/1‑70 (MQ=255)
                    cGGTTCGCGGCGTGGTAACCATGAAGTGATGATGCGCGGCACCTTCACCAATATTCGCATCCGTAATGaaa                                             >  1:1855207/1‑71 (MQ=255)
                    cGGTTCGCGGCGTGGTAACCATGAAGTGATGATGCGCGGCACCTTCACCAATATTCGCATCCGTAATGaaa                                             >  1:1739409/1‑71 (MQ=255)
                        tCGCGGCGTGGTAACCATGAAGTGATGATGCGCGGCACCTTCACCAATATTCGCATCCGTAATGAAATGGt                                         >  1:2044441/1‑71 (MQ=255)
                         cgcgGCGTGGTAACCATGAAGTGATGATGCGCGGCACCTTCACCAATATTCGCATCCGTAATGAAATGGt                                         <  1:559070/70‑1 (MQ=255)
                         cgcgGCGTGGTAACCATGAAGTGATGATGCGCGGCACCTTCACCAATATTCGCATCCGTAATGAAATGGt                                         <  1:2238071/70‑1 (MQ=255)
                                  gTAACCATGAAGTGATGATGCGCGGCACCTTCGCCAATATTCGCATCCGTAATGAAATGGTGCCTGGCGt                                <  1:487627/70‑1 (MQ=255)
                                  gTAACCATGAAGTGATGATGCGCGGCACCTTCACCAATATTCGCATCCGTAATGAAATGGTGCCTGGCGtt                               >  1:2785389/1‑71 (MQ=255)
                                  gTAACCATGAAGTGATGATGCGCGGCACCTTCACCAATATTCGCATCCGTAATGAAATGGTGCCTGGCGtt                               >  1:1150736/1‑71 (MQ=255)
                                   tAACCATGAAGTGATGATGCGCGGCACCTTCACCAATATTCTCATCCGTAATGAAATGGTGCCTGGCg                                 <  1:2069329/68‑1 (MQ=255)
                                   tAACCATGAAGTGATGATGCGCGGCACCTTCACCAATATTCGCATCCGTAATGaa                                              <  1:3523039/55‑1 (MQ=255)
                                   tAACCATGAAGTGATGATGCGCGGCACCTTCACCAATATTCGCATCCGTAATGaa                                              <  1:3518682/55‑1 (MQ=255)
                                   tAACCATGAAGTGATGATGCGCGGCACCTTCACCAATATTCGCATCCGTAATGaa                                              <  1:919780/55‑1 (MQ=255)
                                   tAACCATGAAGTGATGATGCGCGGCACCTTCACCAATATTCGCATCCGTAATGAAATGGTGCCTGGCg                                 <  1:24826/68‑1 (MQ=255)
                                   tAACCATGAAGTGATGATGCGCGGCACCTTCACCAATATTCGCATCCGTAATGAAATGGTGCCTGGCGTTg                              <  1:2425121/71‑1 (MQ=255)
                                   tAACCATGAAGTGATGAAGCGCGGCACCTTCACCaa                                                                 <  1:407838/36‑1 (MQ=37)
                                    aaCCATGAAGTGATGATGCGCGGCACCTTCACCaa                                                                 <  1:3597736/35‑1 (MQ=255)
                                     aCCATGAAGTGATGATGCGCGGCACCTTCACCAATATTCGCATCCGTAATGAAATGGTGCCTGGCGTTg                              >  1:2568674/1‑69 (MQ=255)
                                     aCCATGAAGTGATGATGCGCGGCACCTTCACCAATATTCGCATCCGTAATGAAATGGTGCCTGGCGTTg                              >  1:2841060/1‑69 (MQ=255)
                                     aCCATGAAGTGATGATGCGCGGCACCTTCACCAATATTCGCATCCGTAATGAAATGGTGCCTGGCGTTGaa                            <  1:2024293/71‑1 (MQ=255)
                                     aCCATGAAGTGATGATGCGCGGCACCTTCACCAATATTCGCATCCGTAATGAAATGGTGCCTGGCGTTGaa                            <  1:1058855/71‑1 (MQ=255)
                                     aCCATGAAGTGATGATGCGCGGCACCTTCACCATATTCGCATCCGTAATGAAATGGTGCCTGGCGTTGaa                            <  1:430657/70‑1 (MQ=255)
                                               gatgatGCGCGGCACCTTCACCAATATTCGCATCCGTAATGAAATGGTGCCTGGCGTTGAAGGGGGGATGa                  >  1:2803422/1‑71 (MQ=255)
                                               gatgatGCGCGGCACCTTCACCAATATTCGCATCCGTAATGAAATGGTGCCTGGCGTTGAAGGGGGGATGa                  >  1:566859/1‑71 (MQ=255)
                                                         ggCACCTTCACCAATATTCGCATCCGTAATGAAATGGTGCCTGGCGTTGAAGGGGGGATGACGCGGCAttt        <  1:2500488/71‑1 (MQ=255)
                                                               ttCACCAATATTCGCATCCGTAATGAAATGGTGCCTGGCGTTGAAGGGGGGATGACGCGGCATTTACCTGa  >  1:2151860/1‑71 (MQ=255)
                                                               ttCACCAATATTCGCATCCGTAATGAAATGGTGCCTGGCGTTGAAGGGGGGATGACGCGGCATTTACCTGa  >  1:1926968/1‑71 (MQ=255)
                                                                  |                                                                   
CGAAAAGACTTTAACTCCTACGGTTCGCGGCGTGGTAACCATGAAGTGATGATGCGCGGCACCTTCGCCAATATTCGCATCCGTAATGAAATGGTGCCTGGCGTTGAAGGGGGGATGACGCGGCATTTACCTGA  >  W3110S.gb/1339642‑1339775

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