Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077793,284,7760TG50.0% 40.4 / 14.8 36G79G (GGT→GGGagaCPTS system N‑acetylgalactosamine‑specific transporter subunit IIC
Reads supporting (aligned to +/- strand):  ref base T (2/16);  new base G (18/0);  total (20/16)
Fisher's exact test for biased strand distribution p-value = 4.19e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GACTGGCTTAATTACCGGTGGTCTGACAGAGTTGGCTTTCGCCGGATTAACCCCTGCAGGTGGTGTTCAGCCGCCCAACCCGATTATGGCGGGTCTGATGACCACCGTCATTGCATGGTCTACGGGCGTTGATGCCAAAACAGCAATTGGTCTTGGCCTGCCGTTTAGTTTGTTAATGCAGTACGTCATT  >  NC_007779/3284683‑3284872
                                                                                             |                                                                                                
gACTGGCTTAATTACCGGTGGTCTGACAGAGTTGGCTTTCGCCGGATTAACCCCTGCAGGTGGTGTTCAGCCGCCCAACCCGATTATGGCGGGGCtgatga                                                                                           >  2:85668/1‑101 (MQ=255)
     gCTTAATTACCGGTGGTCTGACAGAGTTGGCTTTCGCCGGATTAACCCCTGCAGGTGGGGTTCAGCCGCCCAACCCGATTATGGCGGGTCtgctgagaacc                                                                                      >  1:514113/1‑96 (MQ=255)
         aaTTACCGGTGGTCTGACAGAGTTGGCTTTCGCCGGATTAACCCCTGCAGGTGGTGTTCAGCCGCCCAACCCGATTATGGCGGGGGTGGGGaccacccccc                                                                                  >  1:17584/1‑97 (MQ=255)
            tACCGGTGGTCTGACAGAGTTGGCTTTCGCCGGATTAACCCCTGCAGGTGGTGTTCAGCCGCCCAACCCGATTATGGCGGGGCCGGTGactaccccccttc                                                                               >  1:585950/1‑94 (MQ=255)
                ggtggtCTGACAGAGTTGGCTTTCGCCGGATTAACCCCTGCAGGTGGTGTTCAGCCGCCCAACCCGATTATGGCGGGTCTGATGACCACCGTCATTGCATg                                                                           <  2:116183/101‑1 (MQ=255)
                  tggtCTGACAGCGTTGGCTTTCGCCGGATTAACCCCTGCAGGTGGTGTTCAGCCGCCCAACCCGATTATGGCGGGGCTGATGACCACCCTCATTGCATggg                                                                         >  4:8994/1‑100 (MQ=255)
                   ggtCTGACAGAGTTGGCTTTCGCCGGATTAACCCCTGCAGGTGGTGTTCAGCCGCCCAACCCGATTATGGCGGGTCTGATGACCACCGTCATTGCATGGTc                                                                        <  1:453105/101‑1 (MQ=255)
                     tCTGACAGAGTTGGCTTTCGCCGGATTAACCCCTGCAGGTGGTGTTCAGCCGCCCAACCCGATTATGGCGGGGGTGATGACCACCGTCATTGCATGGTCTa                                                                      >  2:29041/1‑101 (MQ=255)
                        gACAGAGTTGGCTTTCGCCGGATTAACCCCTGCAGGTGGTGTTCAGCCGCCCAACCCGATTATGGCGGGTCTGATGACCACCGTCATTGCATGGTCTACgg                                                                   <  3:751770/101‑1 (MQ=255)
                        gACAGAGTTGGCTTTCGCCGGATTAACCCCTGCAGGTGGTGTTCAGCCGCCCAACCCGATTATGGCGGGTCTGATGACCACCGTCATTGCATGGTCTACgg                                                                   <  2:4864/101‑1 (MQ=255)
                          cAGAGTTGGCTTTCGCCGGATTAACCCCTGCAGGTGGTGTTCAGCCGCCCAACCCGATTATGGCGGGGCTGGTGACCCCCCGCCTTTCATTGTCTTCCGGc                                                                 >  3:1220870/1‑101 (MQ=255)
                             agTTGGCTTTCGCCGGATTAACCCCTGCAGGTGGTGTTCAGCCGCCCAACCCGATTATGGCGGGGCTGATGAGCACCCGCCTTTCTTTGTCTCCGGGCggt                                                              >  1:179243/1‑99 (MQ=255)
                                tGGCTTTCGCCGGATTAACCCCTGCAGGTGGTGTTCAGCCGCCCAACCCGATTATGGCGGGGGTGATGACCGCCGGCCTTTCATGGTCTTCGGGCggcggt                                                           >  1:216547/1‑96 (MQ=255)
                                  gCTTTCGCCGGATTAACCCCTGCAGGTGGTGTTCAGCCGCCCAACCCGATTATGGCGGGTCTGATGACCACCGTCATTGCATGGTCTACGGGCGTTGATGc                                                         <  3:1594030/101‑1 (MQ=255)
                                  gCTTTCGCCGGATTAACCCCTGCAGGTGGTGTTCAGCCGCCCAACCCGATTATGGCGGGGGTGGTGACCACCCCCCTTTCCTGGTCTACGGGCGgtggtgg                                                         >  3:614878/1‑100 (MQ=255)
                                          cGGATTCACCCATGCAGGTGGTGTTCAGCCGCCCAACCCGATTATGGCGGGGCTGGTGGCCACCCTCCTTGCATGGGCTACGGGCGGTGGTGGCCAAACAg                                                 >  4:120765/1‑101 (MQ=255)
                                                aaCCCCTGCAGGTGGTGTTCAGCCGCCCAACCCGATTATGGCGGGGCTGGTGACCACCCCCCTTGCATGGTCTACGGGCGGGGGTGCCAAAACAACAAATg                                           >  4:3530/1‑101 (MQ=255)
                                                 aCCCCTGCAGGTGGTGTTCAGCCGCCCAACCCGATTATGGCGGGGCTGATGACCACCGGCATTTCATGGGCTACGGGCGGTGGTGGCAAAACAACAAATgg                                          >  1:1525308/1‑101 (MQ=255)
                                                    ccTGCAGGTGGTGTTCAGCCGCCCAACCCGATTATGGCGGGTCTGATGACCACCGTCATTGCATGGTCTACGGGCGTTGATGCCAAAACAGCAATTGGTCt                                       <  4:15914/101‑1 (MQ=255)
                                                     cTGCAGGTGGTGTTCAGCCGCCCAACCCGATTATGGCGGGTCTGATGACCACCGTCATTGCATGGTCTACGGGCGTTGATGCCAAAACAGCAATTGGTCtt                                      <  1:587978/101‑1 (MQ=255)
                                                     cTGCAGGTGGTGTTCAGCCGCCCAACCCGATTATGGCGGGTCTGATGACCACCGTCATTGCATGGTCTACGGGCGTTGATGCCAAAACAGCAATTGGTCtt                                      <  1:791572/101‑1 (MQ=255)
                                                          ggtggtGTTCAGCCGCCCAACCCGATTATGGCGGGGCTGATGACCACCGTCATTTCATGGTCTACGGGCGTTGATGCCAAAACAACAATTGGTCTTGGCCt                                 >  3:503992/1‑101 (MQ=255)
                                                           gtggtgTTCAGCCGCCCAACCCGATTATGGCGGGGCTGATGACCACCCCCCTTTCATGGTCTACGGGCGGTGATGGCAAAAAAGCAAATGGGCTTGGcccg                                >  1:1509880/1‑99 (MQ=255)
                                                                 ttCAGCCGCCCAACCCGATTATGGCGGGTCTGATGACCACCGTCATTGCATGGTCTACGGGCGTTGATGCCAAAACAGCAATTGGTCTTGGCCTGCCGttt                          <  1:265215/101‑1 (MQ=255)
                                                                  tCAGCCGCCCAACCCGATTATGGCGGGTCTGATGACCACCGTCATTGCATGGTCTACGGGCGTTGATGCCAAAACAGCAATTGGTCTTGGCCTGCCGTTTa                         <  3:1126201/101‑1 (MQ=255)
                                                                  tCAGCCGCCCAACCCGATTATGGCGGGTCTGATGACCACCGTCATTGCATGGTCTACGGGCGTTGATGCCAAAACAGCAATTGGTCTTGGCCTGCCGTTTa                         <  1:913085/101‑1 (MQ=255)
                                                                    aGCCGCCCAACCCGATTATGGCGGGTCTGATGACCACCGTCATTGCATGGTCTACGGGCGTTGATGCCAAAACAGCAATTGGTCTTGGCCTGCCGTTTAGt                       <  1:865646/101‑1 (MQ=255)
                                                                       cgccCAACCCGATTATGGCGGGGCTGATGACCACCCTCATTGCATGGTCTACGGGCGGTGATGCCAAAACAACAAATTGTCTTTGCCTGCCGTTTAGtttt                    >  2:153138/1‑100 (MQ=255)
                                                                         cccAACCCGATTATGGCGGGGCTGATGACCACCCCCATTGCATGGGCTTCGGGGGGTGGTGGCAAAACAACAAATGGGCTTGGCCTGCCGTTTtgtttgtt                  >  1:1540246/1‑101 (MQ=255)
                                                                            aaCCCGATTATGGCGGGGTTGATGCCCACCGTGATTTGATGGGCTACCGGGGGTGGTGCCAAAACAGCCAATGGGCTTTGGCTGCCGGTTTGTTTGTTTAt               >  1:1180740/1‑101 (MQ=255)
                                                                              cccGATTATGGCGGGTCTGATGACCACCGTCATTGCATGGTCTACGGGCGTTGATGCCAAAACAGCAATTGGTCTTGGCCTGCCGTTTAGTTTGTTAATGc             <  3:988839/101‑1 (MQ=255)
                                                                                 gATTATGGCGGGGCTGATGACCACCCTCATTGCATGGTCTACGGGGGGTGATGCCAAAACAGCAAATGGTCTTGGCCTGCCGTTTTGTTTTTTTATGCAGt          >  1:1448319/1‑101 (MQ=255)
                                                                                      tGGCGGGTCTGATGACCACCGTCATTGCATGGTCTACGGGCGTTGATGCCAAAACAGCAATTGGTCTTGGCCTGCCGTTTAGTTTGTTAATGCAGTACGTc     >  3:623714/1‑101 (MQ=255)
                                                                                        gCGGGTCTGATGACCACCGTCATTGCATGGTCTACGGGCGTTGATGCCAAAACAGCAATTGGTCTTGGCCTGCCGTTTAGTTTGTTAATGCAGTACGTCAt   <  2:77044/101‑1 (MQ=255)
                                                                                         cGGGTCTGATGACCACCGTCATTGCATGGTCTACGGGCGTTGATGCCAAAACAGCAATTGGTCTTGGCCTGCCGTTTAGTTTGTTAATGCAGTACGTCAtt  <  1:499513/101‑1 (MQ=255)
                                                                                         cGGGTCTGATGACCACCGTCATTGCATGGTCTACGGGCGTTGATGCAAAAACAGCAATTGGTCTTGGCCTGCCGTTTAGTTTGTTAATGCAGTACGTCAtt  <  3:553959/101‑1 (MQ=255)
                                                                                             |                                                                                                
GACTGGCTTAATTACCGGTGGTCTGACAGAGTTGGCTTTCGCCGGATTAACCCCTGCAGGTGGTGTTCAGCCGCCCAACCCGATTATGGCGGGTCTGATGACCACCGTCATTGCATGGTCTACGGGCGTTGATGCCAAAACAGCAATTGGTCTTGGCCTGCCGTTTAGTTTGTTAATGCAGTACGTCATT  >  NC_007779/3284683‑3284872

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