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→CGAcitEcitrate 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

GTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGC  >  NC_000913/650330‑650610
                                                                                                                                       |                                                                                                                                                  
gtgtCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGGGTTTTTCGTTGTTTCAGCg                                                                                                                                       >  1:58668/1‑148 (MQ=255)
gtgtCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGaa                                                                                                                                     <  2:114727/150‑1 (MQ=255)
  gtCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGc                                                                                                                                   <  1:724642/150‑1 (MQ=255)
  gtCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGGGGGTTTTTCGTTTTTGGAGCGAAGCg                                                                                                                                  >  2:669170/1‑151 (MQ=255)
  gtCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGCGGGTGTTTTCGTTCTGAAGCGAAg                                                                                                                                    >  2:549635/1‑149 (MQ=255)
      ttttCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGaa                                                                                                                               <  2:715581/150‑1 (MQ=255)
          cACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGGGAGGTTTTTGGTGGTGCAGGGAAGCGGaaagtat                                                                                                                          >  1:734778/1‑147 (MQ=255)
                aTGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCaa                                                                                                                    <  1:22813/151‑1 (MQ=255)
                     aCGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGCGGGGTTTTCGTTGTTTCAGCGAAGCGGGAAAAAATTGCAAtcc                                                                                                                 >  1:542067/1‑149 (MQ=255)
                        gAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGGGGGGTTTTCGTTGTTGCAGCGAAGCGGGAAATATTGGCAATCCTccc                                                                                                             >  1:437685/1‑150 (MQ=255)
                         aGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGGGCGGTTTTATTTGTTTCAGCGCAGAGGGAATCAATGGCaagacttc                                                                                                              >  1:209952/1‑142 (MQ=255)
                          gTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATg                                                                                                          <  2:295018/151‑1 (MQ=255)
                             ttCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGGGAGTTTTTCGCTTTTGTAGAGCAACGGGAAACATTTGCCATCCTccc                                                                                                             >  2:744281/1‑145 (MQ=255)
                                       aaCATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCgg                                                                                                <  1:56090/148‑1 (MQ=255)
                                        aCATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGGGGGTTTTTCGTTGTTGCAGCGGAACGGGAATCAATGGGAAAcccccccagggcgaggcggg                                                                                               >  1:690397/1‑134 (MQ=255)
                                             aGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGGGGGTTTTTTGTTTTTGCAGCGCAAAGGAAAACATTTGCAAACCCCCCCTCGCCAGGCCGGGGTAccc                                                                                         >  2:215527/1‑148 (MQ=255)
                                                         gggTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCcag                                                                           <  2:188182/151‑1 (MQ=255)
                                                              gATGAAGGTGTTGCTGACCAGCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTTCGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCAGGATACCGCTGGCGCGTGCcagcatg                                                                       <  2:505622/150‑3 (MQ=255)
                                                               aTGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGAGGCGGGGTTTTAGTTGTTTCTGCGGAAGGGAAAACAATGGCAAACCCCCCCAGGCAGAGCCGGGAAACCCCCCGCGCGTGCcaccccgg                                                                      >  1:89748/1‑144 (MQ=255)
                                                                  aaGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGGGTTTTTCGTTGTTTCAGCGAAGCGGAAATCATTTGCAATACTCCCCCGGCAGAGCCCGGATACCGCTGGCGCGTGCCCGCAGGGcct                                                                   >  1:345396/1‑148 (MQ=255)
                                                                    ggAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCa                                                                <  2:168745/151‑1 (MQ=255)
                                                                             gACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCAc                                                       <  1:46485/151‑1 (MQ=255)
                                                                                  tCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGGGGTTTTTGTTTTTGCCGCGAAGCGGGAAACAATTGCAAACCCCCCCTGGCAGAGCCGGGATACCGCCGGCGCGGGCCAGCAGGGGCTCCAGTCGTGCCCGTaa                                                   >  2:497897/1‑150 (MQ=255)
                                                                                   cgcgGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGCGGGGTTTTCGTTGTTGCAGCGGAGCGGGAATCATTGGCAATCCTCCCCTGGCCGAGCCGGGGAACCGCTGGCGCGGGGCCGCAGGGCTTCCAGGCGTGCACGTaaa                                                  >  1:524569/1‑150 (MQ=255)
                                                                                    gcgGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAAt                                                <  2:78449/151‑1 (MQ=255)
                                                                                        cATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGGGGGGTTTTTGGTGTTGGCACGGCACGGAAAACATTGGGAATCCCTCCCCGGCCGAGGCGCGATTACCCTGGTGGGGGCCCGCAGGGCGTCCAGTCGTGGACgggaa                                                  >  2:156288/1‑141 (MQ=255)
                                                                                                ccAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTgc                                    <  1:105722/151‑1 (MQ=255)
                                                                                                      acaaacaaCATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCttt                              <  1:609556/151‑1 (MQ=255)
                                                                                                            aacaTGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATc                        <  1:454492/151‑1 (MQ=255)
                                                                                                              caTGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGa                      <  1:435151/151‑1 (MQ=255)
                                                                                                                      ggcggGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGc              <  1:324659/151‑1 (MQ=255)
                                                                                                                        cggGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTg            >  2:212035/1‑151 (MQ=255)
                                                                                                                           gTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTAc         >  1:75259/1‑151 (MQ=255)
                                                                                                                             gCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTAcg        <  1:111324/150‑1 (MQ=255)
                                                                                                                                 gTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCcgcg   <  1:431351/151‑1 (MQ=255)
                                                                                                                                  ttttACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCgcgc  >  2:285400/1‑151 (MQ=255)
                                                                                                                                       |                                                                                                                                                  
GTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGC  >  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: