Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1435,222,5100GC34.6% 31.6 / 12.3 27intergenic (‑133/+43)cntO_7/CJ019_04917Metal‑pseudopaline receptor CntO/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (1/16);  new base C (7/2);  total (9/18)
Fisher's exact test for biased strand distribution p-value = 3.97e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGC‑GGG‑TAAACCCGCCCC‑TA‑CGGGGGGTTTACCC‑GTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCA  >  ACB143/5222377‑5222629
                                                                                                                                       |                                                                                                                          
tGGAAGCAGACCTGTGATGGGGTGATTATTGCGACTCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGC‑GGG‑TAAACCCg                                                                                                                            <  8:223426/134‑1 (MQ=255)
    cGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGC‑GGG‑TAAACCCGCCCC‑t                                                                                                                      <  5:63280/135‑1 (MQ=255)
    cGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGC‑GGG‑TAAACCCGCCCC‑t                                                                                                                      <  1:210371/135‑1 (MQ=255)
    cGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGC‑GGG‑TAAACCCGCCCC‑t                                                                                                                      <  5:282835/135‑1 (MQ=255)
      cAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGC‑GGG‑TAAACCCGCCCC‑TA‑c                                                                                                                   <  1:210533/135‑1 (MQ=255)
         cccTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGG‑TATACCC‑CCCC‑GA‑Cggg                                                                                                                >  8:1474/2‑135 (MQ=255)
               gAGGGGGTGATTATTGCGCATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGC‑GGG‑TAAACCCGCCCC‑TA‑CGGGGGGtt                                                                                                           <  1:202878/134‑1 (MQ=255)
               gAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCGCTTCGC‑GGG‑TAAACCCGCCCC‑TA‑CGGGGGGtt                                                                                                           <  6:363174/134‑1 (MQ=255)
                 gggggTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGC‑GGGTTAAACCCCCCCC‑AC‑CGGGGGGTTTa                                                                                                         >  2:119542/1‑135 (MQ=255)
                 gggggTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGC‑GGGTTAAACCCCCCCA‑AC‑CGGGGGGTTTa                                                                                                         >  6:280196/1‑135 (MQ=255)
                 gggggTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGC‑GGG‑GTTACACCCCCC‑AA‑Ccggggggggta                                                                                                        >  7:175375/1‑130 (MQ=255)
                      tGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGG‑TAACACCCCCCC‑AC‑CGGGGGGGTTAcagg                                                                                                     >  3:23020/1‑131 (MQ=255)
                      tGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGC‑GGG‑GTTACCCCCCCC‑AACCGGGGGGTTTGccggg                                                                                                    >  1:113326/1‑132 (MQ=255)
                        attattGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGG‑TAAACCCCCCCC‑AC‑CGGGGGGTTTGCCC‑ggtg                                                                                                 >  3:22374/1‑132 (MQ=255)
                        attattGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGG‑TAAAACCCCCCC‑AC‑CCGGGGGTTTGCCC‑g                                                                                                    >  2:30828/1‑132 (MQ=255)
                                    caatcGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGC‑GGG‑TTAACCCGCCCC‑TA‑CGGGGGGGTTTCCC‑GGGGAGAGGCggggccc                                                                                    >  2:251976/1‑131 (MQ=255)
                                                  cTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGC‑GGG‑TAAACCCGCCCC‑TA‑CGGGGGGTTTACCC‑GTGATGTGGCTGGCACAGGCAACTACCGaa                                                                       <  3:107686/134‑1 (MQ=255)
                                                         cagggATGCCACGTGGTGTTACTTCACCTGGGCCTGTAAATATTTTTGCTGTGGCGGCCTCTTCGC‑GGG‑AAAACCCGCCCC‑TA‑CGGGGGGTTAACCC‑GTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTc                                                               <  8:294180/132‑1 (MQ=255)
                                                          acGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGC‑GGG‑TAAACCCGCCCC‑TA‑CGGGGGGTTTACCC‑GTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTc                                                               <  7:109673/134‑1 (MQ=255)
                                                          acGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGC‑GGG‑TAAACCCGCCCC‑TA‑CGGGGGGTTTACCC‑GTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTc                                                               <  2:56798/134‑1 (MQ=255)
                                                                       gcggTTCCCTCACCTGTGCCTGTAAATATTTTTGCTGTGGCGGCCTCTTCGC‑GGG‑AAAACCCGCCCC‑TA‑GGGGGGGTTTACCC‑GTGATGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAacac                                                 <  7:344228/132‑1 (MQ=255)
                                                                                   ccTGTGCCTGTATATATTTTTGCTGGGGGGGCCTCTTCGC‑GGG‑AAACCCCGCCCC‑TA‑CGGGGGTTTTACCC‑GTGATTAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTg                                      <  6:81681/134‑1 (MQ=255)
                                                                                        tccTGTAAATATTTTTGCTGGGGCGGCCTCTTCGC‑GGG‑AACACCCGCCCC‑TC‑CGGGGGGTTTACCC‑GGGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAAc                                <  4:299547/134‑1 (MQ=255)
                                                                                        gCCTGTAAATTTTTTTGGTGGGGGGGCCTCTGCGC‑GGG‑TAACCCCGCCCC‑TC‑GGGGGGTTTTCCCC‑GTTAAGAGGGTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAAc                                <  5:273091/135‑1 (MQ=255)
                                                                                        gCCTGTAAATATTTTTGCTGGGGCGGCCTCTCCGC‑GGG‑TAAACCCGCCCC‑AA‑CAGGGGGTTAACCC‑GTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAAc                                <  4:350251/135‑1 (MQ=255)
                                                                                        gCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGC‑GGG‑TAAACCCGCCCC‑TA‑CGGGGGGTTTACCC‑GTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAAc                                <  5:185959/135‑1 (MQ=255)
                                                                                                            gggcggcCTATCCGC‑‑GG‑TAAACCCCCCCCATA‑GGGGGGTTTTACCC‑GTTAATAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGaa            <  5:99365/134‑1 (MQ=255)
                                                                                                                         cggggg‑gAAACCCCCCCC‑TT‑CGGGGGTTTTACCCGGTTATTCGGCGGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCca  <  5:318881/131‑1 (MQ=255)
                                                                                                                                       |                                                                                                                          
TGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGC‑GGG‑TAAACCCGCCCC‑TA‑CGGGGGGTTTACCC‑GTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCA  >  ACB143/5222377‑5222629

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