Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913650,8090TG34.8% 74.4 / 18.8 46T350P (ACC→CCC) citCcitrate lyase synthetase
Reads supporting (aligned to +/- strand):  ref base T (11/19);  new base G (16/0);  total (27/19)
Fisher's exact test for biased strand distribution p-value = 1.74e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCGGGAAGCGGATAT  >  NC_000913/650684‑650947
                                                                                                                             |                                                                                                                                          
tGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTACAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTg                                                                                                                              <  2:129128/140‑1 (MQ=255)
      gCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGt                                                                                                                        <  1:294579/140‑1 (MQ=255)
      gCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTTACGAGCTGCCGCGt                                                                                                                        >  1:380615/1‑140 (MQ=255)
       cAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTTATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTc                                                                                                                       >  2:799166/1‑140 (MQ=255)
       cAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTc                                                                                                                       >  2:1011381/1‑140 (MQ=255)
                    tGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGGTTATTTTCATTTTTCACCTGTTTTTCATGCGGGGGGCTTTTTACGAGCTGCCGCGTCCTGGCgggggggg                                                                                                          >  2:997398/1‑135 (MQ=255)
                    tGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGGTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTCGCgggggggg                                                                                                          >  2:541005/1‑135 (MQ=255)
                     gCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTTACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCggggggggt                                                                                                         >  2:879611/1‑134 (MQ=255)
                         aTCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTTACGAGCTGCCGCGTCCTGGCGGGGGGGGTCaaa                                                                                                     >  2:1006906/1‑139 (MQ=255)
                            gAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAg                                                                                                  <  1:657060/140‑1 (MQ=255)
                                ggggCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTc                                                                                              >  1:665340/1‑140 (MQ=255)
                                    cGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTTACCTGTTTCTCATGCGGGGGGCTTTTTACGAGCTGCCGCGTCCTGGCGGGGGGGGTCAAACCGGGTCTGCa                                                                                          >  2:1067572/1‑140 (MQ=255)
                                    cGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTTACGAGCTGCCGCGTCCTGGCGGGGGGGGTCCAGCAGGGTCTTCa                                                                                          >  2:657060/1‑140 (MQ=255)
                                      aTGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCCCCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCaaa                                                                                        <  2:673317/140‑1 (MQ=255)
                                      aTGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCaaa                                                                                        <  2:665340/140‑1 (MQ=255)
                                       tGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGTGCTCAAGCAGGGTTTTCaaaa                                                                                       >  2:611734/1‑139 (MQ=255)
                                            atcatcACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGc                                                                                  <  1:541005/140‑1 (MQ=255)
                                              catcaCATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAg                                                                                <  1:879611/140‑1 (MQ=255)
                                                   caTCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGa                                                                           <  1:1067572/140‑1 (MQ=255)
                                                   caTCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGa                                                                           <  2:1004962/140‑1 (MQ=255)
                                                     tCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTTACGAGCTGCCGCGTCCTGGCGGGGGGGCTCAAACAGGGTTTTCAAAAAAAGCAGCGTGGCt                                                                         >  2:874977/1‑140 (MQ=255)
                                                      ccccAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTTATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGGTCAAGCAGGGTCTTCAAATAATGCAGCGGGACTg                                                                        >  2:854007/1‑140 (MQ=255)
                                                      ccccAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGGGCTCAAGCAGGGTTTGCAAATAATGCAGCGGGACTg                                                                        >  2:232486/1‑140 (MQ=255)
                                                        ccAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCa                                                                      <  2:1002372/140‑1 (MQ=255)
                                                                aaGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGGTCCAGCAGGGTTTTCAAATAATGCAGCGGGGCTGCAGGGGCCAg                                                              >  2:1079080/1‑140 (MQ=255)
                                                                 aGGGTGCCTGCAACGGCGGGCTGGTTTATTTTTCTTTTTTACCTGTTTCTCATGCGGGGGGCTTTTTACCAGCTGCCGCGTCCTGGCGGGAGGGCTCCAACAGGTTCTGCAAAAAAATCAGCGTGACTGCAGGGACCaac                                                             >  2:84059/1‑138 (MQ=255)
                                                                 aGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAgc                                                             <  1:1011381/140‑1 (MQ=255)
                                                                 aGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAgc                                                             <  1:799166/140‑1 (MQ=255)
                                                                   ggTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAgcgg                                                           <  1:312422/140‑1 (MQ=255)
                                                                       ccTGCAACGGCGGGCCGGGTTATTTTCATTTTTCACCTGTTTCTTATGGGGGGGGCTTTTTACGAGCTGCCGCGTCCTGGCGGGGGTGCTCTAGCAGGGTCTTCAAATAATACAGCGCGGATGCAGGGACCAGTGGcgcg                                                       >  2:531900/1‑140 (MQ=255)
                                                                        cTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGCTCAAACAGGGTCTGCCAAAAAAGCAGCGTGACTGCAGGGGCCAGCGGcgcgg                                                      >  2:312422/1‑139 (MQ=255)
                                                                           cAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAg                                                   <  1:1006906/140‑1 (MQ=255)
                                                                                gcggGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAAAGCAGCGTGACTGCCGGGGCCAGCGGCGCGGTAGCCGTg                                              >  2:86934/1‑140 (MQ=255)
                                                                                   ggCTGGGTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGCTCAAGCAGGTTTTGCAAATAATGCAGCGGGGCTGCAGGGGCCAGCGGCGCGATAGCCGTgaga                                           >  2:348046/1‑140 (MQ=255)
                                                                                           tATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTTCTGGCGGGACTGCTCCAGCAGGGTCTCCAAAAAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGGTCGttttt                                   >  2:1154508/1‑140 (MQ=255)
                                                                                                tCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCca                              >  1:646557/1‑140 (MQ=255)
                                                                                                  aTTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCcagc                            >  2:149551/1‑140 (MQ=255)
                                                                                                      ttCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGAGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGtt                        >  1:481234/1‑140 (MQ=255)
                                                                                                          ccTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACg                    <  1:874977/140‑1 (MQ=255)
                                                                                                             gTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTAc                 >  1:202399/1‑140 (MQ=255)
                                                                                                              tttCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACg                    >  1:359144/1‑136 (MQ=255)
                                                                                                              tttCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACg                    <  2:359144/136‑1 (MQ=255)
                                                                                                                tctcATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCg              >  2:564239/1‑140 (MQ=255)
                                                                                                                     tGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCGGGAAg         <  1:348046/140‑1 (MQ=255)
                                                                                                                     tGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCGGGAAg         <  1:191843/140‑1 (MQ=255)
                                                                                                                     tGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCGGGAAg         <  2:202399/140‑1 (MQ=255)
                                                                                                                       cGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCGGGAAGCg       >  2:267835/1‑140 (MQ=255)
                                                                                                                         ggggTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCGGGAAGCGGa     >  2:434931/1‑140 (MQ=255)
                                                                                                                            gTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCGGGAAGCGGatat  >  2:435203/1‑140 (MQ=255)
                                                                                                                             |                                                                                                                                          
TGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCGGGAAGCGGATAT  >  NC_000913/650684‑650947

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑