Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24402,950,2010TG59.1% 4.9 / 15.9 23intergenic (+99/‑97)spuC_2/puuC_1Putrescine‑‑pyruvate aminotransferase/NADP/NAD‑dependent aldehyde dehydrogenase PuuC
Reads supporting (aligned to +/- strand):  ref base T (2/7);  new base G (0/13);  total (3/20)
Fisher's exact test for biased strand distribution p-value = 1.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.85e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGTGCGGCGTGCTTGCGGACTGTGTGGTACCTGTTACTTGCGCAAGAAGGACCGCCACGACGGTCGTTTGGAGTACAAATGTACACTCTCGAATTCTGGCAACAACGCGCC  >  KT2440/2950091‑2950330
                                                                                                              |                                                                                                                                 
ggCACGCGGTACCCCGCGTACAAGGGCCCCGCGCTCCCCTGGGGGCGGGGGCACGCCCGCGACCCCCGGCGAACCCGGTGCCATCCCCCGGGGCGCCTGATCCGCGGGCAGGCCCGCGCCCACAGAGGGAGTgcg                                                                                                           <  8:170540/135‑1 (MQ=255)
ggCAAGCGGAACGCCGCCCACAAGGCCACCGCGCTCCCCGGGGGGCGGGGCCCCGCCCGCCACCACCGGCGACGCCGGGGCCAGCCCCCGGGCCCCCTTATGCGCGGCCCGCCCCGCGCCCACAGAGGGAGTgcg                                                                                                           <  4:320209/135‑1 (MQ=255)
          atccggCTTACATGGGCATGGCGATTCCCTGTGGGAGGGGGCAGGGCGGGTATCACGGGCGAAGCCGGTGCCAGACCCGGTGTCGCGTGTTTCGCGGGCATGCCGTCTCCCACAGAGGGAGTGCGGCGTGCTTGc                                                                                                 <  8:250728/130‑1 (MQ=255)
          acgCCGCTGACAAGGGCAACGCGCTGCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCGCCCACAGAGGGAGTGCGGCGTGCTTGc                                                                                                 <  4:284447/133‑1 (MQ=255)
          accccGCGGACACGGGCAACGCGCTGCCCGGGGGGCGCGGGCCCGCCCGCCAACACCGGCGAAGCGGGGGCCATCCACCGTGGCGCCTGATGCGCGGGCAGGCCCGCCCCCACAGAGGGAGTGCGGCGTGCTTGc                                                                                                 <  2:29948/132‑1 (MQ=255)
            gCCGGTGACCCGGGCACCGCGCTGCCCGGGGGGCGCGGGCACGCCCGCGACCACCGGCGACGCGGGTGCCATCCCCCGGGGCGCCGGATTCGCGGGCATGCCGGCTCCCACAGAGGGAGTGCGGCGTGCTTGCgg                                                                                               <  8:33505/135‑1 (MQ=255)
                 tGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGTGCGGCGTGCTTGCGGACtgt                                                                                          <  2:44751/135‑1 (MQ=255)
                 ggaaacgggCACCGCGCTTCCGTGGGGCGTGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGTGCGGCGTGCTTGCGGACtgtg                                                                                         <  6:231235/129‑1 (MQ=25)
                    cAAGGGCACCGCGCTGCCCTGGGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGGCGCCTGATTCGCGGGCAGGCCCGCCCCCACAGAGGGAGTGCGGCGTGCTTGCGGACtgtg                                                                                         <  7:238829/133‑1 (MQ=25)
                            aCGCGCTGCCGTGGGGGCGCGGGCAGGCCCGCGAACACGGGCGACGCGGGGGCCAGCCACCGGGGCCCCTGATGCGCGGGCTGGCCCGCTCCCACAGAGGGAGTGCGGCGTGCTTGCGGACTGTGTGGTACCTGt                                                                               <  2:66899/135‑1 (MQ=255)
                                                   cACGCCCGCGACCACGGGCGCCGCCGGTGCCCCCCCCCGTGGCGCCTGTGGCGCGGGCAGGCCCGCGCCCACAGAGGGAGTGCGGCGTGCTTGCGGACGGTGTGGTACCTGTTACTTGCGCAAGACGGCCCGCCa                                                        <  2:270347/135‑1 (MQ=255)
                                                   cACGCCCGCCAACACCGGCGACGCGGGTGCCAGCCACCGTGGCGCCTGATTCGCGGCCAGGCCCGCTCCCACAGAGGGAGTGCGGCGTGCTTGCGGACTGTGTGGTACCTGTTACTTGCGCAAGAAGGACCGCCa                                                        <  8:71941/135‑1 (MQ=255)
                                                     cGCCCGCTAACCCCGGGGACGCGGGTGCCAGCCCCCGGGGCCCCTGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGTGCGGCGTGCTTGCGGACTGTGTGGTACCTGTTACTTGCGCAAGAAGGACCGCCacg                                                      <  6:299677/135‑1 (MQ=255)
                                                                cccGGCCACGGCGGGGCCAGCCCCGGGGGCGCCTGCTGCGGGGCCAGCCCCGCCCCCACAGAGGGAGTGCGGCGTGCTTGCGGACTGTGTGGTACCTGTTACTTGCGCAAGAAGGACCGCCACGACGGTCGTTTg                                           <  5:104385/134‑1 (MQ=255)
                                                                acCGGCTAAGCCGGTGCCAGCCACCGTGTCGCCTGATTCGCGGGCAGGCCCGCTCCCACAGAGGGAGTGCGGCGTGCTTGCGGACTGTGTGGTACCTGTTACTTGCGCAAGAAGGACCGCCACGACGGTCGTTTg                                           <  6:367005/135‑1 (MQ=255)
                                                                  cGGCGAAGCGGGTGCCATCCCCCGTGTCGCCTGATGCGCGGGCATGCCCGCTCCCACAGAGGGAGTGCGGCGTGCTTGCGGACTGTGTGGTACCTGTTACTTGCGCAAGAAGGACCGCCACGACGGTCGTTTGGa                                         <  7:266263/135‑1 (MQ=255)
                                                                     cGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGTGCGGCGTGCGTGCGGACTGTGTGGTACCTGTTACTTGCGCAAGAAGGACCGCCCCGCCGGTCGTTTGGAGt                                       >  8:37934/1‑134 (MQ=255)
                                                                          cgggTGCCAGCCCCCGGGTCGCCTGATTCGCGGGCAGGCCCGCTCCCACAGAGGGAGTGCGGCGTGCTTGCGGACTGTGTGGTACCTGTTACTTGCGCAAGAAGGACCGCCACGACGGTCGTTTGGAGTACAAAt                                 <  8:343709/133‑1 (MQ=39)
                                                                          ccGTTGCCATCCCCCGGCGCGCCTTTTTGGAGGCTAGGCCGGCGCCCACAGAGGGAGTGCGGCGTGCTTGCGGAGTGTGTGGTACCTGTTACTTGCGCAAGAAGGACCGCCACGACGGTCGTTTGGAGTACAAAt                                 <  7:361047/135‑1 (MQ=255)
                                                                                           gtCGCCTGCTTCGGGTGCAAGCCCGCTCCCTCAGAGGGAGTGCGGCCTTCTGGGGGCCTGTGTGGTTCCAGTTAGTGGCGCAAGAAGGACCGCCAAGGCGGGCCGTTTGCGAAGAAAAGTACCCTCACCAATTCt                >  6:287526/1‑135 (MQ=255)
                                                                                             cGCCTGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGTGCGGCGTGCTTGCGGCCTGTGTGGTACCTGTTACGTGCGCAAGAAGGACCGCCACGACGGTCGTTTTGTGTACAAATGTACCCTCTCGGATTCTgg              >  8:235890/1‑135 (MQ=255)
                                                                                                    tccgcgGGCAGGCCCGCCCCCACAGAGGGAGTGCGGCGTGCTTGCGGACTGTGTGGTACCTGTTACTTGCGCAAGAAGGACCGCCACGACGGTCGTTTGGAGTACAAATGTACACTCTCGAATTCTGGcaaca         <  6:166876/131‑1 (MQ=255)
                                                                                                         gggCAGGCCCGCACCCACAGAGGGAGTGCGGCGTGCTTGCGGACTGTGTGGTACCTGTTACTTGCGCAAGAAGGACCGCCACGACGGTCGTTTGGAGTACAAATGTACACTCTCGAATTCTGGCAACAACGCGcc  <  5:174002/135‑1 (MQ=255)
                                                                                                              |                                                                                                                                 
GGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGTGCGGCGTGCTTGCGGACTGTGTGGTACCTGTTACTTGCGCAAGAAGGACCGCCACGACGGTCGTTTGGAGTACAAATGTACACTCTCGAATTCTGGCAACAACGCGCC  >  KT2440/2950091‑2950330

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: