Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1435,383,7580CG50.0% 22.7 / 15.8 27intergenic (+528/‑330)CJ019_05064/CJ019_05067putative protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/13);  new base G (13/0);  total (13/14)
Fisher's exact test for biased strand distribution p-value = 1.92e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ACACGGCTAAGCATGTAGTACCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTC  >  ACB143/5383658‑5383889
                                                                                                     |                                                                                                                                   
acacGGCTAAGCATGTAGTACCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTGTTGTTTTGGGGCTGGAATTTGGGCGCGGCGGGTCg                                                                                                  >  8:183453/1‑136 (MQ=255)
        aaGCATGTAGTACCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTGGGTTTTTTTTGGGGGGATTTGGTGGCGGCGGGGCGaaaaaaag                                                                                          >  2:27665/1‑133 (MQ=255)
            aTGTCGTACCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTAGGCGCTGCGGGTCGAAAAATCGGAtt                                                                                      <  1:318447/136‑1 (MQ=255)
                 gtccccaCAGCGGAGTACTGGCGGACGTGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGCCGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAAcc                                                                                 <  8:139341/130‑1 (MQ=255)
                 gtaCCTACCGCGGAGTCCGGGCGGACGGGGGTTCAAAACCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAAcc                                                                                 <  6:211110/136‑1 (MQ=255)
                 gtaCCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTGGTTTTTTTTTGCTGGAATTTTGGGGCTGGGGGGGGAAAAATAGGGTTTAAcc                                                                                 >  8:16093/1‑136 (MQ=255)
                 gtaCCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTGGTTTTTTGTGGCTGGGATTTGGGGGCGGCGGGGGGAAAAAAAGGTTTTAAca                                                                                 >  6:337628/1‑135 (MQ=255)
                 gtaCCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAAcc                                                                                 <  8:103552/136‑1 (MQ=255)
                 gtaCCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAAcc                                                                                 <  2:28022/136‑1 (MQ=255)
                 ggaCCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAACAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAAcc                                                                                 <  4:205289/134‑1 (MQ=255)
                   aCCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTGGTTTTTTTTTGGGGGGATTTTGGGGGTGCGGGTCGGAAAATAGGGTTTaaacac                                                                               >  1:2637/1‑132 (MQ=255)
                   aCCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTGTTTTTTTTTGCCTGGAATTTTGGCGCTGCGGGTCGAAAAATAGGATTTAAaccc                                                                               >  7:22014/1‑136 (MQ=255)
                         aGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCACAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCtt                                                                         <  3:255991/136‑1 (MQ=255)
                          gCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTGGTTTTTTTTGGCGGGGATTTTTGGGCGGGGGGGGGCAAAATCGGATTTAACACCGTTCTTg                                                                        >  6:342961/1‑136 (MQ=255)
                          gCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTGGTTTTTTTTGGCGGGAATTTGGGGGCTGCGGGTGGAAAAATAGGATTTTACCCCGTTCTTg                                                                        >  1:235610/1‑136 (MQ=255)
                              aGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCGCGGGCGTGGTTTTTTTTTGCTGGGATTTTGGGGGGGGGGGGGGGAAAAAAGGAGTTAAACCCGTtgttgtttt                                                                    >  6:112486/1‑133 (MQ=255)
                               gTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAggg                                                                   <  6:343800/136‑1 (MQ=255)
                                                ttCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTGGTTTTTTTTGGCTGGGATTTTGGGGGGGGGGGGGGGAAAATAGGTATTAAACCCGGTGTTTTTGGGCGTTTGCCaaaaaaaat                                                  >  2:28280/1‑130 (MQ=255)
                                                     aTCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTGGTTTTTTGTGGCGGGGATTTGGGGGGGGCGGGTCGAAAAATAGGGTTTAAACCCTTTCTTGTATGGCGGTTGCCAAAGATActgcctg                                             >  1:32192/1‑136 (MQ=255)
                                                      tCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTGGGTTTTTTTTGGGGGCATTTTGGGGGTGCGGGGGGGAAAATACGATTTAAACACAATGATGTAGGGCGTGTGCCCAAAATAatggaaaa                                            >  2:97245/1‑131 (MQ=255)
                                                                            aaaGAAGAAGACGCCCTCGGGCGTGGTTTTTTTTGGGGGGGATTTTGGGGCGGGGGGGGGTAAAATAGGATTTTAACCCGTTGTTGTTGGGGGTTTTCGAAATATAATTCCTGCGCCGGTTTTAAGTTGGGGTGAt                      >  4:173363/1‑136 (MQ=255)
                                                                                aagaagaCGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCAttttt                  <  3:187884/136‑1 (MQ=255)
                                                                                aagaagaCGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCAttttt                  <  1:28280/136‑1 (MQ=255)
                                                                                 agatgaCGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCAttttt                  <  2:142702/135‑1 (MQ=255)
                                                                                        gCCCTCGGGCGTCTTTTTTTGTGCCTTGACTTTGGGCGCTGCTGTTCGAAAAATCTGATTTAACCCCTTTATTGTAGGGCGTTTCACAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCATTTTTTCAGGGtt          <  8:310995/136‑1 (MQ=255)
                                                                                                tcGTATTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGTTACTGCCTGCGCCTGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGtctc  <  8:261056/135‑1 (MQ=255)
                                                                                                gCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGTGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGGTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGtctc  <  8:261053/136‑1 (MQ=255)
                                                                                                     |                                                                                                                                   
ACACGGCTAAGCATGTAGTACCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTC  >  ACB143/5383658‑5383889

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: