Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,057,8730CG64.7% 37.0 / 38.2 34W418C (TGG→TGCyihOMFS transporter
Reads supporting (aligned to +/- strand):  ref base C (7/5);  new base G (12/10);  total (19/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.73e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ATTTCCTCAATAATCAATGCAAACCGTTTTTCATTGAGTGTGTAAAAGAATCCCATCGTCAATGCAGCAATAATTGCCAGTGCACAAGGCCAGATAAAAATGAGCTGACGCAAACCCTGTAATGTCGCATCGCTTTGCGCGATGTTGGGAACATAACCAATTTGCGTCAGCATAATGCC  >  NZ_CP009273/4057784‑4057962
                                                                                         |                                                                                         
aTTTCCTCAATAATCAATGCAAACCGTTTTTCATTGAGTGTGTAAAAGAATCCCATCGTCAATGCAGCAATAATTGCCAGTGCACAAGGc                                                                                           <  2:427124/90‑1 (MQ=255)
           aatCAATGCAAACCGTTTTTCATTGAGTGTGTAAAAGAATCCCATCGTCAATGCAGCAATAATTGCCAGTGCACAAGGGCAGATAAAAAt                                                                                <  1:167652/90‑1 (MQ=255)
              cAATGCAAACCGTTTTTCATTGAGTGTGTAAAAGAATCCCATCGTCAATGCAGCAATAATTGCCAGTGCACAAGGGCAGATAAAAATGAg                                                                             >  2:197845/1‑90 (MQ=255)
              cAATGCAAACCGTTTTTCATTGAGTGTGTAAAAGAATCCCATCGTCAATGCAGCAATAATTGCCAGTGCACAAGGGCAGATAAAAATGAg                                                                             >  1:523622/1‑90 (MQ=255)
               aaTGCAAACCGTTTTTCATTGAGTGTGTAAAAGAATCCCATCGTCAATGCAGCAATAATTGCCAGTGCACAAGGGCAGATAAAAATGAGc                                                                            <  1:332064/90‑1 (MQ=255)
                 tGCAAACCGTTTTTCATTGAGTGTGTAAAAGAATCCCATCGTCAATGCAGCAATAATTGCCAGTGCACAAGGCCAGATAAAAATGAGCTg                                                                          >  2:197413/1‑90 (MQ=255)
                     aaCCGTTTTTCATTGAGTGTGTAAAAGAATCCCATCGTCAATGCAGCAATAATTGCCAGTGCACAAGGGCAGATAAAAATGAGCTGACGc                                                                      >  2:368105/1‑90 (MQ=255)
                     aaCCGTTTTTCATTGAGTGTGTAAAAGAATCCCATCGTCAATGCAGCAATAATTGCCAGTGCACAAGGGCAGATAAAAATGAGCTGACGc                                                                      >  1:677183/1‑90 (MQ=255)
                       ccGTTTTTCATTGAGTGTGTAAAAGAATCCCATCGTCAATGCAGCAATAATTGCCAGTGCACAAGGCCAGATAAAAATGAGCTGACGCaa                                                                    <  2:611181/90‑1 (MQ=255)
                         gTTTTTCATTGAGTGTGTAAAAGAATCCCATCGTCAATGCAGCAATAATTGCCAGTGCACAAGGGCAGATAAAAATGAGCTGACGCAAAc                                                                  <  1:774708/90‑1 (MQ=255)
                          tttttCATTGAGTGTGTAAAAGAATCCCATCGTCAATGCAGCAATAATTGCCAGTGCACAAGGGCAGATAAAAATGAGCTGACGCAAAcc                                                                 >  1:59011/1‑90 (MQ=255)
                              tCATTGAGTGTGTAAAAGAATCCCATCGTCAATGCAGCAATAATTGCCAGTGCACAAGGGCAGATAAAAATGAGCTGACGCAAACCCTGt                                                             <  1:710552/90‑1 (MQ=255)
                                   gAGTGTGTAAAAGAATCCCATCGTCAATGCAGCAATAATTGCCAGTGCACAAGGGCAGATAAAAATGAGCTGACGCAAACCCTGTAATGt                                                        >  1:824187/1‑90 (MQ=255)
                                             aaGAATCCCATCGTCAATGCAGCAATAATTGCCAGTGCACAAGGGCAGATAAAAATGAGCTGACGCAAACCCTGTAATGTCGCATCGCtt                                              <  2:59011/90‑1 (MQ=255)
                                               gAATCCCATCGTCAATGCAGCAATAATTGCCAGTGCACAAGGGCAGATAAAAATGAGCTGACGCAAACCCTGTAATGTCGCATCGCTTTg                                            >  1:115646/1‑90 (MQ=255)
                                                 aTCCCATCGTCAATGCAGCAATAATTGCCAGTGCACAAGGGCAGATAAAAATGAGCTGACGCAAACCCTGTAATGTCGCATCGCTTTgcg                                          >  2:570206/1‑90 (MQ=255)
                                                 aTCCCATCGTCAATGCAGCAATAATTGCCAGTGCACAAGGGCAGATAAAAATGAGCTGACGCAAACCCTGTAATGTCGCATCGCTTTgcg                                          >  1:14854/1‑90 (MQ=255)
                                                 aTCCCATCGTCAATGCAGCAATAATTGCCAGTGCACAAGGGCAGATAAAAATGAGCTGACGCAAACCCTGTAATGTCGCATCGCTTTgcg                                          >  1:442986/1‑90 (MQ=255)
                                                 aTCCCATCGTCAATGCAGCAATAATTGCCAGTGCACAAGGCCAGATAAAAATGAGCTGACGCa                                                                     <  1:348629/63‑1 (MQ=255)
                                                 aTCCCATCGTCAATGCAGCAATAATTGCCAGTGCACAAGGCCAGATAAAAATGAGCTGACGCa                                                                     >  2:348629/1‑63 (MQ=255)
                                                 aTCCCATCGTCAATGCAGCAATAATTGCCAGTGCACAAGGCCAGATAAAAATGAGCTGACGCAAACCCTGTAATGTCGCATCGCTTTgcg                                          >  2:102730/1‑90 (MQ=255)
                                                  tCCCATCGTCAATGCAGCAATAATTGCCAGTGCACAAGGGCAGATAAAAATGAGCTGACGCAAACCCTGTAATGTCGCATCGCTTTgcgc                                         <  1:368105/90‑1 (MQ=255)
                                                  tCCCATCGTCAATGCAGCAATAATTGCCAGTGCACAAGGCCAGATAAAAATGAGCTGACGCAAACCCTGTAATGTCGCATCGCTTTgcgc                                         <  1:156702/90‑1 (MQ=255)
                                                    ccATCGTCAATGCAGCAATAATTGCCAGTGCACAAGGGCAGATAAAAATGAGCTGACGCAAACCCTGTAATGTCGCATCGCTTTGCGCGa                                       <  1:288310/90‑1 (MQ=255)
                                                    ccATCGTCAATGCAGCAATAATTGCCAGTGCACAAGGGCAGATAAAAATGAGCTGACGCAAACCCTGTAATGTCGCATCGCTTTGCGCGa                                       <  1:152569/90‑1 (MQ=255)
                                                         gTCAATGCAGCAATAATTGCCAGTGCACAAGGGCAGATAAAAATGAGCTGACg                                                                       >  1:267097/1‑53 (MQ=255)
                                                         gTCAATGCAGCAATAATTGCCAGTGCACAAGGGCAGATAAAAATGAGCTGACg                                                                       <  2:267097/53‑1 (MQ=255)
                                                                    aataatTGCCAGTGCACAAGGGCAGATAAAAATGAGCTGACGCAAACCCTGTAATGTCGCATCGCTTTGCGCGATGTTGGGAACATAAcc                       >  1:523714/1‑90 (MQ=255)
                                                                              aGTGCACAAGGGCAGATAAAAATGAGCTGACGCAAACCCTGTAATGTCGCATCGCTTTGCGCGATGTTGGGAACATAACCAATTTGCGTc             <  2:43106/90‑1 (MQ=255)
                                                                              aGTGCACAAGGCCAGATAAAAATGAGCTGACGCAAACCCTGTAATGTCGCATCGCTTTGCGCGATGTTGGGAACATAACCAATTTGCGTc             <  2:104098/90‑1 (MQ=255)
                                                                              aGTGCACAAGGCCAGATAAAAATGAGCTGACGCAAACCCTGTAATGTCGCATCGCTTTGCGCGATGTTGGGAACATAACCAATTTGCGTc             <  1:223110/90‑1 (MQ=255)
                                                                               gTGCACAAGGCCAGATAAAAATGAGCTGACGCAAACCCTGTAATGTCGCATCGCTTTGCGCGATGTTGGGAACATAACCAATTTGCGTCa            >  1:691555/1‑90 (MQ=255)
                                                                               gTGCACAAGGCCAGATAAAAATGAGCTGACGCAAACCCTGTAATGTCGCATCGCTTTGCGCGATGTTGGGAACATAACCAATTTGCGTCa            >  2:608253/1‑90 (MQ=255)
                                                                               gTGCACAAGGCCAGATAAAAATGAGCTGACGCAAACCCTGTAATGTCGCATCGCTTTGCGCGATGTTGGGAACATAACCAATTTGCGTCa            >  1:130925/1‑90 (MQ=255)
                                                                                      aGGCCAGATAAAAATGAGCTGACGCAAACCCTGTAATGTCGCATCGCTTTGCGCGATGTTGGGAACATAACCAATTTGCGTCAGCATAAt     >  1:84483/1‑90 (MQ=255)
                                                                                         ccAGATAAAAATGAGCTGACGCAAACCCTGTAATGTCGCATCGCTTTGCGCGATGTTGGGAACATAACCAATTTGCGTCAGCATAATGcc  >  1:500046/1‑90 (MQ=255)
                                                                                         |                                                                                         
ATTTCCTCAATAATCAATGCAAACCGTTTTTCATTGAGTGTGTAAAAGAATCCCATCGTCAATGCAGCAATAATTGCCAGTGCACAAGGCCAGATAAAAATGAGCTGACGCAAACCCTGTAATGTCGCATCGCTTTGCGCGATGTTGGGAACATAACCAATTTGCGTCAGCATAATGCC  >  NZ_CP009273/4057784‑4057962

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