Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1225,223,1590GC34.2% 42.6 / 12.3 38intergenic (‑133/+43)cntO_7/CJ019_04917Metal‑pseudopaline receptor CntO/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (8/17);  new base C (13/0);  total (21/17)
Fisher's exact test for biased strand distribution p-value = 6.20e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGC‑GGG‑TAAACCCGCCCCT‑A‑C‑GGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGC  >  ACB122/5223028‑5223291
                                                                                                                                     |                                                                                                                                       
gACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGC‑GGG‑TAAACCCGCCCCt                                                                                                                                    <  8:103012/137‑1 (MQ=255)
gACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATAGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGAGGCGGCCTCTTCCCGGGG‑TAAGcccccccc                                                                                                                                     >  8:8096/1‑137 (MQ=255)
       aCCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGC‑GGG‑TTAAACCCCCCCC‑A‑C‑cgggg                                                                                                                          >  3:173281/1‑137 (MQ=255)
       aCCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGC‑GGG‑TAAACCCGCCCCT‑A‑C‑ggggg                                                                                                                          <  8:196059/137‑1 (MQ=255)
       aCCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGC‑GGG‑TAAACCCGCCCCT‑A‑C‑ggggg                                                                                                                          <  8:245991/137‑1 (MQ=255)
       aCCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTCTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGC‑GGG‑TAAACCCGCCCCT‑A‑C‑ggggg                                                                                                                          <  1:400972/137‑1 (MQ=255)
           gtgAGGGTGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGC‑GGG‑TAAACCCGCCCCC‑C‑C‑GGGGGGttt                                                                                                                      <  6:343326/137‑1 (MQ=255)
              aGGGGGTGATTATTGCGAATCAATCTCACTTTCACTCATTACACGGATGACACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGC‑GGG‑TAAACCCGCCCCA‑A‑C‑GGGGGGTTTAcc                                                                                                                   >  8:299700/1‑137 (MQ=255)
              aGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGC‑GGG‑TAAACCCGCCCCT‑A‑A‑GGGGGGGTTaac                                                                                                                   >  8:305955/1‑135 (MQ=255)
               gggggTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGC‑GGGTTAAACCCCCCCAA‑A‑G‑GGGGGGTTTAcc                                                                                                                   >  5:304948/1‑137 (MQ=255)
                    tGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGC‑GGGTTAAACCCCCCCCT‑C‑C‑GGGGGGTTTAccgggtt                                                                                                              >  2:403225/1‑132 (MQ=255)
                      attattGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGC‑GGG‑TTAACCCCCCCCA‑A‑C‑CGGGGGTTTTAACGTGAagt                                                                                                           >  5:347664/1‑136 (MQ=255)
                      attattGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGC‑GGG‑TAAACCCCCCCCT‑A‑C‑GGGGGGGTTTCCCGTGAagg                                                                                                           >  6:71219/1‑136 (MQ=255)
                         attGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGC‑GGG‑TAAACCCGCCCCt                                                                                                                                    <  8:391034/112‑1 (MQ=255)
                         attGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGC‑GGG‑TAAACCCGCCCCt                                                                                                                                    >  7:391034/1‑112 (MQ=255)
                              gAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGC‑GGGTTAAACCCCCCCCA‑A‑G‑GGGGGGTTTACCCG‑GAAAAGCCTGGcc                                                                                                   >  8:284390/1‑136 (MQ=255)
                              gAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGC‑GGG‑TAAACCCGCCCCT‑A‑A‑GGGGGGTTTACCCGGGAAGAGGCTgggc                                                                                                   >  2:275290/1‑135 (MQ=255)
                                atcaatcCCACTTTCCCTCATTACCCGGATGCCACGGGCTGTTACCTCACCTGTCCCTGGAAATATTTTTGCGGGGGCGGCCTCTTCGC‑GGG‑TAAACCCGCCCCT‑A‑C‑GGGGGGTTTACCCGTTAAGAGGCTGGcaca                                                                                                 <  7:284961/137‑1 (MQ=38)
                                   aatcGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGC‑GGG‑TAAACCCCCCCCT‑ACC‑GGGGGGTTTACCC‑TGTAGAGGCTGGGACAGGc                                                                                              >  7:329446/1‑137 (MQ=255)
                                        cACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGGAATATTTTTGCTGTGGCGGCCTCTTCGC‑GGG‑TAAAACACCCCCA‑A‑C‑GGGGGGTTTAACCGTGAATAGGGTGGCACACGCAACTa                                                                                         >  6:392167/1‑137 (MQ=255)
                                           tttCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGTAAATATTTTTGCTGTGGCGGCCTCTTCGC‑GGT‑TAAACCCGCCCCT‑A‑C‑GGGGGGTTTACCCGTGAATAGGCTGGCACAGGCAACTACCg                                                                                      <  4:160595/137‑1 (MQ=255)
                                                      acacGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGC‑GGG‑TTAACCCCCCCCTAA‑C‑GGGGGGTTTA‑CCGTGTAGAGGCTGGCACAGGCAAATACCCAACCAGGGGCa                                                                           >  6:185208/1‑137 (MQ=255)
                                                                       tGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGC‑GGG‑TTAACCCCCCCCA‑A‑C‑GGGGGGGTTAAACGTGAAGAGGGTGGGCCACGCAACTAACGAACCCCGGGGAAACCCCCCAACCCgggg                                                          >  5:73098/1‑135 (MQ=255)
                                                                         ttACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGC‑GGG‑TTAAACCCCCCCG‑A‑A‑GGGGGGGTTTCCCGGGTAGAGGCTGGCCCCGGCAACTACCGCACCAGGGGCAAACGCCCCAACACGGtgtg                                                        >  1:60419/1‑137 (MQ=255)
                                                                           aCCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGC‑GGG‑TTAACCCCCCCCA‑A‑C‑CGGGGGGTTTCCCCTGAAGAGGCGGGGACCTGGCACAACCGCACCCAGGTGAAACCCCCCAACACGGGGTGcc                                                      >  4:215364/1‑137 (MQ=255)
                                                                             cTCACCTGTGCCTGTAAATATTTTTGCTGTGGCGGCCTCTTCGC‑GGG‑TAAACCCGCCCCT‑A‑C‑GGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGt                                                    <  5:185208/137‑1 (MQ=255)
                                                                              tCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGC‑GGG‑TAAACCCGCCCCT‑A‑C‑GGGGGGTTTACCCGTGAAGAGGCGGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTg                                                   <  5:67636/137‑1 (MQ=255)
                                                                              acACCTGTGCCTGTAAATTTTTTTGCGGTGGGGCCCTCTGCGC‑GGG‑TAACCCCGCCCCT‑C‑C‑GGGGGGTTTACCGGTGAAGAGGCGGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTg                                                   <  3:241809/136‑1 (MQ=255)
                                                                                    ttgCCTGGAAATATTTTTGCGGGGGCGGCCTCTTCGC‑GGG‑TAAACCCGCCCCT‑C‑C‑GGGGGGTTTACCGGTGAAGTGGCGGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAAc                                             <  6:136044/136‑1 (MQ=255)
                                                                                    gtgCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGC‑GGG‑TAAACCCCCCCCA‑A‑CGGGGGGGTTTACCCG‑GAAGAGGCGGGCACACGCACCAACCCAACCACGGGCAAAACCCCCCACACGGTCTGCCCGGAAAAAc                                             >  4:92666/1‑137 (MQ=255)
                                                                                        cTGGAAATATTTTTGCTGTGGCGGCCTCTTCGC‑GGG‑TTAAACCGCCCCA‑A‑A‑GGGGGGGTTTACCGGGGAGCGGCTGGCCCAGGCAACGACCCCACCAGGGGGAAACGCCCCCACCCCGTCTGCCGGGCAAAACCCaa                                         >  2:303566/1‑136 (MQ=255)
                                                                                        cTGGAAATATTTTTGCTGTGGCGGCCTCTTCGC‑GGG‑TAAACCCGCCCCT‑A‑C‑GGGGGG‑TTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCAta                                        >  3:43123/1‑136 (MQ=255)
                                                                                         tgtaaaTATTTTTGCTGTGGCGGCCTCTTCGC‑GGG‑TAAACCCGCCCCT‑C‑C‑GGGGGGTTTACCCGTGAAGAGGCGGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATc                                        <  6:346024/134‑1 (MQ=255)
                                                                                             aaTATTTTTGCTGTGGCGGCCTCTTCGC‑GGG‑TAAACCCGCCCCT‑C‑C‑GGGGGGTTTACCGGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCt                                    <  4:361694/137‑1 (MQ=255)
                                                                                                ttttttGCTGGGGCGGCATCTTCGG‑GGG‑TAAACCCGCCCCG‑A‑G‑GGGGGGTTTACCTGGTAATAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGc                                 <  8:127829/136‑1 (MQ=255)
                                                                                                  ttttGCTGTGGCGGCCTCTTCGC‑GGG‑TTAACCCGCCCCC‑A‑A‑GGGGGGGTTTACCGGGAAGGGGGGGGCCCATGGAACTACCGCACCAGGGGGAAACGCCCCAACCCCGTCAGGCGGGAAAAAACCAACACCATGccc                               >  2:272893/1‑137 (MQ=255)
                                                                                                         gtgGGGGCCTCTTCGC‑GGG‑AAAACCCGCCCCT‑A‑G‑GGGGGGTTTACCCGTGATGTGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAAt                        <  7:135969/137‑1 (MQ=255)
                                                                                                             gggcCTCTGCGC‑GTG‑TACACCCGCCCCT‑C‑G‑GGGGGTTTTACCGGTTATGTGGCGGCCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGGCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAg                    <  8:157998/136‑1 (MQ=255)
                                                                                                                                 aaccGCCCCT‑A‑A‑CGGGGGGTTACCCGTGAAGAGGCTGGCACACGCAACAACCCAACCACGGTGAAAACCCCCAAAACCGACGGCCCTGAAAAACCCCACACCTAGCCCAGCGAAACAACCGCCACCACCCGCGCGGc  >  4:72290/3‑137 (MQ=255)
                                                                                                                                     |                                                                                                                                       
GACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGC‑GGG‑TAAACCCGCCCCT‑A‑C‑GGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGC  >  ACB122/5223028‑5223291

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