Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24402,113,8950TG30.4% 67.4 / 10.8 46R663R (CGT→CGGCJ019_02007hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (8/24);  new base G (14/0);  total (22/24)
Fisher's exact test for biased strand distribution p-value = 1.33e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GAAGGACGGGGCTCGGTCGTCGCTGAGGGCGATGGCAGGCGTTGCACCTATACCGCAGGCAGTAAGGTGGATGCCAACGTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGTTCAATCTATGTGCTGGTACGCCAGCGGCGGCCTGAACTGGTGATCGACGTCAAGCAGTTGCCCGATGGCGACCTGCAGTTGACAGGTATGTTCGATGGGGTAGGCCAGCCAGACACCCA  >  KT2440/2113761‑2114014
                                                                                                                                      |                                                                                                                        
gaagGACGGGGCTCGGTCGTCGCTGAGGGCGATGGCAGGCGTTGCACCTATACCGCAGGCAGTAAGGTGGATGCCAACGTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGTTc                                                                                                                        <  7:114052/137‑1 (MQ=255)
 aagGACGGGGCTCGGTCGTCGCTGAGGGCGATGGCAGGCGTTGCACCTATACCGCAGGCAGTAAGGTGGATGCCAACGTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGTTCa                                                                                                                       >  4:133134/1‑137 (MQ=255)
           cTCGGTCGTCGCTGAGGGCGATGGCAGGCGTTGCACCTATACCGCAGGCAGTAAGGTGGATGCCAACGTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGTTCAATCTATGTGc                                                                                                             <  3:133134/137‑1 (MQ=255)
                  gtcgCTGAGGGCGATGGCAGGCGTTGCACCTATACCGCAGGCAGTAAGGTGGATGCCAACGTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGTTCa                                                                                                                       <  1:99111/120‑1 (MQ=255)
                  gtcgCTGAGGGCGATGGCAGGCGTTGCACCTATACCGCAGGCAGTAAGGTGGATGCCAACGTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGTTCa                                                                                                                       >  2:99111/1‑120 (MQ=255)
                     gCTGAGGGCGATGGCAGGCGTTGCACCTATACCGCAGGCAGTAAGGTGGATGCCAACGTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGTTCAATCTATGTGCTGGTACGCCa                                                                                                   <  5:21893/137‑1 (MQ=255)
                     gCTGAGGGCGATGGCAGGCGTTGCACCTATACCGCAGGCAGTAAGGTGGATGCCAACGTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGTTCAATCTATGTGCTGGTACGCCa                                                                                                   <  4:171943/137‑1 (MQ=255)
                      cTGAGGGCGATGGCAGGCGTTGCACCTATACCGCAGGCAGTAAGGTGGATGCCAACGTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGGTGAATCTATGTGGTGGTACGCCac                                                                                                  >  4:67247/1‑136 (MQ=255)
                                tGGCAGGCGTTGCACCTATACCGCAGGCAGTAAGGTGGATGCCAACGTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGGTAAAAAAATTTTCTTGTACCCCAGCCGCGGcccg                                                                                        >  6:101785/1‑135 (MQ=255)
                                       cGTTGCACCTATACCGCAGGCAGTAAGGTGGATGCCAACGTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGTTCAATCTATGTGCTGGTACGCCAGCGGCGGCCTGAACTGGt                                                                                 <  7:67518/137‑1 (MQ=255)
                                                 ataCCGCAGGCAGTAAGGTGGATGCCAACGTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGTTCAATCTATGTGCTGGTACGCCAGCGGCGGCCTGAACTGGTGATCGACGTc                                                                      <  2:438094/137‑1 (MQ=255)
                                                 ataCCGCAGGCAGTAAGGTGGATGCCAACGTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGTTCAATCTATGTGCTGGTACGCCAGCGGCGGCCTGAACTGGTGATCGACGTc                                                                      <  6:440659/137‑1 (MQ=255)
                                                    ccGCAGGCAGTAAGGTGGATGCCAACGTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGTTCAATCTATGTGCTGGTACGCCAGCGGCGGCCTGAACTGGTGATCGACGTCAAg                                                                   <  1:15861/137‑1 (MQ=255)
                                                    ccGCAGGCAGTAAGGTGGATGCCAACGTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGTTCAATCTATGTGCTGGTAAGCCACCCGCGGCCGGAACTGGGGATCGACGTCAAg                                                                   >  4:349727/1‑137 (MQ=255)
                                                    ccGCAGGCAGTAAGGTGGATGCCAACGTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGGTGAATCTATTTGCTGGTACGCCAGCGGCGGCCTGAACTGGTGATCGACGTCaac                                                                   >  8:67518/1‑136 (MQ=255)
                                                    ccGCAGGCAGTAAGGTGGATGCCAACGTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGGTGAATCGATGTGCTGGTACGCCCGCGGCGGCCTGCACGGGTGGTCGACGTCAAg                                                                   >  8:332435/1‑137 (MQ=255)
                                                    ccGCAGGCAGTAAGGTGGATGCCAACGTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGGTCAATCTATGTGGTTGTAAGCCAGCGGCGGCCTGAACTGGTGCTAGCCGTCaac                                                                   >  2:95838/1‑136 (MQ=255)
                                                    ccGCAGGCAGTAAGGTGGATGCCAACGTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGGTCAATATAATTGCTGGTAAGCCACCGGCGGCCGGAACAGGTGATCGACCTCaaa                                                                   >  4:365933/1‑136 (MQ=255)
                                                         ggcagTAAGGTGGATGCCAACGTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGTTGAATCTATTTGGGGGTAACCCAGCGGCGGGCGGGACTGGGGGTCGGCGTCAAGCAGtt                                                              >  6:106921/1‑137 (MQ=255)
                                                          gcagTAAGGTGGATGCCAACGTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGTTCAATAGATGTGCTGGTACCCCACCGCCCGCCAGAACTGGTGATAGACCTCAAGCACTTg                                                             >  7:48638/1‑137 (MQ=255)
                                                          gcagTAAGGTGGATGCCAACGTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGGTGTATCTATGTTGTGGTACGCCCCCGGCGGGCGGCACAGGTGATCGACGTCAAGCAgatt                                                             >  4:51971/1‑134 (MQ=255)
                                                          gcagTAAGGTGGATGCCAACGTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGGTGAATCTAAGTTGAGGTACGCCCTCGGCGGCCGGAACTTGTGGTAGACGTCAACAAGTTg                                                             >  3:355873/1‑137 (MQ=255)
                                                          gcagTAAGGTGGATGCCAACGTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGGTCAAAATATTTGATGGTACGCCAGCCGCGGCCTGAACTGGTGATCGCCGTCAAGCACttt                                                             >  7:245997/1‑136 (MQ=255)
                                                          gcagTAAGGTGGATGCCAACGTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGGGGAATCTATGTGCTGGTTCGCCCCCGGCGGCCGGAACTGGTGATCGACGTCAAGCAGttt                                                             >  5:73157/1‑136 (MQ=255)
                                                                      aTGCCAACGTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGTTCAATCTATGTGCTGGTACGCCAGCGGCGGCCTGAACTGGTGATCGACGTCAAGCAGTTGCCCGATGGCGAc                                                 <  3:365933/137‑1 (MQ=255)
                                                                         ccAACGTGGCCTACGTGCTAGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCTGTGAATCTATTTGGAGGTTCACCCGCGCCGGCCGCCAATGGTGGTAGACCTCCAGCAGGTGCAACAGGGCGAcccg                                              >  8:136306/1‑135 (MQ=255)
                                                                           aaCGTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGTTGAATCAATGGTCTGGGACGCCCGCGGCGGCCGGCACTGGTGATCGACCGCCAAGATATTGCCCGTGGCGAccgtga                                            >  8:423399/1‑133 (MQ=255)
                                                                           aaCGTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGTTCAATCTATGTGCTGGTACGCCAGCGGCGGCCTGAACTGGTGATCGACGTCAAGCAGTTGCCCGATGGCGACCTGCa                                            <  3:67247/137‑1 (MQ=255)
                                                                           aaCGTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGGTGAATCTATGTGGTTGGACCCCACCCGCGCCCTGCTCTTGGTATCGAAGGGAAGCACTTTCACGATGGCGAACtcac                                            >  8:423394/1‑134 (MQ=255)
                                                                            aCGTGGCCTACGTGCTCGATCAGGTCAGGGGACAGAACCGGGAGACCAGCGAGTCGCGTTCAATCTATGTGCTGGTACGCCAGCGGCGGCCTGAACTGGTGATCGACGTCAAGCAGTTGCCCGATGGCGACCTGCAg                                           <  4:139814/137‑1 (MQ=255)
                                                                             cGTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGTTCAATCTATGTGCTGGTACGCCAGCGGCGGCCTGAACTGGTGATCGACGTCAAGCAGTTGCCCGATGGCGACCTGCAGt                                          <  1:159067/137‑1 (MQ=255)
                                                                             cGTGGCCTAAGTGCTCGATCAGGTCAGGGTACAGAACCGGGATACCAGCGAGTCTCGTTCATTCTATGTGCTGGTACGCCAGCGGCGGCCTGAACTGGTGATCGACGTCAAGCAGTTGCCCGATGGCGACCTGCAGt                                          <  5:182665/137‑1 (MQ=255)
                                                                                   cTACGTGCTCGATCAGGTCAGCGTACAGAACCGGGAGACCAGCGAGTCGCGTTCAATCTATGTGCTGGTACGCCAGCGGCGGCCTGAACTGGTGATCGACGTCAAGCAGTTGCCCGATGGCGACCTGCAGTTGACAg                                    <  8:310586/137‑1 (MQ=255)
                                                                                         gCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGTTCAATCTATGTGCTGGTACGCCAGCGGCGGCCTGAACTGGTGATCGACGTCAAGCAGTTGCCCGATGGCGACCTGCAGTTGACAGGTATGt                              <  3:196811/137‑1 (MQ=255)
                                                                                          cTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGTTCAATCTATGTGCTGGTACGCCAGCGGCGGCCTGAACTGGTGATCGACGTCAAGCAGTTGCCCGATGGCGACCTGCAGTTGACAGGTATGtt                             <  1:99722/137‑1 (MQ=255)
                                                                                          cTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGTTCAATCTATGTGCTGGTACGCCAGCGGCGGCCTGAACTGGTGATCGACGTCAAGCAGTTGCCCGATGGCGACCTGCAGTTGACAGGTATGtt                             <  6:218064/137‑1 (MQ=255)
                                                                                                    tcaggGTACAGAACCGGGAGACCAGCGAGTCGCGTTCAATCTATGTGCTGGTACGCCAGCGGCGGCCTGAACTGGTGATCGACGTCAAGCAGTTGCCCGATGGCGACCTGCAGTTGACAGGTATGTTCGATGGGGTa                   <  8:267369/137‑1 (MQ=255)
                                                                                                    tcaggGTACAGAACCGGGAGACCAGCGAGTCGCGTTCAATCTATGTGCTGGTACGCCAGCGGCGGCCTGAACTGGTGATCGACGTCAAGCAGTTGCCCGATGGCGACCTGCAGTTGACAGGTATGTTCGATGGGGTa                   <  3:51971/137‑1 (MQ=255)
                                                                                                     caggGTACAGAACCGGGAGACCAGCGAGTCGCGGTCAATATATTTGCTGGGAAGCCACCGGCGGGCGGCACTGGTGATCGCCGTCAACCAGTTGCCCCATTGCGACCTGCAGTTGACAGGGATTTTAGATGGGGTAg                  >  6:448099/1‑137 (MQ=255)
                                                                                                      aggGTACAGAACCGGGAGACCAGCGAGTCGCGTTCAATCTATGTGCTGGTACGCCAGCGGCGGCCTGAACTGGTGATCGACGTCAAGCAGTTGcc                                                           >  7:260769/1‑95 (MQ=255)
                                                                                                      aggGTACAGAACCGGGAGACCAGCGAGTCGCGTTCAATCTATGTGCTGGTACGCCAGCGGCGGCCTGAACTGGTGATCGACGTCAAGCAGTTGcc                                                           <  8:260769/95‑1 (MQ=255)
                                                                                                      aggGTACAGAACCGGGAGACCAGCGAGTCGCGTTCAATCTATGTGCTGGTACGCCAGCGGCGGCCTGAACTGGTGATCGACGTCAAGCAGTTGCCCGATGGCGACCTGCAGTTGACAGGTATGTTCGATGGGGTAgg                 <  6:162597/137‑1 (MQ=255)
                                                                                                              caaCCGGGAGACCAGCGAGTCGCGTTCAATCTCTGTGCTGGTACGCCAGCGGCGGCCTGAACTGGTGATCGACGTCAAGCAGTTGCCCGATGGCGACCTGCAGTTTACATGTATTTTCGATGTGGTAGGccagccag         <  1:197999/136‑1 (MQ=255)
                                                                                                                aCCGGGAGACCAGCGAGTCGCGGTGAATATATTTGCTGGTTCGCCAGCGGCGGCCGGCACTGGGGATCGACGGCAAGCAGTTGCCCGATGGCGCCCGGCAGTTGACAGGGTTGTTCGATGGGGTAGgccagccaccc       >  6:136588/1‑134 (MQ=255)
                                                                                                                   gggAGACCAGCGAGTCGCGTTCAATCTATGTGCTGGTACGCCAGCGGCGGCCTGAACTGGTGATCGACGTCAAGCAGTTGCCCGATGGCGACCTGCAGTTGACAGGTATGTTCGATGGGGTAGGCCAGCCAGACAcc    <  6:73157/137‑1 (MQ=255)
                                                                                                                     gagaCCAGCGAGTCGCGTTCAATCAATTTGGTGGTACGCCCCCGGCCGCCGGGACTGGTGATCGACCTCAAGCACTTTGCCGATGGCGCCCAGCAGTTGACAGGTATGTTGGATGGGGTAGGCCCGCCACACCCCCa  >  4:151035/1‑137 (MQ=255)
                                                                                                                                      |                                                                                                                        
GAAGGACGGGGCTCGGTCGTCGCTGAGGGCGATGGCAGGCGTTGCACCTATACCGCAGGCAGTAAGGTGGATGCCAACGTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGTTCAATCTATGTGCTGGTACGCCAGCGGCGGCCTGAACTGGTGATCGACGTCAAGCAGTTGCCCGATGGCGACCTGCAGTTGACAGGTATGTTCGATGGGGTAGGCCAGCCAGACACCCA  >  KT2440/2113761‑2114014

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: