Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913650,8090TG36.7% 36.4 / 13.3 30T350P (ACC→CCC) citCcitrate lyase synthetase
Reads supporting (aligned to +/- strand):  ref base T (6/13);  new base G (11/0);  total (17/13)
Fisher's exact test for biased strand distribution p-value = 3.17e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.57e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CTGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCT‑CAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCGGGAAG  >  NC_000913/650677‑650940
                                                                                                                                    |                                                                                                                                    
ctgctATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTTTTTTTATTTTTTACCTGTTTTTCATGCGGGGGGCttttta                                                                                                                               >  2:395938/1‑138 (MQ=255)
 tgctATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGAc                                                                                                                              >  1:357729/1‑140 (MQ=255)
           agacAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCcgc                                                                                                                    <  1:42966/137‑1 (MQ=255)
                   cGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTTTTTTCCTTTTTTACCTGTTTTTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGc                                                                                                            >  2:652956/1‑140 (MQ=255)
                     atatCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCgg                                                                                                          <  1:57692/140‑1 (MQ=255)
                            gCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGGTTATTTTTATTTTTTACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCgggctgtct                                                                                                   >  2:42966/1‑134 (MQ=255)
                                 tCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCT‑CAAGc                                                                                             <  2:357729/140‑1 (MQ=255)
                                    aGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGGT‑CAAGCAgg                                                                                          >  2:506263/1‑140 (MQ=255)
                                              tGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTTTTTTTATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTTACGAGCTGCCGCGTCCTGGCGGGGGGGCT‑CAAACAGGGTTTGCaaaa                                                                                >  2:173849/1‑139 (MQ=255)
                                               gCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGGTTATTTTCATTTTTTACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGGT‑CCAACAGGGTCTTCCaata                                                                               >  2:616463/1‑140 (MQ=255)
                                                   atcatcACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTTATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTTACGAGCTGCCGCGTCCTGGCGGGGGGGCT‑CAAACAGGGTTTGCAAATAAAGc                                                                           >  2:631683/1‑140 (MQ=255)
                                                        cacaTCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCT‑CAAGCAGGTTCTGCAAATAATGCAGCGt                                                                      <  1:395938/140‑1 (MQ=255)
                                                        cacaTCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCT‑CAAGCAGGTTCTGCAAATAATGCAGCGt                                                                      <  1:656724/140‑1 (MQ=255)
                                                          caTCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCT‑CAAGCAGGTTCTGCAAATAATGCAGCGTGa                                                                    <  1:173849/140‑1 (MQ=255)
                                                          agtCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCT‑CAAGCAGGTTCTGCAAATAATGCAGCGTGa                                                                    <  1:203714/138‑1 (MQ=255)
                                                                cAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGCT‑CAAGCAGGGTCTGCAAAAAAAGCAGCGGGGCTTCAg                                                              >  2:229568/1‑140 (MQ=255)
                                                                       aaGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTTTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGGTCCAACCGGGTTTTCCAAAAAAT‑CAGCGGGGCTGCCGGGGCCaa                                                       >  2:603848/1‑139 (MQ=255)
                                                                       aaGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCT‑CAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAg                                                       <  1:603848/140‑1 (MQ=255)
                                                                       aaGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGGT‑CAAGCAGGGTTTGCAAAAAAAGCAGCCTGGCTGCAGGGGcccg                                                       >  2:524398/1‑138 (MQ=255)
                                                                       aaGGGTGCCTGCAACGGCGGGCTGGGTTATTTTCCTTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGCGCTGCCGCGTCCTCGCGGGGGGGGT‑CAAGAAAGGTCTGGCAATAATGCGGCGCGACTGCAGGGGCCac                                                       >  2:538016/1‑139 (MQ=255)
                                                                             gCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCCTGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGGT‑CCAACAGGGTTTTCAAAAAAAGCAGCGGGACTGCCGGGGCCCGCGgcgc                                                 >  2:101973/1‑140 (MQ=255)
                                                                              ccTGCAACGGCGGGCTGGTTTATTTTCATTTTTTACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGTGCT‑CAAACAGGGTTTTCCAAAAAAGCAGGGTGGCTGCAGGGACCAGCGGcgcg                                                >  2:515831/1‑140 (MQ=255)
                                                                              ccTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCT‑CAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGcgcg                                                <  1:251736/140‑1 (MQ=255)
                                                                                                      ttCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCT‑CAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGcc                        <  2:214856/140‑1 (MQ=255)
                                                                                                              tCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCT‑CAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGt                  <  2:234855/138‑1 (MQ=255)
                                                                                                              tCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCT‑CAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGt                  >  1:234855/1‑138 (MQ=255)
                                                                                                                     tttCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCT‑CAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTAcc         <  1:524398/140‑1 (MQ=255)
                                                                                                                        ctcATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCT‑CAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCgg      <  1:210043/140‑1 (MQ=255)
                                                                                                                          cATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCT‑CAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCGGGa    >  2:15072/1‑140 (MQ=255)
                                                                                                                          cATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCT‑CAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCGGGa    >  1:38251/1‑140 (MQ=255)
                                                                                                                            tGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCT‑CAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCGGGAAg  >  1:193620/1‑140 (MQ=255)
                                                                                                                                    |                                                                                                                                    
CTGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCT‑CAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCGGGAAG  >  NC_000913/650677‑650940

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