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

GGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCGATCAGGATCAGCGTA  >  NC_002947/1544360‑1544645
                                                                                                                                                 |                                                                                                                                             
ggggggggggggggCGTGGGCGCGGGGGCCGCTGTGGTTTGGGGGGGGGGGGCGGCCGGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTa                                                                                                                                           <  7:135619/149‑1 (MQ=38)
  ggggggggggggggTTGGCGGTGGGCCCGCTGTGGTTTGTGGTGGGGGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAgg                                                                                                                                         <  4:296762/149‑1 (MQ=21)
   ggTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCCCATAGCCGGGTTAACCCCACCAACACAAGGCCCCGCCGCGTTGTTggg                                                                                                                                        >  5:293579/1‑148 (MQ=255)
            ggCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTa                                                                                                                               <  2:131221/149‑1 (MQ=255)
                        gtgGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCCCAACAACACACGGCCGCGCCGCTTTGTATGAGCGGGTTTTCACGCGAAGAgg                                                                                                                   >  7:189114/1‑149 (MQ=38)
                             cGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAAAACCCAGCTAAACAGGGCCGCGCCGCGTTGTAGGAGCGGGGTTTCCCCCGCAGAGGCCAg                                                                                                               >  3:63637/1‑148 (MQ=255)
                                   gttttGTGGTGGGTGGGCGGCCTGGGTATGGTGGTGTAGGTCGATAATGTGTATTGCCGGGGCGGGCGTTTTCGGGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGc                                                                                                        <  6:195399/147‑1 (MQ=25)
                                            tgggtgggCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTGCTAAACAGGGCCGCGCCGCGTTGTAGGAGCGGGTTTACACCCGAAGAGGCCAGCAAACGCAATAAAAAg                                                                                               >  7:181283/1‑149 (MQ=37)
                                                tgggGGGCCTGGGTATGGGGCTGTAGCTCGATAATGGGTATTGCCGGGGCCGGCCTCTTCGGGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTc                                                                                           <  6:99228/149‑1 (MQ=37)
                                                 gttgggcGTGGGTAAGGGGTGGTTGTTCTATATTCGGTATGGCGGGGGCGGGCTTGTTCGGGGTTAAACCCGCTCCAACACCGGGCGGGGCGTTTTTTTTGTAGGGGTTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCa                                                                                          <  6:293579/145‑1 (MQ=255)
                                                                       tAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTATACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTgg                                                                    <  2:280455/149‑1 (MQ=255)
                                                                                                         ttCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCt                                  >  7:155698/1‑149 (MQ=255)
                                                                                                           cgcgGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTc                                    >  1:125162/1‑145 (MQ=255)
                                                                                                                    aaCCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTg                        >  1:88113/1‑149 (MQ=255)
                                                                                                                    aaCCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCCGGCAATACCAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTg                        >  1:200313/1‑149 (MQ=18)
                                                                                                                        cGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCCGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAgg                    >  6:4229/1‑149 (MQ=21)
                                                                                                                        cGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAgg                    >  7:141177/1‑149 (MQ=37)
                                                                                                                        cGCTCATACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAgg                    <  4:63637/149‑1 (MQ=255)
                                                                                                                         gCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCCGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAggg                   <  5:68431/149‑1 (MQ=21)
                                                                                                                         gCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCCGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAggg                   >  6:68431/1‑149 (MQ=21)
                                                                                                                           tCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCGCTGCAGGGc                  >  5:301127/1‑148 (MQ=255)
                                                                                                                           tCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCCGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGc                  >  5:277810/1‑148 (MQ=21)
                                                                                                                           tCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCt                                  >  1:146331/1‑132 (MQ=25)
                                                                                                                              tACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCg                                                                                    <  7:127716/79‑1 (MQ=255)
                                                                                                                                          cgccgcTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCGCTGCAGGGCGATCAGGATCAGCGTa  <  6:301127/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                             
GGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCGATCAGGATCAGCGTA  >  NC_002947/1544360‑1544645

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