Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913953,0881.C20.0% 48.9 / 14.1 20coding (467/2283 nt)pflBformate C‑acetyltransferase 1, anaerobic, pyruvate formate‑lyase 1
Reads supporting (aligned to +/- strand):  ref base . (7/9);  new base C (2/2);  total (9/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AACTGTGCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAGTTCGCGGTTGTACGCTTTGCAGGAACCTTCGATCATTTTGATACCACCGAACGGGAT  >  NC_000913/952953‑953224
                                                                                                                                        |                                                                                                                                        
aaCTGTGCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCg                                                                                                                                       >  1:239640/1‑139 (MQ=255)
aaCTGTGCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCg                                                                                                                                       <  1:97693/139‑1 (MQ=255)
                      tcaGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGcc                                                                                                                 <  2:498519/139‑1 (MQ=255)
                                   aCCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGtt                                                                                                    <  2:22778/139‑1 (MQ=255)
                                     cGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGtttt                                                                                                  <  2:180562/139‑1 (MQ=255)
                                       tACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTa                                                                                                 >  1:495781/1‑139 (MQ=255)
                                       tACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGATCAGAACACCAGATTTACGGCAACGCAGGATGTCCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTa                                                                                                 >  1:233853/1‑139 (MQ=255)
                                            cgcAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGccc                                                                                                                >  1:155073/1‑118 (MQ=255)
                                            cgcAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGccc                                                                                                                <  2:155073/118‑1 (MQ=255)
                                            cgcAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTAt                                                                                           >  1:60543/1‑139 (MQ=255)
                                               aaCGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGccc                                                                                                                >  1:204496/1‑115 (MQ=255)
                                               aaCGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGccc                                                                                                                <  2:204496/115‑1 (MQ=255)
                                                                   tgatACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGatcat                                                                    <  2:239640/139‑1 (MQ=255)
                                                                     atACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTAc                                                                                                >  1:62341/1‑109 (MQ=255)
                                                                     atACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTAc                                                                                                <  2:62341/109‑1 (MQ=255)
                                                                         ggccacggccaTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGAt                                                               <  2:495781/139‑1 (MQ=255)
                                                                         ggccacggccaTATGCATCTGGCAGACCGATCAGAACACCAGATTTACGGCAACGCAGGATGTCCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGAt                                                               <  2:233853/139‑1 (MQ=255)
                                                                                           cTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTtgtg                                                                                                        <  2:251701/79‑1 (MQ=255)
                                                                                           cTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTtgtg                                                                                                        >  1:251701/1‑79 (MQ=255)
                                                                                                                                     tGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAGTTCGCGGTTGTACGCTTTGCAGGAACCTTCGATCATTTTGATACCACCGAACGGGAt  >  2:449447/1‑139 (MQ=255)
                                                                                                                                        |                                                                                                                                        
AACTGTGCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAGTTCGCGGTTGTACGCTTTGCAGGAACCTTCGATCATTTTGATACCACCGAACGGGAT  >  NC_000913/952953‑953224

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: