Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1435,768,3320TG45.1% 48.3 / 13.1 53S609R (AGC→CGC) caiTL‑carnitine/gamma‑butyrobetaine antiporter
Reads supporting (aligned to +/- strand):  ref base T (3/25);  new base G (22/1);  total (27/26)
Fisher's exact test for biased strand distribution p-value = 3.19e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TCGAGCACATCACGCACGATCTGGTCTGGGGCAAAGCCCATCAGGTCGTACTCCTGACTGCCGTCATGCAGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGCCTCGCTCGGCGCGGCCATATAGCCGTCAAGGCTCACTTCGTAGACGAACGGGTTGCCCTCTTCCATCATCACCCGCAAGCCCATCAGGTTGCGCGACTGGCCAACGCGGGTTTCCACCTCGAAGCCCTGCGTCTG  >  ACB143/5768199‑5768461
                                                                                                                                     |                                                                                                                                 
tCGAGCACATCACGCACGATCTGGTCTGGGGCAAAGCCCATCAGGTCGTACTCCTGACTGCCGTCATGCAGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGCctcgctcg                                                                                                                                 <  5:320857/136‑1 (MQ=255)
 cGAGCACATCACGCACGATCTGGTCTGGGGCAAAGCCCATCAGGTCGTACTCCTGACTGCCGTCATGCAGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGCctcgctcg                                                                                                                                 <  7:62501/135‑1 (MQ=255)
  gAGCACATCACGCACGATCTGGTCTGGGGCAAAGCCCATCAGGTCGTACTCCTGACTGCCGTCATGCAGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGCCGCGCGCgg                                                                                                                                >  4:334636/1‑135 (MQ=255)
  gAGCACATCACGCACGATCTGGTCTGGGGCAAAGCCCATCAGGTCGTACTCCTGACTGCCGTCATGCAGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGCAGCGCGCgg                                                                                                                                >  3:95856/1‑135 (MQ=255)
  gAGCACATCACGCACGATCTGGTCTGGGGCAAAGCCCATCAGGTCGTACTCCTGACTGCCGTCATGCAGGTAGACCTCGGCGCGGTAGAAGCGCTGGAGCAGTTCAGGCTCACCCTCGACCGGCGCCGCGCCCgg                                                                                                                                >  4:220377/1‑135 (MQ=255)
  gAGCACATCACGCACGATCTGGTCTGGGGCAAAGCCCATCAGGTCGTACACCTGACTGCCGTCATGCAGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGAGCCTCCCTCTgc                                                                                                                               >  6:306150/1‑136 (MQ=255)
   aGCACATCACGCACGATCTGGTCTGGGGCAAAGCCCATCAGGTCGTACTCCTGACTGCCGTCATGCAGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGCCTCGCTCgg                                                                                                                                <  1:10140/134‑1 (MQ=255)
            cgcacgATCTGGTCTGGGGCAAAGCCCATCAGGTCGTACTCCTGACTGCCGTCATGCAGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGCCGCGCGGGGCCCGGcccg                                                                                                                       >  7:294182/1‑132 (MQ=255)
            cgcacgATCTGGTCTGGGGCAAAGCCCATCAGGTCGTACTCCTGACTGCCGTCATGCAGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGCCGCGCGCGGCGCGGccct                                                                                                                       >  8:375332/1‑132 (MQ=255)
            cgcacgATCTGGTCTGGGGCAAAGCCCATCAGGTCGTACTCCTGACTGCCGTCATGCAGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGCCGCGCGCGGCGCGGcccga                                                                                                                      >  5:182792/1‑132 (MQ=255)
            cgcacgATCTGGTCTGGGGCAAAGCCCATCAGGTCGTACTCCTGACTGCCGTCATGCAGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGCCGCGCGCGGCGCGGCCag                                                                                                                       >  7:74141/1‑133 (MQ=255)
                 gATCTGGTCTGGGGCAAAGACCATCAGGTCGTACTCCTGACTGCCGTCATGCAGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGTCGCCGCGCGCGCCGCCGCCCTATAAcc                                                                                                                 >  8:85254/1‑135 (MQ=255)
                     tggtctggGGCAAAGCCCATCAGGTCGTACTCCTGACTGCCGTCATGCAGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGTCTCACCCTCGACCGGCGCCTCGCGCGGCGCGGCCATAAAGACggaaa                                                                                                            >  6:71394/1‑132 (MQ=255)
                     tggtctggGGCAAAGCCCATCAGGTCGTACTCCTGACTGCCGTCATGCAGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGGGGCGCGCGGGGCGAGGCCATATAGCAGTCaa                                                                                                            >  6:167178/1‑136 (MQ=255)
                         ctggGGCAAAGCCCATCAGGTCGTACTCCTGACTGCCGTCATGCAGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGCCTCGCTCGGCGCGGCCATATAGCCGTCAAGGc                                                                                                         <  2:364117/135‑1 (MQ=255)
                          tggGGCAAAGCCCATCAGGTCGTACTCCTGACTGCCGTCATGCAGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGCCTCCCAGGGCGCGGCCATATACACGTCAACGca                                                                                                        >  7:320341/1‑134 (MQ=255)
                          tggGGCAAAGCCCATCAGGTCGTACTCCTGACTGCCGTCATGCAGGTAGACATCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGGGGCGCGCGGCGCGGCCAGATAGCCGGCAggggt                                                                                                        >  8:38634/1‑133 (MQ=255)
                           ggggCAAAGCCCATCAGGTCGTACTCCTGACTGCCGTCATGCAGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGGGCCGCCCGCGGCGCCGCCCTATACCCGTCAAAGCtaa                                                                                                      >  6:72943/1‑134 (MQ=255)
                                aaaGCCCATCAGGTCGTACTCCTGACTGCCGTCATGCAGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGCCGCGCGCGGCGCGGCCCTAGAGCCGGCACAGCGCACTTCg                                                                                                 >  6:59706/1‑136 (MQ=255)
                                    cccATCAGGTCGTACTCCTGACTGCCGTCATGCAGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGGCTCGCGCGGCGCGGGCCGAAAGCCGGCAACGCTCACTTCGTAg                                                                                              >  5:205496/1‑135 (MQ=255)
                                    cccATCAGGTCGTACTCCTGACTGCCGTCATGCAGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGCCGCGCGCGGCGCGGGCAGATAGCCGGCAAGGCGCCCATCGTAg                                                                                              >  2:7856/1‑135 (MQ=255)
                                    cccATCAGGTCGTACTCCTGACTGCCGTCATGCAGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGCCGCCCGCGGCGCGGCCCGCAAGCCGGCAAGGCTCCCTTCGTAg                                                                                              >  8:133727/1‑135 (MQ=255)
                                     ccagccggTCGTAATCCTGACTGCCGTCATGCAGTTAGACCTCGGCGCTTTTGAAGAGCTGGGGCAGTTCAGGCTCACCCTCGACCGGCGGCTCGCGCGGCGCGGCAATATAGCCGTCAAGGCTCACTTCGTCGa                                                                                             <  3:11275/129‑1 (MQ=255)
                                          aGGTCGTACTCCTGACTGCCGTCATGCAGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGCCTCGCTCGGCGCGGCCATATAGCCGTCAAGGCTCACTTCGTAGACGAAc                                                                                        <  5:279415/135‑1 (MQ=255)
                                          aGGTCGTACTCCTGACTGCCGTCATGCAGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGCCGCGCGCGGCGCGGCCCGATACACGACAAGGCTCCCTTAGTAGACgcaa                                                                                        >  8:76500/1‑132 (MQ=255)
                                           ggTCGTACTCCTGACTGCCGTCATGCAGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGCCTCGCTCGGCGCGGCCATATAGCCGTCAAGGCTCACTTCGTAGACGAAc                                                                                        <  3:131481/134‑1 (MQ=255)
                                           ggTCGTACTCCTGACTGCCGTCATGCAGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGCCTCGCTCGGCGCGGCCATATAGCCGTCAAGGCTCACTTCGTAGACGAACg                                                                                       <  7:76500/135‑1 (MQ=255)
                                             tCGTACTCCTGACTGCCGTCATGCAGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGCCTCGCTCGGCGCGGCCATATAGCCGTCAAGGCTCACTTCGTAGACGAACggg                                                                                     <  2:240228/135‑1 (MQ=255)
                                               ggaCTCCTGACTGCCGTCATGCAGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGCCTCGCTCGGCGCGGCCATATAGCCGTCAAGGCTCACTTCGTAGACGAACGGGtt                                                                                   <  6:359788/133‑1 (MQ=255)
                                               gTACTCCTGACTGCCGTCATGCAGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGCCTCGCTCGGCGCGGCCATATAGCCGTCAAGGCTCACTTCGTAGACGAACGGGtt                                                                                   <  5:92550/135‑1 (MQ=255)
                                               gTACTCCTGACTGCCGTCATGCAGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGCCTCGCTCGGCGCGGCCATATAGCCGTCAAGGCTCACTTCGTAGACGAACGGGtt                                                                                   <  6:182695/135‑1 (MQ=255)
                                                tACTCCTGACTGCCGTCATGCAGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGCCTCGCTCGGCGCGGCCATATAGCCGTCAAGGCTCACTTCGTAGACGAACGGGtt                                                                                   <  4:21340/134‑1 (MQ=255)
                                                tACTCCTGACTGCCGTCATGCAGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGCCTCGCTCGGCGCGGCCATATAGCCGTCAAGGCTCACTTCGTAGACGAACGGGtt                                                                                   <  8:30941/134‑1 (MQ=255)
                                                 aCTCCTGACTGCCGTCATGCAGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGGGGCGCGCGCGGCGCGGCCCGCTAGCCGGCAACGCTGACATCGGAGGCGAACGGGGTgg                                                                                 >  1:341286/1‑134 (MQ=255)
                                                 aCTCCTGACTGCCGTCATGCAGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGCCTCGCTCGGCGCGGCCATATAGCCGTCAAGGCTCACTTCGTAGACGAACGGGTTGc                                                                                 <  7:322102/135‑1 (MQ=255)
                                                    ccTGACTGCCGTCATGCAGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGCCGCGCGCGGCGCGGCCCGATAGCCGTCAAGGCTCACTACGTAGACGAACGGGGTGgccgc                                                                             >  5:158375/1‑134 (MQ=255)
                                                        aCTGCCGTCATGCAGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGCCTCGCTCGGCGCGGCCATATAGCCGTCAAGGCTCAAGTCGTAGACGAACGGGTTGCCcgctgc                                                                          >  1:344670/1‑133 (MQ=255)
                                                             cGTCATGCAGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGCCTCGCTCGGCGCGGCCATATAGCCGTCAAGGCTCACTTCGTAGACGAACGGGTTGCCCTCTTCcat                                                                       <  4:119141/133‑1 (MQ=255)
                                                                  tGCAGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGCCGCGCGCGGCGCGGCCAAATAGCCGTCAAGGCTCACGTCGTAGACGCACGGGTGGCGCTCGTCCAGCATCAccc                                                               >  6:72541/1‑136 (MQ=255)
                                                                     aGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGCCTCGCTCGGCGCGGCCATATAGCCGTCAAGGCTCACTTCGTAGACGAACGGGTTGCCCTCTTCCATCATCACCCg                                                              <  2:354739/134‑1 (MQ=255)
                                                                       gTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGGCGCGCGCGGCGCGGCCCGAAAGCCGTAAAGGCGCACGGCGGAGAAGAACCGGGTGCGCGCGGCCCTCCTCAccccccc                                                           >  2:239213/1‑131 (MQ=255)
                                                                                 gcgcGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGCCTCGCTCGGCGCGGCCATATAGCCGTCAAGGCTCACTTCGTAGACGAACGGGTTGCCCTCTTCCATCATCACCCGCAAGCCCATCAgg                                                 <  8:191599/135‑1 (MQ=255)
                                                                                   gcgGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGCCTCGCTCGGCGCGGCCATATAGCCGTCAAGGCTCACTTCGTAGACGAACGGGTTGCCCTCTTCCATCATCACCCGCAAGCCCATCAGGtt                                               <  1:239213/135‑1 (MQ=255)
                                                                                                    gcAGTTCAGGCTCACCCTCGACCGGCGCCGCGCGCGGCGCGGCCCTATATCCGTGAAAGCGCACATAGTTGACGAACCGGGTGGCCTCGTCCATCATAACCCCCCACCCCAACACGTTGCGCGGCTGGCCCAcac                              >  4:286237/1‑133 (MQ=255)
                                                                                                        ttCAGGCTCACCCTCGACCGGCGCCTCGCTCGGCGCGGCCATATAGCCGTCAAGGCTCACTTCGTAGACGAACGGGTTGCCCTCTTCCATCATCACCCGCAAGCCCATCAGGTTGCGCGACTGGCCAACGCGGGt                          <  7:282861/135‑1 (MQ=255)
                                                                                                        ttCAGGCTCACCCTCGACCGGCGCCTCGCTCGGCGCGGCCATATAGCCGTCAAGGCTCACTTCGTAGACGAACGGGTTGCCCTCTTCCATCATCACCCGCAAGCCCATCAGGTTGCGCGACTGGCCAACGCGGGt                          <  7:362492/135‑1 (MQ=255)
                                                                                                        ttCAGGCTCACCCTCGACCGGCGCCTCGCTCGGCGCGGCCATATAGCCGTCAAGGCTCACTTCGTAGACGAACGGGTTGCCCTCTTCCATCATCACCCGCAAGCCCATCAGGTTGCGCGACTGGCCAACGCGGGt                          <  2:229025/135‑1 (MQ=255)
                                                                                                              cTCACCCTCGACCGGCGCCTCGCTCGGCGCGGCCATATAGCCGTCAAGGCTCACTTCGTAGACGAACGGGTTGCCCTCTTCCATCATCACCCGCAAGCCCATCAGGTTGCGCGACTGGCCAACGCGGGTTTCCAcc                   <  5:216045/136‑1 (MQ=255)
                                                                                                                     tCGACCGGCGCCGCGCGCGGCGCGGCCCTATAGCCGTCAACGCTCACTTCGTAGACGAACGGGGTGGGCGCGTGCATCATCACCCCCAACCCCCACAGGTTGCGCGCCTGGGCCACCCCGGGTTGCACCACGAAGc            >  6:72900/1‑136 (MQ=255)
                                                                                                                        aCCGGCGCCTCGCTCGGCGCGGCCATATAGCCGTCAAGGCTCACTTCGTAGACGAACGGGTTGCCCTCTTCCATCATCACCCGCAAGCCCATCAGGTTGCGCGACTGGCCAACGCGGGTTTCCACCTCGAAGc            <  2:341286/133‑1 (MQ=255)
                                                                                                                        aCCGGCGCCTCGCTCGGCGCGGCCATATAGCCGTCAAGGCTCACTTCGTAGACGAACGGGTTGCCCTCTTCCATCATCACCCGCAAGCCCATCAGGTTGCGCGACTGGCCAACGCGGGTTTCCACCTCGAAGc            <  3:10291/133‑1 (MQ=255)
                                                                                                                        aCCGGCGCCTCGCTCGGCGCGGCCATATAGCCGTCAAGGCTCACTTCGTAGACGAACGGGTTGCCCTCTTCCATCATCACCCGCAAGCCCATCAGGTTGCGCGACTGGCCAACGCGGGTTTCCACCTCGAAGc            <  1:249717/133‑1 (MQ=255)
                                                                                                                             cgcCTCGCTCGGCGCGGCCATATAGCCGTCAAGGCTCACTTCGTAGACGAACGGGTTGCCCTCTTCCATCATCACCCGCAAGCCCATCAGGTTGCGCGACTGGCCAACGCGGGTTTCCACCTCGAAGCCCTGCGt     <  1:48164/135‑1 (MQ=255)
                                                                                                                               ccgcgcgcgGCGCGGCCAGATAGCCGTCAAGGCTCACTTCGTAGACGAACGGGTTGCCCTCGTCCATCATAACCCCCCAGCCCAGCAGGTTGCGCGCCTGGCCCACGCGGGGGTCCACCACTCAGCCCGGCGTCt   >  2:125513/4‑135 (MQ=255)
                                                                                                                               ccTCGCTCGGCGCGGCCATATAGCCGTCAAGGCTCACTTCGTAGACGAACGGGTTGCCCTCTTCCATCATCACCCGCAAGCCCATCAGGTTGCGCGACTGGCCAACGCGGGTTTCCACCTCGAAGCCCTGCGTCTg  <  5:218218/136‑1 (MQ=255)
                                                                                                                               ccTCGCTCGGCGCGGCCATATAGCCGTCAAGGCTCACGTCGTAGACGACCGGGGTGCCCGCGTCCATCATCACCCGCACGCCCAACACGTTGCGCGCCTGGGCAACGCGGGGGTCCACCTCGAAGCCCTGCGGCt   >  8:223434/1‑135 (MQ=255)
                                                                                                                                ctcgctcgGCGCGGCCATATAGCCGTCAAGGCTCACTTCGTAGACGACCGGGTTGGCCGCTTCCATCATAACCCCCCAGCCCAAAAGGTTGCGCGCCTGGCCAACGCGGGTTTCCCCCTCGAAGCCCTGCGTCt   >  8:287052/1‑134 (MQ=255)
                                                                                                                                ctcgctcgGCGCGGCCATATAGCCGTCAAGGCTCACTTCGTAGACGAACGGGTTGCCCTCTTCCATCATCACCCGCAAGCCCATCAGGTTGCGCGACTGGCCAACGCGGGTTTCCACCTCGAAGCCCTGCGTCTg  <  4:305595/135‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                 
TCGAGCACATCACGCACGATCTGGTCTGGGGCAAAGCCCATCAGGTCGTACTCCTGACTGCCGTCATGCAGGTAGACCTCGGCGCGGTAGAAGCGCTGGCGCAGTTCAGGCTCACCCTCGACCGGCGCCTCGCTCGGCGCGGCCATATAGCCGTCAAGGCTCACTTCGTAGACGAACGGGTTGCCCTCTTCCATCATCACCCGCAAGCCCATCAGGTTGCGCGACTGGCCAACGCGGGTTTCCACCTCGAAGCCCTGCGTCTG  >  ACB143/5768199‑5768461

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: ‑