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,7080GA80.7% 64.0 / 10.5 31A722T (GCC→ACC) acnAaconitate hydratase 1
Reads supporting (aligned to +/- strand):  ref base G (5/1);  new base A (9/16);  total (14/17)
Fisher's exact test for biased strand distribution p-value = 6.71e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AAAGACTTTAACTCCTACGGTTCGCGGCGTGGTAACCATGAAGTGATGATGCGCGGCACCTTCGCCAATATTCGCATCCGTAATGAAATGGTGCCTGGCGTTGAAGGGGGGATGACGCGGCATTTACCT  >  W3110S.gb/1339645‑1339773
                                                               |                                                                 
aaaGACTTTAACTCCTACGGTTCGCGGCGTGGTAACCATGAAGTGATGATGCGCGGCACCTTCACCAat                                                              <  1:1190418/69‑1 (MQ=255)
 aaGACTTTAACTCCTACGGTTCGCGGCGTGGTAACCATGAAGTGATGATGCGCGGCACCTTCGCCAata                                                             >  1:467718/1‑69 (MQ=255)
 aaGACTTTAACTCCTACGGTTCGCGGCGTGGTAACCATGAAGTGATGATGCGCGGCACCTTCGCCAata                                                             >  1:1216260/1‑69 (MQ=255)
    aCTTTAACTCCTACGGTTCGCGGCGTGGTAACCATGAAGTGATGATGCGCGGCACCTTCACCAATAtt                                                           >  1:3283124/1‑68 (MQ=255)
        tAACTCCTACGGTTCGCGGCGTGGTAACCATGAAGTGATGATGCGCGGCACCTTCACCAATATTCGCAt                                                      >  1:555543/1‑69 (MQ=255)
         aaCTCCTACGGTTCGCGGCGTGGTAACCATGAAGTGATGATGCGCGGCACCTTCACCAATATTCGCATc                                                     >  1:1605085/1‑69 (MQ=255)
            tCCTACGGTTCGCGGCGTGGTAACCATGAAGTGATGATGCGCGGCACCTTCACCAATATTCGCATCCGt                                                  <  1:234908/69‑1 (MQ=255)
            tCCTACGGTTCGCGGCGTGGTAACCATGAAGTGATGATGCGCGGCACCTTCACCAATATTCGCATCCGt                                                  <  1:2787329/69‑1 (MQ=255)
            tCCTACGGTTCGCGGCGTGGTAACCATGAAGTGATGATGCGCGGCACCTTCACCAATATTCGCATCCGt                                                  <  1:1156283/69‑1 (MQ=255)
                   gTTCGCGGCGTGGTAACCATTAAGTGATGATGCGCGGCACCTTCACCAATATTCGCATCCGTAATGa                                             <  1:2132885/67‑1 (MQ=255)
                   gTTCGCGGCGTGGTAACCATGAAGTGATGATGCGCGGCACCTTCACCAATATTCGCATCCGTAATGa                                             <  1:2094993/67‑1 (MQ=255)
                   gTTCGCGGCGTGGTAACCATGAAGTGATGATGCGCGGCACCTTCACCAATATTCGCATCCGTAATGa                                             <  1:2285279/67‑1 (MQ=255)
                    ttCGCGGCGTGGTAACCATGAAGTGATGATGCGCGGCACCTTCACCAATATTCGCATCCGTAATGAAAt                                          <  1:1955030/69‑1 (MQ=255)
                    ttCGCGGCGTGGTAACCATGAAGTGATGATGCGCGGCACCTTCACCAATATTCGCATCCGTAATGAAAt                                          >  1:1538307/1‑69 (MQ=255)
                    ttCGCGGCGTGGTAACCATGAAGTGATGATGCGCGGCACCTTCACCAATATTCGCATCCGTAATGAAAt                                          <  1:268527/69‑1 (MQ=255)
                    ttCGCGGCGTGGTAACCATGAAGTGATGATGCGCGGCACCTTCACCAATATTCGCATCCGTAATGAAAt                                          <  1:3655370/69‑1 (MQ=255)
                     tCGCGGCGTGGTAACCATGAAGTGATGATGCGCGGCACCTTCACCAATATTCGCATCCGt                                                  >  1:2026376/1‑60 (MQ=255)
                         ggcgTGGTAACCATGAAGTGATGATGCGCGGCACCTTCACCAATATTCGCATCCGTAATGAAATGGt                                       <  1:1267015/67‑1 (MQ=255)
                                  aCCATGAAGTGATGATGCGCGGCACCTTCACCAATATTCGCATCCGTAATGAAATGGt                                       <  1:1140800/58‑1 (MQ=255)
                                  aCCATGAAGTGATGATGCGCGGCACCTTCACCAATATTCGCATCCGTAATGAAATGGTGCCTGGCGt                              >  1:810653/1‑67 (MQ=255)
                                    cATGAAGTGATGATGCGCGGCACCTTCACCAATATTCGCATCCGTAATGAAATGGTGCCTGGCGTTGa                           <  1:1387190/68‑1 (MQ=255)
                                        aaGTGATGATGCGCGGCACCTTCACCAATATTCGCATCCGTAATGAAATGGTGCCTGGCGt                              <  1:2742093/61‑1 (MQ=255)
                                            gatgatGCGCGGCACCTTCGCCAATATTCGCATCCGTAATGAAATGGTGCCTGGCGTTGAAGGGGGGAt                  >  1:2950792/1‑69 (MQ=255)
                                            gatgatGCGCGGCACCTTCGCCAATATTCGCATCCGTAATGAAATGGTGCCTGGCGTTGAAGGGGGGAt                  >  1:3538093/1‑69 (MQ=255)
                                            gatgatGCGCGGCACCTTCGCCAATATTCGCATCCGTAATGAAATGGTGCCTGGCGTTGAAGGGGGGAt                  >  1:459211/1‑69 (MQ=255)
                                               gatgCGCGGCACCTTCACCAATATTCGCATCCGTAATGAAATGGTGCCTGGCGTTGAAGGGGGGATGAc               >  1:2866239/1‑69 (MQ=255)
                                                 tgCGCGGCACCTTCACCAATATTCGCATCCGTAATGAAAt                                          <  1:3231536/40‑1 (MQ=255)
                                                         aCCTTCGCCAATATTCGCATCCGTAATGAAATGGTGCCTGGCGTTGAAGGGGGGATGACGCGGCAttt      <  1:1130948/68‑1 (MQ=255)
                                                            ttCACCAATATTCGCATCCGTAATGAAATGGTGCCTGGCGTTGAAGGGGGGATGACGCGGCATTTACCt  >  1:841741/1‑69 (MQ=255)
                                                            ttCACCAATATTCGCATCCGTAATGAAATGGTGCCTGGCGTTGAAGGGGGGATGACGCGGCATTTACCt  <  1:24813/69‑1 (MQ=255)
                                                            ttCACCAATATTCGCATCCGTAATGAAATGGTGCCTGGCGTGGAAGGGGGGATGACGCGGCATTTAc    >  1:3620491/1‑67 (MQ=255)
                                                               |                                                                 
AAAGACTTTAACTCCTACGGTTCGCGGCGTGGTAACCATGAAGTGATGATGCGCGGCACCTTCGCCAATATTCGCATCCGTAATGAAATGGTGCCTGGCGTTGAAGGGGGGATGACGCGGCATTTACCT  >  W3110S.gb/1339645‑1339773

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: