Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1225,223,1640TC38.1% ‑1.1 / 11.3 27intergenic (‑138/+38)cntO_7/CJ019_04917Metal‑pseudopaline receptor CntO/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (4/9);  new base C (6/2);  total (16/11)
Fisher's exact test for biased strand distribution p-value = 8.05e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.57e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

AGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGC‑GGGT‑AAACCCGCCCCTA‑C‑GGGGGGTTTACCCGTG‑AAGAGGC‑TGGCACAGGCAACTACCGAACC‑AGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGG  >  ACB122/5223033‑5223297
                                                                                                                                     |                                                                                                                                          
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)
                                                                                                                                     tA‑C‑GGGGGGTTTACCCGTG‑AAGAGGC‑TGGCACAGGCAACTACCGAACC‑AGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCgg  >  8:17751/1‑134 (MQ=255)
                                                                                                                                     |                                                                                                                                          
AGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGC‑GGGT‑AAACCCGCCCCTA‑C‑GGGGGGTTTACCCGTG‑AAGAGGC‑TGGCACAGGCAACTACCGAACC‑AGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGG  >  ACB122/5223033‑5223297

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