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

TCAACGCTGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCGGGAAGCGG  >  NC_000913/650671‑650943
                                                                                                                                          |                                                                                                                                      
tcAACGCTGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCTATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTc                                                                                                                                       <  1:1170919/140‑1 (MQ=255)
tcAACGCTGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTc                                                                                                                                       <  1:81344/140‑1 (MQ=255)
 cAACGCTGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTTACCTGTTTCTCATGCGGGGGGCt                                                                                                                                      >  2:614068/1‑140 (MQ=255)
   aCGCTGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTTACCTGTTTTTCATGCGGGGGGCttt                                                                                                                                    >  2:81344/1‑140 (MQ=255)
   aCGCTGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCttt                                                                                                                                    <  1:704110/140‑1 (MQ=255)
     gctgctATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTg                                                                                                                                  <  2:1060003/140‑1 (MQ=255)
      ctgctATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGTGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTTACCTGTTTTTCCTGCGGGGGTCTTTTgt                                                                                                                                 >  2:949737/1‑139 (MQ=255)
      ctgctATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGa                                                                                                                                 >  1:70563/1‑140 (MQ=255)
      ctgctATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGGTTATTTTCATTTTTTACCCGTTTCTCATGCGGGGGGCtttttt                                                                                                                                 >  2:321483/1‑138 (MQ=255)
         ctATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGGTTATTTTCATTTTTTACCTGTTTCTCATGCGGGGGGCTTTTGACGa                                                                                                                              >  2:1049641/1‑140 (MQ=255)
          tATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAg                                                                                                                             <  2:1025585/140‑1 (MQ=255)
                 ttGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGGTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCcgc                                                                                                                      >  1:234308/1‑140 (MQ=255)
                   gCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGt                                                                                                                    <  1:59535/140‑1 (MQ=255)
                   gCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTTACGAGCTGCCGCGt                                                                                                                    >  2:354719/1‑140 (MQ=255)
                                 tGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTGTCATTTTTCACCTGTTTCTCATGCGGGGGGATTTTTACGAGCTGCCGCGTCCTGGCgggggggg                                                                                                      >  2:1008457/1‑135 (MQ=255)
                                          aGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGGTCAAACAgg                                                                                             >  2:4770/1‑140 (MQ=255)
                                                         atcatcACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGGTCTGCAAATAATGc                                                                              >  1:1178788/1‑140 (MQ=255)
                                                          tcatcaCATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGGTCAAGCAGGGTTTGCAaaaaattgc                                                                             >  2:745832/1‑137 (MQ=255)
                                                              cacaTCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGGTTTTTTTCATTTTTCACCTGTTTCTCCTGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGGTCAAACAGGGTCTTCAAAAAAAGCAGCgg                                                                         >  2:1189758/1‑139 (MQ=255)
                                                                caTCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGa                                                                       <  2:70563/140‑1 (MQ=255)
                                                                caTCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGa                                                                       <  1:401460/140‑1 (MQ=255)
                                                                  tCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACt                                                                     <  1:417734/140‑1 (MQ=255)
                                                                   ccccAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTTATTTTTCACCTGGTTCTCATGCGGGGGGCTTTTGTCGCGCTTCCGCGTCCTGGCGGGGGGGCGCAAGCACGGTTTTCAAAAAATGCCGCCGGACTg                                                                    >  2:336492/1‑140 (MQ=255)
                                                                      cAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGGGCTCAAGCAGGGTTTGCAAAAAAAGCAGCGTGACTGCAg                                                                 >  1:1148455/1‑140 (MQ=255)
                                                                             aaGGGTGCCTGCAACGGCGGGCTGGTTTTTTTTTTTTTTTTACCCGTTTCTCATGCGGGGGGCTTTTTTCGAGCTGCCGCGTCCTGGCGGGGGGGCCCCAGGAGGTTCTGTCAATAATACAGAGGGGCCGCAGGGGCCAg                                                          >  2:1263061/1‑140 (MQ=255)
                                                                             aaGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGCTCAAACAGGGTCTGCAAATAATGCAGCGTGACTGCAGGGGCCAg                                                          >  1:1249934/1‑140 (MQ=255)
                                                                             aaGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAg                                                          <  1:639551/140‑1 (MQ=255)
                                                                             aaGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGCTCAAGCAGGGTCTGCAAATAATGCAGCGTGACTGCAGGGGCCAg                                                          >  1:82171/1‑140 (MQ=255)
                                                                             aaGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGCTCAAACAGGGTCTGCCAATAAAGCAGCGTGACTGCAGGGAcccg                                                          >  2:417734/1‑138 (MQ=255)
                                                                             aaGGGTGCCTGCAACGGCGGGCTGGGTTATTTTCATTTTTTACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGCTCAAGCAGGGTTTTCAAATAAATCAGCGTGGCTGCCGGGGCCAg                                                          >  2:184777/1‑140 (MQ=255)
                                                                             aaGGGTGCCTGCAACGGCGGGCTGGGTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGTGGTCAAACAGGGTCTGCAAAAAATGCAGCGGGACTGCAGGGGCCAg                                                          >  2:807408/1‑140 (MQ=255)
                                                                                ggTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAgcgg                                                       <  1:1189758/140‑1 (MQ=255)
                                                                                ggTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAgcgg                                                       <  1:479031/140‑1 (MQ=255)
                                                                                    ccTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGcgcg                                                   <  1:745832/140‑1 (MQ=255)
                                                                                    ccTGCAACGGCGGGCTGGGTTTTTTTTATTTTTTACCTGTTTTTCATGCGGGGGGCTTTTTACGAGCTGCCGCGTCCTGGCGGGGGGGGTCCAACAGGGTTTGCCAAAAAATCAGCGGGGCTGCCGGGGGCAGCGGggcg                                                   >  2:622155/1‑140 (MQ=255)
                                                                                      tGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGAt                                                 <  1:1211520/140‑1 (MQ=255)
                                                                                                            ttCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGcc                           <  1:807408/140‑1 (MQ=255)
                                                                                                               aTTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGTTCGTTTTTCGCcagc                        <  1:1263061/140‑1 (MQ=255)
                                                                                                               aTTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCcagc                        >  1:1051152/1‑140 (MQ=255)
                                                                                                               aTTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCcagc                        >  1:93588/1‑140 (MQ=255)
                                                                                                                           tttCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTAcc            >  1:682671/1‑140 (MQ=255)
                                                                                                                             tctcATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCg          >  2:5551/1‑140 (MQ=255)
                                                                                                                                  tGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCGGGAAg     >  2:144853/1‑140 (MQ=255)
                                                                                                                                     gggggTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCGGGAAGCgg  >  2:390061/1‑140 (MQ=255)
                                                                                                                                          |                                                                                                                                      
TCAACGCTGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCGGGAAGCGG  >  NC_000913/650671‑650943

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: ‑