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.

TGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCA  >  ACB143/5222377‑5222629
                                                                                                                                       |                                                                                                                          
tGGAAGCAGACCTGTGATGGGGTGATTATTGCGACTCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCg                                                                                                                            <  8:223426/134‑1 (MQ=255)
    cGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCt                                                                                                                      <  5:63280/135‑1 (MQ=255)
    cGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCt                                                                                                                      <  1:210371/135‑1 (MQ=255)
    cGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCt                                                                                                                      <  5:282835/135‑1 (MQ=255)
      cAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTAc                                                                                                                   <  1:210533/135‑1 (MQ=255)
         cccTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGGTATACCCCCCCGACggg                                                                                                                >  8:1474/2‑135 (MQ=255)
               gAGGGGGTGATTATTGCGCATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGtt                                                                                                           <  1:202878/134‑1 (MQ=255)
               gAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCGCTTCGCGGGTAAACCCGCCCCTACGGGGGGtt                                                                                                           <  6:363174/134‑1 (MQ=255)
                 gggggTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTTAAACCCCCCCCACCGGGGGGTTTa                                                                                                         >  2:119542/1‑135 (MQ=255)
                 gggggTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTTAAACCCCCCCAACCGGGGGGTTTa                                                                                                         >  6:280196/1‑135 (MQ=255)
                 gggggTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGGTTACACCCCCCAACcggggggggta                                                                                                        >  7:175375/1‑130 (MQ=255)
                      tGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGGTAACACCCCCCCACCGGGGGGGTTAcagg                                                                                                     >  3:23020/1‑131 (MQ=255)
                      tGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGGTTACCCCCCCCAACCGGGGGGTTTGccggg                                                                                                    >  1:113326/1‑132 (MQ=255)
                        attattGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGGTAAACCCCCCCCACCGGGGGGTTTGCCCggtg                                                                                                 >  3:22374/1‑132 (MQ=255)
                        attattGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGGTAAAACCCCCCCACCCGGGGGTTTGCCCg                                                                                                    >  2:30828/1‑132 (MQ=255)
                                    caatcGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTTAACCCGCCCCTACGGGGGGGTTTCCCGGGGAGAGGCggggccc                                                                                    >  2:251976/1‑131 (MQ=255)
                                                  cTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGATGTGGCTGGCACAGGCAACTACCGaa                                                                       <  3:107686/134‑1 (MQ=255)
                                                         cagggATGCCACGTGGTGTTACTTCACCTGGGCCTGTAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGAAAACCCGCCCCTACGGGGGGTTAACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTc                                                               <  8:294180/132‑1 (MQ=255)
                                                          acGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTc                                                               <  7:109673/134‑1 (MQ=255)
                                                          acGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTc                                                               <  2:56798/134‑1 (MQ=255)
                                                                       gcggTTCCCTCACCTGTGCCTGTAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGAAAACCCGCCCCTAGGGGGGGTTTACCCGTGATGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAacac                                                 <  7:344228/132‑1 (MQ=255)
                                                                                   ccTGTGCCTGTATATATTTTTGCTGGGGGGGCCTCTTCGCGGGAAACCCCGCCCCTACGGGGGTTTTACCCGTGATTAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTg                                      <  6:81681/134‑1 (MQ=255)
                                                                                        tccTGTAAATATTTTTGCTGGGGCGGCCTCTTCGCGGGAACACCCGCCCCTCCGGGGGGTTTACCCGGGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAAc                                <  4:299547/134‑1 (MQ=255)
                                                                                        gCCTGTAAATTTTTTTGGTGGGGGGGCCTCTGCGCGGGTAACCCCGCCCCTCGGGGGGTTTTCCCCGTTAAGAGGGTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAAc                                <  5:273091/135‑1 (MQ=255)
                                                                                        gCCTGTAAATATTTTTGCTGGGGCGGCCTCTCCGCGGGTAAACCCGCCCCAACAGGGGGTTAACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAAc                                <  4:350251/135‑1 (MQ=255)
                                                                                        gCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAAc                                <  5:185959/135‑1 (MQ=255)
                                                                                                            gggcggcCTATCCGCGGTAAACCCCCCCCATAGGGGGGTTTTACCCGTTAATAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGaa            <  5:99365/134‑1 (MQ=255)
                                                                                                                         cggggggAAACCCCCCCCTTCGGGGGTTTTACCCGGTTATTCGGCGGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCca  <  5:318881/131‑1 (MQ=255)
                                                                                                                                       |                                                                                                                          
TGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCA  >  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: