Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273341,0150GA36.0% 23.5 / 10.2 25intergenic (+313/‑100)BW25113_RS01680/BW25113_RS01685hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (13/3);  new base A (0/9);  total (13/12)
Fisher's exact test for biased strand distribution p-value = 1.08e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGG  >  NZ_CP009273/340935‑341100
                                                                                |                                                                                     
ttCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGc                                                                              <  2:279573/90‑1 (MQ=255)
     atGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTaa                                                                         <  1:191263/90‑1 (MQ=255)
        aaaaaaaTCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGc                                                                      >  1:103841/1‑90 (MQ=255)
                             tGTTTGTGGGGGGTTTTTTCTGTATCTTACGCACCGCCCTCAACCCTGAAAAAAAATATGCTGAAAGGCTCATATAAAAAACCGCCATTa                                                 <  2:27204/90‑1 (MQ=255)
                              gTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCaaa                                                            >  1:116116/1‑78 (MQ=255)
                              gTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCaaa                                                            <  2:116116/78‑1 (MQ=255)
                              gTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAg                                                >  1:27204/1‑90 (MQ=255)
                                ttGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCACATCAAAAACCGCCATTAGCt                                              >  2:149503/1‑90 (MQ=255)
                                ttGTGGGGGTTTTTTTCTTTATCTTCCCCACCCCCCCCAAGCCTAAAAAAAAATATGCTGAAAGGCTCATAAAAAAAACCGCCATTAGCt                                              <  1:109076/90‑1 (MQ=255)
                                ttGTGGGGGGTTTTTTCTGTATCTTCCGCATCGCACTAAACCCTGAAAAAAAATATGCTGTAAGGCTCATATAAAAAACCGCCATTAGCt                                              <  2:214037/90‑1 (MQ=255)
                                ttGTGGGGGGTTTTTTCTGTATCTTACGCATCCCACTCAAGCCTGACAAAAAATATGCTGTAAGGCTCATATAAAAACCCGCCATTAGCt                                              <  2:310126/90‑1 (MQ=255)
                                 tgtgtgGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCtc                                             >  1:328912/1‑90 (MQ=255)
                                     tgGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCtcatca                                         >  2:105552/1‑90 (MQ=255)
                                       gggTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAgg                                       >  2:305289/1‑90 (MQ=255)
                                        ggTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGa                                      >  1:281594/1‑90 (MQ=255)
                                                    aTCTTACGCATCGCACTCAAGCCTGACAAAAAATATGCTGTAAGGCTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAgacgac                          <  2:54622/90‑1 (MQ=255)
                                                       ttACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAAc                       >  2:44887/1‑90 (MQ=255)
                                                       ttACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAAc                       >  2:49409/1‑90 (MQ=255)
                                                            cATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGa                  >  1:274861/1‑90 (MQ=255)
                                                              tCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAtt                >  1:324291/1‑90 (MQ=255)
                                                              tCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAtt                >  2:144764/1‑90 (MQ=255)
                                                                  ccTCAAGCCTGAAAAAAAATATGTTGTAGGGCTAATAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAttgttg            <  2:182818/89‑1 (MQ=255)
                                                                  aCTCAAGCCTAAAAAAAAATATGCTGAAAGGCTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAttgttg            <  2:31056/90‑1 (MQ=255)
                                                                      aaGCCTAAAAAAAAATATGCTGAAAGCCTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATg        <  1:78705/90‑1 (MQ=255)
                                                                            gACAAAAAATGTGCTGTAAGGCTCATAAAAAAACCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACgg  <  1:170124/90‑1 (MQ=255)
                                                                                |                                                                                     
TTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGG  >  NZ_CP009273/340935‑341100

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