Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,018,3930TC53.8% 7.3 / 16.4 13E319G (GAG→GGG) yccSYccS family putative transporter
Reads supporting (aligned to +/- strand):  ref base T (1/5);  new base C (7/0);  total (8/5)
Fisher's exact test for biased strand distribution p-value = 4.66e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.56e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TTATTATGGGGTAGTGCCTGGGCCTGTTCTGATTCAATTGTTGCCAGTTGGGCATCAATGGCGCGTAAATTGTTCAGCAAAAATCCCAGTGTTTTGAGTAAATCGGCGGGTGCGCCGTTATCGCGCATCCGCTCCAGCGCAGCATCAATATGCGTAAAAGCGCGCTCAAAATGCGGATCATGTTGATAAGGCTGACGCAACAAAATACAGCGTGACAGTTGCTGGCACGCC  >  NZ_CP009273/1018261‑1018491
                                                                                                                                    |                                                                                                  
ttattaTGGGGTAGTGCCTGGGCCTGTTCTGATTCAATTGTTGCCAGTTGGGCATCAATGGCGCGTAAATTGTTCAGCAAAAATCCCAGTGTTTTGAGTAAATCGGCGGGTGCGCCGTTATCGCGCATCCGCTCCAgcg                                                                                              <  1:467995/139‑1 (MQ=255)
    agacaggTAGTGCCTGGGCCTGTTCTGATTCAATTGTTGCCAGTTGGGCATCAATGGCGCGTAAATTGTTCAGCAAAAATCCCAGTGTTTTGAGTAAATCGGCGGGTGCGCCGTTATCGCGCATCCGCTCCAGCgcagc                                                                                          <  1:79974/134‑1 (MQ=255)
        gggTAGTGCCTGGGCCTGTTCTGATTCAATTGTTGCCAGTTGGGCATCAATGGCGCGTAAATTGTTCAGCAAAAATCCCAGTGTTTTGAGTAAATCGGCGGGTGCGCCGTTATCGCGCATCCGCTCCAGCGCAGCATCa                                                                                      <  1:388910/139‑1 (MQ=255)
           tAGTGCCTGGGCCTGTTCTGATTCAATTGTTGCCAGTTGGGCATCAATGGCGCGTAAATTGTTCAGCAAAAATCCCAGTGTTTTGAGTAAATCGGCGGGGGGGCCGGTATCGCGCATCCCCCCCACGGCCGCATCCata                                                                                   >  2:467995/1‑139 (MQ=255)
              tGCCTGGGCCTGTTCTGATTCAATTGTTGCCAGTTGGGCATCAATGGCGCGTAAATTGTTCAGCAAAAATCCCAGTGTTTTGAGTAAATCGGCGGGGGGGCCGTTTTCGCGCATCCCCCCCCGCGCAGCATCAatattc                                                                                >  2:79181/1‑137 (MQ=255)
               gCCTGGGCCTGTTCTGATTCAATTGTTGCCAGTTGGGCATCAATGGCGCGTAAATTGTTCAGCAAAAATCCCAGTGTTTTGAGTAAATCGGCGGGGGGGCCGTTATTGCCCATCCGCCCCAGCGCAGCATCAataatcg                                                                               >  1:136526/1‑135 (MQ=255)
                     gCCTGTTCTGATTCAATTGTTGCCAGTTGGGCATCAATGGCGCGTAAATTGTTCAGCAAAAATCCCAGTGTTTTGAGTAAATCGGCGGGTGGGCCGTTTTCGCGCATCCGCCCCAGCGCAGCATCAATATGCgctgtct                                                                         >  2:79974/1‑133 (MQ=255)
                        tGTTCTGATTCAATTGTTGCCAGTTGGGCATCAATGGCGCGTAAATTGTTCAGCAAAAATCCCAGTGTTTTGAGTAAATCGGCGGGTGGGCCGTTATCGCGCCAACGCCCCAGCGCAGCATCAATATATGTTAAAgagc                                                                      >  1:11122/1‑139 (MQ=255)
                               aTTCAATTGTTGCCAGTTGGGCATCAATGGCGCGTAAATTGTTCAGCAAAAATCCCAGTGTTTTGAGTAAATCGGCGGGGGCGCCGTTTTCGCGCCTCCGCCCCCGCGCCGCATCAATATGCGTAAAAACGCGCTcca                                                                >  2:98291/1‑136 (MQ=255)
                                                                 tAAATTGTTCAGCAAAAATCCCAGTGTTTTGAGTAAATCGGCGGGTGGGCCGGTATCGCGCATCCGCCCCAGCGCAGCATCCATATGCGTAAAAACGCGCTCAAAATGCGGATCATGTTGATAAGGCTGACGcaacaa                              >  1:167118/1‑138 (MQ=255)
                                                                                             ttGAGTAAATCGGCGGGTGCGCCGTTATCGCGCATCCGCTCCAGCGCAGCATCAATATGCGTAAAAGCGCGCTCAAAATGCGGATCATGTTGATAAGGCTGACGCAACAAAATACAGCGTGACAGTTGCTGGCACGcc  <  2:167118/138‑1 (MQ=255)
                                                                                                                             cATCCGCTCCAGCGCAGCATCAATATGCGTAAAAGCGCGCTCAAAATGCGGATCAt                                                    <  1:121953/56‑1 (MQ=255)
                                                                                                                             cATCCGCTCCAGCGCAGCATCAATATGCGTAAAAGCGCGCTCAAAATGCGGATCAt                                                    >  2:121953/1‑56 (MQ=255)
                                                                                                                                    |                                                                                                  
TTATTATGGGGTAGTGCCTGGGCCTGTTCTGATTCAATTGTTGCCAGTTGGGCATCAATGGCGCGTAAATTGTTCAGCAAAAATCCCAGTGTTTTGAGTAAATCGGCGGGTGCGCCGTTATCGCGCATCCGCTCCAGCGCAGCATCAATATGCGTAAAAGCGCGCTCAAAATGCGGATCATGTTGATAAGGCTGACGCAACAAAATACAGCGTGACAGTTGCTGGCACGCC  >  NZ_CP009273/1018261‑1018491

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