Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24404,135,6070AC46.0% 7.4 / 16.9 51Q32P (CAG→CCG) tgnC(Z)‑2‑((N‑methylformamido)methylene)‑5‑hydroxybu tyrolactone dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (13/14);  new base C (18/5);  total (31/20)
Fisher's exact test for biased strand distribution p-value = 4.17e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.18e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CGCACCCTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCCGCGGCAAGCTGCTGCGTCGTCTGG  >  KT2440/4135478‑4135741
                                                                                                                                 |                                                                                                                                      
cGCACCCTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCACCCCGGCCACCGCCCCCGCCt                                                                                                                                   >  7:55318/1‑135 (MQ=255)
cGCACCCTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCACCCCGGCCACCGCCCAGGCCt                                                                                                                                   >  2:161959/1‑135 (MQ=255)
cGCACCCTCACCGAATCCGACAAGGTGATTGTTGCATGCAACTCGTTCGTTTCCATATGTGCATCGACGTCCAATTGCGCGATGAGAAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGACCACGCCt                                                                                                                                   >  5:163146/1‑135 (MQ=255)
     ccTCACCGATCCCGCAATGGGGATTGTTGACTGACCCTGTTTCGTTTCCAGATGTCCATCGACGGCCAATGGCGCTATGCCCAGCGGGCCAAGACCTTCGACTGCATCAACCCGGCCACCGCCCAGGCCTGGgcg                                                                                                                              <  8:124467/135‑1 (MQ=255)
           cGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCGGGGCGCACCTg                                                                                                                        >  3:235576/1‑135 (MQ=255)
           cGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCAAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCGGCCTGGGCGCAACTg                                                                                                                        >  3:243154/1‑135 (MQ=255)
            gAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGCCCTTCGACAGCCTCAACCCGGCCCCCGCCCAGGCCTGGGCGCGCCTGc                                                                                                                       >  1:288905/1‑135 (MQ=255)
            gAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAAAGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCGCCCAGGGCGCAACGGc                                                                                                                       >  3:300587/1‑135 (MQ=255)
                      aGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCCGCCAGGGCGCACCAGCCCGACCCCGa                                                                                                             >  8:5166/1‑135 (MQ=255)
                        ggtaTTGTTGCATGACCCTCGTTCGTTTCCAGTTGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATg                                                                                                           <  3:87871/132‑1 (MQ=255)
                        gggATTGTTGCATGACCCTCGTTCGTTTCCAGATGTCCAGCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCCCCCAGGCCTGGGCGCAACTGCCCGACGCCGATg                                                                                                           <  6:66649/135‑1 (MQ=255)
                          gATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCa                                                                                                                                        <  5:162978/104‑1 (MQ=255)
                           aTTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAg                                                                                                        <  6:78578/135‑1 (MQ=255)
                               ttgCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGCGCGGCAAGACCTTCGACAGCCTCAACCCGCCCACCGCCCCCGCCTCGGCGCCACGGCCCCCCGCCGCGGAAGcc                                                                                                      >  7:121646/1‑133 (MQ=255)
                                tgCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCGGCCAGGGCGGCCCTGCCCGCCGCCGCGGCAGCCgcg                                                                                                   >  3:153067/1‑133 (MQ=255)
                                ggcttGACCCTCGTTCTTTTCCAGATGTGCATCGACGGCAAATGGGGCGATGCCCAGAGCGCCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGAt                                                                                                   <  2:301136/131‑1 (MQ=255)
                                     gACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCACCCCGGCCACCGCCCCCCCCAGGGCGCAACAGCCCCCCGCCGCTGAAGCCGATGTGGa                                                                                              >  1:84185/1‑135 (MQ=255)
                                     gACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCACGCCTGGGCGCCCCAGCCCGACGCCGCGGCACCCGATGTggc                                                                                              >  6:198755/1‑134 (MQ=255)
                                      aCCCTCGTTCGTTACCAGATGTGCAGCCCCGCCCACTGGCGCGAGGCCCAGCGCGCCAACACCTGCGACACCCTCAACCCGCCCACCCCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAg                                                                                             <  1:24325/135‑1 (MQ=255)
                                       cccTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCGGCCGGGGCGGAACTGCCCGCCGCCGATGAAGCCCATGTGGAGc                                                                                            >  6:309625/1‑135 (MQ=255)
                                       cccTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCCGCCTGGGCGGCACGGCGCGACCCCGAGGCAGACGAGGTGGAGc                                                                                            >  1:232577/1‑135 (MQ=255)
                                       cccTCGTGCGTTCCCAGATGTGCATCGACGCCCAATGGCGCTATGCCAAGAGCGGAAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGc                                                                                            <  8:202713/135‑1 (MQ=255)
                                            gttcgttTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCGGGGCGCCACTGCCCCCCGCCCCTGAAGCCGCTGTGGAGCTGGcc                                                                                       >  6:333324/1‑135 (MQ=255)
                                            gttcgttTCCAGATGTGCATCGACGGCCAATGGCACGATGCCCAGAGCGGCAAGCCCTTCGACAGCCTCAACCCGGCCACCGCCCCGGCCGGGGCCGAACGGCCCGACCCCGCAGAAGACGAAGTGGAGCTGGcc                                                                                       >  6:147785/1‑135 (MQ=255)
                                             ttcgttTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCGCGGCCCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCg                                                                                      >  8:147626/1‑135 (MQ=255)
                                              tcgttTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCCACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGt                                                                                     >  6:255324/1‑135 (MQ=255)
                                                gttTCCTGATGGGCTGCGCCGGCCAAGGGCGCGATGCCCAGAGCGGAAAGACATTCGCCAGCCTCACCCCGGCCACCGCCCCGGCGGGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGtt                                                                                    <  1:253602/134‑1 (MQ=255)
                                                gttTCCAGATGTCCTGCGCCGCCCATTGGCGCGATGCCCAGAGCGCCAACACCTTCCACACCCCCAACCCGCCCCCCGCCCCGGCCGGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGtt                                                                                    <  2:5425/134‑1 (MQ=255)
                                                    ccAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCCCCGCCCCGGCCTGTGCCCCACTACCCCACGCCGCAGAAGCCGAGGGGGAGCAGGACCTTCCGgc                                                                               >  2:229898/1‑135 (MQ=255)
                                                    ccAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCGGGGCGCAACTGCCCGACGCCGCTGAAGCCGATGTGGAGCTGGCCGTTCAGgc                                                                               >  8:111318/1‑135 (MQ=255)
                                                     cAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCG‑CCAGGCCTGGGCGCAACTGCCCGACGCCGATGAATCCGATGTGGAGCTGGCCGTTCAGgccc                                                                             >  3:30874/1‑134 (MQ=255)
                                                      aGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCCGCCCGGGCGCCCCTGCCCGCCGCCGCGGAAGCCGCGGTGGCGCTGGCCGGTGAGgccc                                                                             >  5:170705/1‑134 (MQ=255)
                                                      aGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCCCCCCCGCCGGGGCGCCACAGCCCGCCGCCCCGGAAGCCCCAGGGGAGCGGGCCGGTCAGgccc                                                                             >  1:312000/1‑134 (MQ=255)
                                                            gCATCGACGGCCACGGGCGCGAGGCCCAGAGCGGCAAGACCTTCGACCGCCTCACCCCGGCCACCGCCCCGGCCGGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGc                                                                       <  1:242351/135‑1 (MQ=255)
                                                            gCAGCGACGGCCACGGGCGCGATGCCCAGCGCGCCAAGACCTGCCACACCCTCACCCCGCCCCCCCCCCCGGCCGGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGc                                                                       <  7:152084/135‑1 (MQ=255)
                                                              aTCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCGGCCGGGGCGGCCCGGCCCGCCGCCGAGGGAGCCGCTGTGGAGCTGGCCGGTCAGGCCGCCCCGCGt                                                                     >  1:4048/1‑135 (MQ=255)
                                                              aTCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCCGGCTGGGGGCCCCTGCCCGCGGCCGCTGAAGCCGCTGTGGCGCTGGGCGGTACGGCCGCCCcgcgg                                                                     >  6:6066/1‑134 (MQ=255)
                                                                cGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCGGCCTGGGCGCCACCGCCCGCCGCCGCTGAAGCCGATGTGGAGCGGGGCGTTCAGGCCGCCCCCCGTGc                                                                   >  3:317517/1‑135 (MQ=255)
                                                                cGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCGGGGCGCAACTGGCCGACGCCGATGAACACGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCCTGc                                                                   >  8:30538/1‑135 (MQ=255)
                                                                                   gCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCGGCCTGGGCGCAACTGCCCGCCGCCGCTGAAGCCGATGTGGATCTGGCCGTTCAGGCCGCCCCGCGTGCCGTCGCCAGCAAGGCGTgg                                                >  2:293729/1‑135 (MQ=255)
                                                                                                 gACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGCCCCCCATGAAGCCGCTGTGGAGCTGGCCGTTCAGGCCGGCCCCCCTGCCGTCGACAGCAAGGCGGGGCGCCGCCTCACCGc                                  >  1:111538/1‑135 (MQ=255)
                                                                                                     ttCGCCAGCCTCACCCCGGCCCCGCCCCAGCCCGGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTAcc                              <  6:8956/135‑1 (MQ=255)
                                                                                                     ttCGACAGCCTCAACCCGGCCCCCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTAcc                              <  4:218641/135‑1 (MQ=255)
                                                                                                     ttCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTAcc                              <  4:94532/135‑1 (MQ=255)
                                                                                                        caccgCCCCACCCCCGCCCCCGCCCTGGCCGGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTAccgc                            <  5:99776/130‑1 (MQ=255)
                                                                                                         aCAGCCTCAACCCGGCCACCCCCCCGGCCGGGGCGCCACTGCCCGCCGCCGCTGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTGGCCAGCAAGGCGTGGCGCCGGCTCACCGCTAccgccc                          >  8:133726/1‑135 (MQ=255)
                                                                                                         aCAGCCACAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTAccgccc                          <  2:82801/135‑1 (MQ=255)
                                                                                                            gCCTCACCCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCcgcg                       <  8:140815/135‑1 (MQ=255)
                                                                                                            gCCTCAACCCGCCCACCGCCCCGGCCGGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCcgcg                       <  4:320100/135‑1 (MQ=255)
                                                                                                                cAACCCGGCCACCGCCCACGCCAGGGCGCAACAGCCCCCCGCCGAATAAGCCGATGTGGAGCTGGCCGTGAAGGCCGCCCCCCGTGCCGTCGACAGCAAGGCGGGGCGCAGGAACCCCCCGACCCCCCGCCGCac                   >  2:84462/1‑134 (MQ=255)
                                                                                                                    ccGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCCGCGGCAAgctg               <  6:81820/135‑1 (MQ=255)
                                                                                                                          aCCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCCGCGGCAAGCTGCTGcgt         <  2:84185/135‑1 (MQ=255)
                                                                                                                              ccccGGCCGGGGCGCCACTGCCCGCCGCCGCTGAAGCCGATGTGGAGCTGGCCCGTCAGGCCGCCCCGCCTGCCGTCGCCCGCCACGCGTGGCGCCGCCTCACCGCGACCCCCCGCGGCACCCTGCTgcggcggc     >  2:74851/1‑133 (MQ=255)
                                                                                                                                 aGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCGTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCCGCGGCAAGCTGCTGCGTCGTCTgg  >  6:326657/1‑135 (MQ=255)
                                                                                                                                 |                                                                                                                                      
CGCACCCTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCCGCGGCAAGCTGCTGCGTCGTCTGG  >  KT2440/4135478‑4135741

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