Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913650,4640AT43.6% 76.6 / 25.8 56R9R (CGT→CGAcitEcitrate lyase, citryl‑ACP lyase (beta) subunit
Reads supporting (aligned to +/- strand):  ref base A (1/30);  new base T (24/0);  total (26/30)
Fisher's exact test for biased strand distribution p-value = 1.00e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATGCGGCGGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGC  >  NC_000913/650317‑650610
                                                                                                                                                   |                                                                                                                                                    
cATGCGGCGGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACgt                                                                                                                                                   <  1:997777/151‑1 (MQ=255)
    cggcggGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACgttgtt                                                                                                                                               <  2:799766/151‑1 (MQ=255)
            ggTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGaa                                                                                                                                       <  2:312727/151‑1 (MQ=255)
            ggTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGGGGGGTTTTCgttttttgagc                                                                                                                                          >  1:742720/1‑143 (MQ=255)
             gtgtCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAg                                                                                                                                      <  2:199903/151‑1 (MQ=255)
               gtCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGc                                                                                                                                     <  1:1021234/150‑1 (MQ=255)
                   ttttCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGaa                                                                                                                                 <  1:515013/150‑1 (MQ=255)
                   ttttCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGa                                                                                                                                  <  1:812587/149‑1 (MQ=255)
                    tttCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGaaa                                                                                                                                <  2:857561/150‑1 (MQ=255)
                      tCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGCGGGTTTTTCGGTTTTGTAGCGaaaagg                                                                                                                                   >  2:683639/1‑141 (MQ=255)
                       cACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGAGAAGCGGAAATCAt                                                                                                                          <  2:882820/151‑1 (MQ=255)
                        aCGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGCGGGGTTTTCGTTGTTTGCGCGAAACGGaaa                                                                                                                                >  2:690270/1‑146 (MQ=255)
                             aTGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCaa                                                                                                                    <  1:834591/151‑1 (MQ=255)
                               gCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAAt                                                                                                                   <  2:197563/150‑1 (MQ=255)
                                 tACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTCGGAGTTTTTTGAAGGTTCAGCGAAAAGGGAAACATTGGGaa                                                                                                                    >  2:418369/1‑147 (MQ=255)
                                    ggAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTccc                                                                                                             <  2:795921/151‑1 (MQ=255)
                                    ggAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGGGTTTAGTTTTTGCAGCGGAACGGAGATAAATGGGaaaaccccc                                                                                                             >  1:585812/1‑144 (MQ=255)
                                    ggAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGGGGGGTTTTTGTTTTTGCAGCGAAACGGGAAACATTTGGCATCCTccc                                                                                                             >  1:2495/1‑151 (MQ=255)
                                      aGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCAt                                                                                                           <  2:581554/151‑1 (MQ=255)
                                       gTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATg                                                                                                          <  2:205044/151‑1 (MQ=255)
                                       gTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGCGGGGTTTTCGTTGTTTCAGCGGAGCGGAAAACATTTGCCATCCTcccc                                                                                                            >  2:78701/1‑148 (MQ=255)
                                          ttCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGCGGGATTTTCGTTGTTGCAGAGCAGAGGAAATAAATggtgaacccccc                                                                                                             >  2:785965/1‑135 (MQ=255)
                                            cGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCaga                                                                                                     <  2:269206/151‑1 (MQ=255)
                                             gAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCagag                                                                                                    <  1:58179/151‑1 (MQ=255)
                                                    aaCATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGGGGGGTTTTCGTTGTTGGTGCGGAACGGGAAACAATGGGCAACCCCCCCCGGCAGAGGCggg                                                                                               >  1:275776/1‑149 (MQ=255)
                                                    aaCATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGCGAGTTTTTCGTTGTTGCAGCGAAGCGGAAAACAATGGCAATCCTCCCCTGGCAGAGCCggg                                                                                               >  2:231153/1‑149 (MQ=255)
                                                       aTCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGCGCGGTTTGTGGTTTTTGTGAGACACGGGAATCATTTTGAATCCCCCCCCGGCGGGGCCGGGATaa                                                                                           >  2:1009008/1‑149 (MQ=255)
                                                        tCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGCGGGGTTTTCGTTGTTGCGGCGGAGCGGGAAACATTTGGAAACCCCCCCTGGCAGAGCCgggggac                                                                                           >  1:70644/1‑145 (MQ=255)
                                                         cAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGGGGGTTTTTCGTTGTTGCAGCGAAGCGGGAAACAATTGCAATCCTTCCACGGCGGAGCCGGGAGACCGc                                                                                        >  1:757645/1‑151 (MQ=255)
                                                            ggCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGGGTTTTTAGTTGTTGCAGCGAAGCGGGAATCATTGGCAATCCTCCCACGGCAGAGCCGGGAGAcc                                                                                          >  1:338775/1‑146 (MQ=255)
                                                                  aGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTg                                                                               <  1:616399/151‑1 (MQ=255)
                                                                    ccGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGGGGGTTTTTCTTTTTTTCGCCGAAACGGGAAACAATGGGCAACCCCCCATGGCAGGGCCGGGGTACCCCCCGCGCGCGcc                                                                             >  1:21608/1‑151 (MQ=255)
                                                                        gTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCcagca                                                                         <  1:678556/151‑1 (MQ=255)
                                                                          aGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGCGGGGTTTTCGGTGTTGCAGCGAAGCGGGAAACATTGGGAATCCTCCCCTGGCAGAGCCGGGGTAACGCTGGCGCGGGCCCGCAgg                                                                       >  2:248285/1‑151 (MQ=255)
                                                                            aTGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGGGGTTTTCGTTTTTTCAGCGAAGCGGGAAACAATTGCAATCCCCCCCCGGCAGGGCCGGGGTAACGCCGGCGCGGGCCCGCAggg                                                                      >  1:705684/1‑150 (MQ=255)
                                                                              gAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCt                                                                    <  1:353546/150‑1 (MQ=255)
                                                                                     ttGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGGGGTTTTCGTTGTTGCAGCGAAGCGGAAAACAATGGCAATCCTCCCCTGGCAGAGCCGGGGGACCGCTGGCGCGTGCCAGCCGGGCTTcccg                                                               >  2:114234/1‑146 (MQ=255)
                                                                                      tGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTca                                                            <  1:708791/150‑2 (MQ=255)
                                                                                         tGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTATGTTTTTGCAGGGAAGAGGAAATCAATTGGCAACCTCCCATCGCAGAGGCCGGGTAACGCCGGCGCGGGGCCGCAGGGGTTTCaat                                                              >  2:555449/1‑143 (MQ=255)
                                                                                            ccATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGGGTTTTTCGTTGTTGCAGCGAAGCGGGAAACATTGGCAATCCTCCCATGGCCGAGCCGGGGTACCGCCGGCGCGGGCCAGCAGGGCTTCCAGTCGTTCACGt                                                     >  2:319229/1‑151 (MQ=255)
                                                                                                cgcgGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTTCGGTTTTGCAGCGAAACGGGAAACATTTGCAAACCTCCCCTGGCAGAGCCGGGGGACCGCTGGCGCGGGCCAGCAGGGCTTCCCGTCGTGCACGTaaa                                                  >  1:387239/1‑150 (MQ=255)
                                                                                                    gCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCa                                             <  2:484314/151‑1 (MQ=255)
                                                                                                        tGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGGTTTTCGTTGTTTCAGCGAAGCGGAAAACAATGGCAATCCCCCCCTGGCAGAGCCCGGGTAACCCCGGCGCGGGCCCGCAGGGCTTCCAGTCGGGCACGGAAAAAGCAGGcc                                         >  2:913422/1‑151 (MQ=255)
                                                                                                        tGGCACCAGGCACAAACAACATGCTGCGGCGGGGGGGGGGTTTTCGTTGTTGCAGCGAAGAGGGAAACATTTGGCAACCTCCCCTGGCAGAGGCGGGGTACCGCTGGCGCGGGCCCGCAGGGCCTCCCGTCGTGCCCGTAAAAAGCAGGcc                                         >  2:964116/1‑151 (MQ=255)
                                                                                                         ggCACCAGGCACAAACAACATGCTGCGGCGGGGGGGGGTTTTTCCATGTTGATACGAAGCGGAAAACATTTGCAATTCCCCCCTGGCAGGGCCGGGGTACCCCGCGGGCGGGCCAGCCGGGGCTCCCCTCGTGCCCggacaa                                                 >  1:183061/1‑137 (MQ=255)
                                                                                                           cACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGt                                      <  2:714318/151‑1 (MQ=255)
                                                                                                               aGGCACAAACAACATGCTGCGGCGGGTGGGGGGTTTTCGCTTTTGGAGGGAAGAGGGAAACCTTGGGAATCCCCCCATGTCAGAGCCGGGGAATCGCCGTGGCGTGCCCACACGGCTTCCCCTGGGGCACGAAAAAAGTAGTcccctg                                     >  1:629432/1‑144 (MQ=255)
                                                                                                                ggCACAAACAACATGCTGCGGCGGGGGGGAGTTTTTAGGTGTTGCAGCGAAGGGGGAAACATTAGCCATCCTCCCCTGGCAGAAACGGGGAACCCCCGGCGCGTGCCAGCCGGGGTTCCCGTCGTGCACGTAAAAAGCAGTcccg                                       >  2:433678/1‑143 (MQ=255)
                                                                                                                 gCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCt                                <  2:58291/151‑1 (MQ=255)
                                                                                                                      aacaacATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTc                           <  2:70652/151‑1 (MQ=255)
                                                                                                                        caacaTGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCAt                         <  2:2495/151‑1 (MQ=255)
                                                                                                                              gctgcGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACAt                   <  2:210600/151‑1 (MQ=255)
                                                                                                                                     cggGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCta            <  1:634073/151‑2 (MQ=255)
                                                                                                                                        gTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTAc         <  2:884381/151‑1 (MQ=255)
                                                                                                                                              gTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCcgcg   <  1:690217/151‑1 (MQ=255)
                                                                                                                                                tttACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCgcgc  <  1:948570/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CATGCGGCGGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGC  >  NC_000913/650317‑650610

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: