Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1433,989,9750GC33.3% 50.4 / 12.3 36intergenic (+25/+20)CJ019_03768/CJ019_03769hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (1/23);  new base C (12/0);  total (13/23)
Fisher's exact test for biased strand distribution p-value = 1.04e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CCATGACGCGCTCAAGGAACTGCGCGTACTGTCCACCGTGAAGTTCTACGGCTGCCTGATGGTGGACGGGCTGGGTCCTGCCATGCGTTTCGACATCCTGCCCAAATGACCCTCGACGGCCTCAGGCGCCCCCGCGCCCGAG‑GCCG‑CAAACC‑GTCATTTTTCCTGCAACACCGCCCTGCCCTTCACCGCCCTCGGACCGCGCCAGGCCCAGAACAGCAAACCGGGGAAGGGTGCGTAGACCACGAAGACGATCCACAGCATCTTGCGCT  >  ACB143/3989842‑3990110
                                                                                                                                     |                                                                                                                                          
caaTGACGCGCTAAAGTAACTGCGCTTAGTGTCCACCGTTAAGTTCTACGGCTGCCTGAGGGGGGAGGGGCTGGGTCCTGCCATGCGTTTCGAAATCCTGCCCAAATGACCCTCGACGGCCTCCGGCGCCCCCgcgc                                                                                                                                         <  4:191026/135‑1 (MQ=255)
     aCGCGCTCAAGGAACTGCGCGTACTGTCCACCGTGAAGTTCTACGGCTGCCTGATGGTGGACGGGCTGGGTCCTGCCATGCGTTTCGACATCCTGCCCAAATGACCCTCGACGGCCGCAGGCGCccccccccccccc                                                                                                                                    >  8:191880/1‑134 (MQ=255)
             aaGGAACTGCGCGTACTGTCCACCGTGAAGTTCTACGGCTGCCTGATGGTGGACGGGCTGGGTCCTGCCATGCGTTTCGACATCCTGCCCAAATGACCCTCGACGGCCTCCGGCGCCCCCCCGCCCGCG‑GGCG‑Caaa                                                                                                                          >  4:164821/1‑137 (MQ=255)
             aaGGAACTGCGCGTACTGTCCACCGTGAAGTTCTACGGCTGCCTGATGGTGGACGGGCTGGGTCCTGCCATGCGTTTCGACATCCTGCCCAAATGACCCTCGACCGCCCGCGGCCCCCCCCCGCCCCCC‑GCCG‑Caaa                                                                                                                          >  4:220434/1‑137 (MQ=255)
                    tGCGCGTACTGTCCACCGTGAAGTTCTACGGCTGCCTGATGGTGGACGGGCTGGGTCCTGCCATGCGTTTCTACATCCTGCCCAAATGACCCTCGACGGCCTCAGGCGCCCCCGCGCCCGAG‑GCCG‑CAAACC‑GTCAt                                                                                                                  <  7:54092/137‑1 (MQ=255)
                    tGCGCGTACTGTCCACCGTGAAGTTCTACGGCTGCCTGATGGGGGAGGGGCTGGGGCCTGCCATGCGTTTCGACATCCTGCCCAAATGACCCTCGACGGCCTCAGGCGCCCCCGCGCCCGAG‑GCCG‑CAAACC‑GTCAt                                                                                                                  <  1:303236/137‑1 (MQ=255)
                    tGCGCGTACTGTCAACCGTGAAGTTCTAGGGCTGCCTTATGGTGGACGGGCTGGGTCCTGCCATGCGTTTCTACATCCTCCCCAATTGACCCTCGACGGCCTCAGGCGCCCCCGCGCCCGAG‑GCCG‑CAAACC‑GTCAt                                                                                                                  <  5:327914/137‑1 (MQ=255)
                    ggcggGTACTGGCCACCGTGAAGTGCTCGGGGTGCCTGATGGTGGACGGGGTGGGTCCTCCCATGGTTTTCGACATCCTGCCCAAATGACCCTCGACGGCCTCAGGCGCCCCCGCGCCCGAG‑GCCG‑CAAACC‑GTCAt                                                                                                                  <  6:58104/136‑1 (MQ=255)
                    ggcgcGTACTGTCCACCGTGAAGTTCTCCGGCGGCCTGATGGGGGACGGGCTGGGTCCTGCCATGCGTTTCCACATCCTGCCCAAATGACCCGCGACGGCCTCCGGCGCCCCCGCGCCCGAG‑GCCG‑CAAACC‑GTCAt                                                                                                                  <  8:6713/136‑1 (MQ=255)
                    ggcccgTACTGTCCACCGTGAATGTCATCGGCTGGCCGCTGGTGGACGGGCGGGGTCCTGCCATGCGTTTCGACATCCTGCCCAAATGACCCTCGACGGCCTCAGGCGCCCCCGCGCCCGAG‑GCCG‑CAAACC‑GTCAt                                                                                                                  <  7:373560/133‑1 (MQ=255)
                       gcgTACTGTCCACCGTGAAGTTCTACGGCTGCCTGATGGTGGACGGGCTGGGTCCTGCCATGCGTTTCGACATCCTGCCCAAATGACCCTCGACGGCCTCAGGCGCCCCCGCGCCCGAG‑GCCG‑CAAACC‑GTCAtttt                                                                                                               <  7:359633/137‑1 (MQ=255)
                       gcgTACTGTCCACCGTGAAGTTCTACGGCTGCCTGATGGTGGACGGGCTGGGTCCTGCCATGCGTTTCGACATCCTGCCCAAATGACCCTCGACCGCCACACGCGCCCCCCCCCCCCAGCGCGC‑CAAACC‑GCAAttt                                                                                                                >  5:135152/1‑137 (MQ=255)
                        cgTACTGTCCACCGTGAAGTTCTACGGCTGCCTGATGGTGGACGGGCTGGGTCCTGCCATGCGTTTCGACATCCTGCCCAAATGACCCTCGACGGCCTCCCGCGCCCCCCCGCCCCCG‑GCCG‑CCAACC‑GGCAttttt                                                                                                              >  6:160169/1‑137 (MQ=255)
                        cgTACTGTCCACCGTGAAGTTCTACGGCTGCCTGATGGTGGACGGGCTGGGTCCCGCCATGCGTTTCGACATCCTGCCCAAATGACCCTCGAGGGCCTCCCGCGCCCCCCCCCCCGGG‑GCCC‑GCAACC‑GGCCttttt                                                                                                              >  7:254301/1‑137 (MQ=255)
                         gTACTGTCCACCGTGAAGTTCTACGGCTGCCTGATGGTGGACGGGCTGGGTCCTGCCATGCGTTTCGACATCCTGCCCAAATGACCCTCGACGGCCTCAGGCGCCCCCGCGCCCGAG‑GCCG‑CAAACC‑GTCATTTTTc                                                                                                             <  5:61218/137‑1 (MQ=255)
                          tACTGTCCACCGTGAAGTTCTACGGCTGCCTGATGGTGGACGGGCTGGGTCCTGCCATGCGTTTCGACATCCTGCCCAAATGACCCTCGACCGCCGCAGGCCCCCCCCCGCCCGAG‑GCCC‑CCAACC‑CGCAattttga                                                                                                            >  4:180150/1‑135 (MQ=255)
                              gTCCACCGTGAAGTTCTACGGCTGCCTGATGGTGGACGGGCTGGGTCCTGCCATGCGTTTCGACATCCTGCCCAAATGACCCTCGACGGCCTCAGGCGCCCCCGCGCCCGAG‑GCCG‑CAAACC‑GTCATTTTTCCTGCa                                                                                                        <  5:189328/137‑1 (MQ=255)
                                 cACCGTGAAGTTCTACGGCTGCCTGATGGTGGACGGGCTGGGTCCTGCCATGCGTTTCGACATCCTGCCCAAATGACCCTCGACGGCCTCCCGCGCCCCCCCGCCCCCG‑GCCC‑CCAACC‑CTCATTTTTACTGCAaca                                                                                                     >  8:359633/1‑137 (MQ=255)
                                       gAAGTTCTACGGCTGCCTGATGGTGGACGGGCTGGGTCCTGCCATGCGTTTCGACATCCTGCCCAAATGACCCTCGACCGCCACACCCGCCCCCCCCCCCCCG‑CCCG‑CCCACC‑GCCAATTTTAATACAACAACcccc                                                                                               >  7:3352/1‑137 (MQ=255)
                                         aGTTCTACGGCTGCCTGATGGTGGACGGGCTGGGTCCTGCCATGCGTTTCGACATCCTGCCCAAATGACCCTCGACGGCCTCAGGCGCCCCCGCGCCCGAG‑GCCG‑CAAACC‑CTCATTTTTCCTGCCACAccccccgg                                                                                             >  2:21202/1‑135 (MQ=255)
                                            tCTACGGCTGCCTGATGGTGGACGGGCTGGGTCCTGCCATGCGTTTCGACATCCTGCCCAAATGACCCTCGACGGCCTCCGGCGCCCCCCCCCCCCCG‑GCCG‑CCAACC‑CGGAATTTTACTGCAACCCCccccagccc                                                                                          >  5:325309/1‑137 (MQ=255)
                                              aaCGGCTGCCTGATGGTGGACGGGCTGGGTCCTGCCATGCGTTTCGACATCCTGCCCAAATGACCCTCGACGGCCTCAGGCGCCCCCGCGCCCGAG‑GCCG‑CAAACC‑GTCATTTTTCCTGCAACACCGCCCTGCCCtt                                                                                        <  5:330045/136‑1 (MQ=255)
                                                        tGATGGTGGACGGGCTGGGTCCTGCCATGCGTTTCGACATCCTGCCCAAATGACCCTCGAGCGCCCCAGGCCCCCCCCCCCCCGAC‑GGCG‑CACACC‑GGCCTTTTTACTGAAAAAACCCCCCTCCCGTCACCAcccac                                                                              >  1:71768/1‑135 (MQ=255)
                                                               ggACGGGCTGGGTCCTGCCATGCGTTTCGACATCCTGCCCAAATGACCCTCGACGGCCTCAGGCGCCCCCGCGCCCGAG‑GCCG‑CAAACC‑GTCATTTTTCCTGCAACACCGCCCTGCCCTTCACCGCCCTCGGACcgc                                                                       <  7:191880/137‑1 (MQ=255)
                                                               ggACGGGCTGGGTCCTGCCATGCGTTTCGACATCCTGCCCAAATGACCCTCGACGGCCTCAGGCGCCCCCGCGCCCGAG‑GCCG‑CAAACC‑GTCATTTTTCCTGCAACACCGCCCTGCCCTTCACCGCCCTCGGACcgc                                                                       <  4:96322/137‑1 (MQ=255)
                                                               ggACGGGCTGGGTCCTGCCATGCGTTTCGACATCCTGCCCAAATGACCCTCGACGGCCTCAGGCGCCCCCGCGCCCGAG‑GCCG‑CAAACC‑GTCATTTTTCCTGCAACACCGCCCTGCCCTTCACCGCCCTCGGACcgc                                                                       <  8:212914/137‑1 (MQ=255)
                                                                  cGGGCTGGGTCCTGCCATGCGTTTCGACATCCTGCCCAAATGACCCTCGACGGCCTCAGGCGCCCCCGCGCCCGAG‑GCCG‑CAAACC‑GTCATTTTTCCTGCAACACCGCCCTGCCCTTCACCGCCCTCGGACCGCGcc                                                                    <  6:135629/137‑1 (MQ=255)
                                                                            ccTGCCATGCGTTTCGACATCCTGCCCAAATGACCCTCGACGGCCTCAGGCGCCCCCGCGCCCGAG‑GCCG‑CAAACC‑GTCATTTTTCCTGCAACACCGCCCTGCCCTTCACCGCCCTCGGACCGCGCCAGGCCCAGaa                                                          <  6:135152/137‑1 (MQ=255)
                                                                              tGCCATGCGTTTCGACATCCTGCCCAAATGACCCTCGCCGGCCTCAGGCGCCCCCGCGCCCGAG‑GCCG‑CAAACC‑GTCATTTTTCCTGCAACACCGCCCTGCCCTTCACCGCCCTCGGACCGCGCCAGGCCCAGAACa                                                        <  5:160169/137‑1 (MQ=255)
                                                                                  aTGCGTTTCGACATCCTGCCCAAATGACCCTCGACGGGCGCCGGCGCCCCCCCCCCCGCG‑GCCGCCAAACC‑G‑CCTTTTTGCTGACACACCCCCCGGCCCGTCCCCCCCCACGGTCCCCCCCAGGCCCCCAACAGCaa                                                    >  5:331840/1‑137 (MQ=255)
                                                                                        ttCGACATCCTGCCCAAATTACCCTCGACGGCCTCATGCGCCCCCGCGCCCGAG‑GCCG‑CAAACC‑GTCATTTTTCCTGCAACACCGCCCTGCCCTTCACCGCCCTCGGACCGCGCCAGGCCCAGAACAGCAAACCggg                                              <  8:310513/137‑1 (MQ=255)
                                                                                          cGACATCCTGCCCAAATGACCCTCGACGGCCTCAGGCGCCCCCGCGCCCGAG‑GCCG‑CAAACC‑GTCATTTTTCCTGCAACACCGCCCTGCCCTTCACCGCCCTCGGACCGCGCCAGGCCCAGAACAGCAAACCGGGGa                                            <  8:171840/137‑1 (MQ=255)
                                                                                             cATCCTGCCCAAATGACCCTCGACGGCCGCAGGCGCCCCCCCGCCCCCG‑GCCG‑CCAACCGGGCATTTTT‑ATGCAAAAACCCCCAGCCCATGCCCCCCCCCGGCCCCCGCCCCGGCCAGAACCGCAAACCCGGGGAgg                                         >  4:169403/1‑137 (MQ=255)
                                                                                                      cAAATGACCGTCGACGGCCTCAGGCGCCCCCGCGCCCGAG‑GCCG‑CAAACC‑GTCATTTTTCCTGCAACACCGCCCTGCCCTTCACCGCCCTCGGACCGCGCCAGGCCCAGAACAGCAAACCGGGGAAGGGTGCGTAGa                                <  8:3352/137‑1 (MQ=255)
                                                                                                                gcGACGGCCTCAGGCGCCCCCGCGCCCGAG‑GCCG‑CAAACC‑GTCATTTTTCCTGCAACACCGCCCTGCCCTTCACCGCCCTCGGACCGCGCCAGGCCCAGAACAGCAAACCGGGGAAGGGTGCGTAGACCACGAAGAc                      <  3:164821/136‑1 (MQ=255)
                                                                                                                  gccGGCCTCAGGCGCCCCCGCGCCCGAG‑GCCG‑CAAACC‑GTCATTTTTCCTGCAACACCGCCCTGCCCTTCACCGCCCTCGGACCGCGCCAGGCCGAGAACAGCAAACCGGGGAAGGGTGCGTAGACCACGAAGACGa                    <  8:299362/135‑1 (MQ=255)
                                                                                                                     ggCGTCAGGCGCCCGGGCGCGCGAG‑GCCG‑CAAACC‑GTCATTTTTCCTGCAACACCGCGGTGCCCTTCACCTCCCTCTGAGCGCGGCAGGCCCAGAACTGGAAACAGGGGAAGGGTGCGTAGACGACGAAGACGATcc                 <  7:294833/137‑1 (MQ=255)
                                                                                                                                 ccccGCGCCCGAG‑GCCG‑CAAACC‑GTCATTTTTCCTGCAACACCGCCCTGCCCTTCACCGCCCTCGGACCGCGCCAGGCCCAGAACAGCAAACCGGGGAAGGGTGCGTAGACCACGAAGACGATCCACAGCATCTTgc     <  2:357540/137‑1 (MQ=255)
                                                                                                                                    cgcgCCCGAG‑GCCG‑CAATCC‑GTCATTTTTCCTCCAACCCCGCCCTGCCCTTCCCCGCCCTCGGACCGCGCCAGGCCCACAACAGCAAACCGGGGAAGGGTGCGTAGACCACGAAGACGAACCCCAGCATCTTgcgcg  >  8:45650/1‑136 (MQ=255)
                                                                                                                                    cgcgCCCGAG‑GCCG‑CAAACC‑GTCATTTTTCCTGCAACACCGCCCTGCCCTTCACCGCCCTCGGACCGCGCCAGGCCCCGAACAGCAACCAGGCGAAGGCTGCGTAGACCACGAAGACGATCCCCAGCATCTTGAGCt  >  2:183054/1‑137 (MQ=255)
                                                                                                                                     |                                                                                                                                          
CCATGACGCGCTCAAGGAACTGCGCGTACTGTCCACCGTGAAGTTCTACGGCTGCCTGATGGTGGACGGGCTGGGTCCTGCCATGCGTTTCGACATCCTGCCCAAATGACCCTCGACGGCCTCAGGCGCCCCCGCGCCCGAG‑GCCG‑CAAACC‑GTCATTTTTCCTGCAACACCGCCCTGCCCTTCACCGCCCTCGGACCGCGCCAGGCCCAGAACAGCAAACCGGGGAAGGGTGCGTAGACCACGAAGACGATCCACAGCATCTTGCGCT  >  ACB143/3989842‑3990110

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