Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913650,8090TG30.9% 55.9 / 18.1 42T350P (ACC→CCC) citCcitrate lyase synthetase
Reads supporting (aligned to +/- strand):  ref base T (10/19);  new base G (13/0);  total (23/19)
Fisher's exact test for biased strand distribution p-value = 4.59e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CGCTGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCGGGAAGCGGA  >  NC_000913/650675‑650944
                                                                                                                                      |                                                                                                                                        
cGCTGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCtttt                                                                                                                                     >  1:1162327/1‑140 (MQ=255)
   tgctATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTTACCTGTTTTTCCTGCGGGGGGCTTTTgtc                                                                                                                                  >  2:405673/1‑138 (MQ=255)
      tATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAg                                                                                                                               <  1:552342/140‑1 (MQ=255)
             ttGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTTACCTGTTTTTCTTGCGGGGGGGTTTTTACGAGCTGCcgc                                                                                                                        >  2:552342/1‑140 (MQ=255)
             ttGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCcgc                                                                                                                        <  2:1195754/140‑1 (MQ=255)
                     cGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGc                                                                                                                >  1:531147/1‑140 (MQ=255)
                       atatCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTTACGAGCTGCCGCGTCCTGGCgg                                                                                                              >  2:1197226/1‑140 (MQ=255)
                            cTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTg                                                                                                         <  2:193882/140‑1 (MQ=255)
                              gCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGGTTATTTTTATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTTACGAGCTGCCGCGTCCTGGCggggggggt                                                                                                       >  2:1185920/1‑134 (MQ=255)
                               cGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTTATTTTTTACCTGTTTCTCCTGCGGGGGGCTTTTGACGAGCTTCCGCGTCCTGGCggggggggtc                                                                                                      >  2:614944/1‑133 (MQ=255)
                               cGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTTACGAGCTGCCGCGTCCTGGCGGGGGTGGTc                                                                                                      >  2:819129/1‑140 (MQ=255)
                               cGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGGTTATTTTTATTTTTTACCTGGTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGCTc                                                                                                      >  2:1235240/1‑140 (MQ=255)
                                     gAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAg                                                                                               <  1:614944/140‑1 (MQ=255)
                                      aGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTTACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGGCCTGGCGGGGGGGCTCAAGCAgg                                                                                              >  2:337972/1‑140 (MQ=255)
                                         ggggCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTc                                                                                           <  1:337972/140‑1 (MQ=255)
                                          gggCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCt                                                                                          <  1:980851/140‑1 (MQ=255)
                                              gATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCaa                                                                                      <  2:1162327/140‑1 (MQ=255)
                                                       catcaCATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAg                                                                             <  1:1235240/140‑1 (MQ=255)
                                                          cacaTCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGt                                                                          <  1:405673/140‑1 (MQ=255)
                                                          cacaTCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGt                                                                          <  2:531147/140‑1 (MQ=255)
                                                          cacaTCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCCGTTTCTCCTGCGGGGGGCTTTTTACGAGCTGCCGCGTCCTGGCGGGGGGGCTCCAACAGGGTCTGCCAAAAAATCAGCgg                                                                          >  2:1177343/1‑139 (MQ=255)
                                                            caTCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAACAGGTTCTGCAAATAATGCAGCGGGa                                                                        >  1:19656/1‑140 (MQ=255)
                                                            caTCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGCTCAAACAGGGTCTGCCAAAAAAGCAGCGTgg                                                                        >  1:1234696/1‑139 (MQ=255)
                                                               ccccAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTg                                                                     <  1:1185920/140‑1 (MQ=255)
                                                               ccccAGACTCAAGGGTGCCTGCAACGGCGGGCCGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGGTCAAGCAGGGTCTGCAAATAATACAGCGGGGCTg                                                                     >  2:69325/1‑140 (MQ=255)
                                                                    gACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGGTCAAGCAGGGTCTTCAAATAATGCAGCGGGGCTGCAggg                                                                >  2:991639/1‑140 (MQ=255)
                                                                    gACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGGGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAggg                                                                >  1:718187/1‑140 (MQ=255)
                                                                         aaGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCCGGCGGGGGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGCGGCTGCAGGGGcccg                                                           >  2:1065096/1‑138 (MQ=255)
                                                                          aGGGTGCCTGCAACGGCGGGCTGGTTTTTTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTTACGAGCTGCCGCGTCCTGGCGGGGGGGGTCAAGCAGGGTCTGCAAATAATGCAGCGGGACTGCAGGGGCCCgc                                                          >  2:975708/1‑140 (MQ=255)
                                                                          aGGGTGCCTGCAACGGCGGGCTGGTTTTTTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGCTCAAACAGGGTTTGCAAATAAATCAGCGTGGCTGCCGGGGCCCgc                                                          >  2:361352/1‑140 (MQ=255)
                                                                          aGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGCTCAAGCAGGGTCTGCAAATAATGCAGCGTGACTGCAGGGACCAgc                                                          >  1:1031513/1‑140 (MQ=255)
                                                                              tGCCTGCAACGGCGGGCTGGTTTATTTTTATTTTTCACCTGTTTTTCACGAGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTTCAGCAAGGTTTTTGAATAATGCCGCGTGTATGCAGGGGACCgcggcg                                                      >  2:1225644/1‑140 (MQ=255)
                                                                               gCCTGCAACGGCGGGCTGGTTTTTTTTTATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGCTCCAGGAGGGTTTTCCAAAAATTGCGCGGGGCTGCCGGGGGCCGCGgcgc                                                     >  2:1262159/1‑140 (MQ=255)
                                                                                ccTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGcgcg                                                    <  1:991639/140‑1 (MQ=255)
                                                                                ccTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGcgcg                                                    <  1:975708/140‑1 (MQ=255)
                                                                                ccTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGcgcg                                                    <  1:739891/140‑1 (MQ=255)
                                                                                 cTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGa                                                   <  1:1209888/140‑1 (MQ=255)
                                                                                     aaCGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGc                                               <  1:894452/140‑1 (MQ=255)
                                                                                          cggGCTTGTTTATTTTCATTTTTTCCCTGTTTTTCCTGCGGGGGGCTTTTGCCGGGAGGCCGCGTCCTGGCGGGAGTGCTTAAGAAGGGTCTGCAAAAAATGCAGTGTTTCTGCACGGGCCCGGGGGGCGCTAGCCCTga                                          >  2:480450/1‑140 (MQ=255)
                                                                                          cggGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGCTCAAACAGGGTCTGCAAATAATGCAGCGGGGCTGCCGGGGCCCACGGCGCGAGAACCGTgg                                          >  2:894452/1‑139 (MQ=255)
                                                                                          cggGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGCTCAAACCGGGTCTGCAAATAATGCAGCGTGACTGCAGGGGCCAGCGGCGCGGTAGCCgggg                                          >  2:1209888/1‑137 (MQ=255)
                                                                                                   ttATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGtttt                                 >  1:1001039/1‑140 (MQ=255)
                                                                                                        ttCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGcc                            <  1:1229985/140‑1 (MQ=255)
                                                                                                        ttCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGcc                            <  2:159008/140‑1 (MQ=255)
                                                                                                           aTTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCcagc                         >  1:792120/1‑140 (MQ=255)
                                                                                                             ttttCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCcagcag                       <  1:69325/140‑1 (MQ=255)
                                                                                                                tCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTg                    >  2:83047/1‑140 (MQ=255)
                                                                                                                       tttCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTAcc             >  2:1137434/1‑140 (MQ=255)
                                                                                                                                  ggggTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCGGGAAGCGGa  >  1:925048/1‑140 (MQ=255)
                                                                                                                                      |                                                                                                                                        
CGCTGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCGGGAAGCGGA  >  NC_000913/650675‑650944

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: