Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,544,5040TG29.6% 41.4 / 18.0 27intergenic (+63/+51)mdfA/ampGmultidrug MFS transporter/muropeptide permease AmpG
Reads supporting (aligned to +/- strand):  ref base T (10/9);  new base G (2/6);  total (12/15)
Fisher's exact test for biased strand distribution p-value = 2.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGT‑AAACCCGCTCCTACACAGGGCCGCGCCGC‑TTTGTAGGAGCGGGTTTACCCGCGAAG‑AGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCGATCAGGATCAGCGTAGGCACC  >  NC_002947/1544361‑1544651
                                                                                                                                                 |                                                                                                                                                    
ggggggggggggGCGTTGGGGCTGGGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGT‑AAACCCGCTCCTACAAAGCGGCGCGGCCC‑TGTGTAg                                                                                                                                                 <  8:20156/149‑1 (MQ=25)
ggggggggggggGCGTGGGCGCTGGGGCCGCTGTGGTTGGTGGGGGGTGGGCGGCCTGGGGAACGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAAACCCGCTCCTACAAAGCGGCGCGGCCC‑TGTGTa                                                                                                                                                  <  1:38711/149‑1 (MQ=21)
  gggggggggggCGTGGGCGCGGGGGCCGCTGGGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGT‑AAACCCGCTCCTACACAGGGCCGCGCCGC‑TTTGTAGGa                                                                                                                                               <  6:80046/143‑1 (MQ=38)
    gggggggggCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGGGGGGGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGT‑AAACCCGCTCCTACAAAGCGGCGCGGCCC‑TGTGTAGGAGc                                                                                                                                             <  4:196883/145‑1 (MQ=25)
      gggggggCGTGGGCGCTGTGGCCCCTGGGGTGTGGGGTGGGGGGCCGGCCTGGGGAAGGGTCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGT‑AAACCCGCTCCTACAAAGCGGCGCGGCCC‑TGTGTAGGAGCgg                                                                                                                                           <  7:63100/147‑1 (MQ=18)
      gCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGT‑AAACCCGCTCCTACACAGGGCCGCGCCGC‑TTTGTa                                                                                                                                                  >  5:215112/1‑142 (MQ=255)
      gCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGT‑AAACCCGCTCCTACACAGGGCCGCGCCGC‑TTTGTa                                                                                                                                                  <  6:215112/142‑1 (MQ=255)
      cggggggCGGTGGCGCTGTGGCCGCTGGGGTTTGGGGGGGGGGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGT‑AAACCCGCTCCTACAAAGCGGCGCGGCCC‑TGTGTAGGAGCgg                                                                                                                                           <  6:110393/147‑1 (MQ=21)
              gtgggCGCTGTGGCTGCTGTGGTTGGTGGTGGGGGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGT‑AAACCCGCTCCTACACAGGGCCGCGCCGC‑TTTGTAGGAGCGGGTTTAccc                                                                                                                                   <  3:36215/146‑1 (MQ=255)
                          gCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTAGCGGGT‑TAAACCCCTCCGACACAGGGCCGCGCCGCGTTTTTTGAGGGGGGTTTACCCCCGAA‑‑AGGccc                                                                                                                      >  3:121109/1‑148 (MQ=255)
                                                  gcggcCTGGGTAAGGGGGTGTAGCTCTATAATCGGTATTGCCGGGGCGGGCCTCTTCGCGGTT‑AAACCCGCTCCTACACAGGGCCGCGCGGC‑TTTGTAGGAGGGGGTTTACCCGCGAAG‑AGGCCAGCACAGGCAATACAAAGCTTCAgg                                                                                              <  4:121109/149‑1 (MQ=37)
                                                             aaGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCACTTAGCGGGT‑GAACCCGCTACAACACAGGGCCGCGCCGC‑TTTGTAGGCGCGGGGTTACACGCGAAC‑AGGCCAGCACAGGCAAAACAAACCATAAAGGCCGCTCAga                                                                                    >  6:228845/1‑147 (MQ=255)
                                                                   cggaaGCTCGAAAACCGGTATTGCCGGGGCCGGCCTCTTCGCGGGA‑AACCCCGCTCCACCACCGGCCCGCGCCGC‑TTTGTAGCAGCGGGTTAACCCGCGAAGCAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGa                                                                              <  1:236583/145‑1 (MQ=255)
                                                                                                   gccTCTTCGCGGGT‑AAACCCGCTCCTACACAGGGCCGCGCCGC‑GTTGTAGGAGCGGGGTTACCCGCGAAG‑AGGCCAGCCCCGGCAATACAAAGCTTCAGGGCCGCTCAGCCCCCGAGTGGTCGTCCGCCACCGGCGGCCcccgccccgc                                             >  4:283138/1‑149 (MQ=25)
                                                                                                      tctTCGCGGGT‑AAACCCGCTCCTACACAGGGCCGCGCCGC‑TTTGTAGGAGCGGGTTTACCCGCGAAG‑AGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTc                                          >  5:334017/1‑149 (MQ=255)
                                                                                                         tCGCGGGT‑AAACCCGCTCCTACAAAGCGGCGCGGCCC‑TGTGTAGGAGCGGGTTTACCCGCGAAG‑AGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCt                                        >  3:76475/1‑148 (MQ=21)
                                                                                                                      ccGCTCCTACACAGGGCCGCGCCGC‑TTTGTAGGAGCGGGTTTACCCGCGAAG‑AGGCCAGCACAGGCAATACAAAGCTTCa                                                                                                <  2:17575/80‑1 (MQ=255)
                                                                                                                      ccGCTCCTACACAGGGCCGCGCCGC‑TTTGTAGGAGCGGGTTTACCCGCGAAG‑AGGCCAGCACAGGCAATACAAAGCTTCa                                                                                                >  1:17575/1‑80 (MQ=255)
                                                                                                                        gCTCCTACAAAGCGGCGCGGCCC‑TGTGTAGGAGCGGGTTTACCCGCGAAG‑AGGCCGGCCCCGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAggg                         >  8:300695/1‑149 (MQ=21)
                                                                                                                          tCCTACACAGGGCCGCGCCGC‑TTTGTAGGAGCGGGTTTACCCGCGAAG‑AGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGc                        >  2:124110/1‑148 (MQ=255)
                                                                                                                          tCCTACACAGGGCCGCGCCGC‑TTTGTAGGAGCGGGTTTACCCGCGAAG‑AGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGc                        >  1:36938/1‑148 (MQ=255)
                                                                                                                            cTACACAGGGCCGCGCCGC‑TTTGTAGGAGCGGGTTTACCCGCGAAG‑AGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgat                     >  5:130653/1‑149 (MQ=255)
                                                                                                                               cacaGGGCCGCGCCGC‑TTTGTAGGAGCGGGTTTACCCGCGAAG‑AGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatcag                  <  5:228845/149‑1 (MQ=255)
                                                                                                                               cacaGGGCCGCGCCGC‑TTTGTAGGAGCGGGTTTACCCGCGAAG‑AGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatca                   <  6:334017/148‑1 (MQ=255)
                                                                                                                                      gcgcggccc‑TGTGTAGGAGCGGGTTTACCCGCGAAG‑AGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCGATCAGGATCAGc           <  4:76475/148‑1 (MQ=255)
                                                                                                                                          gccgc‑TTTGTAGGAGCGGGTTTACCCGCGAAG‑AGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCGATCAGGATCAGCGTAg       <  1:124110/149‑1 (MQ=255)
                                                                                                                                           ccgc‑TTTGTAGGAGCGGGTTTACCCGCGAAG‑AGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCGATCAGGATCAGCGTAgg      >  6:52727/1‑149 (MQ=255)
                                                                                                                                                tttGTAGGAGCGGGTTTACCCGCGAAG‑AGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCGATCAGGATCAGCGTAGGCAcc  >  1:87555/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGT‑AAACCCGCTCCTACACAGGGCCGCGCCGC‑TTTGTAGGAGCGGGTTTACCCGCGAAG‑AGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCGATCAGGATCAGCGTAGGCACC  >  NC_002947/1544361‑1544651

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑