Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1225,223,1590GC55.6% ‑4.2 / 16.2 27intergenic (‑133/+43)cntO_7/CJ019_04917Metal‑pseudopaline receptor CntO/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (3/9);  new base C (13/2);  total (16/11)
Fisher's exact test for biased strand distribution p-value = 1.95e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.43e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGC‑GGGT‑AAACCCGCCCCTA‑C‑GGGGGGTTTACCCGTG‑AAGAGGC‑TGGCACAGGCAACTACCGAACC‑AGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCG  >  ACB122/5223028‑5223284
                                                                                                                                     |                                                                                                                                  
gACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGC‑GGGT‑AAAcccccc                                                                                                                                  >  4:21963/1‑134 (MQ=255)
     aGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGC‑GGGT‑TAACACCCCCCAA‑C‑g                                                                                                                          >  6:167930/1‑135 (MQ=255)
       aCCTGTGAGGGGGTGATTATTGCGGATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGC‑GGGT‑AAACCCCCCCCTA‑C‑cgg                                                                                                                        >  8:126611/1‑135 (MQ=255)
       aCCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCACATCGC‑GGGT‑TAAAACCCCCCCA‑C‑gcg                                                                                                                        >  5:133534/1‑133 (MQ=255)
          tgtgAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGC‑GGGT‑TAAACCCCCCCCA‑A‑gggggg                                                                                                                     >  2:48494/1‑135 (MQ=255)
               gggggTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGC‑GGGT‑TAAACCCCCCCCA‑C‑GGGGGGgttta                                                                                                                >  4:150406/1‑133 (MQ=255)
               gggggTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGC‑GGGG‑TAACCCCCCCCCA‑C‑cgggggggtta                                                                                                                >  2:123417/1‑130 (MQ=255)
                    tGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGC‑GGGTAAAACCCCCCCCAA‑C‑GGGGGGTTTAccgt                                                                                                             >  8:214462/1‑132 (MQ=255)
                    tGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGC‑GGGT‑TAACCCGCCCCGA‑C‑GGGGGGGTTTCCCg                                                                                                             >  4:36743/1‑133 (MQ=255)
                    tGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGC‑GGGT‑AAACCCGCCCCGA‑C‑GGAGGGTTTACCCGt                                                                                                            >  4:45214/1‑134 (MQ=255)
                    tGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGC‑GGGT‑AAAACCCCCCCCA‑C‑GGGGGGGTTTCCCgg                                                                                                            >  1:49587/1‑133 (MQ=255)
                      attattGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGC‑GGGT‑TAACCCGCCCCAA‑C‑GGGGGGGTTACCCCTgca                                                                                                         >  8:304819/1‑133 (MQ=255)
                      attattGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGC‑GGGT‑AAACCCCCCCCTC‑C‑GGGGGGGTTTCCCGTgta                                                                                                         >  8:320099/1‑133 (MQ=255)
                      attattGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTCCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGGT‑AAACCCCCCCCAC‑C‑GGGGGGTTTGCCCgggt                                                                                                          >  7:53770/1‑132 (MQ=255)
                             ccaatcaagcGCACTTACACTCATAACCCGCAGGCCAGGGGCTGTACCCTCACCGGTGCAGGTAAATATTTTGGCGGTGGCGGCCGCTCCGC‑GGGA‑AAACCCGCCCCTA‑C‑GGGGGGTTAACCGGGT‑AAGTGGC‑Tg                                                                                                <  6:377430/133‑1 (MQ=255)
                                               aCTCATTACACGGATGCCCCGGTCTTTTACCGCACCTGTGCCTGGAAATTTTTTTGCGGGGGCGGCCTGTTCGC‑GGGT‑AAACCCGCCCCTA‑C‑GGGGGGTTTACCCGTG‑AAGAGGC‑TGGCACAGGCAACTACCGaa                                                                              <  8:282776/135‑1 (MQ=38)
                                                                cACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGC‑GGGT‑TAACCCCCCCCCA‑CAGGGGGGTTTACCCG‑G‑GACAGGC‑CTGGACAGGCAACAACCGCACCAAGGGCCAAACcccc                                                             >  6:358868/1‑135 (MQ=255)
                                                                cACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGC‑GGGT‑TAACACCCCCCTA‑C‑GGGGGGGGTACCCGGG‑AAGAGGC‑TGGCCCCGGCAACTACCCAACC‑AGGGGCACACcccc                                                             >  8:51016/1‑134 (MQ=255)
                                                                                cccggggCCTGCAAATATTTTTGCTGGGGCGGCCGCTTCGC‑GGGA‑AAACCCGCCCCTC‑G‑GGGGGGTTAACCCGTG‑AAGTGGC‑TGGCACAGGCAACTACCGAACC‑AGGGTCAAACGCCCCAACACGGTCTGCCGTg                                            <  8:126624/129‑1 (MQ=255)
                                                                                aCCTGTGCCTGGAAATATTTTTGCTGGGGCGGCCTCTTTGC‑GGGA‑AAACCCGCCCCTC‑C‑GGGGGGTTTACCCGTG‑AAGAGGC‑TGGCACAGGCAACTACCGAACC‑AGGGTCAAACGCCCCAACACGGTCTGCCGTg                                            <  7:356796/135‑1 (MQ=255)
                                                                                 cgtgtgCCTGCAAATATTTTTGCGGGGGCGGCCTGTTCGC‑GGGA‑ACACCCGCCCCTC‑C‑GGGGGGTTAACCGGTG‑AAGAGGC‑TGGCACAGGCAACTACCGAACC‑AGGGTCAAACGCCCCAACACGGTCTGCCGt                                             <  8:225535/131‑1 (MQ=255)
                                                                                      ccctccaaaTATTTTGGGGGGGGGGGCCTC‑TCCC‑GGGA‑AAACCCCCCCCCACC‑GGGGGGTTAACCGGGG‑AAGTGGC‑GGGCACAGGCAACTACCGAACC‑AGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAAc                                      <  5:376399/129‑1 (MQ=255)
                                                                                        cTGGAAATATTTTTGCTGTGGCGGCCTCTTCGC‑GGGT‑‑AACCCCCCCCTA‑C‑CGGGGGTTTAACCGGG‑AAGCGGC‑GGGCACCGGCAACTACCGAACC‑AGGGTCAAACGCCCCAACACGGTCTGCCGGGCAAAACCCa                                   >  8:27627/1‑135 (MQ=255)
                                                                                               tatTTTTGCTGGGGGGGCCGCTTCGC‑GGGT‑AACCCCGCCCATA‑C‑GGGGGGTTACCCCGGT‑AAGTGGC‑TGGCACAGGCAACTACCGAACC‑AGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCt                             <  8:76059/135‑1 (MQ=255)
                                                                                               tatTTTTGCTGGGGGGGCCGCTGCGC‑GGGT‑AAACCCGCCCCTA‑G‑GGGGGGTTTACCGGGT‑AAGAGGC‑TGGCACAGGCAACTACCGAACC‑AGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCt                             <  3:172380/135‑1 (MQ=255)
                                                                                                        tgggggggcCTCTTCGG‑GGGT‑AACCCCGCCCCTA‑G‑GGGGGGTTAACCCGTG‑AAGAGGC‑TGGCACAGGCAACTACCGAACC‑AGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCg                    <  2:49587/129‑1 (MQ=255)
                                                                                                                    tgcgc‑gGGT‑ACACCGGCCCCCT‑G‑GGGGGGTGTACCGGGT‑ACGTGGC‑GGGCACAGGCAACTACCGAACC‑AGGGTCAACCGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCca        <  7:182369/133‑1 (MQ=255)
                                                                                                                             gccaaCCCCGCCCTA‑C‑GGGGGGTTTACCCGGGCACTCGGCGTGCCAC‑GGCAACTACCGAAGC‑AGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTCCCCAGCGAATCCAGCGCCAGCAGCg  <  8:204195/132‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                  
GACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGC‑GGGT‑AAACCCGCCCCTA‑C‑GGGGGGTTTACCCGTG‑AAGAGGC‑TGGCACAGGCAACTACCGAACC‑AGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCG  >  ACB122/5223028‑5223284

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: ‑