Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB143188,0210GC44.4% 19.0 / 14.1 36E960D (GAG→GAC) CJ019_00185hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (11/9);  new base C (0/16);  total (11/25)
Fisher's exact test for biased strand distribution p-value = 5.53e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGT‑AATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAG  >  ACB143/187887‑188149
                                                                                                                                       |                                                                                                                                
aaCGGCCAAACGCTTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACATGCGGT‑AATTTCGAg                                                                                                                                  >  1:113127/1‑135 (MQ=255)
 aCGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGT‑AATTTCGAGCa                                                                                                                                <  5:82894/136‑1 (MQ=37)
                  gTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGT‑AATTTCGAGCAACTGACCCCGAACCCGa                                                                                                               >  7:204889/1‑136 (MQ=34)
                       aGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGT‑AATTTCGAGCAACTGACCCCGAACCCGACCCCg                                                                                                          <  7:145903/136‑1 (MQ=37)
                              aaaaCCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGT‑AATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAg                                                                                                    >  8:214809/1‑135 (MQ=37)
                              aaaaCCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACAGTCAGCGTCACCATCGAGAACGCCACAGGCGGT‑AATTTCGAGCAACTGCCCCCGAACCCGACCCCGGCTCAGa                                                                                                   >  4:233649/1‑136 (MQ=255)
                                               cGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGT‑AATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGa                                                                                                   <  8:72399/119‑1 (MQ=35)
                                               cGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGT‑AATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGa                                                                                                   >  7:72399/1‑119 (MQ=35)
                                               cGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGT‑AATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCg                                                                                  >  1:280763/1‑136 (MQ=37)
                                               cGACTTCCAGACCCCGGCGAATGACGTCGACAACAACGGTTCGACAGTCAGCGTCACCATCGAGAACGCCAAAGGCGGT‑AATTTCGAGCAACAGACCCCGAAACCGACACCGGCTCAGACCACGATCAACGACTCg                                                                                  >  1:184452/1‑136 (MQ=11)
                                                 aattCCCGACCCCGGCGAATGACGTCTACAACACCGGTTCGACTGTCAGCGTCACCATCGAGACCGCCACAGGCGGT‑AATTTCGAGCAACTGACCCCGAAACCCACCCCGGCTCAGACCACGATCAACGACTCGGt                                                                                <  6:160474/134‑1 (MQ=11)
                                                       aGACCCCGGCTAATCACTTATACAACAAGGGTTCGACTGTCAGCGTCACAAGCGATACCGCCACGGGCGTT‑AATTTCGAGCAATTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGacac                                                                          <  2:280763/136‑4 (MQ=1)
                                                            cgggCTAATTAGGTATAAAACATCGGTTCGACGGTCTGTGTAACCTGCGACACCGCCAGTGGCGTA‑AATTTCCACCATCTCACCCCTAACCCGCCCCCGGCTCAGACCACGATCAACGACTCGGTCGATGccacca                                                                     <  6:99517/134‑1 (MQ=2)
                                                             cGGCCATTGACGTCAAAAACACCGGTTCGCGTGTCAGCTCCACAATCGACACCCCCAC‑TGCGGTAAATTTCGAGCAACTCCCCCCGACCCCGACCCCGGCTCACACCACTATCAACGACTCGGTCGATGCcaccac                                                                    <  8:5999/136‑1 (MQ=2)
                                                                                      ttCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGT‑AATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACAACCGCGACCCTGACGGCGAGCCCGTCg                                           >  7:293876/1‑136 (MQ=21)
                                                                                      ttCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGT‑AATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATACCACCACCGCGACCCTGACGGCGAGCCCGTCg                                           >  7:188894/1‑136 (MQ=21)
                                                                                        cGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGT‑AATTTCGAGCAACTGACCCCGAAACCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGAGCCCACCACCGCGACCCTGACGGCGAGCCCGTCGGt                                         >  6:248773/1‑136 (MQ=12)
                                                                                                   gcaacaaTCGACACCCCCACAGGCGTT‑AATTTCCACCAACTCCCCCCGAACCCCCCCCCGGCTCACCCCACGACCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcg                              <  3:123044/130‑1 (MQ=11)
                                                                                                     caccacCGACAACGCCACCGGCGGT‑AATTTCGACCAACTGCCCCCTAACCCGCCCCCGGCTCACACCCCCACCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcggc                            <  2:82472/136‑1 (MQ=11)
                                                                                                     cACCATCGATAACGCACTCGGCGGT‑AATTTCGACCACCTCACCCCGACCCGCACCGCGGATCAGACCACTACCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcggc                            <  2:18991/136‑1 (MQ=11)
                                                                                                     cACCATCGAGAACGCCACCGGCGTT‑AATTTCGCCCACCTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcggc                            <  7:86216/136‑1 (MQ=11)
                                                                                                     cACCATCGAGAACGCCACCGGCGGT‑AATTTCGCGCAACTGACCCCGACCCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcggc                            <  1:282652/136‑1 (MQ=12)
                                                                                                     cACCATCGAGAACGCCACCGGCGGT‑AATTTCGACCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcggc                            <  3:238462/136‑1 (MQ=17)
                                                                                                     cACCATCGAGAACGCCACCGGCGGA‑AATTTCGACCACCTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcggc                            <  8:104265/136‑1 (MQ=255)
                                                                                                     cACCATCGAGAACGCCACAGGCGTT‑AATTTCCAGCAACTGACCCCCAACCCGCCCCCGGCTCAGACCACGACCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcggc                            <  4:121311/136‑1 (MQ=11)
                                                                                                     cACCATCGAGAACCCCACTGGCGGA‑AATTGCCACCACCTCACCCCGAACCCGACCCCGGCTCACACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcggc                            <  3:33941/136‑1 (MQ=11)
                                                                                                     cACCATCGACACCGCCACAGGCGGT‑AATTTCGACCAACTGACCCCGAACCCGCCCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcggc                            <  3:248448/136‑1 (MQ=255)
                                                                                                      aCCATCTACAACGCCACAGGCGTT‑AATTTCGACCACCTGCCCCCGACCCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACTGATAGCCCGTCGGTCACCGAAggcgga                            <  6:148926/135‑2 (MQ=11)
                                                                                                       ccATCGACACCGCCACTGGCGGA‑AAATCCTACCAACTGACCCCTAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGt                          <  6:50805/136‑1 (MQ=11)
                                                                                                         aTCGAGAACGCCACAGGCGGT‑AATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGtga                        >  6:160204/1‑134 (MQ=11)
                                                                                                             agaACGCCACAGGCGGT‑AATTTCGACAAACTCACCCCGAACCCGACCCCGGCGCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCAc                    <  6:311071/136‑1 (MQ=255)
                                                                                                             agaACGCAACGGGCGTT‑AATTTCTAGAAACTGACCCCTAACCCGACCCCGGCTAAGACCACGATCCACGACTCGGTCGAATCCACAACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCAc                    <  4:77314/136‑1 (MQ=11)
                                                                                                               caCGCCACGGGGGGT‑AATTTCGACCACCGGACCCCGAACCCGCCCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCt                  <  6:106333/135‑1 (MQ=11)
                                                                                                                  gCCACGGGCGGA‑AATTTCTACCACCGGACCCCGACCCCGACCCCGGCTCAGACCACTAACAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTaca               <  2:127557/136‑1 (MQ=11)
                                                                                                                  cccACCGGGGTT‑AATTACCACCAACTCCCCCCGAACCCGCCCCCGGCTCACACCACTAACAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTaca               <  2:358227/135‑1 (MQ=11)
                                                                                                                            gT‑AAATTCCACCACCTGACCCCGAACCCGACCCCGGCTAAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCt     <  2:62577/136‑1 (MQ=11)
                                                                                                                                aTTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGAACCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACGTACACCGTGACCCTGAg  >  2:253069/1‑136 (MQ=12)
                                                                                                                                       |                                                                                                                                
AACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGT‑AATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAG  >  ACB143/187887‑188149

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