Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,434,4380TC72.7% 24.3 / 15.5 22E216G (GAA→GGA) fabBbeta‑ketoacyl‑ACP synthase I
Reads supporting (aligned to +/- strand):  ref base T (2/4);  new base C (7/9);  total (9/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.36e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAAACACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGTTA  >  NZ_CP009273/2434307‑2434572
                                                                                                                                   |                                                                                                                                      
aGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGgt                                                                                                                                 >  1:557259/1‑139 (MQ=255)
 gCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGgtg                                                                                                                                <  1:326761/139‑1 (MQ=255)
        aTTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGt                                                                                                                         <  1:353646/139‑1 (MQ=255)
           tCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGTAt                                                                                                                       <  2:557259/138‑1 (MQ=255)
                     tgtgAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTCCCGGGGTGTCGTTGTATTTAGTAGACAg                                                                                                            <  2:452355/139‑1 (MQ=255)
                       tgAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTCCCGGGGTGTCGTTGTATTTAGTAGACAgcg                                                                                                          <  2:353009/139‑1 (MQ=255)
                              cacGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTCCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCAt                                                                                                   >  1:492145/1‑139 (MQ=255)
                                                       ttCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTCCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTCGAATTCGCAAGCCATTTcc                                                                          >  1:345294/1‑139 (MQ=255)
                                                                        cgccgccgcCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTcgcc                                                         >  2:112951/1‑139 (MQ=255)
                                                                          ccgccgccAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTCCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTCGAATTCGCAAGCCa                                                                               <  2:67284/115‑1 (MQ=255)
                                                                          ccgccgccAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTCCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTCGAATTCGCAAGCCa                                                                               <  2:506135/115‑1 (MQ=255)
                                                                          ccgccgccAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTCCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTCGAATTCGCAAGCCa                                                                               >  1:67284/1‑115 (MQ=255)
                                                                          ccgccgccAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTCCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTCGAATTCGCAAGCCa                                                                               >  1:506135/1‑115 (MQ=255)
                                                                              cgccAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTCCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCgccgccgcc                                                   <  1:15737/139‑1 (MQ=255)
                                                                                            aaaCCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAAACACGATGt                                     <  1:112951/139‑1 (MQ=255)
                                                                                            aaaCCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTCCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAAACACGATGt                                     <  1:332427/139‑1 (MQ=255)
                                                                                            aaaCCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTCCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAAACACGATGt                                     >  1:224453/1‑139 (MQ=255)
                                                                                                     cGGTGAGCGTCGTAAGTACGGGAGGCTTTTCCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAAACACGATGTCCTGTTTGc                            <  2:224453/139‑1 (MQ=255)
                                                                                                            cgtcgtAAGTACGGGAGGCTTTTCCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAAACACGATGTCCTGTTTGCCCAGTTg                     >  2:113869/1‑139 (MQ=255)
                                                                                                             gtcgtAAGTACGGGAGGCTTTTCCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAAACACGATGTCCTGTTTGCCCAGTTgg                    <  2:345294/139‑1 (MQ=255)
                                                                                                              tcgtAAGTACGGGAGGCTTTTCCCGGGGTGTCGtt                                                                                                                           >  2:27228/1‑35 (MQ=255)
                                                                                                              tcgtAAGTACGGGAGGCTTTTCCCGGGGTGTCGtt                                                                                                                           <  1:27228/35‑1 (MQ=255)
                                                                                                                               ttttcCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAAACACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGTTa  >  1:479531/1‑139 (MQ=255)
                                                                                                                                   |                                                                                                                                      
AGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAAACACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGTTA  >  NZ_CP009273/2434307‑2434572

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