Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1222,950,1660CG56.2% 4.7 / 16.6 16intergenic (+64/‑132)spuC_2/puuC_1Putrescine‑‑pyruvate aminotransferase/NADP/NAD‑dependent aldehyde dehydrogenase PuuC
Reads supporting (aligned to +/- strand):  ref base C (2/5);  new base G (0/9);  total (2/14)
Fisher's exact test for biased strand distribution p-value = 1.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.80e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGTGCGGCGTGCTTGCGGACTGTGTGGTACCTGTTACTTGCGCAAGAAGGACCGCCACGACGGTCGTTTGGA  >  ACB122/2950062‑2950291
                                                                                                        |                                                                                                                              
tcgCCGTAGCCCGCACGGCGGGGCTGGGCGGCAAGCTGTACGCCGCTGACAAGGGCAACGCGCTTCCCTGGGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTcgcg                                                                                                  <  6:94654/133‑1 (MQ=37)
tAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTcgcg                                                                                                  <  3:660/135‑1 (MQ=255)
  gCCGTAGACCGCAGGGGGGGGGGGGTCGGAAAGCTTTAAGCGGGTGACTTGGGAAACGCGCTTCCCTGTGGGTGGGGGCACGCCGGCGAACACGGGCGAAGCGGGTGCAATCCACCGTGTGGCCTGATTCGCgg                                                                                                 <  4:43438/134‑1 (MQ=255)
                             ggCACGCTGAACGCCGCGCACAAGGCCACCGCGCTGCCCTGGGGGAGGGGGCACCCCCGCCAACACCGGCGACGCGGGTGCCACCCCCCGTGCCGCCTGATGCGCGGGCAGGCCCGCCCCCACAGAGGGAGTgcg                                                                    <  4:112818/135‑1 (MQ=255)
                             gccccgCTGTCCGCCGCGCACACGGGCACCGCGCTCCCCGGGGGGCGCGGGCACCCCCGCGAACCCCGGCCACGCGGGTGCCACCCCCCGGGGCGCCGGATTCGCGGGCATGCCCGCCCCCACAGAGGGAGTgcg                                                                    <  2:44997/130‑1 (MQ=12)
                                       acgCCGCTGACTCGGGACACGGGCTGCGGTGTGGGTGCGGGCAGGCCCGCGTACACCGGCCACGCGGGTGCCTTCCACGTTGGCGCCTTTTTGTGTGTATTGCCCGTTCCCACAGAGGGAGTGCGGCGTGCtttg                                                          <  8:176602/133‑3 (MQ=255)
                                        aGCCGCTGACAAGGGCACCGCGCTTCCCTGTGGGTGCGGGCACGCCCGCGAACACCGGCTAAGCGGGTGCCATCCACGGTGTCGCCTGATTCGCGGGCATGCCCGCCCCCACAGAGGGAGTGCGGCGTGCTTGc                                                          <  8:113664/134‑1 (MQ=37)
                                                 cAAGGGCAACGCGGTGCCGTGGGGGAGGGGGCACGCCCGGTAACACCGGCGAAGCGGGTGCCATCCACCGGGTCGCCTGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGTGCGGCGTGCTTGCGGACtgtg                                                  <  6:329607/133‑1 (MQ=25)
                                                         aCGCGCGCCCGTGGGGGCGCGGGCAGGCCCGCTAACACGGGCCACGCGGGTGCCATCCACCGGGCCGCCTGATGCGCGGGCAGGCCCGCTCCCACAGAGGGAGTGCGGCGTGCTTGCGGACTGTGTGGTACCTGt                                        <  4:359804/135‑1 (MQ=255)
                                                              cTTCCCGGGGGGCGGGGGCACGCCCGCGAACACCGGCGACGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGTGCGGCGTGCTTGCGGACTGTGTGGTACCTGTTACtt                                   <  1:27701/135‑1 (MQ=37)
                                                                   cTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGcc                                                                                                           <  7:247391/59‑1 (MQ=255)
                                                                   cTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGcc                                                                                                           >  8:247391/1‑59 (MQ=255)
                                                                                               cGGCTAAGCGGGTGCCATCCACCGGGCCGCCTGATGCGCGGGCTGCCCCGCCCCCACAGAGGGAGTGCGGCGTGCTTGCGGACTGTGTGGTACCTGTTACTTGCGCAAGAAGGACCGCCACGACGGTCGTTTGGa  <  3:182478/135‑1 (MQ=255)
                                                                                               cGGCGACGCGGGTGCCACCCACGGTGTGGCCTTATTGGCGGCCTTGCCCGCGCCCACAGAGGGAGTGCGGCGTGCTTGCGGACTGTGTGGTACCTGTTACTTGCGCAAGAAGGACCGCCACGACGGTCGTTTGGa  <  3:167601/135‑1 (MQ=255)
                                                                                                  cGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGTGCGGCGTGCTTGCGGACTGTGTGGTACCTGTTACTTGCGCAAGAAGGACCGCCACGACGGTCGt       <  7:85270/127‑1 (MQ=255)
                                                                                                  cGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGTGCGGCGTGCTTGCGGACTGTGTGGTACCTGTTACTTGCGCAAGAAGGACCGCCACGACGGTCGt       >  8:85270/1‑127 (MQ=255)
                                                                                                        |                                                                                                                              
TAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGTGCGGCGTGCTTGCGGACTGTGTGGTACCTGTTACTTGCGCAAGAAGGACCGCCACGACGGTCGTTTGGA  >  ACB122/2950062‑2950291

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