Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,032,3980AC62.5% 8.9 / 19.4 24G235G (GGT→GGGxerDsite‑specific tyrosine recombinase XerD
Reads supporting (aligned to +/- strand):  ref base A (6/3);  new base C (0/15);  total (6/18)
Fisher's exact test for biased strand distribution p-value = 6.24e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.11e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

ACGCGTAAATCCGCACCATGATTTAATAAATGAGTGGCAAAAGCGTGACGCAACACATGCGGTGACAGCTTTTCGCTGTCGATACCCGCCAGCACAGCATAATGTTTAATACGGTGCCAGAAGGTCTGTCGCGTCATCTGCTGCGCACGCTGGCTGGGAAACAACACGT  >  NZ_CP009273/3032315‑3032483
                                                                                   |                                                                                     
ccgcgTAAACCCGCACCATGATTTAATAAAGGAGGGGAAAAAGCGTGACGCAACACATGGGGTGACAGCTTTTCGCTGTCGATCCCCGcc                                                                                 <  2:678055/89‑1 (MQ=255)
acGCGTAAATCCGCACCATGATTTAATAAATGAGTGGCAAAAGCGTGACGCAACACATGCGGTGACAGCTTTTCGCTGTCGATACCCGcc                                                                                 >  2:265927/1‑90 (MQ=255)
   cgTAAATCCGCACCATGATTTAATAAATGAGTGGCAAAAGCGTGACGCAACACAGGGGGTGACAGCTTTTCGCTGTCGATACCCGCcagc                                                                              <  2:721511/90‑1 (MQ=255)
    gTAAATCGGCCCTATGATTTAATAAATGGGTGGAAAAAGCGTGACGAAACAAGGGCGGAAAGAGCTTTTCGGTGTCGATCCCCGCcagca                                                                             <  2:545817/90‑1 (MQ=255)
     aaaaCCCGCACCATGATTTAATAAATGGGGGGAAAAGGCGTGACCCAACACATGGGGTGAAAGCTTTTCGCTGTCGATCCCCGCcagcac                                                                            <  2:351938/89‑1 (MQ=255)
          ccccACCTTAATTAAATAACAGGGTGCAAAACGCTTGACCAACCACAGGCGGTGACACTTTTCCGCTGCCGATCCCCGCCAGCACAGCAt                                                                       <  2:121123/87‑1 (MQ=255)
           cGCACCATGATTTAATAAATGAGTGGCAAAAGCGTGACGCAACACATGCGGTGACAGCTTTTCGCTGTCGATACCCGCCAGCACAGCATa                                                                      >  1:545817/1‑90 (MQ=255)
            gCACCATGATTTAATAAATGAGTGGAAAAAGCGTGACGCAACACATGCGGTGACAGCTTTTCGCTGTCGATCCCCGCCAGCACAGCATaa                                                                     <  1:38114/90‑1 (MQ=255)
            gCACCATGATTTAATAAATGAGGGGCAAAAGCGTGACGCACCACATGCGGTGACATTTTTTCGCTGTCGATCCCCGCCAGCACAGCATaa                                                                     <  1:265927/90‑1 (MQ=255)
            ccacaATGATTAAATAAATGAGTGGAAAAAGCCGGACGCAACACATGCGGTGACAGCTTTCCGCTGCCGATCCCCGCCAGCACAGCATaa                                                                     <  2:628886/89‑1 (MQ=255)
            ccACCATGTTTAAATAAATGGGTGAAAAAAGCGGGCCGCAACACATGGGGTGACAGCTTTTCGCTGTCGATCCCCGCCAGCACAGCATaa                                                                     <  2:401655/89‑1 (MQ=255)
            ccACCATGATTTAATAAATGAGTGGAAAAAGCGTGACGCAACACATGCGGTGACAGCTTTTCGCTGTCGATACCCGCCAGCACAGCATaa                                                                     <  1:390361/89‑1 (MQ=255)
            ccACCATAATTAAATAAAGGAGTGGAAAAAGCGTGACCCAACAAATGCGGTGAAAGCTTTCCGTTGTCGATCCCCGCCAGCACAGCATaa                                                                     <  2:372369/89‑1 (MQ=255)
            ccACCATAATTAAAAAAATGAGTGACAAAAGCGTGCCCCAACACATGCGGTGACAGTTTTCCGCTGTCGATCCCCGCCAGCACAGCATaa                                                                     <  1:199328/89‑1 (MQ=255)
                 aTGATTTAATAAATGAGTGGCAAAAGCGTGACGCAACACATGCGGTGACAGCTTTTCGCTGTCGATACCCGCCAGCACAGCATAATGttt                                                                >  1:164069/1‑90 (MQ=255)
                 aTGATTTAATAAATGAGTGGAAAAAGCGTGACCCAAAACATGTGGTGACAGTTTTTCGCTGTCTATCCCCGCCAGCACAGCATAATGttt                                                                <  2:358924/90‑1 (MQ=255)
                      ttAAAAAATGAGTGGCAAAAGCGCGACGCAACACATGCGGTAACAGCTTTTCGTTGTCGATCCCCGCCAGCACAGCATAATGTTTAATAc                                                           <  2:130566/90‑1 (MQ=255)
                       taaaaaAGGGGTGGCAAAAGGGTTACCCAACCCATGGGGTGCCAGTTTTTCGCTTTCAATCCCCGCCAGCACAGCATAATGTTTAATACg                                                          <  1:603940/90‑1 (MQ=255)
                        aataaATGGGTGGAAAAGGCGTGACAAAACACCGGCGGTGACAGTTTTTCGCGGTCGATCCCCGCCAGCACAGCATAATGTTTAATACgg                                                         <  2:567911/90‑1 (MQ=255)
                         ataaATGAGTGGCAAAAGCGTGACGCAACACATGCGGTGACAGCTTTTCGCTGTCGATACCCGCCAGCACAGCATAATGTTTAATACGGt                                                        >  1:130566/1‑90 (MQ=255)
                                 gTGGCAAAAGCGTGACGCAACACATGCGGTGACAGCTTTTCGCTGTCGATACCCGCCAGCACAGCATAATGTTTAATACGGTGCCAGAAg                                                >  1:289367/1‑90 (MQ=255)
                                                            ggTGACAGCTTTTCGCTGTCGATACCCGCCAGCACAGCATaa                                                                     <  1:618103/42‑1 (MQ=255)
                                                            ggTGACAGCTTTTCGCTGTCGATACCCGCCAGCACAGCATaa                                                                     >  2:618103/1‑42 (MQ=255)
                                                                               cGATCCCCGCCAGCACAGCATAATGTTTAATACGGTGCCAGAAGGTCTGTCGCGTCATCTGCTGCGCACGCTGGCTGGGAAACAACACGt  <  1:734769/90‑1 (MQ=255)
                                                                                   |                                                                                     
ACGCGTAAATCCGCACCATGATTTAATAAATGAGTGGCAAAAGCGTGACGCAACACATGCGGTGACAGCTTTTCGCTGTCGATACCCGCCAGCACAGCATAATGTTTAATACGGTGCCAGAAGGTCTGTCGCGTCATCTGCTGCGCACGCTGGCTGGGAAACAACACGT  >  NZ_CP009273/3032315‑3032483

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