Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1431,188,7210AC56.7% 9.9 / 24.6 34Q52P (CAG→CCG) CJ019_01117hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (9/4);  new base C (5/12);  total (14/20)
Fisher's exact test for biased strand distribution p-value = 6.34e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.40e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TCACCGGTGCTGGCCTGACCCTGGCCGGCTGGCTTTCTGCGCAGTGTGTGCTGCAGTTGAGCCGGCAGCCACACCCGGCCACGGCGCCTGCGGTGGCCGACAGCCCGCTGCCCGGGCTGATGGTTGGTCACTGGCAGGCGCCGCCCGACGACGGCGCCATCACCATCACTCGTCTGCCGCTGCAATACCTCGGCGGCCTCAAGGCCCAGCCATTGGCCGCCAGTGTGGTGGTGCTGCGCTACGGCGAGCAAGTGCGCACCCTGGC  >  ACB143/1188586‑1188850
                                                                                                                                        |                                                                                                                                 
tCACCGGTGCTGGCCTGACCCTGGCCGGCTGGCTTTCTGCGCACTGTGTGCTGCAGTTGAGCCGGCAGCCACACCCGGCCACGGCGCCTGCGGTGGCCGACCGACCCCAGCCCCCGCTGATGGTTGGTCACAGGCa                                                                                                                                   >  1:202224/1‑136 (MQ=255)
            gCCTGACCCTGGCCGGCTGGCTTTCTGCGCAGTGTGTGCTGCAGTTGAGCCGGCAGCCACACCCGGCCACGGCGCCTGCGGTGGCCGACCGCCCGCGGCCCGGGCTGCTGGTTGGTCACGGGCAGGCGCCCCCcgc                                                                                                                       >  7:217231/1‑135 (MQ=255)
            gCCTGACCCTGGCCGGCTGGCTTTCTGCGCAGTGTGTGCTGCAGTTGAGCCGGCAGCCACACCCGGCCACGGCGCCTGCGGTGGCCGACAGCCCGCAGCCCGGGCTGATGGTTGGGCACTGGCCGGCGcccccccc                                                                                                                       >  2:284804/1‑134 (MQ=255)
            gCCTGACCCTGGCCGGCTGGCTTTCTGCGCAGTGTGTGCTGCAGTTGAGCCGGCAGCCACACCCGGCCACGGCGCCTGCGGTGGCCGACACCCCCCGGCCCCGGCAGCAGGGTGGGCACTGGCCGGCGcccccccc                                                                                                                       >  1:122809/1‑134 (MQ=255)
               tGACCCTGGCCGGCTGGCTTTCTGCGCAGTGTGTGCTGCAGTTGAGCCGGCAGCCACACCCGGCCACGGCGCCTGCGGTGGCCGACAGCCCGCGGCCCGGGCGGATGGGTGGTGCCTGGCCCGCGccccccccccc                                                                                                                    >  1:263148/1‑131 (MQ=255)
                         cGGCTGGCTTTCTGCGCAGTGTGTGCTGCAGTTGAGCCGGCAGCCACACCCGGCCACGGCGCCTGCGGTGGCCGACAGCCCGCTGCCCGGGCGGATGGTTGGTGACTGGCCGGCGCCGCCCGCCCCCGGCGcccga                                                                                                          >  7:269180/1‑133 (MQ=255)
                         cGGCTGGCTTTCTGCGCAGTGTGTGCTGCAGTTGAGCCGGCAGCCACACCCGGCCACGGCGCCTGCGGTGGCCGACAGCCCGCGGCCCGGGCTGATGGTTGGTCACTGGCAGGCGCCCCCCCCCCCCGGCGccct                                                                                                           >  1:282836/1‑133 (MQ=255)
                         cGGCTGGCTTTCTGCGCAGTGTGTGCTGCAGTTGAGCCGGCAGCCACACCCGGCCACGGCGCCTGCGGTGGCCGACAGCCCGCGGCCCGGGCGGATGGGTGGGGACTGGCCGGGGCCCCCCCCCCCCGGCGcccg                                                                                                           >  6:112791/1‑133 (MQ=255)
                                                      tgTTGCGCCGGCAGCCACCCCCGCCCACGGCGCCGGCGGGGGCCGACAGCCCGCTGCCCGGGCTGATGGTTGGTCACTGGCCGGCGCCGCCCGACGACGGCGCCATCACCATCACTCGTCTGCCGCTGCAATACCt                                                                             <  1:55013/135‑1 (MQ=255)
                                                       gTTGTGCCGGCAGCAACACCGGGCCAGGGCGCCGGGGGTGGCGGACAGCCCGCTGCCCGGGGTGATGGTTGGTCACTGGCTGGGGCCGCCCGACGACGGCGCCATCACCATCACTCGTCTGCCGCTGCAATACCt                                                                             <  5:259402/135‑1 (MQ=255)
                                                             ccGGCAGCCACACCCGGCCACGGCGCCTGCGGTGGCCGACAGCCCGCTGCCCGGGCTGATGGTTGGTCACTGGCAGGCGCCCCCCCCCCCCGGGGCCCTCACCATCCAGCGTCTGCCGCAGGCATACCAAggcggc                                                                      >  4:278139/1‑136 (MQ=255)
                                                             ccGGCACCCAACCCCGGCCCGGCCCCCGGGGGGGCCCTACACCCCGCTGCCCGGGCTGAGGGTTGGTCAGGGGCCGGGGCCGCCCGACGACGGCGCCATCACCATCACTCGTCTGCCGCTGCAATACCTCtgcggc                                                                      <  1:116764/136‑1 (MQ=255)
                                                                   gCCACACCGGGCCCCGGCGCGGGGGGGGGCCGACCGCCCGCTGCCGGGGGTGCGGGTTGGGCCGGGGCTGGGGCGGCCCGACGACGGCGCCATCACCATCACTCGTCTGCCGCTGCAATACCTCGGCGGCCTCAAg                                                                <  7:122715/136‑1 (MQ=255)
                                                                              ccccGGCGCGTGGGGGGGCCGACAGCCCGCGGCCCGGGCTGAGGGTGGGGCACTGCCCGGCGCCGCCCGACGACGGCGCCATCACCATCACTCGTCTGCCGCTGCAATACCTCGGCGGCCTCAAGGCCCAGCCAtt                                                     <  6:96198/133‑1 (MQ=255)
                                                                                 gggCGCCGGGGGTGGCCGACAGCCCGCTGCCCGGGCTGAGGGTGGGTCACTGGCCGGCGCCGCCCGACGACGGCGCCATCACCATCACTCGTCTGCCGCTGCAATACCTCGGCGGCCTCAAGGCCCAGCCATTGgc                                                  <  2:295050/135‑1 (MQ=255)
                                                                                         gCGGGGGCCGACAGCCCGCTGCCCGGGCTGATGGTGGGTCACTGGCAGGCGCCGCCCGACGACGGCGCCATCACCATCACTCGTCTGCCGCTGCAATACCTCGGCGGCCTCAAGGCCCAGCCATTGGCCGCCAgtg                                          <  7:323707/136‑1 (MQ=255)
                                                                                              ggCCGACAGCCCGCTGCCCGGGTTTTTGTTTGGTCCCTGCCAGGCGCCGCCCGACGACGGCGCCATCACCATCACTCGTCTGCCGCTGCAATACCTCGGCGGCCTCAAGGCCCAGCCATTGGCCGCCAGTgtggtg                                     <  1:112627/136‑1 (MQ=255)
                                                                                              ggCCCACCGCCCGCGCCCCGGGCTTATGGTTGGCCACTGCCCGGCGCCGCCCGACGACGGCGCCATCACCATCACTCGTCTGCCGCTGCAATACCTCGGCGGCCTCAAGGCCCAGCCATTGGCCGCCAGTgtggtg                                     <  5:317743/136‑1 (MQ=255)
                                                                                                         cggcggCCGGGGCTGTTGGGTGGGCACTGGCTGGCGCCGCCCGACGACGGCGCCATCACCATCACTCGTCTGCCGCTGCAATACCTCGGCGGCCTCAAGGCCCAGCCATTGGCCGCCAGTGTGGTGGTGCTGCGCt                           <  2:232490/131‑1 (MQ=255)
                                                                                                                cggggCTTTTGGTGGGTCACTGGCAGGCGCCGCCCGACGACGGCGCCATCACCATCACTCGTCTGCCGCTGCAATACCTCGGCGGCCTCAAGGCCCAGCCATTGGCCGCCAGTGTGGTGGTGCTGCGCTACGGCGa                    <  4:311194/134‑1 (MQ=255)
                                                                                                                 ggggCTGTTGGGTGGCCACGGCCCGCCCCCCCCCGACGACGGCGCCATCACCATCACTCGTCTGCCGCTGCAATACCTCGGCGGCCTCAAGGCCCAGCCATTGGCCGCCAGTGTGGTGGTGCTGCGCTACGGCGAg                   <  3:343679/135‑1 (MQ=255)
                                                                                                                     cTGATGGTGGGCCAGGGGACGGCCGCGCCCGACGACGGCGCCATCACCATCACTAGTCTGCCGCTGCAATACCTCGGCGGCCTCAAGGCCCAGCCATTGGCCGCCAGTGTGGTGGTTCTGCGCTACGGCGAGCAAg               <  4:84534/136‑1 (MQ=255)
                                                                                                                      tgttggttggtCACTGGCCGGCGCCGCCCGACGACGGCGCCATCACCATCACTCGTCTGCCGCTGCAATACCTCGGCGGCCTCAAGGCCCAGCCATTGGCCGCCAGTGTGGTGGTGCTGCGCTACGGCGAGCAAGt              <  5:323469/133‑1 (MQ=255)
                                                                                                                       gATGGTGGGTCACTGGCCGGCGCCGCCCGACGACGGCGCCATCACCATCACTCGTCTGCCGCTGCAATACCTCGGCGGCCTCAAGGCCCAGCCATTGGCCGCCAGTGTGGTGGTGCTGCGCTACGGCGAGCAAGt              <  2:326577/135‑1 (MQ=255)
                                                                                                                        aTGGTTGGTCACTGGCAGGCGCCGCCCGACGACGGCGCCATCACCATCACTCGTCTGCCGCTGCAATACCTCGGCGGCCTCAAGGCCCAGCCATTGGCCGCCAGTGTGGTGGTGCTGCGCTACGGCGAGCAAGt              <  7:175217/134‑1 (MQ=255)
                                                                                                                          ggttggtCACTGGCTGGGGCGGCCCGACGACGGCGCCATCACCATCACTCGTCTGCCGCTGCAATACCTCGGCGGCCTCAAGGCCCAGCCATTGGCCGCCAGTGTGGTGGTGCTGCGCTACGGCGAGCAAGTgcgc          <  2:202224/136‑1 (MQ=255)
                                                                                                                           ggtggtCACTGGCCGGCGCGGCCCGACGACGGCGCCATCACCATCACTCGTCTGCCGCTGCAATACCTCGGCGGCCTCAAGGCCCAGCCATTGGCCGCCAGTGTGGTGGTGCTGCGCTACGGCGAGCAAGTgcgc          <  7:21788/133‑1 (MQ=255)
                                                                                                                            ttggtCACTGGCCGGCGCCGCCCGACGACGGCGCCATCACCATCACTCGTCTGCCGCTGCAATACCTCGGCGGCCTCAAGGCCCAGCCATTGGCCGCCAGTGTGGTGGTGCTGCGCTACGGCGAGCAAGTGCGCAc        <  2:49983/136‑1 (MQ=255)
                                                                                                                            ttggtCACTGGCAGGCGCCGCCCGACGACGGCGCCATCACCATCACTAGGCTGCCGCTGCAATACCTCGGCGGCCGCAAGGCCCCGCCAATGGCCCCCAGTGTGGTGGTGGTGCGCTACGGCGCGGAAGTGCGCAc        >  2:18960/1‑136 (MQ=255)
                                                                                                                            ttggtCACTGGCAGGCGCCGCCCGACCACCGCGCCATACCCAACACTCGTCAGGCGCTGCAAAAACTCGGCGGCCTCAAGGCCCAGCCCTTTGCCGCCAGGGTGGTGGTTCTGCGCTAAGGGGAGCAAGTgcgcca        >  6:102905/1‑134 (MQ=255)
                                                                                                                                 cACTGGCAGGCGCCGCCCGCCGACGGCGCCATCACCATCACTCGTCTGCCGCTGCAATACCTCGGCGGCCTCAAGGCCCAGCCATTGGCCGCCAGTGTGGTGGTGCTGCGCTACGGCGAGCAAGTGCGCACCCTgg   >  1:282340/1‑136 (MQ=255)
                                                                                                                                 cACTGGCAGGCGCCGCCCGACGCCGGCGCCATCACCATCACTCGGCTGCCGCGGCAATACCACGGCGGCCTCAACGCCCCGCCAATTGCCCCCACTGTGGTGGTGCTGCGCTAAGGCGAGCACTGGCGCAcccat    >  2:118209/1‑133 (MQ=255)
                                                                                                                                 cACTGGCAGGCGCCGCCCGACGACGGCGCCATCACCATCACTCGTCTGCCGCTGCAATACCTCGGCGGCCTCAAGGCCCCGCCAATGGCCGCCAGTGTGGTGGTGCTGCGCTACGGCGAGCAAGTGCGCACCCTgg   >  2:261004/1‑136 (MQ=255)
                                                                                                                                  aCTGGCCGGCGCCGCCCGACGACGGCGCCATCACCATCACTCGTCTGCCGCTGCAATACCTCGGCGGCCTCAAGGCCCAGCCATTGGCCGCCAGTGTGGTGGTGCTGCGCTACGGCGAGCAAGTGCGCACCCTGGc  <  8:345150/136‑1 (MQ=255)
                                                                                                                                  aCTGGCAGGCGCCGCCCGCCGCCGGCGCCCTCACCATAACTCGTCTGCCGCTGCAATAACTCGGCGGCCGCAAGGCCCCGCCATTGGCCCCCACTGTGGTGGTGCTGCGCTACGGCGCGCAAGTGCGCCCCCTGGc  >  8:175217/1‑136 (MQ=255)
                                                                                                                                        |                                                                                                                                 
TCACCGGTGCTGGCCTGACCCTGGCCGGCTGGCTTTCTGCGCAGTGTGTGCTGCAGTTGAGCCGGCAGCCACACCCGGCCACGGCGCCTGCGGTGGCCGACAGCCCGCTGCCCGGGCTGATGGTTGGTCACTGGCAGGCGCCGCCCGACGACGGCGCCATCACCATCACTCGTCTGCCGCTGCAATACCTCGGCGGCCTCAAGGCCCAGCCATTGGCCGCCAGTGTGGTGGTGCTGCGCTACGGCGAGCAAGTGCGCACCCTGGC  >  ACB143/1188586‑1188850

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: