Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,544,2650AT70.0% 28.8 / 17.0 20I104N (ATC→AAC) yfeYRpoE‑regulated lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (3/3);  new base T (9/5);  total (12/8)
Fisher's exact test for biased strand distribution p-value = 6.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TTATCATCACCATCAGCTTTTTGGCAATTGCCAAATGCTTTGCTGTAAAGGTCGCTAAACGGTGTACCGATTTTAACACCGGTGTCAGCAGGAATATCGCTATCCAGCACGTCAATGCGACTGATCGTGCCTTGATCGCCGTTAATCACCATCGCCACGTTGTCGCCCTTCATCACTTCAAAAAAGCGCACCACGTTGCCGTTCGCGGTTTTCATTCCGCTGCGCAGGCGATAATCGCCA  >  NZ_CP009273/2544142‑2544381
                                                                                                                           |                                                                                                                    
ttATCATCACCATCAGCTTTTTGGCAATTGCCAAATGCTTTGCTGTAAAGGTCGCTAAACGGTGTACCGATTTTAACACCGGTGTCAGCAGGAATATCGCTATCCAGCACGTCAATGCGACTGATCGTGCCTTGATCGc                                                                                                       <  2:158390/139‑1 (MQ=255)
     atcaCCATCAGCTTTTTGGCAATTGCCAAATGCTTTGCTGTAAAGGTCGCTAAACGGTGTACCGATTTTAACACCGGTGTCAGCAGGAATATCGCTATCCAGCACGTCAATGCGACTGTTCGTGCCTTGATCGCCGtt                                                                                                   >  2:520721/1‑138 (MQ=255)
          cATCAGCTTTTTGGCAATTGCCAAATGCTTTGCTGTAAAGGTCGCTAAACGGTGTACCGATTTTAACACCGGTGTCAGCAGGAATATCGCTATCCAGCACGTCAATGCGACTGTTCGTGCCTTGATCGCCGTTAATCAc                                                                                             >  2:201308/1‑139 (MQ=255)
                                  aTGCTTTGCTGTAAAGGTCGCTAAACGGTGTACCGATTTTAACACCGGTGTCAGCAGGAATATCGCTATCCAGCACGTCAATGCGACTGTTCGTGCCTTGATCGCCGTTAATCACCATCGCCACGTTGTCGCCCTtcat                                                                     <  1:296373/139‑1 (MQ=255)
                                  aTGCTTTGCTGTAAAGGTCGCTAAACGGTGTACCGATTTTAACACCGGTGTCAGCAGGAATATCGCTATCCAGCACGTCAATGCGACTGATCGTGCCTTGATCGCCGTTAATCACCATCGCCACGTTGTCGCCCTtcat                                                                     <  2:492683/139‑1 (MQ=255)
                                                    cGCTAAACGGTGTACCGATTTTAACACCGGTGTCAGCAGGAATATCGCTATCCAGCACGTCAATGCGACTGTTCGTGCCTTGATCGCCGTTAATCACCATCGCCACGTTGTCGCCCTTCATCACTTCAAAAAAGCGcac                                                   >  1:178441/1‑139 (MQ=255)
                                                     gCTAAACGGTGTACCGATTTTAACACCGGTGTCAGCAGGAATATCGCTATCCAGCACGTCAATGCGACTGTTCGTGCCTTGATCGCCGTTAATCACCATCGCCACGTTGTCGCCCTTCATCACTTCAAAAAAGCGcac                                                   >  2:199828/1‑138 (MQ=255)
                                                     gCTAAACGGTGTACCGATTTTAACACCGGTGTCAGCAGGAATATCGCTATCCAGCACGTCAATGCGACTGTTCGTGCCTTGATCGCCGTTAATCACCATCGCCACGTTGTCGCCCTTCATCACTTCAAAAAAGCGcac                                                   >  2:135297/1‑138 (MQ=255)
                                                       tAAACGGTGTACCGATTTTAACACCGGTGTCAGCAGGAATATCGCTATCCAGCACGTCAATGCGACTGATCGTGCCTTGATCGCCGTTAATCACCATCGCCACGTTGTCGCCCTTCATCACTTCAAAAAAGCGcaccac                                                >  1:265703/1‑139 (MQ=255)
                                                            ggTGTACCGATTTTAACACCGGTGTCAGCAGGAATATCGCTATCCAGCACGTCAATGCGACTGTTCGTGCCTTGAt                                                                                                          >  1:18856/1‑76 (MQ=255)
                                                            ggTGTACCGATTTTAACACCGGTGTCAGCAGGAATATCGCTATCCAGCACGTCAATGCGACTGTTCGTGCCTTGAt                                                                                                          <  2:18856/76‑1 (MQ=255)
                                                                      ttttAACACCGGTGTCAGCAGGAATATCGCTATCCAGCACGTCAATGCGACTGTTCGTGCCTTGATCGCCGTTAATCACCATCGCCACGTTGTCGCCCTTCATCACTTCAAAAAAGCGCACCACGTTGCCGTTCGCGGt                                 <  2:178441/139‑1 (MQ=255)
                                                                                      agcagGAATATCGCTATCCAGCACGTCAATGCGACTGATCGTGCCTTGATCGCCGTTAATCACCATCGCCACGTTGTCGCCCTTCATCACTTCAAAAAAGCGCACCACGTTGCCGTTCGCGGTTTTCATTCCGCTgcgc                 <  2:265703/139‑1 (MQ=255)
                                                                                         agGAATATCGCTATCCAGCACGTCAATGCGACTGATCGTGCCTTGATCGCCGTTAATCACCATCGCCACGTTGTCGCCCTTCATCACTTCAAAAAAGCGCACCACGTTGCCGTTCGCGGTTTTCATTCCGCTGCGCAgg              >  1:268550/1‑139 (MQ=255)
                                                                                                     aTCCAGCACGTCAATGCGACTGTTCGTGCCTTGATCGCCGTTAATCACCATCGCCACGTTGTCGCCCTTCATCACTTCAAAAAAGCGCACCACGTTGCCGTTCGCGGTTTTCATTCCGCTGCGCAGGCGATAATCGCCa  >  2:467510/1‑139 (MQ=255)
                                                                                                     aTCCAGCACGTCAATGCGACTGATCGTGCCTTGATCGCCGTTAATCACCATCGCCACGTTGTCGCCCTTCATCACTTCAAAAAAGCGCACCACGTTGCCGTTCGCGGTTTTCATTCCGCTGCGCAGGCGATAATCGCCa  >  1:224589/1‑139 (MQ=255)
                                                                                                      tCCAGCACGTCAATGCGACTGTTCGTGCCTTGATCGCCGTTAATCACCATCGCCACGTTGt                                                                               <  1:566978/61‑1 (MQ=255)
                                                                                                      tCCAGCACGTCAATGCGACTGTTCGTGCCTTGATCGCCGTTAATCACCATCGCCACGTTGt                                                                               >  2:566978/1‑61 (MQ=255)
                                                                                                           cACGTCAATGCGACTGTTCGTGCCTTGATCGCCGTTAATCACCATCGCCACGTTGTCGCCCTTCATCACTTCAAAAAAGCGCACCACGTTGCCGTTCGCGGt                                 >  1:326616/1‑102 (MQ=255)
                                                                                                           cACGTCAATGCGACTGTTCGTGCCTTGATCGCCGTTAATCACCATCGCCACGTTGTCGCCCTTCATCACTTCAAAAAAGCGCACCACGTTGCCGTTCGCGGt                                 <  2:326616/102‑1 (MQ=255)
                                                                                                                           |                                                                                                                    
TTATCATCACCATCAGCTTTTTGGCAATTGCCAAATGCTTTGCTGTAAAGGTCGCTAAACGGTGTACCGATTTTAACACCGGTGTCAGCAGGAATATCGCTATCCAGCACGTCAATGCGACTGATCGTGCCTTGATCGCCGTTAATCACCATCGCCACGTTGTCGCCCTTCATCACTTCAAAAAAGCGCACCACGTTGCCGTTCGCGGTTTTCATTCCGCTGCGCAGGCGATAATCGCCA  >  NZ_CP009273/2544142‑2544381

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