Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913650,4640AT44.4% 48.6 / 15.6 36R9R (CGT→CGA) citEcitrate lyase, citryl‑ACP lyase (beta) subunit
Reads supporting (aligned to +/- strand):  ref base A (3/17);  new base T (15/1);  total (18/18)
Fisher's exact test for biased strand distribution p-value = 4.06e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGT‑TTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGC  >  NC_000913/650330‑650610
                                                                                                                                       |                                                                                                                                                  
gtgtCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGGGT‑TTTTCGTTGTTTCAGCg                                                                                                                                       >  1:58668/1‑148 (MQ=255)
gtgtCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGT‑TTTACGTTGTTGCAGCGaa                                                                                                                                     <  2:114727/150‑1 (MQ=255)
  gtCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGT‑TTTACGTTGTTGCAGCGAAGc                                                                                                                                   <  1:724642/150‑1 (MQ=255)
  gtCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGGGGGT‑TTTTCGTTTTTGGAGCGAAGCg                                                                                                                                  >  2:669170/1‑151 (MQ=255)
  gtCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGCGGGTGTTTTCGTT‑CTGAAGCGAAg                                                                                                                                    >  2:549635/1‑149 (MQ=255)
      ttttCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGT‑TTTACGTTGTTGCAGCGAAGCGGaa                                                                                                                               <  2:715581/150‑1 (MQ=255)
          cACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGGGAGG‑TTTTTGGTGGTGCAGGGAAGCGGaaagtat                                                                                                                          >  1:734778/1‑147 (MQ=255)
                aTGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGT‑TTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCaa                                                                                                                    <  1:22813/151‑1 (MQ=255)
                     aCGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGCGGGG‑TTTTCGTTGTTTCAGCGAAGCGGGAAAAAATTGCAAtcc                                                                                                                 >  1:542067/1‑149 (MQ=255)
                        gAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGGGGGG‑TTTTCGTTGTTGCAGCGAAGCGGGAAATATTGGCAATCCTccc                                                                                                             >  1:437685/1‑150 (MQ=255)
                         aGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGGGCGG‑TTTTATTTGTTTCAGCGCAGAGGGAATCAATGGCaagacttc                                                                                                              >  1:209952/1‑142 (MQ=255)
                          gTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGT‑TTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATg                                                                                                          <  2:295018/151‑1 (MQ=255)
                             ttCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGGGAGT‑TTTTCGCTTTTGTAGAGCAACGGGAAACATTTGCCATCCTccc                                                                                                             >  2:744281/1‑145 (MQ=255)
                                       aaCATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGT‑TTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCgg                                                                                                <  1:56090/148‑1 (MQ=255)
                                        aCATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGGGGGT‑TTTTCGTTGTTGCAGCGGAACGGGAATCAATGGGAAAcccccccagggcgaggcggg                                                                                               >  1:690397/1‑134 (MQ=255)
                                             aGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGGGGGT‑TTTTTGTTTTTGCAGCGCAAAGGAAAACATTTGCAAACCCCCCCTCGCCAGGCCGGGGTAccc                                                                                         >  2:215527/1‑148 (MQ=255)
                                                         gggTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGT‑TTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCcag                                                                           <  2:188182/151‑1 (MQ=255)
                                                              gATGAAGGTGTTGCTGACCAGCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGT‑TTTTCGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCAGGATACCGCTGGCGCGTGCcagcatg                                                                       <  2:505622/150‑3 (MQ=255)
                                                               aTGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGAGGCGGGG‑TTTTAGTTGTTTCTGCGGAAGGGAAAACAATGGCAAACCCCCCCAGGCAGAGCCGGGAAACCCCCCGCGCGTGCcaccccgg                                                                      >  1:89748/1‑144 (MQ=255)
                                                                  aaGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGGGT‑TTTTCGTTGTTTCAGCGAAGCGGAAATCATTTGCAATACTCCCCCGGCAGAGCCCGGATACCGCTGGCGCGTGCCCGCAGGGcct                                                                   >  1:345396/1‑148 (MQ=255)
                                                                    ggAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGT‑TTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCa                                                                <  2:168745/151‑1 (MQ=255)
                                                                             gACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGT‑TTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCAc                                                       <  1:46485/151‑1 (MQ=255)
                                                                                  tCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGGGG‑TTTTTGTTTTTGCCGCGAAGCGGGAAACAATTGCAAACCCCCCCTGGCAGAGCCGGGATACCGCCGGCGCGGGCCAGCAGGGGCTCCAGTCGTGCCCGTaa                                                   >  2:497897/1‑150 (MQ=255)
                                                                                   cgcgGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGCGGGG‑TTTTCGTTGTTGCAGCGGAGCGGGAATCATTGGCAATCCTCCCCTGGCCGAGCCGGGGAACCGCTGGCGCGGGGCCGCAGGGCTTCCAGGCGTGCACGTaaa                                                  >  1:524569/1‑150 (MQ=255)
                                                                                    gcgGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGT‑TTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAAt                                                <  2:78449/151‑1 (MQ=255)
                                                                                        cATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGGGGGG‑TTTTTGGTGTTGGCACGGCACGGAAAACATTGGGAATCCCTCCCCGGCCGAGGCGCGATTACCCTGGTGGGGGCCCGCAGGGCGTCCAGTCGTGGACgggaa                                                  >  2:156288/1‑141 (MQ=255)
                                                                                                ccAGGCACAAACAACATGCTGCGGCGGGTGCGAGT‑TTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTgc                                    <  1:105722/151‑1 (MQ=255)
                                                                                                      acaaacaaCATGCTGCGGCGGGTGCGAGT‑TTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCttt                              <  1:609556/151‑1 (MQ=255)
                                                                                                            aacaTGCTGCGGCGGGTGCGAGT‑TTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATc                        <  1:454492/151‑1 (MQ=255)
                                                                                                              caTGCTGCGGCGGGTGCGAGT‑TTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGa                      <  1:435151/151‑1 (MQ=255)
                                                                                                                      ggcggGTGCGAGT‑TTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGc              <  1:324659/151‑1 (MQ=255)
                                                                                                                        cggGTGCGAGT‑TTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTg            >  2:212035/1‑151 (MQ=255)
                                                                                                                           gTGCGAGT‑TTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTAc         >  1:75259/1‑151 (MQ=255)
                                                                                                                             gCGAGT‑TTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTAcg        <  1:111324/150‑1 (MQ=255)
                                                                                                                                 gT‑TTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCcgcg   <  1:431351/151‑1 (MQ=255)
                                                                                                                                  t‑tttACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCgcgc  >  2:285400/1‑151 (MQ=255)
                                                                                                                                       |                                                                                                                                                  
GTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGT‑TTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGC  >  NC_000913/650330‑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: ‑