Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913753,6910GT41.4% 21.8 / 32.1 29L260M (CTG→ATG) gltAcitrate synthase
Reads supporting (aligned to +/- strand):  ref base G (9/8);  new base T (9/3);  total (18/11)
Fisher's exact test for biased strand distribution p-value = 2.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ACGCGGTGACCGAAGCCCATCAGGCGGAAAGAATCATTTTTGTCTTTCGCACGACGAACAAATTCCGGAATGTGTTTAACGGAGCTGATTTCTTCCAGCATTTTCAGCGCCGCTTCGTTAGCACCGCCGTGCGCAGGTCCCCACAGTGAAGCAATACCTGCTGCGATACAGGCAAACGGGTTCGCACCCGAAGAGCCAGCGGTACGCACGGTGGAGGTAGAGGCGTTCTGTTCATGGTCAGCGTGCAGGATCAGAATACGGTCCATAGCACGTTCCAGAATCGGATT  >  NC_000913/753546‑753832
                                                                                                                                                 |                                                                                                                                             
acGCGGTGACCGAAGCCCATCAGGCGGAAAGAATCATTTTTGTCTTTCGCACGACGAACAAATTCCGGAATGTGTTTAACGGAGCTGATTTCTTCCAGCATTTTCAGCGCCGCTTCGTTAGCACCGCCGTGCGCAGGTCCCCACATTg                                                                                                                                             >  1:167903/1‑148 (MQ=255)
                 cATCAGGCGGAAAGAATCATTTTTGTCTTTCGCACGACGAACAAATTCCGGAATGTGTTTAACGGAGCTGATTTCTTCCAGCATTTTCAGCGCCGCTTCGTTAGCACCGCCGTGCGCAGGTCCCCACATTGAAGCAATACCTGCTGCGa                                                                                                                           <  2:167903/149‑1 (MQ=255)
                     aGGCGGAAAGAATCATTTTTGTCTTTCGCACGACGAACAAATTCCGGAATGTGTTTAACGGAGCTGATTTCTTCCAGCATTTTCAGCGCCGCTTCGTTAGCACCGCCGTGCGCAGGTCCCCACAGTGAAGCAATACCTGCTGCGATa                                                                                                                         >  5:30992/1‑147 (MQ=255)
                       gCGGAAAGAATCATTTTTGTCTTTCGCACGACGAACAAATTCCGGAATGTGTTTAACGGAGCTGATTTCTTCCAGCATTTTCAGCGCCGCTTCGTTAGCACCGCCGTGCGCAGGTCCCCACAGTGAAGCAATACCTGCTGCGATACAgg                                                                                                                     <  6:30992/149‑1 (MQ=255)
                                aTCATTTTTGTCTTTCGCACGACGAACAAATTCCGGAATGTGTTTAACGGAGCTGATTTCTTCCAGCATTTTCAGCGCCGCTTCGTTAGCACCGCCGTGCGCAGGTCCCCACATTGAAGCAATACCTGCTGCGATACAGGCAAACGGGt                                                                                                            >  4:27898/1‑149 (MQ=255)
                                aTCATTTTTGTCTTTCGCACGACGAACAAATTCCGGAATGTGTTTAACGGAGCTGATTTCTTCCAGCATTTTCAGCGCCGCTTCGTTAGCACCGCCGTGCGCAGGTCCCCACAGTGAAGCAATACCTGCTGCGATACAGGCAAACGGGt                                                                                                            <  8:57826/149‑1 (MQ=255)
                                       ttGTCTTTCGCACGACGAACAAATTCCGGAATGTGTTTAACGGAGCTGATTTCTTCCAGCATTTTCAGCGCCGCTTCGTTAGCACCGCCGTGCGCAGGTCCCCACATTGAAGCAATACCTGCTGCGATACAGGCAAACGGGTTCGCAcc                                                                                                     >  7:13785/1‑149 (MQ=255)
                                          tCTTTCGCACGACGAACAAATTCCGGAATGTGTTTAACGGAGCTGATTTCTTCCAGCATTTTCAGCGCCGCTTCGTTAGCACCGCCGTGCGCAGGTCCCCACAGTGAAGCAATACCTGCTGCGATACAGGCAAACGGGTTCGCACCCg                                                                                                   >  1:119680/1‑148 (MQ=255)
                                            tttCGCACGACGAACAAATTCCGGAATGTGTTTAACGGAGCTGATTTCTTCCAGCATTTTCAGCGCCGCTTCGTTAGCACCGCCGTGCGCAGGTCCCCACAGTGAAGCAATACCTGCTGCGATACAGGCAAACGGGTTCGCACCCg                                                                                                   >  6:112303/1‑146 (MQ=255)
                                             ttCGCACGACGAACAAATTCCGGAATGTGTTTAACGGAGCTGATTTCTTCCAGCATTTTCAGCGCCGCTTCGTTAGCACCGCCGTGCGCAGGTCCCCACAGTGAAGCAATACCTGCTGCGATACAGGAAAACGGGTTCGCACCCGaa                                                                                                 >  1:151270/1‑147 (MQ=255)
                                                    gacgaACAAATTCCGGAATGTGTTTAACGGAGCTGATTTCTTCCAGCATTTTCAGCGCCGCTTCGTTAGCACCGCCGTGCGCAGGTCCCCACATTGCAGCAATACCTGCTGCGATACAGGCAAACGGGTTCGCACCCGAAGAGCCAGCg                                                                                        >  3:82237/1‑149 (MQ=255)
                                                                ccGGAATGTGTTTAACGGAGCTGATTTCTTCCAGCATTTTCAGCGCCGCTTCGTTAGCACCGCCGTGCGCAGGTCCCCACATTGAAGCAATACCTGCTGCGATACAGGCAAACGGGTTCGCACCCGAAGAGCCAGCGGTACGCACGGt                                                                             <  8:13785/148‑1 (MQ=255)
                                                                   gAATGTGTTTAACGGAGCTGATTTCTTCCAGCATTTTCAGCGCCGCTTCGTTAGCACCGCCGTGCGCAGGTCCCCACAGTGAAGCAATACCTGCTGCGATACAGGCAAACGGGTTCGCACCCGAAGAGCCAGCGGTACGCACGGTGGAg                                                                         <  5:112303/149‑1 (MQ=255)
                                                                                  aGCTGATTTCTTCCAGCATTTTCAGCGCCGCTTCGTTAGCACCGCCGGGCGCAGGTCCCCACAGTGAAGCAATACCTGCTGCGATACAGGCAAACGGGTTCGCACCCGAAGAGCCAGCGGTACGCACGGTGGAGGTAGAGGCGTTCTGt                                                          <  8:97680/149‑1 (MQ=255)
                                                                                      gATTTCTTCCAGCATTTTCAGCGCCGCTTCGTTAGCACCGCCGTGCGCAGGTCCCCACAGTGAAGCAATACCTGCTGCGATACAGGCAAACGGGTTCGCACCCGAAGAGCCAGCGGTACGCACGGTGGAGGTAGAGGCGTTCTGTTCAt                                                      >  1:150093/1‑149 (MQ=255)
                                                                                               cAGCATTTTCAGCGCCGCTTCGTTAGCACCGCCGTGCGCAGGTCCCCACATTGAAGCAATACCTGCTGCGATACAGGCAAACGGGTTCGCACCCGAAGAGCCAGCGGTACGCACGGTGGAGGTAGAGGCGGTCTGTTCATGGTCAGCGt                                             >  7:109532/1‑149 (MQ=255)
                                                                                                 gCATTTTCAGCGCCGCTTCGTTAGCACCGCCGTGCGCAGGTCCCCACAGTGAAGCAATACCTGCTGCGATACAGGCAAACGGGTTCGCACCCGAAGAGCCAGCGGTACGCACGGTgg                                                                           >  8:71709/1‑117 (MQ=255)
                                                                                                 gCATTTTCAGCGCCGCTTCGTTAGCACCGCCGTGCGCAGGTCCCCACAGTGAAGCAATACCTGCTGCGATACAGGCAAACGGGTTCGCACCCGAAGAGCCAGCGGTACGCACGGTgg                                                                           <  7:71709/117‑1 (MQ=255)
                                                                                                 gCATTTTCAGCGCCGCTTCGTTAGCACCGCCGTGCGCAGGTCCCCACAGTGAAGCAATACCTGCTGCGATACAGGCAAACGGGTTCGCACCCGAAGAGCCAGCGGTACGCACGGTGGAGGTAGAGGCGTTCTGTTCATGGTCAGCGTg                                            >  8:79736/1‑148 (MQ=255)
                                                                                                      ttCAGCGCCGCTTCGTTAGCACCGCCGTGCGCAGGTCCCCACAGTGAAGCAATACCTGCTGCGATACAGGCAAACGGGTTCGCACCCGAAGAGCCAGCGGTACTCACGGTGGAGGTAGAGGCGTTCTGTTCATGGTCAGCGTGCAGGa                                       <  2:119680/148‑1 (MQ=255)
                                                                                                                cTTCGTTAGCACCGCCGTGCGCAGGTCCCCACATTGAAGCAATACCTGCTGCGATACAGGCAAACGGGTTCGCACCCGAAGAGCCAGCGGTACGCACGGTGGAGGTAGAGGCGTTCTGTTCATGGTCAGCGTGCAGGATCAGAATACgg                            >  6:120968/1‑149 (MQ=255)
                                                                                                                   cGTTAGCACCGCCGTGCGCAGGTCCCCACAGTGAAGCAATACCTGCTGCGATACAGGCAAACGGGTTCGCACCCGAAGAGCCAGCGGTACGCACGGTGGAGGTAGAGGCGTTCTGTTCATGGTCAGCGTGCAGGATCAGAATACGGTc                          <  2:151270/148‑1 (MQ=255)
                                                                                                                     ttAGCACCGCCGTGCGCAGGTCCCCACATTGAAGCAATACCTGCTGCGATACAGGCAAACGTGTTCGCACCCGAAGAGCCAGCGGTACGCACGGTGGAGGTAGAGGCGTTCTGTTCATGGTCAGCGTGCAGGATCAGAATACGGTCCa                        >  1:88056/1‑148 (MQ=255)
                                                                                                                        gCACCGCCGTGCGCAGGTCCCCACATTGAAGCAATACCTGCTGCGATACAGGCAAACGGGTTCGCACCCGAAGAGCCAGCGGTACGCACGGTGGAGGTAGAGGCGTTCTGTTCATGGTCAGCGTGCAGGATCAGAATACGGTCCATAg                     <  8:109532/148‑1 (MQ=255)
                                                                                                                         cACCGCCGTGCGCAGGTCCCCACAGTGACGCACTACCTGCTGCGATAAAGGCGACCGGGTTCTCACCCGAAGAGCCCGCGGTACGCACGGTGGTGGTAGCGGCGTTCTGTTCATGGTCAGCGTGCAGGATCAGAATACGGTCCATCGCa                   >  4:92949/1‑149 (MQ=255)
                                                                                                                                    gcAGGTCCCCACAGTGAAGCAATACCTGCTGCGATACAGGCAAACGGGTTCGCACCCGAAGAGCCAGCGGTACGCACGGTGGAGGTAGAGGCGTTCTGTTCATGGTCAGCGTGCAGGATCAGAATACGGTCCATAGCACGTTc              >  2:126312/1‑143 (MQ=255)
                                                                                                                                    gcAGGTCCCCACAGTGAAGCAATACCTGCTGCGATACAGGCAAACGGGTTCGCACCCGAAGAGCCAGCGGTACGCACGGTGGAGGTAGAGGCGTTCTGTTCATGGTCAGCGTGCAGGATCAGAATACGGTCCATAGCACGTTc              <  1:126312/143‑1 (MQ=255)
                                                                                                                                     cAGGTCCCCACATTGAAGCAATACCTGCTGCGATACAGGCAAACGGGTTCGCACCCGAAGAGCCAGCGGTACGCACGGTGGAGGTAGAGGCGTTCTGTTCATGGTCAGCGTGCAGGATCAGAATACGGTCCATAGCACGTTCCAGAATc       >  4:115360/1‑149 (MQ=255)
                                                                                                                                           cccACATTGAAGCAATACCTGCTGCGATACAGGCAAACGGGTTCGCACCCGAAGAGCCAGCGGTACGCACGGTGGAGGTAGAGGCGTTCTGTTCATGGTCAGCGTGCAGGATCAGAATACGGTCCATAGCACGTTCCAGAATCGGAtt  >  3:21544/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                             
ACGCGGTGACCGAAGCCCATCAGGCGGAAAGAATCATTTTTGTCTTTCGCACGACGAACAAATTCCGGAATGTGTTTAACGGAGCTGATTTCTTCCAGCATTTTCAGCGCCGCTTCGTTAGCACCGCCGTGCGCAGGTCCCCACAGTGAAGCAATACCTGCTGCGATACAGGCAAACGGGTTCGCACCCGAAGAGCCAGCGGTACGCACGGTGGAGGTAGAGGCGTTCTGTTCATGGTCAGCGTGCAGGATCAGAATACGGTCCATAGCACGTTCCAGAATCGGATT  >  NC_000913/753546‑753832

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑