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.

GACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGC  >  ACB122/5223028‑5223291
                                                                                                                                     |                                                                                                                                       
gACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCt                                                                                                                                    <  8:103012/137‑1 (MQ=255)
gACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATAGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGAGGCGGCCTCTTCCCGGGGTAAGcccccccc                                                                                                                                     >  8:8096/1‑137 (MQ=255)
       aCCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTTAAACCCCCCCCACcgggg                                                                                                                          >  3:173281/1‑137 (MQ=255)
       aCCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACggggg                                                                                                                          <  8:196059/137‑1 (MQ=255)
       aCCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACggggg                                                                                                                          <  8:245991/137‑1 (MQ=255)
       aCCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTCTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACggggg                                                                                                                          <  1:400972/137‑1 (MQ=255)
           gtgAGGGTGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCCCCGGGGGGttt                                                                                                                      <  6:343326/137‑1 (MQ=255)
              aGGGGGTGATTATTGCGAATCAATCTCACTTTCACTCATTACACGGATGACACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCAACGGGGGGTTTAcc                                                                                                                   >  8:299700/1‑137 (MQ=255)
              aGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTAAGGGGGGGTTaac                                                                                                                   >  8:305955/1‑135 (MQ=255)
               gggggTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTTAAACCCCCCCAAAGGGGGGGTTTAcc                                                                                                                   >  5:304948/1‑137 (MQ=255)
                    tGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTTAAACCCCCCCCTCCGGGGGGTTTAccgggtt                                                                                                              >  2:403225/1‑132 (MQ=255)
                      attattGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTTAACCCCCCCCAACCGGGGGTTTTAACGTGAagt                                                                                                           >  5:347664/1‑136 (MQ=255)
                      attattGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCCCCCCTACGGGGGGGTTTCCCGTGAagg                                                                                                           >  6:71219/1‑136 (MQ=255)
                         attGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCt                                                                                                                                    <  8:391034/112‑1 (MQ=255)
                         attGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCt                                                                                                                                    >  7:391034/1‑112 (MQ=255)
                              gAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTTAAACCCCCCCCAAGGGGGGGTTTACCCGGAAAAGCCTGGcc                                                                                                   >  8:284390/1‑136 (MQ=255)
                              gAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTAAGGGGGGTTTACCCGGGAAGAGGCTgggc                                                                                                   >  2:275290/1‑135 (MQ=255)
                                atcaatcCCACTTTCCCTCATTACCCGGATGCCACGGGCTGTTACCTCACCTGTCCCTGGAAATATTTTTGCGGGGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTTAAGAGGCTGGcaca                                                                                                 <  7:284961/137‑1 (MQ=38)
                                   aatcGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCCCCCCTACCGGGGGGTTTACCCTGTAGAGGCTGGGACAGGc                                                                                              >  7:329446/1‑137 (MQ=255)
                                        cACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGGAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAAACACCCCCAACGGGGGGTTTAACCGTGAATAGGGTGGCACACGCAACTa                                                                                         >  6:392167/1‑137 (MQ=255)
                                           tttCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGTAAATATTTTTGCTGTGGCGGCCTCTTCGCGGTTAAACCCGCCCCTACGGGGGGTTTACCCGTGAATAGGCTGGCACAGGCAACTACCg                                                                                      <  4:160595/137‑1 (MQ=255)
                                                      acacGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTTAACCCCCCCCTAACGGGGGGTTTACCGTGTAGAGGCTGGCACAGGCAAATACCCAACCAGGGGCa                                                                           >  6:185208/1‑137 (MQ=255)
                                                                       tGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTTAACCCCCCCCAACGGGGGGGTTAAACGTGAAGAGGGTGGGCCACGCAACTAACGAACCCCGGGGAAACCCCCCAACCCgggg                                                          >  5:73098/1‑135 (MQ=255)
                                                                         ttACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTTAAACCCCCCCGAAGGGGGGGTTTCCCGGGTAGAGGCTGGCCCCGGCAACTACCGCACCAGGGGCAAACGCCCCAACACGGtgtg                                                        >  1:60419/1‑137 (MQ=255)
                                                                           aCCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTTAACCCCCCCCAACCGGGGGGTTTCCCCTGAAGAGGCGGGGACCTGGCACAACCGCACCCAGGTGAAACCCCCCAACACGGGGTGcc                                                      >  4:215364/1‑137 (MQ=255)
                                                                             cTCACCTGTGCCTGTAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGt                                                    <  5:185208/137‑1 (MQ=255)
                                                                              tCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCGGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTg                                                   <  5:67636/137‑1 (MQ=255)
                                                                              acACCTGTGCCTGTAAATTTTTTTGCGGTGGGGCCCTCTGCGCGGGTAACCCCGCCCCTCCGGGGGGTTTACCGGTGAAGAGGCGGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTg                                                   <  3:241809/136‑1 (MQ=255)
                                                                                    ttgCCTGGAAATATTTTTGCGGGGGCGGCCTCTTCGCGGGTAAACCCGCCCCTCCGGGGGGTTTACCGGTGAAGTGGCGGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAAc                                             <  6:136044/136‑1 (MQ=255)
                                                                                    gtgCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCCCCCCAACGGGGGGGTTTACCCGGAAGAGGCGGGCACACGCACCAACCCAACCACGGGCAAAACCCCCCACACGGTCTGCCCGGAAAAAc                                             >  4:92666/1‑137 (MQ=255)
                                                                                        cTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTTAAACCGCCCCAAAGGGGGGGTTTACCGGGGAGCGGCTGGCCCAGGCAACGACCCCACCAGGGGGAAACGCCCCCACCCCGTCTGCCGGGCAAAACCCaa                                         >  2:303566/1‑136 (MQ=255)
                                                                                        cTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCAta                                        >  3:43123/1‑136 (MQ=255)
                                                                                         tgtaaaTATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTCCGGGGGGTTTACCCGTGAAGAGGCGGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATc                                        <  6:346024/134‑1 (MQ=255)
                                                                                             aaTATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTCCGGGGGGTTTACCGGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCt                                    <  4:361694/137‑1 (MQ=255)
                                                                                                ttttttGCTGGGGCGGCATCTTCGGGGGTAAACCCGCCCCGAGGGGGGGTTTACCTGGTAATAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGc                                 <  8:127829/136‑1 (MQ=255)
                                                                                                  ttttGCTGTGGCGGCCTCTTCGCGGGTTAACCCGCCCCCAAGGGGGGGTTTACCGGGAAGGGGGGGGCCCATGGAACTACCGCACCAGGGGGAAACGCCCCAACCCCGTCAGGCGGGAAAAAACCAACACCATGccc                               >  2:272893/1‑137 (MQ=255)
                                                                                                         gtgGGGGCCTCTTCGCGGGAAAACCCGCCCCTAGGGGGGGTTTACCCGTGATGTGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAAt                        <  7:135969/137‑1 (MQ=255)
                                                                                                             gggcCTCTGCGCGTGTACACCCGCCCCTCGGGGGGTTTTACCGGTTATGTGGCGGCCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGGCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAg                    <  8:157998/136‑1 (MQ=255)
                                                                                                                                 aaccGCCCCTAACGGGGGGTTACCCGTGAAGAGGCTGGCACACGCAACAACCCAACCACGGTGAAAACCCCCAAAACCGACGGCCCTGAAAAACCCCACACCTAGCCCAGCGAAACAACCGCCACCACCCGCGCGGc  >  4:72290/3‑137 (MQ=255)
                                                                                                                                     |                                                                                                                                       
GACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGC  >  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: