Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1435,096,4150CG32.6% 58.6 / 10.6 47intergenic (+28/‑162)CJ019_04792/CJ019_04793Putative pterin‑4‑alpha‑carbinolaminedehydratase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (13/18);  new base G (12/3);  total (25/22)
Fisher's exact test for biased strand distribution p-value = 2.60e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCC‑TTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGC  >  ACB143/5096282‑5096540
                                                                                                                                     |                                                                                                                              
cACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCG‑tt                                                                                                                             >  5:150413/1‑136 (MQ=255)
cACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCC‑tt                                                                                                                             >  1:245490/1‑136 (MQ=255)
   cgggACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCC‑ttttt                                                                                                                          <  8:309359/133‑1 (MQ=255)
   cgggACCTGGGGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCC‑ttttt                                                                                                                          <  2:231267/133‑1 (MQ=255)
   cGTGACGTGGTGGAGCCACTCGATCAAAGGCCTGCCCCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGACGGCCGTAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCC‑ttttt                                                                                                                          <  6:192931/136‑1 (MQ=255)
   cGTGACCTGGTGGAGCCAGTAGATCAAAGGCATGCACCGCAACGAATTGATCATTTTTTCACGCACTGACAAGGGGGCGGAAACGGCTGAAGGCAGGATGTAAACACAAGAAGGAAGTTGCCAGAGGCCCG‑ttttt                                                                                                                          <  2:92887/136‑1 (MQ=255)
   cGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCTCCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGCAAGGGGCCAGAGGCCCC‑ttttt                                                                                                                          <  1:302926/136‑1 (MQ=255)
   cGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCC‑ttttt                                                                                                                          <  8:297957/136‑1 (MQ=255)
   cGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCC‑ttttt                                                                                                                          <  6:230510/136‑1 (MQ=255)
   cGTGACCGGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCTTGTGTGCACGCACTTACAAGGTGGCGGAAACGGCTGAAGGCCGTAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCC‑ttttt                                                                                                                          <  8:180581/136‑1 (MQ=255)
     tGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCC‑tttttt                                                                                                                         >  7:149680/1‑135 (MQ=255)
     tGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCC‑TTTTTTc                                                                                                                        <  4:203198/136‑1 (MQ=255)
     tGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCC‑TTTTTTc                                                                                                                        <  6:100482/136‑1 (MQ=255)
      gACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCC‑GTTTTTCa                                                                                                                       >  5:114556/1‑136 (MQ=255)
          tggtggAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCATTTTTTCATgg                                                                                                                    >  1:9846/1‑136 (MQ=255)
            gtggAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCC‑TTTTTTCATggg                                                                                                                   <  7:302371/134‑1 (MQ=255)
            gtggAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCC‑GTTGTTAATGGGGc                                                                                                                 >  1:327751/1‑136 (MQ=255)
               gAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGCGGCCCG‑GTTTTTCATGGGGGata                                                                                                              >  3:89663/1‑136 (MQ=255)
                aGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCC‑TTTTTTCATGGGGCatat                                                                                                             <  4:346915/136‑1 (MQ=255)
                    aCTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCC‑TTTTTTCATGGGGCATATGTgc                                                                                                         >  1:2568/1‑136 (MQ=255)
                         aTCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCG‑GTTTTTCATGGGGCATATGTGCGCTGt                                                                                                    >  3:211665/1‑136 (MQ=255)
                                       ccGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCG‑GGTTTTATTGGGGGATATTTTGGCTGTTGGCCCGCCAAAcc                                                                                      >  6:40263/1‑136 (MQ=255)
                                       ccGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCC‑GTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAAcc                                                                                      >  1:276247/1‑136 (MQ=255)
                                             cGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCG‑GTTTTTGATGGGGCATATTTTCGCTGTTCGCCCGCAAAACCGACaaa                                                                                >  7:294880/1‑136 (MQ=255)
                                             cGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCG‑GTTTTTCATGGGGCATATTTGCGGTGTTAGGCCGCAAAACCGACaa                                                                                 >  1:258057/1‑135 (MQ=255)
                                             cGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCC‑TTTTTTCATGGGGCATATGTGCGCTGTTAGCCCGCAAAACCCACaa                                                                                 >  4:327260/1‑135 (MQ=255)
                                             cGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCC‑GTTTTTGATTGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACaaa                                                                                >  1:210073/1‑136 (MQ=255)
                                                       atGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAATTAAACACCAGAAGGAAGGGGCCAGAGGCCCC‑TTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGa                                                                       <  3:266569/135‑1 (MQ=255)
                                                            tgCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCG‑GTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGCCAAACCCCCGGGcaaaaa                                                                  >  4:96042/1‑135 (MQ=255)
                                                             gcacgcacTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCC‑GTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAAc                                                                >  3:66010/1‑136 (MQ=255)
                                                                 gcacTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCG‑GTTTTTAATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCaa                                                             >  4:68756/1‑135 (MQ=255)
                                                                 gcacTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACAAGAAGGAAGGGGCCAGAGGCCCG‑TATTTTCAAGGGGCTATTGTGCGCTGTTCGCCCGCAAAACCGACAAACCCCCGGGAAAAAAACCaaa                                                            >  2:344169/1‑136 (MQ=255)
                                                                      gACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCG‑GTTTTTAATGGGGCATATTTGAGCTGTTAGCCCGCAAAACCGACAAACGACCGGGAAAAAAACAAAACTGt                                                        >  4:189502/1‑135 (MQ=255)
                                                                       aCATGGGGGCGTAAAGGGCTGACGGCCGTATGTAAACACCTGACGGACGGGGCCCTAGGCCCT‑TTTTTTCATGGGGAATTTGTGCGCTGTGCGCCCGAAAACCCCACAAACGACCGGGAAAAAAACCAAACTGTc                                                       <  7:291567/135‑1 (MQ=255)
                                                                          tggTGGCGGAAACGGCTGATGGCGGTATGAAAACACCAGACGGACGGGGCCAGCGGCCCC‑TTTTTTCATGGGGAATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATcc                                                   <  3:327260/135‑1 (MQ=255)
                                                                          tggTGCCGGAAACGGCTGAAGGCCGTAAGAAAACACCAGAAGGAAGGGCCCAGAGGCCCG‑TTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATcc                                                   <  5:176980/135‑1 (MQ=255)
                                                                          aGGTGGCGGAAACGGCTGAAGGCCGGATGAAAACACAACAAGCAAGGGGCAAGAGGCCCC‑TTTTTTCATGGGGAATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATcc                                                   <  4:256631/136‑1 (MQ=255)
                                                                          aGGTGGCGGAAACGCGTGAAGGCCGTAAGTAAACACCAGAAGGACGGGGCCCGAGGCCCC‑TTTTTTCATGGGGCATATGTGCGCGGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATcc                                                   <  3:211235/136‑1 (MQ=255)
                                                                          aGGGGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCC‑TTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATcc                                                   <  6:321503/136‑1 (MQ=255)
                                                                              ggCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCG‑GTTTTTTATGGGGGATATGTGGGCTGTTAGCCCGCAAAACCGCCCAAAGACCGGGGAAAAAACCAAACATTCATCCAGa                                                >  2:323719/1‑135 (MQ=255)
                                                                                cGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCG‑TTTTTTAATTGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCCGGAAAAAAACCAAACAGTCATCCAGAcgt                                             >  6:169888/1‑134 (MQ=255)
                                                                                                               aGAAGGAAGGGGCCAGAGGCCCC‑TTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACa              <  2:245490/136‑1 (MQ=255)
                                                                                                                    tacgggcCCACAGGCCCG‑TTTTTTCATGGGGCATTTGGGCGGGGTCCGCCCCAAAACCCGAAAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAAcgc         <  1:309511/133‑1 (MQ=255)
                                                                                                                    gaaggGGCCAGAGGCCCC‑TTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGCCCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCGGCCCGCGCTTTTTGAAGCTACAAAcgc         >  4:79616/1‑136 (MQ=255)
                                                                                                                         ggCCAGAGGCCCC‑GTTTTTGATGGGGCATTTGTTGGCTGTTCGCCCGCAAAACCGACAAAAGACCGGGAAAAAAACCAAACTGTCATCCAGGCTTGGGTATCCTGCCCGTGCTTTTTGAAGCTACAAACGCGCGTg    >  6:129961/1‑136 (MQ=255)
                                                                                                                           cccgcggCCCC‑TTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCCACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGgc  <  1:2489/131‑1 (MQ=255)
                                                                                                                           ccAGAGGCCGC‑GTTTTTAAAGGTGGATATTTGCGGTGTTAGCACGCAAAACAGCCAACCGGCCGCGAAAAAAAACAAACAGGCAACCAGAAATCTGTTTCCGGACACCCCTTTTTGAAGCTACAAACGCCCCAGgc  >  6:122011/1‑136 (MQ=255)
                                                                                                                                     |                                                                                                                              
CACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCC‑TTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGC  >  ACB143/5096282‑5096540

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