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

CAACGCTGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCGGGAAG  >  NC_000913/650672‑650940
                                                                                                                                         |                                                                                                                                    
cAACGCTGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCt                                                                                                                                    >  1:432032/1‑140 (MQ=255)
            tGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTg                                                                                                                        <  1:1069048/140‑1 (MQ=255)
               tttGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGCATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCcg                                                                                                                     <  2:400644/140‑1 (MQ=255)
                  gCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGt                                                                                                                  <  1:559231/140‑1 (MQ=255)
                   cAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGCGCCTGCAACGGCGGGCTGGTTTATTTTTCTTTTTCTACTGTTTCTCATGCGGGGGGCTTTTGTCGCACTGCCGCGTc                                                                                                                 >  2:1041566/1‑140 (MQ=255)
                           tatCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCggg                                                                                                         <  1:460692/140‑1 (MQ=255)
                                tGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCgggggggg                                                                                                    >  2:866450/1‑135 (MQ=255)
                                 gCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGTGCt                                                                                                   >  1:890117/1‑140 (MQ=255)
                                  cGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTc                                                                                                  >  1:24815/1‑140 (MQ=255)
                                       cGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTTACCTGTTTCTCCTGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGGTCAAGCa                                                                                            >  2:366116/1‑140 (MQ=255)
                                          gTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGTGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGt                                                                                         <  1:1041566/140‑1 (MQ=255)
                                             gggCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCt                                                                                      <  1:447827/140‑1 (MQ=255)
                                              ggCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTTATTTTTTACCTGTTTTTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGCTCAAGCAGGGTCTg                                                                                     >  1:876629/1‑140 (MQ=255)
                                               gCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGc                                                                                    <  2:1067622/140‑1 (MQ=255)
                                                   tGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGTGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAaat                                                                                >  1:1090478/1‑140 (MQ=255)
                                                          catcaCATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAg                                                                         <  1:847245/140‑1 (MQ=255)
                                                             cacaTCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGCTCCAGGAGGGTCTGCCAAAAATTCAGCgg                                                                      >  2:892785/1‑139 (MQ=255)
                                                             cacaTCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGGTTATTTTCATTTTTTACCTGTTTTTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGCTCAAGCAAGGTGTGCAAAAAATGCAGCgg                                                                      >  2:870159/1‑139 (MQ=255)
                                                               caTCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGGTTTTTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGa                                                                    >  1:1023944/1‑140 (MQ=255)
                                                                 tCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGGTCAGGCAGGTTCTGCAAAAAATGCAGCGGGAcc                                                                  >  2:802557/1‑139 (MQ=255)
                                                                  ccccAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTTTTTTTATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGGTCAAGCAGGGTTTTCAAAAAAAGCAGCGGGGCTg                                                                 >  2:567483/1‑140 (MQ=255)
                                                                    ccAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCa                                                               <  1:892785/140‑1 (MQ=255)
                                                                     cAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTTTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGGTCAAGCCGGGTCTGCCAAAAATTCAGCGGGACTGccg                                                              >  2:665995/1‑138 (MQ=255)
                                                                       gACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAggg                                                            >  1:584726/1‑140 (MQ=255)
                                                                               ggTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAgcgg                                                    <  1:567483/140‑1 (MQ=255)
                                                                               ggTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAgcgg                                                    <  2:24815/140‑1 (MQ=255)
                                                                                 tGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAgcggcg                                                  <  1:1088583/140‑1 (MQ=255)
                                                                                 tGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAgcggcg                                                  <  1:866450/140‑1 (MQ=255)
                                                                                   ccTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAGGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGcgcg                                                <  1:802557/140‑1 (MQ=255)
                                                                                   ccTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGcgcg                                                <  2:606590/140‑1 (MQ=255)
                                                                                   ccTGCAACGGCGGGCTGGGTTATTTTTATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGTGCTCAAGCAGGGTTTTCCAATAATGCAGCGCGGCTGCCGGGGCCAGCGGggcg                                                >  2:1005014/1‑140 (MQ=255)
                                                                                    cTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGa                                               <  2:584726/140‑1 (MQ=255)
                                                                                    cTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGa                                               <  1:870159/140‑1 (MQ=255)
                                                                                             cggGCTGGGTTATTTTCATTTTTTACCTGTTTCTCATTGGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGGTCAAACCAGGTTTGCAAAAAATGGAGCGTGACTCGAGGGGCCACCGGCGCGGGAGGCGTgg                                      >  2:756132/1‑139 (MQ=255)
                                                                                             cggGCTGGGTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGTGGTCCAACAGGTTCTGCAAATAATGCAGCGGGACTGCAGGGGCCAGCGGCGCGATAGCCgggg                                      >  2:37552/1‑137 (MQ=255)
                                                                                                          tttCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGc                         <  2:1090478/140‑1 (MQ=255)
                                                                                                          tttCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGcc                        >  2:76436/1‑140 (MQ=255)
                                                                                                           ttCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGcc                        <  2:884348/140‑1 (MQ=255)
                                                                                                              aTTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCcagc                     >  2:237012/1‑140 (MQ=255)
                                                                                                              aTTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCcagc                     <  2:1040610/140‑1 (MQ=255)
                                                                                                                      ccTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACg             <  2:432032/140‑1 (MQ=255)
                                                                                                                                 tGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCGGGAAg  <  1:1005014/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                    
CAACGCTGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCGGGAAG  >  NC_000913/650672‑650940

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: