Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24404,446,4610AC44.4% 28.1 / 12.6 45E807A (GAG→GCG) nicBNicotinate dehydrogenase subunit B
Reads supporting (aligned to +/- strand):  ref base A (18/7);  new base C (1/19);  total (19/26)
Fisher's exact test for biased strand distribution p-value = 4.16e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TCGAAGTGGTGGTTCGGCTCGCTGGCCGGGGTGTTCGGCGCCGCATTGGGCATGCTTGCCACCGCCTTGCCCTGGCGCGCCGAAATCGCCCCGGTAACGCCACCTGGTGTCGGCAGCTGGAGCGCCGCCATGCTCGAGCGCGGGCGCCAGGTGGCGGCGGCCGGCGACTGCGCCGTGTGCCATACGGTCAGCGGCGGCAAGGCCAATGCCGGGGGGCTGGCAATGGATACGCCTTTCGGCACGTTGTACAGCACCAATATCACC  >  KT2440/4446325‑4446588
                                                                                                                                        |                                                                                                                               
tCGAAGTGGTGGTTCGGCTCGCTGGCCGGGGTGTTCGGCGCCGCATTGGGCATGCTTGCCACGGCCTTGCCCTGGCGCGCCGAAATCGCCCCGGTAACGCCACCTGGTGTCGGCAGCTGGAGCGCCGCCATGCTCGa                                                                                                                                 <  1:427945/137‑1 (MQ=255)
tCGAAGTGGTGGTTCGGCTCGCTGGCCGGGGTGTTCGGCGCCGCATTGGGCATGCTTGCCACCGCCTTGCCCTGGCGCGCCGAAATCGCCCCGGTAACGCCACCTGGTGTCGGCAGCTGGAGCGCCGCCATGCTCGa                                                                                                                                 >  6:472474/1‑137 (MQ=255)
   aaGTGGTGGTTCGGCTCGCTGGCCGGGGTGTTCGGCGCCGCATTGGGCATGCTTGCCACCGCCTTGCCCTGGCGCGCCGAAATCGCCCCGGTAACGCCACCTGGTGTCGGCAGCTGGAGCGCCGCCATGCTCGAgcg                                                                                                                              >  4:377265/1‑137 (MQ=255)
           gtTCGGCTCGCTGGCCGGGGTGTTCGGCGCCGCATTGGGCATGCTTGCCACCGCCTTGCCCTGGCGCGCCGAAATCGCCCCGGTAACGCCACCTGGTGTCGGCAGCTGGAGCGCCGCCATGCTCGCGCGCGGGCGcc                                                                                                                      >  7:315094/1‑137 (MQ=255)
           gtTCGGCTCGCTGGCCGGGGTGTTCGGCGCCGCATTGGGCATGCTTGCCACCGCCTTGCCCTGGCGCGCCGAAATCGCCCCGGTAACGCCACCTGGTGTCGGCAGCTGGAGCGCCGCCATGCTCGAGCGCGGGCGcc                                                                                                                      >  6:228941/1‑137 (MQ=255)
           gtTCGGCTCGCTGGCCGGGGTGTTCGGCGCCGCATTGGGCATGCTTGCCACCGCCTTGCCCTGGCGCGCCGAAATCGCCCCGGTAACGCCACCTGGTGTCGGCAGCTGGAGCGCCGCCATGCTCGAGCGCGGGCGcc                                                                                                                      >  3:69793/1‑137 (MQ=255)
           gtTCGGCTCGCAGGCCGCCGTGTTCGGCGCCGCATTGGCCATGCTTGCCACCGCCTTGCCCTGGCGCGCCGAAATCGCCCCGGAAACGCCACCTGCTGTCGGCAGCTGGAGCGCCCCCATCCTCGAGCGCGCGCGcc                                                                                                                      >  1:325066/1‑137 (MQ=255)
                 cTCGCTGGCCGGGGTGTTCGGCGCCGCATTGGGCATGCTTGCCACCGCCTTGCCCTGGCGCGCCGAAATCGCCCCGGTAACGCCCCCTGGTGTCGGCAGCTGGAGCGCCGCCATGCTCGAGCGCGGGCGCCAGGTgg                                                                                                                >  7:114896/1‑137 (MQ=255)
                 cTCGCTGGCCGGGGTGTTCGGCGCCGCATTGGGCATGCTTGCCACCGCCTTGCCCTGGCGCGCCGAAATCGCCCCGGTAACGCCACCTGGTGTCGGCAGCTGGAGCGCCGCCATGCTCGAGCGCGGGCGCCAGGTgg                                                                                                                >  1:497068/1‑137 (MQ=255)
                      gggCCGGGGGGTTCGGCGCCGCATGGGGCATGCTTGCCACCGCCTTGCCCGGGCGCGCCTAAATCGCCCCGGTAACGCCACCTGGTGGCGGCCGCTGGAGCGCCGCCAGGCGCGAGCGCGGGCGCCAGGTggcggcg                                                                                                           <  7:463831/136‑1 (MQ=255)
                                  gcgcccccgcATGGGGCATGCTGCCCCCCGCCTGCCCCGGGCGCGCCTAATGCGCCCGGGTACCGCCCCCGGGGGGCGGCAGGTGGAGCGCCGCCAGCGCCGCGCGCGGGCGCCTGGGGGCGGCGGCCGGCGACTgc                                                                                               <  4:144685/131‑1 (MQ=255)
                                  gcGGCGCCGCATTGGCCAGGCGTGCCACCCCCGTGCCCGGCCGCGCCCAAATCGCCCCGGTACCGCCCCCGGGGGGGGGCAGCTGGAGCGCCTCCATGCTCGCGCGCGGGCGCCAGGGGGCGGCGGCCGGCGACTgc                                                                                               <  4:336946/136‑1 (MQ=255)
                                     gcgcCGCATTGGGCATGCTTGCCACCGCCTTGCCCTGGCGCGCCGAAATCGCCCCGGTAACGCCACCTGGTGTCGGCAGCTGGAGCGCCGCCATGCTCGAGCGCGGGCGCCAGGTGGCGGCGGCCGGCGACTGCGcc                                                                                            >  1:48034/1‑137 (MQ=255)
                                     gcgcCGCATTGGGCATGCTTGCCACCGCCTTGCCCTGGCGCGCCGAAATCGCCCCGGTAACGCCACCTGGTGTCGGCAGCTGGAGCGCCGCCATGCTCGAGCGCGGGCGCCAGGTGGCGGCGGCCGGCGACTGCGcc                                                                                            >  3:208950/1‑137 (MQ=255)
                                           cATTGGGCATGCTTGCCACCGCCTTGCCCTGGCGCGCCGAAATCGCCCCGGTAACGCCACCTGGTGTCGGCAGCTGGAGCGCCGCCATGCGCGCGCGCGGGCGCCAGGTGGCGGCGGCCGGCGACTGCGCCGTGTGc                                                                                      <  7:194938/137‑1 (MQ=255)
                                           cATGGGGCATGCTGCCCCCCGCCTTGCCCTGGCGCGCCTAAATCGCCCCGTTAACCCCACCTGGTGGCGGCAGCGGGAGCGCCGCCATGCTCGCGCGCGGGCGCCTGGTGGCGGCGGCCGGCGACTGCGCCGTGTGc                                                                                      <  6:397019/137‑1 (MQ=255)
                                              tGGGCATGCTTGCCACCGCCTTGCCCTGGCGCGCCGAAATCGCCCCGGTAACGCCACCTGGTGTCGGCAGCTGGAGCGCCGCCATGCTCGAGCGCGGGCGCCAGGTGGCGGCGGCCGGCGACTGCGCCGTGTGCCAt                                                                                   >  8:253520/1‑137 (MQ=255)
                                                        tGCCACCGCCTTGCCCTGGCGCGCCGAAATCGCCCCGGTAACGCCACCTGGTGTCGGCAGCTGGAGCGCCGCCATGCTCGAGCGCGGGCGCCAGGTGGCGGCGGCCGGCGACTGCGCCGTGTGCCATACGGTCAgcg                                                                         >  8:495788/1‑137 (MQ=255)
                                                            aCCGCCTTGCCCTGGCGCGCCGAAATCGCCCCGGTAACGCCACCTGGTGTCGGCAGCTGGAGCGCCGCCATGCTCGAGCGCGGGCGCCAGGTGGCGGCGGCCGGCGACTGCGCCGTGTGCCATACGGTCAGcggcgg                                                                     >  6:191523/1‑137 (MQ=255)
                                                             ccGCCTTGCCCTGGCGCGCCGAAATCGCCCCGGTAACGCCACCTGGTGTCGGCAGCTGGAGCGCCGCCATGCTCGAGCGCGGGCGCCAGGTGGCGGCGGCCGGCGACTGCGCCGTGTGCCATACGGTCAGCggcggc                                                                    >  8:174889/1‑137 (MQ=255)
                                                             ccGCCTTGCCCTGGCGCGCCGAAATCGCCCCGGTAACGCCACCTGGTGTCGGCAGCTGGAGCGCCGCCATGCTCGAGCGCGGGCGCCAGGTGGCGGCGGCCGGCGACTGCGCCGTGTGCCATACGGTCAGCggcggc                                                                    >  5:1690/1‑137 (MQ=255)
                                                                cgttGCCGTGGCGCGCCGAATGCGCCGGGGTAACGGCAGCTGGTGGCGGCAGCTGGAGCGCCGCTAGGCTGGCGGGCGGGGGCCTGGTGGCGGCGGCCGGCGACTGCGCCGGGTGCGATACGGTCAGCGGCGGCAAg                                                                 <  8:360475/135‑1 (MQ=255)
                                                                      ccTGGCGCGCCGAAATCGCCCCGGTAACGCCACCTGGTGTCGGCAGCTGGAGCGCCGCCATGCTCGAGCGCGGGCGCCAGGTGGCGGCGGCCGGCGACTGCGCCGTGTGCCATACGGTCAGCGGCGGCAAGGCCAAt                                                           >  3:179704/1‑137 (MQ=255)
                                                                             ggcCTAAATCGCCCCGGAACCCCAACCTGGTGCCGGCAGCTGGCGCGCCGCCAGGCGCGAGCGGGGGCCCCAGGGGGCGGCGGCCGGCGACTGCGCCGTGTGCCATACGGTCAGCGGCGGCAAGGCCAATGCCgggg                                                    <  8:242396/136‑1 (MQ=255)
                                                                             cgcCTAAATCGCCCCGGTAACGCCCCCTGGTGGCGGCAGCGGGCGCGCCGCCATGCTCGAGCGCGGGCGCCAGGTGGCGGCGGCCGGCGACTGCGCCGTGTGCCATACGGTCAGCGGCGGCAAGGCCAATGCCgggg                                                    <  8:223236/137‑1 (MQ=255)
                                                                             cgcCGAAATCGCCCCGGTAACGCCACCTGGTGTCGGCAGCTGGAGCGCCGCCATGCTCGAGCGCGGGCGCCAGGTGGCGGCGGCCGGCGACTGCGCCGTGTGCCATACGGTCAGCGGCGGCAAGGCCAATGCCgggg                                                    >  1:47025/1‑137 (MQ=255)
                                                                             cgcCGAAATCGCCCCGGTAACGCCACCTGGTGTCGGCAGCTGGAGCGCCGCCATGCTCGAGCGCGGGCGCCAGGTGGCGGCGGCCGGCGACTGCGCCGTGTGCCATACGGTCAGCGGCGGCAAGGCCAATGCCgggg                                                    >  6:29985/1‑137 (MQ=255)
                                                                             cgcCGAAACCGCCCGGGAAACCCCCCCTGGGGGCGGCAGCGGGAGCGCCGCCATGCTCGAGCGCGGGCGCCAGGGGGCGGCGGCCGGCGACTGCGCCGTGTGCCATACGGTCAGCGGCGGCAAGGCCAATGCCgggg                                                    <  4:140195/137‑1 (MQ=255)
                                                                             cgcCCAAATCGCCCCGGAACCGCCACGTGGTGTCGGCAGGTGGCGCGCCGCCAGGCTCGCGCGCGGGCGCCAGGTGGCGGCGGCCGGCGACTGCGCCGTGTGCCATACGGTCAGCGGCGGCAAGGCCAATGCCgggg                                                    <  7:257117/137‑1 (MQ=255)
                                                                                  aaaTCGCCCGGGTAACGCCACCTGGTGTCGGCTGCGGGAGCGCCGCCATGCTCGCGCGCGGGCGCCTGGTGGCGGCGGCCGGCGACTGCGCCGTGTGCCATACGGTCAGCGGCGGCAAGGCCAATGCCGGGGGGCTg                                               <  7:216520/137‑1 (MQ=255)
                                                                                  aaaTCGCCCCGGTAACGCCACCTGGTGTCGGCTGGTGGAGCGCCGCCATGCTCGCGCGCGGGCGCCAGGGGGCGGCGGCCGGCGACTGCGCCGTGTGCCATACGGTCAGCGGCGGCAAGGCCAATGCCGGGGGGCTg                                               <  5:474862/137‑1 (MQ=255)
                                                                                  aaaTCGCCCCGGAACCGCCACCTGGTGCCGGCAGCGGGCGCGCCGCCATGCGCGCGCGGGGGCGCCTGGGGGCGGCGGCCGGCGACTGCGCCGTGTGCCATACGGTCAGCGGCGGCAAGGCCAATGCCGGGGGGCTg                                               <  3:377265/137‑1 (MQ=255)
                                                                                        ccccGGTAACGCCCCCTGGGGGGGGCTGCTGGCGCGCCGCCATCGGCGAGCGCGGGCGCCTGGGGGCGGCGGCCGGCGACTGCGCCGTGTGCCATACGGTCAGCGGCGGCAAGGCCAATGCCGGGGGGCTGGCAATg                                         <  7:446330/137‑1 (MQ=255)
                                                                                                  ccccccTGGGGGCGGCAGCGGGCGCGCCGCCATGCGCGCGCGGGGGCGCCAGGGGGCGGCGGCCGGCGACTGCGCCGTGTGCCATACGGTCAGCGGCGGCAAGGCCAATGCCGGGGGGCTGGCAATGGATACGCCtt                               <  5:266802/133‑1 (MQ=255)
                                                                                                    ccccTGGTGGCGGCCGCTGGAGCGCCGCCATGCTCGCGCGCGGGCGCCAGGTGGCGGCGGCCGGCGACTGCGCCGTGTGCCATACGGTCAGCGGCGGCAAGGCCAATGCCGGGGGGCTGGCAATGGATACGCCTTTc                             <  4:215569/135‑1 (MQ=255)
                                                                                                      ccTGGTGTCGGCAGCTGGAGCGCCGCCATGCTCGAGCGCGGGCGCCAGGTGGCGGCGGCCGGCGACTGCGCCGTGTGCCATACGGTCAGCGGCGGCAAGGCCAATGCCGGGGGGCTGGCAATGGATACGCCTTTAgg                           >  2:250021/1‑137 (MQ=255)
                                                                                                       gtGGTGTCGGCAGCGGGAGCGCCGCCATGCTCGCGCGGGGGCGCCAGGGGGCGGCGGCCGGCGACTGCGCCGTGTGCCATACGGTCAGCGGCGGCAAGGCCAATGCCGGGGGGCTGGCAATGGATACGCCTTTCGGc                          <  2:128283/136‑1 (MQ=255)
                                                                                                       cTGGTGTCGGCAGCTGGCGCGCGGCCATGCTCGAGCGGGGGGGCCTGGTGGCGGCGGCCGGCGACTGCGCCGTGTGCCATACGGTCAGCGGCGGCAAGGCCAATGCCGGGGGGCTGGCAATGGATACGCCTTTCGGc                          <  2:48034/137‑1 (MQ=255)
                                                                                                            gtCGGCAGCGGCCGCGCCGCCTGGCGCGCGCGGGGGCGCCCGGGGGCGGCGGCCGGCGACTGCGCCGGGTGCCATACGGTCAGCGGCGGCAAGGCCAATGCCGGGGGGCTGGCAATGGATACGCCTTTCGGCACGtt                     <  1:315389/137‑1 (MQ=255)
                                                                                                            ggcGGCCGCTGGTGCGCCCGCAGGCGCGCGCGGGGGCGCCTGGGGGCGGCGGCCGGCGACTGCGCCGTGTGCCATACGGTCAGCGGCGGCAAGGCCAATGCCGGGGGGCTGGCAATGGATACGCCTTTCGGCACGtt                     <  7:206729/135‑1 (MQ=255)
                                                                                                                gccgCTGGCGCGCCGCCTTGCTCGCGCGCGGGCGCCTGGGGGCGGCGGCCGGCGACTGCGCCGTGTGCCATACGGTCAGCGGCGGCAAGGCCAATGCCGGGGGGCTGGCAATGGATACGCCTTTCGGCACGTTGTAc                 <  7:253520/134‑1 (MQ=255)
                                                                                                                gCAGCTGGAGCGCGGCCATGCTCGCGCGGGGGGGCCAGGTGGCGGCGGCCGGCGACTGCGCCGTGTGCCATACGGTCAGCGGCGGCAAGGCCAATGCCGGGGGGCTGGCAATGGATACGCCTTTCGGCACGTTGTAc                 <  8:114896/137‑1 (MQ=255)
                                                                                                                gCAGCTGGAGCGCCGCCATGCTCGAGCGCGGGCGCCAGGTGGCGGCGGCCGGCGACTGCGCCGTGTGCCATACGGTCAGCGGCGGCAAGGCCAATGCCGGGGGGCTGGCAATGGATACGCCTTTCGGCACGTTGTAc                 <  6:232261/137‑1 (MQ=255)
                                                                                                                   gggggcgcgcCGCCTTGCTCGCGGGGGGGCGCCTGGGGGCGGCGGCCGGCGACTGCGCCGTGTGCCATACGGTCAGCGGCGGCAAGGCCAATGCCGGGGGGCTGGCAATGGATACGCCTTTCGGCACGTTGTACAGc              <  2:424259/131‑1 (MQ=255)
                                                                                                                       gAGCGCCGCCATGCTCGAGCGCGGGCGCCAGGTGGCGGCGGCCGGCGACTGCGCCGTGTGCCATACGGTCAGCGGCGGCAAGGCCAATGCCGGGGGGCTGGCAATGGATACGCCTTTCGGCACGTTGTACAGCACCa          >  8:22821/1‑137 (MQ=255)
                                                                                                                        cgcgcCGCCAGGCGCGCGCGGGGGCGCCTGGGGGCGGCGGCCGGCGACTGCGCCGTGTGCCATACGGTCAGCGGCGGCAAGGCCAATGCCGGGGGGCTGGCAATGGATACGCCTTTCGGCACGTTGTACAGCACCaa         <  7:200008/136‑1 (MQ=255)
                                                                                                                               ccttGCTGGCGCGGGGGCGCCCGGGGGCGGCGGCCGGCGACTGCGCCGTGTGCCATACGGTCAGCGGCGGCAAGGCCAATGCCGGGGGGCTGGCAATGGATACGCCTTTCGGCACGTTGTACAGCACCAATATCAcc  <  8:315094/134‑1 (MQ=255)
                                                                                                                                        |                                                                                                                               
TCGAAGTGGTGGTTCGGCTCGCTGGCCGGGGTGTTCGGCGCCGCATTGGGCATGCTTGCCACCGCCTTGCCCTGGCGCGCCGAAATCGCCCCGGTAACGCCACCTGGTGTCGGCAGCTGGAGCGCCGCCATGCTCGAGCGCGGGCGCCAGGTGGCGGCGGCCGGCGACTGCGCCGTGTGCCATACGGTCAGCGGCGGCAAGGCCAATGCCGGGGGGCTGGCAATGGATACGCCTTTCGGCACGTTGTACAGCACCAATATCACC  >  KT2440/4446325‑4446588

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