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.

GACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCG  >  ACB122/5223028‑5223284
                                                                                                                                     |                                                                                                                                  
gACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAAcccccc                                                                                                                                  >  4:21963/1‑134 (MQ=255)
     aGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTTAACACCCCCCAACg                                                                                                                          >  6:167930/1‑135 (MQ=255)
       aCCTGTGAGGGGGTGATTATTGCGGATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCCCCCCTACcgg                                                                                                                        >  8:126611/1‑135 (MQ=255)
       aCCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCACATCGCGGGTTAAAACCCCCCCACgcg                                                                                                                        >  5:133534/1‑133 (MQ=255)
          tgtgAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTTAAACCCCCCCCAAgggggg                                                                                                                     >  2:48494/1‑135 (MQ=255)
               gggggTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTTAAACCCCCCCCACGGGGGGgttta                                                                                                                >  4:150406/1‑133 (MQ=255)
               gggggTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGGTAACCCCCCCCCACcgggggggtta                                                                                                                >  2:123417/1‑130 (MQ=255)
                    tGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAAACCCCCCCCAACGGGGGGTTTAccgt                                                                                                             >  8:214462/1‑132 (MQ=255)
                    tGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTTAACCCGCCCCGACGGGGGGGTTTCCCg                                                                                                             >  4:36743/1‑133 (MQ=255)
                    tGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCGACGGAGGGTTTACCCGt                                                                                                            >  4:45214/1‑134 (MQ=255)
                    tGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAAACCCCCCCCACGGGGGGGTTTCCCgg                                                                                                            >  1:49587/1‑133 (MQ=255)
                      attattGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTTAACCCGCCCCAACGGGGGGGTTACCCCTgca                                                                                                         >  8:304819/1‑133 (MQ=255)
                      attattGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCCCCCCTCCGGGGGGGTTTCCCGTgta                                                                                                         >  8:320099/1‑133 (MQ=255)
                      attattGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTCCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGGTAAACCCCCCCCACCGGGGGGTTTGCCCgggt                                                                                                          >  7:53770/1‑132 (MQ=255)
                             ccaatcaagcGCACTTACACTCATAACCCGCAGGCCAGGGGCTGTACCCTCACCGGTGCAGGTAAATATTTTGGCGGTGGCGGCCGCTCCGCGGGAAAACCCGCCCCTACGGGGGGTTAACCGGGTAAGTGGCTg                                                                                                <  6:377430/133‑1 (MQ=255)
                                               aCTCATTACACGGATGCCCCGGTCTTTTACCGCACCTGTGCCTGGAAATTTTTTTGCGGGGGCGGCCTGTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGaa                                                                              <  8:282776/135‑1 (MQ=38)
                                                                cACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTTAACCCCCCCCCACAGGGGGGTTTACCCGGGACAGGCCTGGACAGGCAACAACCGCACCAAGGGCCAAACcccc                                                             >  6:358868/1‑135 (MQ=255)
                                                                cACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTTAACACCCCCCTACGGGGGGGGTACCCGGGAAGAGGCTGGCCCCGGCAACTACCCAACCAGGGGCACACcccc                                                             >  8:51016/1‑134 (MQ=255)
                                                                                cccggggCCTGCAAATATTTTTGCTGGGGCGGCCGCTTCGCGGGAAAACCCGCCCCTCGGGGGGGTTAACCCGTGAAGTGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTg                                            <  8:126624/129‑1 (MQ=255)
                                                                                aCCTGTGCCTGGAAATATTTTTGCTGGGGCGGCCTCTTTGCGGGAAAACCCGCCCCTCCGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTg                                            <  7:356796/135‑1 (MQ=255)
                                                                                 cgtgtgCCTGCAAATATTTTTGCGGGGGCGGCCTGTTCGCGGGAACACCCGCCCCTCCGGGGGGTTAACCGGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGt                                             <  8:225535/131‑1 (MQ=255)
                                                                                      ccctccaaaTATTTTGGGGGGGGGGGCCTCTCCCGGGAAAACCCCCCCCCACCGGGGGGTTAACCGGGGAAGTGGCGGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAAc                                      <  5:376399/129‑1 (MQ=255)
                                                                                        cTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAACCCCCCCCTACCGGGGGTTTAACCGGGAAGCGGCGGGCACCGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGGGCAAAACCCa                                   >  8:27627/1‑135 (MQ=255)
                                                                                               tatTTTTGCTGGGGGGGCCGCTTCGCGGGTAACCCCGCCCATACGGGGGGTTACCCCGGTAAGTGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCt                             <  8:76059/135‑1 (MQ=255)
                                                                                               tatTTTTGCTGGGGGGGCCGCTGCGCGGGTAAACCCGCCCCTAGGGGGGGTTTACCGGGTAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCt                             <  3:172380/135‑1 (MQ=255)
                                                                                                        tgggggggcCTCTTCGGGGGTAACCCCGCCCCTAGGGGGGGTTAACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCg                    <  2:49587/129‑1 (MQ=255)
                                                                                                                    tgcgcgGGTACACCGGCCCCCTGGGGGGGTGTACCGGGTACGTGGCGGGCACAGGCAACTACCGAACCAGGGTCAACCGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCca        <  7:182369/133‑1 (MQ=255)
                                                                                                                             gccaaCCCCGCCCTACGGGGGGTTTACCCGGGCACTCGGCGTGCCACGGCAACTACCGAAGCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTCCCCAGCGAATCCAGCGCCAGCAGCg  <  8:204195/132‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                  
GACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCG  >  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: