Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1225,384,4070CG30.0% 43.2 / 12.4 30intergenic (+528/‑330)CJ019_05064/CJ019_05067putative protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (5/16);  new base G (9/0);  total (14/16)
Fisher's exact test for biased strand distribution p-value = 1.40e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TGCCCGTCGGGCGGCTGCGGATTAAATCAGTAGACACGGCTAAGCATGTAGTACCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTT  >  ACB122/5384274‑5384532
                                                                                                                                     |                                                                                                                             
tGCCCGTCGGGCGGCTGCGGATTAAATCAGTAGACACGGCTAAGCATGTAGTACCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCgtt                                                                                                                            >  4:129457/1‑137 (MQ=255)
        gggCGGCTGCGGATTAAATCAGTAGACACGGCTAAGCATGTAGTACCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCACCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTGTTGTTTtgttg                                                                                                                    >  8:154451/1‑135 (MQ=255)
         ggcggcTGCGGATTAAATCAGTACACCCGGCTACCCCTGTAGTACCGACAGCGGCGAACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGcc                                                                                                                   <  2:124883/137‑1 (MQ=255)
                                    cGGCTAAGCATGTAGTACCGACAGCGGAGTACTGGCGGACGGGGGGTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGaaaa                                                                                        <  8:46252/137‑1 (MQ=255)
                                      gCTAAGCATGTAGTATCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTGTTTTTGTGCGTGGAATTTTGGGGCGGGGGGTTGaaaaaa                                                                                      >  6:72985/1‑136 (MQ=255)
                                          aGCATGTAGTACCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGa                                                                                  <  5:99721/137‑1 (MQ=255)
                                             aTGTAGTACCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCGTTTTTTTTGGCTGGAATTTTGGGCCGGCGGGTCGAAAAATCGGAttt                                                                               >  5:125360/1‑137 (MQ=255)
                                               ttagtaCCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTaa                                                                             <  3:216026/136‑1 (MQ=255)
                                               gtagtaCCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCTGGTCGAAAAATCGGATTTaa                                                                             <  5:229333/137‑1 (MQ=255)
                                               gtagtaCCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTaa                                                                             <  5:47481/137‑1 (MQ=255)
                                                 tgtaCCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAAcc                                                                           <  5:265317/136‑1 (MQ=255)
                                                 agtaCCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAAcc                                                                           <  4:286931/137‑1 (MQ=255)
                                                    aCCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTGGTTTTTTTTGGCGGGTATTTTGGGGCTGGGGGGGGAAAAAAAGGATTTTAaccca                                                                        >  4:216026/1‑136 (MQ=255)
                                                    aCCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGGGGTTTTTTTTGGGTGGTATTTTGGGGGTGGGGGGGGGAAAATATGATTTaaaaccc                                                                        >  3:143954/1‑132 (MQ=255)
                                                               aGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAggg                                                             <  8:248334/137‑1 (MQ=255)
                                                                 tACTGGCGGAGGGGGGTTCAAATCCCCCCGGCTCCAGCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGGCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCg                                                           <  2:184409/137‑1 (MQ=255)
                                                                    tGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCAGGCGTCGTTTTTTGTTGCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGttt                                                        >  1:13908/1‑137 (MQ=255)
                                                                      gCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTGTTTTTTTTTTGCTGGTATTTGGGGGGGGCGGGGGGAAAAATAGGATTTAAACCCATTGTTGTAGGGCGgttgg                                                      >  8:10228/1‑135 (MQ=255)
                                                                            gggggTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTGGTTTTTTTTGGCGGGAATTTGGGGGCTGCGGGGGGAAAAATCGGATTTTACCCCGTTATTGTTGGGGGGTTCCCAAAGa                                                >  7:243354/1‑137 (MQ=255)
                                                                            cggggTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCGTTTTTTTTTGCGGGGATTTGGGCGCGGCGGGGCGAAAAATAGGATTTAAACCCTTTCTTTGTGGGGGGTTGCCAAAGa                                                >  2:55665/2‑137 (MQ=255)
                                                                              gggTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTGGTTTTTTGTGCCTGGAATTTGGGGGGGGCGGGTGGCAAAATCGGATTTAACCCCGTTGGTGTAGGGCGTTGGCCAAAGATa                                              >  8:303882/1‑137 (MQ=255)
                                                                                     aaTCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATActgcctg                                       <  3:261130/137‑1 (MQ=255)
                                                                                     aaTCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATActgcctg                                       <  3:256316/137‑1 (MQ=255)
                                                                                                      caAATCCAAAGAAGAAGACGCCCTCGGGCGTGGTTTTTTTTGGCGGGGATTTTGGCGCGGGGGGGGGAAAAAAAGGATTTAAACCCGTTATTTTTGGGGGGTTGACAAAAATAAAGACAGCGCCGGTTTTACCTTgg                      >  4:137901/1‑136 (MQ=255)
                                                                                                                taagaagaCGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCAttttt            <  2:9320/136‑1 (MQ=255)
                                                                                                                gaagaagaCGCCCTCGGGCGTGGTTTTTTGTTGCTGGTATTTTGGGGCGGCGGGTGGAAAAAAAGGATTTTAACCCATTCTTGTAGGGCGGTTACCAAATAAAATGCCTGCGCCCGTTTTACGTTGGGGTCAttttt            >  2:73836/1‑137 (MQ=255)
                                                                                                                gaagaagaCGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCAttttt            <  7:164975/137‑1 (MQ=255)
                                                                                                                     agaCGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCATTTTTTCAgg       <  1:199796/137‑1 (MQ=255)
                                                                                                                       aCGCCCTCGGGCGTGTTTTTTTTTGGCTGGAATTTTGGCGCTGCGGGGGGAAAAAACGGCTTTTAACCCTTTCTTGTAGGGCGTTTGCCAAAGATAATGCCTGCGGCGGTTTTCCGTTGCGGGGAGTTTTTTAgggg     >  8:303929/1‑136 (MQ=255)
                                                                                                                          cccTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCATTTTTTCAGGGtttt  <  7:37289/137‑1 (MQ=255)
                                                                                                                                     |                                                                                                                             
TGCCCGTCGGGCGGCTGCGGATTAAATCAGTAGACACGGCTAAGCATGTAGTACCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTT  >  ACB122/5384274‑5384532

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