Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24405,764,9690TG42.5% 31.7 / 11.9 40E610A (GAG→GCG) caiTL‑carnitine/gamma‑butyrobetaine antiporter
Reads supporting (aligned to +/- strand):  ref base T (3/20);  new base G (17/0);  total (20/20)
Fisher's exact test for biased strand distribution p-value = 2.57e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ATCGAGCACATCACGCACGATCTGGTCTGGGGCAAAGCCCATCAGGTCGTACTCCTGACTGCCGTCATGCAGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGCCTCGCTCGGCGCGGCCATATAGCCGTCAAGGCTCACTTCGTAGACGAACGGGTTGCCCTCTTCCATCATCACCCGCAAGCCCATCAGGTTGCGCGACTGGCCAACGCGGGTTTCCACCTCGAAGCCCTGCGTCTG  >  KT2440/5764839‑5765102
                                                                                                                                  |                                                                                                                                     
aTCGAGCACATCACGCACGATCTGGTCTGGGGCAAAGCCCATCAGGTCGTACTCCTGACTGCCGTCATGCAGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGCctcgctcg                                                                                                                                 <  8:197808/137‑1 (MQ=255)
aTCGAGCACATCACGCACGATCTGGTCTGGGGCAAAGCCCATCAGGTCGTACTCCTGACTGCCGTCATGCAGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGCctcgctcg                                                                                                                                 <  1:329156/137‑1 (MQ=255)
aTCGAGCACATCACGCACGATCTGGTCTGGGGCAAAGCCCATCAGGTCGTACTCCTGACTGCCGTCATGCAGGTAGACCTCGGCGCGGTAGAAGCACTGGCGCAGTTCAGGCTCACCCTCGACCGGCGCctcgctcg                                                                                                                                 <  7:160932/137‑1 (MQ=255)
 tCGAGCACATCACTCACGATCTGGTCTGGGGCAAAGCCCATCAGGTCGTACTCCTGACTGCCGTCATGCAGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGCCTCGCTCgg                                                                                                                                <  2:144010/137‑1 (MQ=255)
   gAGCACATCACGCACTATCTGGTCTGGGGCAAAGCCCATCAGGTCGTACTCCTGACTGCCGTCATGCAGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGCCTCGCTCGgcg                                                                                                                              <  3:428125/137‑1 (MQ=255)
   gAGCACATCACGCACTATCTGGTCTGGGGCAAAGCCCATCAGGTCGTACTCCTGACTGCCGTCATGCAGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGCCTCGCTCGgcg                                                                                                                              <  7:284185/137‑1 (MQ=255)
   gAGCACATCACGCACGATCTGGTCTGGGGCAAAGCCCATCAGGTCGTACTCCTGACTGCCGTCATGCAGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGCCGCGCgcggcg                                                                                                                              >  5:158651/1‑137 (MQ=255)
   gAGCACATCACGCACGATCTGGTCTGGGGCAAAGCCCATCAGGTCGTACTCCTGACTGCCGTCATGCAGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGCCGCGCGGGgcg                                                                                                                              >  3:215463/1‑137 (MQ=255)
     gCACATCACGCACGATCTGGTCTGGGGCAAAGCCCATCAGGTCGTACTCCTGACTGCCGTCATGCAGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGCCTCGCTCGGcgcg                                                                                                                            <  7:23791/137‑1 (MQ=255)
             cgcacgATCTGGTCTGGGGCAAAGCCCATCAGGTCGTACTCCTGACTGCCGTCATGCAGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGGCGCGGTGGGCGCGGccggatt                                                                                                                    >  1:109953/1‑132 (MQ=255)
              gcacgATCTGGTCTGGGGCAAAGCCCATCAGGTCGTACTCCTGACTGCCGTCATGCAGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGATCAGGCTCACCCTCGACCGGCGCCGCGCTCCCCGCGCCCaccgac                                                                                                                   >  1:193246/1‑132 (MQ=255)
                      tggtctggGGCAAAGCCCATCAGGTCGTACTCCTGACTGCCGTCATGCAGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGCCTCGCTCGGCGCGGCCATATAGCCGTCAAg                                                                                                           <  2:220372/137‑1 (MQ=255)
                      tggtctggGGCAAAGCCCATCAGGTCGTACTCCTGACTGCCGTCATGCAGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGCCTCGCTCGGCGCGGCCATATAGCCGTCAAg                                                                                                           <  2:330192/137‑1 (MQ=255)
                      tggtctggGGCAAAGCCCATCAGGTCGTACTCCTGACTGCCGTCATGCAGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGCCTCGCTCGGCGCGGCCATATAGCCGTCAAg                                                                                                           <  2:423151/137‑1 (MQ=255)
                      tggtctggGGCAAAGCCCATCAGGTCGTACTCCTGACTGCCGTCATGCAGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGCCTCGCTCGGCGCGGCCATATAGCCGTCAAg                                                                                                           <  4:358275/137‑1 (MQ=255)
                        gtctggGGCAAAGCCCATCAGGTCGTACTCCTGACTGCCGTCATGCAGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGCCGCGCGCGGCGCGGCCATATAGACGTCAAGGc                                                                                                         >  6:60388/1‑137 (MQ=255)
                                 aaaGCCCATCAGGTCGTACTCCTGACTGCCGTCATGCAGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGCCGCGGGCGGCGCGGCCCAATAGCCGTCCAGGCTCACTTCgg                                                                                                >  4:134001/1‑136 (MQ=255)
                                    gCCCATCAGGTCGTACTCCTGACTGCCGTCATGCAGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGCCTCGCTCGGCGCGGCCATATAGCCGTCAAGGCTCACTTCGTAGa                                                                                             <  7:397241/137‑1 (MQ=255)
                                     cccATCAGGTCGTACTCCTGACTGCCGTCATGCAGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGCCGCGCACGGCGCCGCCAGATATACCTCAACGCTCACTTCGTAGAc                                                                                            >  2:129761/1‑137 (MQ=255)
                                      ccATCAGGTCGTACTCCTGACTGCCGTCATGCAGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGGGGCGCGGGCCGGGCGGCGGGATGGCCGTCAACGCAGACTGCGTAGACg                                                                                           >  3:444528/1‑137 (MQ=255)
                                         tCAGGTCGTACTCCTGACTGCCGTCATGCAGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGCCTCGCTCGGCGCGGCCATATAGCCGTCAAGGCTCACTTCGTAGACGAAc                                                                                        <  3:136761/137‑1 (MQ=255)
                                         tCAGGTCGTACTCCTGACTGCCGTCATGCAGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGCCTCGCTCGGCGCGGCCATATAGCCGTCAAGGCTCACTTCGTAGACGAAc                                                                                        <  5:106426/137‑1 (MQ=255)
                                          cAGGTCGTACTCCTGACTGCCGTCATGCAGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGCCGCGCTCGGCGCGGCCAGATAGACGTCAAGGCTGACATCGTAGACGaccc                                                                                       >  2:160749/1‑134 (MQ=255)
                                          cAGGTCGTACTCCTGACTGCCGTCATGCAGGTACACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGGCGCGCGCGCCGCCCCCAGATAGACGTCAAGCCTCACAACGAAGActagac                                                                                       >  8:260537/1‑132 (MQ=255)
                                              tCGTACTCCTGACTGCCGTCATGCAGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGCCTCGCTCGGCGCGGCCATATAGCCGTCAAGGCTCACTTCGTAGACGAACGGGtt                                                                                   <  1:227036/137‑1 (MQ=255)
                                                     ccTGACTGCCGTCATGCAGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGCCGCGCTCGGCGCGGCCCTATAGCCGGCAAGGCGCACGTCGTAGACGAACGGGGTGGCcgct                                                                            >  4:346458/1‑137 (MQ=255)
                                                     ccTGACTGCCGTAATGCAGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGCCTCGCTCGGCGCGGCCATATAGCCGTCAAGGCTCACTTCGTAGACGAACGGGTTGCCctct                                                                            <  8:127248/137‑1 (MQ=255)
                                                           tGCCGTCATGCAGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGCCGCGCTCGGCGCGGCCATATAGCAGTCAAGGCTCACATCGTAGACGAACGGGGTGCCCTCGTAcata                                                                      >  6:1756/1‑136 (MQ=255)
                                                           tGCCGTCATGCAGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGCCGCGCGCGGCGCGGCCCGATAGCCGTCAAGGCTCACTTCGTAGACGAACGGGTTGCCCGCGTCcata                                                                      >  2:312022/1‑136 (MQ=255)
                                                              cGTCATGCAGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGCCGCGCGCGGCGCGGCGCTAGAGACGACAAGGCTCAATTCGGAGACGAACGGGGTGCCCTCTTGcagcctc                                                                   >  8:307586/1‑132 (MQ=255)
                                                                 cATGCAGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGCCGCGCGCGGCGCGGCCCTATAGCCGGCAAGGCGCACTTCGGAGGCGAACGGGTTGGCCGCTTGCATAATCAcc                                                                >  8:383536/1‑137 (MQ=255)
                                                                     cAGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGCCGCGCGCGGCGCGGCCAGATAGCCGGCAAGGCTCACATGGGAGTCGCACGGGGTGCGCTGGTCCAAGATCACCCGcc                                                            >  8:163418/1‑136 (MQ=255)
                                                                        gTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGCCTCGCTCGGCGCGGCCATATAGCCGTCAAGGCTCACTTCGTAGACGAACGGGTTGCCCTCTTCCATCATCACCCGCAAGc                                                         <  5:167526/137‑1 (MQ=255)
                                                                           gACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCCCCTCGCTCGGCGCGGCCATATAGCCGTCAAGGCTCACTTCGTAGACGAACGGGTTGCCCTCTTCCATCATCACCCGCAAGCCCa                                                      <  6:5103/137‑1 (MQ=255)
                                                                                cGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGCCTCGCTCGGCGCGGCCATATAGCCGTCAAGGCTCACTTCGTAGACGAACGGGTTGCCCTCTTCCATCATCACCCGCAAGCCCATCAgg                                                 <  2:71125/137‑1 (MQ=255)
                                                                                   cgcgGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGGCTCGCTCGGCGCGGCCAAAAAGCCGTCAACGCGCACTTCATTGACGAAAGGGTTGGCCGCTTGCATCATAACCCCCAACCCCATCACGTTg                                              >  8:343924/1‑137 (MQ=255)
                                                                                   cgcgGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGCCGCGCGCGGCGCGGCCCAAGAGCCGGCAAGGCGCACTTCGGAGACGAACGGGGTGCGCGCGTCCCTCATCACCCCCACGCCCATCAAGttt                                              >  6:417020/1‑136 (MQ=255)
                                                                                                   gcgcAGTTCAGGCTCACCCTCGACCGGCGCCTCGCTCGGCGCGGCCATATAGCCGTCAAGGCTCACTGCGTAGACGAACGGGTTGCCCTCTGGCATCATCACCCGCCAGCCCATCAGGTTGGGGGACTGGGCCAcgc                              >  8:186967/1‑137 (MQ=255)
                                                                                                                     cTCGACCGGCGCCTCGCTCGGCGCGGCCATATAGCCGTCAAGGCTCACTTCGTAGACGAACGGGTTGCCCTCTTCCATCATCACCCGCAAGCCCATCAGGTTGCGCGACTGGCCAACGCGGGTTTCCACCTCGAAGc            <  4:144097/137‑1 (MQ=255)
                                                                                                                     cTCGACCGGCGCCTCGCTCGGCGCGGCCATATAGCCGTCAAGGCTCACTTCGTAGACGAACGGGTTGCCCTCTTCCATCATCACCCGCAAGCCCATCAGGTTGCGCGACTGGCCAACGCGGGTTTCCACCTCGAAGc            <  1:129761/137‑1 (MQ=255)
                                                                                                                      tCGACCGGCGCCTCGCTAGGCGCGGCCATATAGACGGCAAGGCTCACTTCGTAGACGAAAGGGGTGGGCTCGTACATCATAACCCGCACGAACAGCACATTTCGCGTCTGGCGAAAGAGGGTTTGCACAACGaacac           >  6:274700/1‑134 (MQ=255)
                                                                                                                            ggCGCCTCGCTCGGCGCGGCCATATAGCCGTCAAGGCTCACTTCGTAGACGAACGGGTTGCCCTCTTCCATCATCACCCGCAAGCCCATCAGGTTGCGCGACTGGCCAACGCGGGTTTCCACCTCGAAGCCCTGCGt     <  8:27611/137‑1 (MQ=255)
                                                                                                                               gcCTCGCTCGGCGCGGCCATATAGCCGTCAAGGCTCACTTCGTAGACGAACGGGTTGCCCTCTTCCATCATCACCCGCAAGCCCATCAGGTTGCGCGACTGGCCAACGCGGGTTTCCACCTCGAAGCCCTGCGTCTg  <  5:1756/137‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                     
ATCGAGCACATCACGCACGATCTGGTCTGGGGCAAAGCCCATCAGGTCGTACTCCTGACTGCCGTCATGCAGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGCCTCGCTCGGCGCGGCCATATAGCCGTCAAGGCTCACTTCGTAGACGAACGGGTTGCCCTCTTCCATCATCACCCGCAAGCCCATCAGGTTGCGCGACTGGCCAACGCGGGTTTCCACCTCGAAGCCCTGCGTCTG  >  KT2440/5764839‑5765102

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