Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273982,4700AC52.9% 12.6 / 20.1 17intergenic (‑32/+571)ompF/asnSporin OmpF/asparagine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base A (6/2);  new base C (0/9);  total (6/11)
Fisher's exact test for biased strand distribution p-value = 2.26e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.53e-03
Rejected as consensus: Frequency below/above cutoff threshold.

TACCTGCTACTAACAGAGCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCAT  >  NZ_CP009273/982384‑982556
                                                                                      |                                                                                      
tACCTGCTACTAACAGAGCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCt                                                                                     >  2:23603/1‑90 (MQ=255)
 aCCTGCTACTAACAGAGCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTg                                                                                    >  1:1677/1‑90 (MQ=255)
         ctaACAGAGCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGcc                                                                            >  2:199168/1‑90 (MQ=255)
         caaaCAGAGCAGGGACGATCCCTCCAAAAATATTCCCCTTCATCATTATTTATTCCCCTCATGGTTTTTTTTATGACCCCTGCCACTGcc                                                                            <  2:697911/88‑1 (MQ=255)
                         gATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTc                                                            >  2:424457/1‑90 (MQ=255)
                                         ttGCGCTTCATAATTATTTATTCCCCTCATGTTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGa                                            <  2:848493/90‑1 (MQ=255)
                                           gcgcTTCATCTTTATTTATTCCCCTCATGTTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAAc                                          <  2:850591/90‑1 (MQ=255)
                                           gcgcTTAATCATTATTTATTCCCCTCATGTTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAAc                                          <  2:397192/90‑1 (MQ=255)
                                                    catTATTTTTTCCCCTCATGTTTTTTTTTATGCCCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATtgatga                                 <  2:37963/90‑1 (MQ=255)
                                                              tccccTCATGTTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGgtg                       <  2:207466/88‑1 (MQ=255)
                                                              tccccTCATGGTTTTTTTTATGCCCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGgtg                       <  1:686163/88‑1 (MQ=255)
                                                              tccccTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGgtg                       <  1:77706/88‑1 (MQ=255)
                                                                   tCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCttt                  >  2:286612/1‑90 (MQ=255)
                                                                    cATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTgg                         >  2:335763/1‑82 (MQ=255)
                                                                    cATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTgg                         <  1:335763/82‑1 (MQ=255)
                                                                         ttttttttATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATgtgt            <  1:176348/90‑1 (MQ=255)
                                                                                   gACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCAt  <  1:199168/90‑1 (MQ=255)
                                                                                      |                                                                                      
TACCTGCTACTAACAGAGCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCAT  >  NZ_CP009273/982384‑982556

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