Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913424,5180AC33.3% 37.7 / 13.4 27V134G (GTG→GGG) acpHacyl carrier protein phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base A (18/0);  new base C (0/9);  total (18/9)
Fisher's exact test for biased strand distribution p-value = 2.13e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CAGACGTGGGCGGCGGCTTGCCATGCCGTTTAACACGTTCTGGATGAAATCCATATCGCGATAGCGCACCAGCCACTGCTCTGACCACAAGTAATTGTTCAGATTGATAAAACGTGGCGGTGAGTCCGGCAAAAT  >  NC_000913/424451‑424585
                                                                   |                                                                   
cAGACGTGGGCGGCGGCTTGCCATGCCGTTTAACACGTTCTGGATGAAATCCATATCGCGATAGCGCa                                                                     >  1:544119/1‑68 (MQ=255)
cAGACGTGGGCGGCGGCTTGCCATGCCGTTTAACACGTTCTGGATGAAATCCATATCGCGATAGCGCa                                                                     >  1:2050180/1‑68 (MQ=255)
  gACGTGGGCGGCGGCTTGCCATGCCGTTTAACACGTTCTGGATGAAATCCATATCGCGATAGCGCAcc                                                                   >  1:1413631/1‑68 (MQ=255)
    cGTGGGCGGCGGCTTGCCATGCCGTTTAACACGTTCTGGATGAAATCCATATCGCGATAGCGCACCac                                                                 >  1:3780700/1‑67 (MQ=255)
    cGTGGGCGGCGGCTTGCCATGCCGTTTAACACGTTCTGGATGAAATCCATATCGCGATAGCGCACCAg                                                                 >  1:905286/1‑68 (MQ=255)
    cGTGGGCGGCGGCTTGCCATGCCGTTTAACACGTTCTGGATGAAATCCATATCGCGATAGCGCACCAg                                                                 >  1:2177206/1‑68 (MQ=255)
    cGTGGGCGGCGGCTTGCCATGCCGTTTAACACGTTCTGGATGAAATCCATATCGCGATAGCGCACCAg                                                                 >  1:207050/1‑68 (MQ=255)
    cGTGGGCGGCGGCTTGCCATGCCGTTTAACACGTTCTGGATGAAATCCATATCGCGATAGCGCACCAg                                                                 >  1:1563578/1‑68 (MQ=255)
         gcggcggcTTGCCATGCCGTTTAACACGTTCTGGATGAAATCCATATCGCGATAGCGCACCAGCCACt                                                            >  1:2156996/1‑68 (MQ=255)
         gcggcggcTTGCCATGCCGTTTAACACGTTCTGGATGAAATCCATATCGCGATAGCGCACCAGCCACt                                                            >  1:839762/1‑68 (MQ=255)
         gcggcggcTTGCCATGCCGTTTAACACGTTCTGGATGAAATACATATCGCGATAGCGCACCAGCCACt                                                            >  1:2678483/1‑68 (MQ=255)
                cTTGCGGTGCGGTTTAACGCGTTCTGGATAAAATCCATATCGGGATAGCCCCCCAGCCACTGCTCTGa                                                     <  1:523196/68‑1 (MQ=255)
                cTTGCCATGCCTTTTAACCCGTTCTGGATAAAACCCATATCGCGATACCGCCCCAGCCACTGCTCTGa                                                     <  1:3160583/68‑1 (MQ=255)
                cTTGCCATGCCGTTTAACCCGTTCTGGATGAAATCCATATCGCGATAGCGCCCCAGCCACTGCTCTGa                                                     <  1:1628058/68‑1 (MQ=255)
                cTTGCCATGCCGTTTAACACGTTCTGGATGAAATCCATATCGCGATAGCGCCCCAGCCACTGCTCTGa                                                     <  1:699678/68‑1 (MQ=255)
                cTTGCCATGCCGTTTAACACGTTCTGGATGAAATCCATATCGCGATAGCCCCCCAGCCACTGCTCTGa                                                     <  1:3007997/68‑1 (MQ=255)
                cTTGCCATGCCGTTTAACACGTTCTGGATAAAATCCATACCGCGATAGCGCCCCAGCCACTGCTCTGa                                                     <  1:1240997/68‑1 (MQ=255)
                      atcccGTTTAACCCGTTCTGGATGAAACCCATATCGCGATAGCGCCCCAGCCACTGCTCTGACCACaa                                               <  1:2622897/65‑1 (MQ=255)
                         ccGTTTAACACGTTCTGGATAAAACCCATACCGCGATGGCGCCCCAGCCACTGCTCTGACCACAAGTa                                            <  1:1714570/68‑1 (MQ=255)
                          cGTTTAACACGTTCTGGATGAAATCCACATCGCGCTAGCGCACCAGCCACTGCTCTGCCCACAAGTaa                                           >  1:2853685/1‑68 (MQ=255)
                              tAACAGTTTCTGGATGAAATCCATTCCGCGATAGCGCACCAGCCACTGCTCTGACCACAAGTAATTGt                                       <  1:1895424/68‑1 (MQ=255)
                              tAACACGTTCTGGATGAAATCCATATCGCGATAGCGCCCCAGCCACTGCTCTGACCACAAGTAATTGt                                       <  1:2669800/68‑1 (MQ=255)
                                acacGTTCTGGATGAAATCCATATCGCGATAGCGCACCAGCCACTGCTCTGACCACAAGTAATTGTTc                                     >  1:2240436/1‑68 (MQ=255)
                                  acGTTCTGGATGAAATCCATATCGCGATAGCGCACCAGCCACTGCTCTGACCACAAGTAATTGTTCAg                                   >  1:3762478/1‑68 (MQ=255)
                                   cGTTCTGGATGAAGTCCATCTCGCGATAGCGCCCCAGCCACTGCTCTGACCACAAGTAATTGTTCAGa                                  <  1:11211/68‑1 (MQ=255)
                                      tCTGGATGAAATCCATATCGCGATAGCGCACCAGCCACTGCTCTGACCACAAGTAATTGTTCAGATTg                               >  1:2911962/1‑68 (MQ=255)
                                      tCTGGATGAAATCCATATCGCGATAGCGCACCAGCCACTGCTCTGACCACAAGTAATTGTTCAGATTg                               >  1:3026253/1‑68 (MQ=255)
                                       cTGGATGAAATCCATATCGCGATAGCGCCCCAGCCACTGCTCTGACCACAAGTAATTGTTCAGATTGa                              <  1:826198/68‑1 (MQ=255)
                                       agggATGAAATCCATATCGCGATAGCGCCCCAGCCACTGCTCTGACCACAAGTAATTGTTCAGATTGa                              <  1:3841953/66‑1 (MQ=255)
                                       agggATAAAAGCCATATCGCGATAGCCCCCCAGCCACTGCTCTGACCACAAGTAATTGTTCAGATTGa                              <  1:1001491/66‑1 (MQ=255)
                                          gATGAAATCCATATCGCGATAGCGCACCAGCCACTGCTCTGACCACAAGTAATTGTTCAGATTGATaa                           >  1:3573675/1‑68 (MQ=255)
                                          gATGAAATCCATATCGCGATAGCGCACCAGCCACTGCTCTGACCACAAGTAATTGTTCAGATTGATaa                           >  1:1222352/1‑68 (MQ=255)
                                            taaaaTCCATATCGCGATAGCGCCCCAGCCACTGTTCTGACCACAAGTAATTGTTCAGATTGATaaaa                         <  1:2677449/66‑1 (MQ=255)
                                                 tCCATATCGCGATAGCGCACCAGCCACTGCTCTGACCACAAGTAATTGTTCAGATTGATAAAACGTgg                    >  1:2121667/1‑68 (MQ=255)
                                                         gcgATACCGCCCCAGCCACTGCTCTGACCACAAGTAATTGTTCAGATTGATAAAACGTGGCGGTGAGt            <  1:50139/68‑1 (MQ=255)
                                                              aGCGCACCAGCCACTGCTCTGACCACCAGTAATTGTTCAGATTGATAAACCGCGGCGGTGAGTCCGGc       >  1:2890411/1‑68 (MQ=255)
                                                              aGCGCACCAGCCACTGCTCTGACCACAAGTAATTGTTCAGATTGATAAAACGTGGCGGTGAGTCCGGc       >  1:2288687/1‑68 (MQ=255)
                                                                 gcACCAGCCACTGCTCTGACCACAAGTAATTGTTCAGATTGATAAAACGTGGCGGTGAGTCCGGCaaa    >  1:3478082/1‑68 (MQ=255)
                                                                 gcACCAGCCACTGCTCTGACCACAAGTAATTGTTCAGATTGATAAAACGTGGCGGTGAGTCCGGCaaa    >  1:1247337/1‑68 (MQ=255)
                                                                   aCCAGCCACTGCTCTGACCACAAGTAATTGTTCAGCTTGATAAAACGTGGCGGTGAGTCCGGCAAAAt  >  1:3191135/1‑68 (MQ=255)
                                                                   aCCAGCCACTGCTCTGACCACAAGTAATTGTTCAGATTGATAAAACGTGGCGGTGAGTCCGGCAAAAt  >  1:2079918/1‑68 (MQ=255)
                                                                   |                                                                   
CAGACGTGGGCGGCGGCTTGCCATGCCGTTTAACACGTTCTGGATGAAATCCATATCGCGATAGCGCACCAGCCACTGCTCTGACCACAAGTAATTGTTCAGATTGATAAAACGTGGCGGTGAGTCCGGCAAAAT  >  NC_000913/424451‑424585

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