Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,026,4520AC60.0% 9.9 / 16.3 20G141G (GGT→GGGyccAFtsH protease modulator YccA
Reads supporting (aligned to +/- strand):  ref base A (5/3);  new base C (0/12);  total (5/15)
Fisher's exact test for biased strand distribution p-value = 3.61e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.09e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

AAGACCGCGCTGATCGCCAGATGCAGAGCAGGCAGCTGCAGGAAGATATTCGCAACCATACCAATCAGCACCACCACAATACCCGCCATCAGCATACCGCCGAGGAACGACATATCTTTGCGGGTGGTCAGCACATATGCAGAGCAGCAGAAGAACACTAACGC  >  NZ_CP009273/1026372‑1026535
                                                                                |                                                                                   
aaGACCGCGCTGATCGCCAGATGCAGAGCAGGCAGCTGCAGGAAGATATTCGCAACCATACCAATCAGCACCACCACAATACCCGCCATc                                                                            >  2:111136/1‑90 (MQ=255)
aaGACCGCGCTGACCGCCAGATGCAGAGCAGGCAGCTGCAGGAAGATATTCGCAACCATACCAACCAGCCCCCCCACAATCCCCGCCATc                                                                            <  1:445922/90‑1 (MQ=255)
     cgcgCTGATCGCCAGATGCAGAGCAGGCAGCTGCAGGAAGATATTCGCAACCATCCCAACCAGCCCCACCACAATACCCGCCATCAGCAt                                                                       <  2:476693/90‑1 (MQ=255)
                  aGATGCAGAGCAGGCAGCTGCAGGAAGATATTCGCAACCATACCAATCAGCACCACCACAATACCCGCCATCAGCATACCGCCGAGGAAc                                                          <  2:232303/90‑1 (MQ=255)
                  aGATCCAGAGCAGGCAGCTGCGGGAAGATATTCGCACCCATCCCAATCAGCACCACCCCAATCCCCGCCATCAGCATACCGCCGAGGAAc                                                          <  2:444292/90‑1 (MQ=255)
                        agagCAGGCAGCTGCAGGAAGATATTCGCAACCATCCCAATCAGCCCCCCCACAATCCCCGCCATCAGCATACCGCCGAGGAACGACata                                                    <  1:46942/90‑1 (MQ=255)
                                  gCTGCAGGAAGATATTCGCAACCATACCAATCAGCCCCCCCACAATCCCCGCCATCAGCATACCGCCGAGGAACGACATATCTTTGCggg                                          <  2:569518/90‑1 (MQ=255)
                                   cTGCAGGAAGATATTCGCAACCATACCAATCAGCACCACCACAATCCCCGCCATCAGCATACCGCCGAGGAACGACATATCTTTGCGggt                                         <  2:255483/90‑1 (MQ=255)
                                        ggAAGATATTCGCAACCATACCAATCAGCCCCCCCACAATCCCCGCCATCAGCATACCGCCGAGGAACGACATATCTTTGCGGGTGGTCa                                    <  1:192514/90‑1 (MQ=255)
                                        ggAAGATATTCCCAACCATACCAATCAGCACCACCACAATACCCGCCATCAGCATACCGCCGAGGAACGACATATCTTTGCGGGTGGTCa                                    <  2:1372/90‑1 (MQ=255)
                                        gaaaGATATCCGCAACCATACCAATCAGCCCCCCCACAATCCCCGCCATCAGCATACCGCCGAGGAACGACATATCTTTGCGGGTGGTCa                                    <  2:571372/88‑1 (MQ=255)
                                            gATATTCGCAACCATACCAATCAGCACCACCACAATACCCGCCATCAGCATACCGCCGAGGAACGACATATCTTTGCGGGTCGTCAGcac                                >  2:534081/1‑90 (MQ=255)
                                             atatCCGCACCCATACCAACCAGCCCCCCCACAATCCCCGCCATCAGCATACCGCCGAGGAACGACATATCTTTGCGGGTGGTCAGcaca                               <  1:119578/90‑1 (MQ=255)
                                                     aaccaaaccaATCAGCACCACCACAATCCCCGCCATCAGCATACCGCCGAGGAACGACATATCTTTGCGGGTGGTCAGCACATATGCaga                       <  2:437349/90‑1 (MQ=255)
                                                          tACCAATCAGCCCCCCCACAATCCCCGCCATCAGCATACCGCCGAGGAACGACATATCTTTGCGGGTGGTCAGCACATATGCAGagcagc                  <  1:37264/90‑1 (MQ=255)
                                                          tACAAACCAGCCCCACCACAATCCCCGCCATCAGCATACCGCCGAGGAACGACATATCTTTGCGGGTGGTCAGCACATATGCAGagcagc                  <  1:129091/90‑1 (MQ=255)
                                                                   gCACCACCACAATACCCGCCATCAGCATACCGCCGAGGAACGACATATCTTTGCGGGTGGTCAGCACATATGCAGAGCAGCAGAAGAaca         >  2:389728/1‑90 (MQ=255)
                                                                   gCACCACCACAATACCCGCCATCAGCATACCGCCGAGGAACGACATATCTTTGCGGGTGGTCAGCACATATGCAGAGCAGCAGAAGAaca         >  1:3927/1‑90 (MQ=255)
                                                                        accacAATACCCGCCATCAGCATACCGCCGAGGAACGACATATCTTTGCGGGTGGTCAGCACATATGCAGAGCAGCAGAAGAACACTAAc    >  2:577976/1‑90 (MQ=255)
                                                                          cacaATCCCCGCCATCAGCATACCGCCGAGGAACGACATATCTTTGCGGGTGGTCAGCACATATGCAGAGCAGCAGAAGAACACTAACGc  <  1:118943/90‑1 (MQ=255)
                                                                                |                                                                                   
AAGACCGCGCTGATCGCCAGATGCAGAGCAGGCAGCTGCAGGAAGATATTCGCAACCATACCAATCAGCACCACCACAATACCCGCCATCAGCATACCGCCGAGGAACGACATATCTTTGCGGGTGGTCAGCACATATGCAGAGCAGCAGAAGAACACTAACGC  >  NZ_CP009273/1026372‑1026535

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: