Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913650,4640AT40.8% 69.6 / 19.6 49R9R (CGT→CGA) citEcitrate lyase, citryl‑ACP lyase (beta) subunit
Reads supporting (aligned to +/- strand):  ref base A (6/23);  new base T (20/0);  total (26/23)
Fisher's exact test for biased strand distribution p-value = 1.63e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCATGCGGCGGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGC‑GAGTTTTACGTTGTTGCAGCGAAGCGGAAATCA‑TTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACG  >  NC_000913/650315‑650613
                                                                                                                                                      |                                                                                                                                                      
aCCATGCGGCGGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGC‑GAGTTTTAc                                                                                                                                                       <  1:876855/151‑1 (MQ=255)
         cgggcggTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGG‑GGGTTTTTCgttgtt                                                                                                                                                 >  2:755935/1‑148 (MQ=255)
               gtgtCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGC‑GAGTTTTACGTTGTTGCAGCGAAg                                                                                                                                        <  1:145829/151‑1 (MQ=255)
                      tttCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGC‑GAGTTTTACGTTGTTGCAGCGAAGCGGaaa                                                                                                                                  <  2:868114/150‑1 (MQ=255)
                        tCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGGGGTTTTTCG‑TGTTGCCGCGAAGCGGgaaac                                                                                                                                >  1:696762/1‑148 (MQ=255)
                                tGCTACGGAGTCTTCGCGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGCGC‑GCGCTTTTCGTTGTTGCAGCGAAGAGGAAAGCA‑CTggt                                                                                                                         >  1:275203/1‑147 (MQ=255)
                                   tACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGG‑GGGGTTTTTGTTGTTTCCGCGAAACGGGAATCATTTGGCAAccc                                                                                                                    >  1:396515/1‑151 (MQ=255)
                                    aCGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGC‑GGGGTTTTCGTTTTTTCCGCGAAGCGGAAAACA‑TTTGCAAtcc                                                                                                                    >  1:942348/1‑149 (MQ=255)
                                      ggAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGC‑GAGTTTTACGTTGTTGCAGCGAAGCGGAAATCA‑TTGGCAATCCTccc                                                                                                                <  1:583996/151‑1 (MQ=255)
                                            ttCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGG‑GGGTTTTTTGTTGTTGCCGCGGAACGGGAAACA‑CTTGGAAACCTccc                                                                                                                >  2:866481/1‑145 (MQ=255)
                                              cGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGC‑GAGTTTTACGTTGTTGCAGCGAAGCGGAAATCA‑TTGGCAATCCTCCCATGGCaga                                                                                                        <  2:824609/151‑1 (MQ=255)
                                              cGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGC‑GAGGTTTAACTTGTTTCAGAGAAGTGGGGATCT‑TTGGGAAACCTcccccggc                                                                                                           >  1:785391/1‑143 (MQ=255)
                                              cGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGC‑GCGGTTTTTGTTGTTGGAGCGAAGCGGAAAACA‑TTGGCAATCCTcccctcgc                                                                                                           >  1:452438/1‑143 (MQ=255)
                                                 ggTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGC‑GAGTTTTACGTTGTTGCAGCGAAGCGGAAATCA‑TTGGCAATCCTCCCATGGCAGAGcc                                                                                                     <  1:846433/151‑1 (MQ=255)
                                                  gTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGC‑GAGTTTTACGTTGTTGCAGCGAAGCGGAAATCA‑TTGGCAATCCTCCCATGGCAGAGCCg                                                                                                    <  1:439009/151‑1 (MQ=255)
                                                     aaaCATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGC‑GAGTTTTACGTTGTTGCAGCGAAGCGGAAATCA‑TTGGCAATCCTCCCATGGCAGAGCCGGGa                                                                                                 <  2:24778/151‑1 (MQ=255)
                                                                cATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGC‑GGGTTTTTCGTTGTTTGAGCGGAACGGAAATCA‑TTTGCAATCCCCCCCTGGCAGAGCCCGGGTAccccccg                                                                                        >  1:958836/1‑144 (MQ=255)
                                                                cATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGC‑GGGATTTTTGGTGTTGCAGGGTAATGGAGATCA‑TTTGCAATCCTCCCCTGGGAGCGACGGGGAATCGCTg                                                                                         >  2:275826/1‑148 (MQ=255)
                                                                   cAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGC‑GAGTTTTACGTTGTTGCAGCGAAGCGGAAATCA‑TTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGt                                                                                   <  1:127835/151‑1 (MQ=255)
                                                                   cAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGC‑GGGGTTTTTGTTGTTTCAGCGAAGCGGGAAACATTTGGCAA‑ACTCCCCTTGCAGAGGCGGGGTACCGCCgg                                                                                        >  1:707103/1‑146 (MQ=255)
                                                                     gCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGC‑GAGTTTTACGTTGTTGCAGCGAAGCGGAAATCA‑TTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGc                                                                                 <  1:660573/151‑1 (MQ=255)
                                                                        gggTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGC‑GAGTTTTACGTTGTTGCAGCGAAGCGGAAATCA‑TTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCcag                                                                              <  2:741701/151‑1 (MQ=255)
                                                                          gTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGC‑GGGGTTTTTGTTGTTTGCGCGAAGCGGAAAACA‑TTGGCAAACCTCCCCCGGCAGAGCCGGGGGAACGCCCGCGCGGGCcag                                                                              >  1:699594/1‑149 (MQ=255)
                                                                            aGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGG‑GGGGTTTTCGTTGTTGCAGCGAAACGGAAAACA‑TTGGGAATCCTCCCCTGGCAGAGGCGGGGTACCGCTGGcgcggg                                                                                  >  2:897010/1‑141 (MQ=255)
                                                                            aGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGC‑GGGGTTTTCGTTGTTGCAGCGAAGCGGGAAACA‑TTTGCAATCCTCCCATGGCAGAGCCGGGGTACCGCTGGCGCGGGCCAGCAgg                                                                          >  2:877748/1‑151 (MQ=255)
                                                                                     aGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGC‑GAGTTTTACGTTGTTGCAGCGAAGCGGAAATCA‑TTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGt                                                                 <  2:483708/151‑1 (MQ=255)
                                                                                            gACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGC‑GGGGTTTTCGTTGTTGGAGCGAAAAGGAAAAAA‑TTTGCAAACCCCCCCTGGCAGAGCCGGGGTACCCCCGGCGCGGGCCCGCAGGGGTTCCCGTc                                                                >  2:726204/1‑145 (MQ=255)
                                                                                            gACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGC‑GAGGGTTTTGGTGTTGCAGCGGAGCGGGAAAAA‑ATAGCAAAACCCCCCCCGCGGAGCCGGGGGACCCCCGGCGCGGGCCCGCAGCGCGTCCCGCCGTg                                                             >  1:672412/1‑148 (MQ=255)
                                                                                              ccATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGC‑GAGTTTTACGTTGTTGCAGCGAAGCGGAAATCA‑TTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGAGCGTGCCAGCAGGGCTTCCAGTCGTGCACGt                                                        <  2:920571/151‑1 (MQ=255)
                                                                                               cATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGC‑GAGTTTTACGTTTTTGCAGCGAAGCGGAAATCA‑TTGGCAATCCTCCCATGGCAGAGCCGGGATTCCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTa                                                       <  1:409894/151‑1 (MQ=255)
                                                                                                aTCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGC‑GGGGTTTTCTTTGGTTCCGCGAAGGCGAAAACA‑TTTGCCATTCCCCCATGGGAAAGGCCGGAGACCCCCCGCGCGGGCCCGCCAGGCCTcccatcgg                                                              >  1:907337/1‑137 (MQ=255)
                                                                                                  cgcgGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGC‑GGGGTTTTCGTTTTTGCAGCGAAGCGGGAATCA‑TTTGGAATCCTCCCCTGGCAGAGCCCGGGGACCGCCGGCGCGTGGCAGCAGGGCCTCCCGTCGTGCACGGaaa                                                     >  2:805033/1‑150 (MQ=255)
                                                                                                        aTTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGC‑GGGTTTTTCGTTGGTGCAGCGAAGCGGAAAACA‑TTGGCAATCCTCCCCAGGGAGAGACGGGAGAACGCTGGCGCGGGGCAACAGGGGTTCCACGCGTGCACGCAAAATGccg                                               >  2:161960/1‑148 (MQ=255)
                                                                                                           gacACCAGGCACAAACAACATGCTGCGGCGGGTGC‑GAGTTTTACGTTGTTGCAGCGAATCGGAAATCA‑TTGGCAATCCTCCAATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGATTCCAGTCGTGCACGTAAAATGCAGTCCa                                           <  1:433521/149‑1 (MQ=255)
                                                                                                            gCACCAGGCACAAACAACATGCTGCGGCGGGTGC‑GAGTTTTACGTTGTTGCAGCGAAGCGGAAATCA‑TTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCa                                           <  2:634241/150‑1 (MQ=255)
                                                                                                            gCACCAGGCACAAACAACATGCTGCGGCGGGTGC‑GAGTTTTACCTTGTTGCAGCGAAAAGGAAAACA‑ATGGGAATACTCCCATGGCAGAGCCGGGGGACCGCCGGCGCGGGCCCCCAGGGCCTCCCGGCGTGCACGTAAAAAGCAGt                                              >  1:172081/1‑147 (MQ=255)
                                                                                                                cAGGCACAAACAACATGCTGCGGCGGGGGG‑GGGTTTTTCGTTTTTGTCGGGGAGCGGGAAACA‑TTTGCAAACCTCCCCAGGGCGAGCCGGGGGACCCCCGGGGCGTGCCAGCAGGGCCTTCCGTCGTGGACGTAAAATGCAGGcc                                            >  2:257221/1‑145 (MQ=255)
                                                                                                                     acaaacaaCATGCTGCGGCGGGGGC‑GGGTTTTTTGTTTTTTCAGCGGAACGGGAAACA‑TTGGCAATCCTCCCCTGGCAGAGCCGGGGTACCCCTGGGGCGTGCCCCCCAGGCTTCCCGTCGCGCACGGAAAATGCAGTCCAGGGCGcc                                    >  2:40541/1‑148 (MQ=255)
                                                                                                                        aaggacaTGCCGCTGCGGGTGC‑GAGTTTTACGTTGTTGCAGCGAAGCGGAAATCA‑TTTGCAATCCTCCCATTGCAGAGCCGGGACACCGCTGGCGCGTGCCTTCAGGGCTTCCAGTCGTGCACGTAAAATGCACTCCAGTGCGGCTTTGt                               <  1:689983/146‑1 (MQ=255)
                                                                                                                              aTGCTGCGGCGGGTGC‑GAGTTTTACGTTGTTGCAGCGAAGCGGAAATCA‑TTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGAc                        <  1:536013/151‑1 (MQ=255)
                                                                                                                                gctgcGGCGGGTGC‑GAGTTTTACGTTGTTGCAGCGAAGCGGAAATCA‑TTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACa                       <  1:214351/150‑1 (MQ=255)
                                                                                                                                   gcggcgGGTGC‑GAGTTTTACGTTGTTGCAGCGAAGCGGAAATCA‑TTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCa                   <  2:414633/151‑1 (MQ=255)
                                                                                                                                   gcggcgGGTGC‑GAGTTGTACGTTGGTGCAGCGAAGCGGAAACCA‑TTGGCAATCCTCCCAGGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCCGTGCGCCTTTGTCATCGACAtt                     <  1:126022/149‑1 (MQ=255)
                                                                                                                                        gggTGC‑GAGTTTTACGTTGTTGCAGCGAAGCGGAAATCA‑TTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGt              <  2:962810/151‑1 (MQ=255)
                                                                                                                                            ttccaGTCTTACGTGGTTGCAGCAATCCGGAAAACA‑TTGGCAATCCTTTCATGGCCCAGCAGGGGTACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGTACGTAAAACGCAGTCTTGTGCCCCTTTGTCACCGACATTCGGCTg               <  1:625511/143‑1 (MQ=255)
                                                                                                                                               gAGTTTTACGTTGTTGCAGCGAAGCGGAAATCA‑TTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTAcgccg        >  1:868950/1‑151 (MQ=255)
                                                                                                                                               gAGTTTTACGTTGTTGCAGCGAAGCGGAAATCA‑TTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTAcgccg        >  2:384059/1‑151 (MQ=255)
                                                                                                                                                   tttACGTTGTTGCAGCGAAGCGGGAATCA‑TTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCgcgc     >  2:254314/1‑150 (MQ=255)
                                                                                                                                                     tACGTTGTTGCAGCGAAGCGGAAATCA‑TTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACg  >  2:326374/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                      
ACCATGCGGCGGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGC‑GAGTTTTACGTTGTTGCAGCGAAGCGGAAATCA‑TTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACG  >  NC_000913/650315‑650613

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