Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604792,137,0550TA30.0% 30.8 / 14.4 20L241Q (CTG→CAG) catB3muconate cyclo‑isomerase
Reads supporting (aligned to +/- strand):  ref base T (8/6);  new base A (3/3);  total (11/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.62e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CAACCAGGCCTGGAACGAAGCGGATGCCGCCACCGGCATTGCCATGCTGGAGACCGCCGGCATCGACCTGATCGAGCAGCCCACGCCACGCGAGCAACGCGCCGCGCTGGCGCGCCTGGCGGCGCGTTTCGTGGTGCCCATCATGGCCGACGAAGCGGTGTGCGGGCCCGAGGACGCGATGGAACTGGCCCGCATCGGCGGCGCCGACGTGTTCGCCCTGAAGATCGCCAAGTCGGGC  >  AM260479/2136948‑2137185
                                                                                                           |                                                                                                                                   
cAACCAGGCCTGGAACGAAGCGGATGCCGCCACCGGCATTGCCATGCTGGAGACCGCCGGCATCGACCTGATCGAGCAGCCCACGCCACGCGAGCAACGCGCCGCGCTGGCGCGCCTGGCGGCGCGt                                                                                                                 >  2:954736/1‑127 (MQ=255)
cAACCAGGCCTGGAACGAAGCGGATGCCGCCACCGGCATTGCCATGCTGGAGACCGCCGGCATCGACCTGATCGAGCAGCCCACGCCACGCGAGCAACGCGCCGCGCTGGCGCGCCTGGCGGCGCGt                                                                                                                 <  1:954736/127‑1 (MQ=255)
                     ggATGCCGCCACCGGCATTGCCATGCTGGAGACCGCCGGCATCGACCTGATCGAGCAGCCCACGCCACGCGAGCAACGCGCCGCGCTGGCGCGCCTGGCGGCGCGTTTCGTGGTGCCcat                                                                                                   <  2:945116/120‑1 (MQ=255)
                     ggATGCCGCCACCGGCATTGCCATGCTGGAGACCGCCGGCATCGACCTGATCGAGCAGCCCACGCCACGCGAGCAACGCGCCGCGCTGGCGCGCCTGGCGGCGCGTTTCGTGGTGCCcat                                                                                                   >  1:945116/1‑120 (MQ=255)
                     ggATGCCGCCACCGGCATTGCCATGCTGGAGACCGCCGGCATCGACCTGATCGAGCAGCCCACGCCACGCGAGCAACGCGCCGCGCTGGCGCGCCTGGCGGCGCGTTTCGTGGTGCCCATCATGGCcgacga                                                                                       <  2:825333/132‑1 (MQ=255)
                     ggATGCCGCCACCGGCATTGCCATGCTGGAGACCGCCGGCATCGACCTGATCGAGCAGCCCACGCCACGCGAGCAACGCGCCGCGCTGGCGCGCCTGGCGGCGCGTTTCGTGGTGCCCATCATGGCcgacga                                                                                       >  1:825333/1‑132 (MQ=255)
                     ggATGCCGCCACCGGCATTGCCATGCTGGAGACCGCCGGCATCGACCTGATCGAGCAGCCCACGCCACGCGAGCAACGCGCCGCGCTGGCGCGCCTGGCGGCGCGTTTCGTGGTGCCCATCATGGCCGACGAAg                                                                                     >  2:105366/1‑134 (MQ=255)
                     ggATGCCGCCACCGGCATTGCCATGCTGGAGACCGCCGGCATCGACCTGATCGAGCAGCCCACGCCACGCGAGCAACGCGCCGCCAGGCGCGCCAGCGCGGCGCGTTTCGTGGTGCCCATCATGGCcgac                                                                                         >  2:433124/1‑130 (MQ=37)
                     ggATGCCGCCACCGGCATTGCCATGCTGGAGACCGCCGGCATCGACCTGATCGAGCAGCCCACGCCACGCGAGCAACGCGCCGCCAGGCGCGCCAGCGCGGCGCGTTTCGTGGTGCCCATCATGGCcgac                                                                                         <  1:433124/130‑1 (MQ=37)
                                     aTTGCCATGCTGGAGACCGCCGGCATCGACCTGATCGAGCAGCCCACGCCACGCGAGCAACGCGCCGCCAGGCGCGccagcg                                                                                                                         <  1:550714/82‑5 (MQ=25)
                                     aTTGCCATGCTGGAGACCGCCGGCATCGACCTGATCGAGCAGCCCACGCCACGCGAGCAACGCGCCGCCAGGCGCGccagcg                                                                                                                         >  2:550714/1‑78 (MQ=25)
                                                       gccgGCATCGACCTGATCGAGCAGCCCACGCCACGCGAGCAACGCGCCGCGCTGGCGCGCCTGGCGGCGCGTTTCGTGGTGCCCATCATGGCcg                                                                                           >  1:599974/1‑94 (MQ=255)
                                                       gccgGCATCGACCTGATCGAGCAGCCCACGCCACGCGAGCAACGCGCCGCGCTGGCGCGCCTGGCGGCGCGTTTCGTGGTGCCCATCATGGCcg                                                                                           <  2:599974/94‑1 (MQ=255)
                                                            cATCGACCTGATCGAGCAGCCCACGCCACGCGAGCAACGCGCCGCCAGGCGCGCCAGCGCGGCGCGt                                                                                                                 <  2:250605/67‑1 (MQ=18)
                                                            cATCGACCTGATCGAGCAGCCCACGCCACGCGAGCAACGCGCCGCCAGGCGCGCCAGCGCGGCGCGt                                                                                                                 >  1:250605/1‑67 (MQ=18)
                                                                      aTCGAGCAGCCCACGCCACGCGAGCAACGCGCCGCGCTGGCGCGCCTGGCGGCGCGTTTCGTGGTGCCCATCATGGCCGACGAAGCGGTGTGCGGGCCCGAGGACGCGATGGAACTGGCCCGCATCGgcggcg                                     <  1:105366/133‑1 (MQ=255)
                                                                                              aaaCGCGCCGCGCTGGCGCGCCTGGCGGCGCGTTTCGTGGTGCCCATCATGGCCGACGAAGCGGTGTGCTGGCCCGAGGACGCGATGGAACt                                                      <  1:183698/91‑1 (MQ=255)
                                                                                              aaaCGCGCCGCGCTGGCGCGCCTGGCGGCGCGTTTCGTGGTGCCCATCATGGCCGACGAAGCGGTGTGCGGGCCCGAGGACGCGATGGAACt                                                      >  2:183698/2‑92 (MQ=255)
                                                                                                   cgccgcgcTGGCGCGCCTGGCGGCGCGTTTCGTGGTGCCCATCATGGCCGACGAAGCGGTGTGCGGGCCCGAGGACGCGATGGAACTGGCCCGCATCGGCGGCGCCGACGTGTTCGCCCTGAAGATCGCCAAGt       >  1:708192/1‑134 (MQ=255)
                                                                                                        cgcTGGCGCGCCTGGCGGCGCGTTTCGTGGTGCCCATCATGGCCGACGAAGCGGTGTGCGGGCCCGAGGACGCGATGGAACTGGCCCGCATCGGCGGCGCCGACGTGTTCGCCCTGAAGATCGCCAAGTCGggc  >  2:458267/1‑134 (MQ=255)
                                                                                                           |                                                                                                                                   
CAACCAGGCCTGGAACGAAGCGGATGCCGCCACCGGCATTGCCATGCTGGAGACCGCCGGCATCGACCTGATCGAGCAGCCCACGCCACGCGAGCAACGCGCCGCGCTGGCGCGCCTGGCGGCGCGTTTCGTGGTGCCCATCATGGCCGACGAAGCGGTGTGCGGGCCCGAGGACGCGATGGAACTGGCCCGCATCGGCGGCGCCGACGTGTTCGCCCTGAAGATCGCCAAGTCGGGC  >  AM260479/2136948‑2137185

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