Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT2440828,8690AC54.5% 6.6 / 18.6 22L87R (CTG→CGG) CJ019_00791putative sulfate transporter
Reads supporting (aligned to +/- strand):  ref base A (7/3);  new base C (1/11);  total (8/14)
Fisher's exact test for biased strand distribution p-value = 6.19e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.97e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GCGCAGCAGCCCGGCAATCACACAGAACGCCCCGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAGCCAGGATAGGCGCGGCCAGCGCCGAGTCCGGGCCCAGCACCAGGATACGGCTGGGACCGAACAGGGCATAGGCCAGCAGGGGGATGATGGTAGCGTACAGGCCATAGATACCGGGGACGCCGGAAGCCTCGGCGTAGGCGATGC  >  KT2440/828745‑828995
                                                                                                                            |                                                                                                                              
gcgcAGCAGCCCGGCAATCACACAGAACGCCCCGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAGCCAGGATAGGCGCGGCCAGCGCCGAGTcc                                                                                                                     >  8:183085/1‑136 (MQ=255)
        gcCCGGCAATCACACAGAACGCCCCGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAGCCAGGATAGGCGCGGCCAGCGCCGAGTCCGGGCCCAg                                                                                                             <  7:183085/136‑1 (MQ=255)
              cAAGCACACAGAACGCCCCGGCGACCAGCGCCATCATGGTGGCGATGGCGTTCGCCCGGTGCGGGGGGCGGGCGGCGTATTGCACCACCACAGCCCGGATAGGCGCGGCCAGCGCCGAGTCCGGGCCCAGCACCAg                                                                                                       <  8:200916/136‑1 (MQ=255)
                        cacccccccGCCGACCCGCGCCATCCTGGGGGGGAGGGCGAGCGCCCGGGGGGGGGCGCGGGGGGCTTATTCCACCCCCCCCGCCCGGATCGGCGCGCCCCGCCCCGAGTCCGGGCCCAGCACCAGGATACGGCTg                                                                                             <  4:32657/131‑1 (MQ=255)
                            gCCCCGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAGCCAGGATAGGCGCGGCCCGCGCCGAGTCCGGGCCCAGCCCCAGGATACGGCTGGGAc                                                                                         >  1:140494/1‑136 (MQ=255)
                            gCCCCGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAGCCAGGATAGGCGCGGCCAGCGCCGAGTCCGGGCCCAGCACCAGGATACGGCTGGGAc                                                                                         >  7:131562/1‑136 (MQ=255)
                             ccccGGCCACCCGCCCCACCTGGCTGGCGTGGGCGACCGCCGGCTGGGGGGCGCGGGCGGCTTATTCCCCCCCCACCCCCCGTATCGGCGCGCCCCGCCCCGAGTCCGGGCCCAGCACCAGGATACGGCTGGGAcc                                                                                        <  4:104431/136‑1 (MQ=255)
                             ccccGCCCACCCGCGCCAACAGGCTGGCGTGGCCGAGCGCCCGCGGGGGTGCGGGGGCGGCGTATTTCCCCACCCCCGCCAGTACAGCCGGGCCCATCGCCGAGTCCGGGCCCAGCACCAGGATACGGCTGGGAcc                                                                                        <  6:184084/136‑1 (MQ=255)
                              cccGGCGACCCGCGCCATCAGCCTGGCGTGGGCGAGCGCCCGCGGCGGGGCGCGGGCGGCGAATTGCACCCCCACCCCCCGCAACGGCGCGCCCCGCGCCGAGTCCGGGCCCAGCACCAGGATACGGCTGGGACCg                                                                                       <  8:131562/136‑1 (MQ=255)
                               cgggCGACCAGCCCCATCATCCGGGCGAGGGCGAGCGCCCGCTGCGGGGCGGGGGCGGCTTATGCCCCCCCCCCCCCCAGCATCGCCGCGGCCCGCGCCGAGTCCGGGCCCAGCACCAGGATACGGCTGGGACCGa                                                                                      <  4:230638/134‑1 (MQ=255)
                               ccGGCGCCCCGCGCCAGCAGGCGGGCGATGGCGAGCGCCCGCGGCGGGGCGCGGGCGGCTTATTCCCCCCCCCCCCCCAGCATCGGCGCGCCCCGCGCCGAGTCCGGGCCCAGCACCAGGATACGGCTGGGACCGa                                                                                      <  8:68390/136‑1 (MQ=255)
                                    cccccgcgcCACCAGGCGGGCGTGGGCGAGCGCCCGCGGCGGGGCGCGGGGGGCGTATGGCCCCCCCCCCCCCAGGACCGGCGCGCCCCGCGCCGAGTCCGGGCCCAGCACCAGGATACGGCTGGGACCGAACAgg                                                                                 <  3:43495/131‑1 (MQ=255)
                                      cccgcgcCTTCGTGGGGGCGTTGGGGAGCGCCCGGGGCGGGTCGGGGGCGGCTTTTTGCCCCCCCCCCGCGCGCCGCGGGCCGCCCCGCGCCGAGTCCGGGCCCAGCACCAGGATACGGCTGGGACCGAACAggg                                                                                <  2:140494/132‑1 (MQ=255)
                                                          atggcgatCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAGCCAGGATAGGCGCGGCCAGCGCCGAGTCCGGGCCCAGCACCAGGATACGGCGGGGACCGAACAGGGCATAGGCCAGCAGGGGgatga                                                           >  8:116523/1‑136 (MQ=255)
                                                                                     gcggagAATTCCACCCCCACCGCCTGGTTCGCCGCGGCGCGCGCCGAGTCCGGGCCCAGCACCAGGATACGGCTGGGACCGAACAGGGCATAGGCCAGCAGGGGGATGATGGTAGCGTACAGGCCATAGATACCgg                                <  2:16495/136‑1 (MQ=255)
                                                                                                 cccaccacACCCTGCATAGGCGGGGCCCCCGCCGAGTCCGGGCCCAGCACCAGGATACGGCTGGGACCGAACAGGGCATAGGCCAGCAGGGGGATGATGGTAGCGTACAGGCCATAGATACCGGGGACGCCGGAAg                    <  8:280/135‑1 (MQ=255)
                                                                                                        caGCCAGGATAGGCGCGGCCAGCGCCGAGTCCGGGCCCAGCACCAGGATACGGCTGGGACCGAACAGGGCATAGGCCAGCAGGGGGATGATGGTAGCGTACAGGCCATAGATACCGGGGACGCCGGAAGCCTCGGc             >  5:80298/1‑136 (MQ=255)
                                                                                                        caGCCAGGATAGGCGCGGCCAGCGCCGAGTCCGGGCCCAGCACCAGGATACGGCTGGGACCGAACAGGGCATAGGCCAGCAGGGGGATGATGGTAGCGTACAGGCCATAGATACCGGGGACGCCGGAAGCCTCGGc             >  7:177219/1‑136 (MQ=255)
                                                                                                           ccAGGATAGGCGTGGCCAGCGCCGAGTCCGGGCCCAGCACCAGGATACGGCTGGGACCGAACAGGGCATAGGCCAGCAGGGGGATGATGGTAGCGTACAGGCCATAGATACCGGGGACGCCGGAAGCCTCGGCg            >  3:141840/1‑134 (MQ=255)
                                                                                                            cAGGATCGGCGCGGCCCGCGCCGAGTCCGGGCCCAGCACCAGGATACGGCTGGGACCGAACAGGGCATAGGCCAGCAGGGGGATGATGGTAGCGTACAGGCCATAGATACCGGGGACGCCGGAAGCCTCGGCGTAg         <  4:110131/136‑1 (MQ=255)
                                                                                                                  aGGCGCGGCCAGCGCCGAGTCCGGGCCCAGCACCAGGATACGGCTGGGACCGAACAGGGCATAGGCCAGCAGGGGGATGATGGTAGCGTACAGGCAATAGATACCGTGGACGCCGGAAGCCTCGGCGAAGGCGATg   >  1:127499/1‑136 (MQ=255)
                                                                                                                   gccgcgCCCCGCCCCGAGTCCGGGCCCAGCACCCGGATACGGCGGGGACCGAACAGGGCATAGGCCAGCAGGGGGATGATGGTAGCGTACAGGCCATAGATACCGGGGACGCCGGAAGCCTCGGCGTAGGCGATGc  <  7:119749/134‑1 (MQ=255)
                                                                                                                            |                                                                                                                              
GCGCAGCAGCCCGGCAATCACACAGAACGCCCCGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAGCCAGGATAGGCGCGGCCAGCGCCGAGTCCGGGCCCAGCACCAGGATACGGCTGGGACCGAACAGGGCATAGGCCAGCAGGGGGATGATGGTAGCGTACAGGCCATAGATACCGGGGACGCCGGAAGCCTCGGCGTAGGCGATGC  >  KT2440/828745‑828995

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑