Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB122994,2290TG50.0% 32.6 / 12.3 46D52E (GAT→GAG) CJ019_00938hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (11/12);  new base G (0/23);  total (11/35)
Fisher's exact test for biased strand distribution p-value = 2.03e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GAGGAAGTGGTGGATGTGGCCGGCTTGCAGGCCATGCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATGGCGCACAACATGCTCGGGCGTTGCCTGGAGCATGGCTGGGGTGGCGA  >  ACB122/994095‑994351
                                                                                                                                      |                                                                                                                          
gAGGAAGTGGTGGATGTGGCCGGCTTGCAGGCCATGCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGAt                                                                                                                            >  1:123555/1‑135 (MQ=255)
  ggAAGTGGTGGATGTGGCCGGCTTGCAGGCCATGCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATgg                                                                                                                          <  3:125059/135‑1 (MQ=255)
  ggAAGGGGTGGATGTGGCCGGCTGGCAGGCCATGCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATgg                                                                                                                          <  5:63719/135‑1 (MQ=255)
      gtggtgGATGTGGCCGGCTTGCAGGCCATGCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGt                                                                                                                      <  4:203923/135‑1 (MQ=255)
            gATGTGGCCGGCTTGCAGGCCATGCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTggggtg                                                                                                                >  1:9157/1‑131 (MQ=255)
            gATGTGGCCGGCTTGCAGGCCATGCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTgggcgg                                                                                                                >  5:296910/1‑131 (MQ=255)
            gATGTGGCCGGCTTGCAGGCCATGCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTgggcgg                                                                                                                >  1:35039/1‑131 (MQ=255)
            gATGTGGCCGGCTTGCAGGCCATGCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCag                                                                                                                 >  8:42328/1‑133 (MQ=255)
            gATGTGGCCGGCTTGCAGGCCATGCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCaaa                                                                                                                >  1:383234/1‑133 (MQ=255)
              gggggCCGGCTTGCAGGCAAGGCGCGAGGAAACCCCCGGCACGGCCGCCCAGGCGATGCTGGCGGCAGCGGGGCAGGGTGGGGGCGAGGCGCAGTGGCTGCTGGGGCAGACCCTGCGGGAGGGGCGTGCCACCca                                                                                                              <  8:345398/132‑1 (MQ=255)
                 ggCCGGCTTGCAGGCCATGCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGGCTCccgca                                                                                                           >  7:345398/1‑135 (MQ=255)
                    cGGCTTGCAGGCCATGCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCCGCAGGa                                                                                                        >  8:384033/1‑135 (MQ=255)
                    cGGCTTGCAGGCCATGCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGa                                                                                                        >  8:296069/1‑135 (MQ=255)
                                         ggAACCTCCCGGCACGGCCGCCCCGGCGATTCTGGCGGCAGCGGGGCTGGGTGGGGTCGAGGCGCTGTTGCTGCGGGGGCAGATCCTGCTGGAGGGGCGGGGCATCCAGCTGGATGCCACGGTTGCCCGGCGTTg                                                                                   <  8:271501/135‑1 (MQ=255)
                                           aaaCTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATGCAGCAGGGTTCCACCGTTGGCCGGCggttgg                                                                                 >  8:281448/1‑130 (MQ=255)
                                            accgcccGGAACGGCGGCCCAGCGTATGCTGGGGCCAGGGGGGCTGGGGGTGGTCGAGGCGCAGTTGGGGCGGGTGAACATCATGCTGGATGGGGTTGGCATCCTGCAGCATGCCACGGTTGCCCGGCGTTGGtt                                                                                <  6:532432/131‑1 (MQ=255)
                                            accgcccGCCACGGCCGCCCCGGCGATGCTGGCGGCAGGGGGGCAGGGTGGGGGCGAGGCGCAGTGGCTGCTGGGGCAGAGCGTGCTGGAGGGGCGTGGCATCCAGCTGGATGCCACGGTTGCCCGGCGTTGGtt                                                                                <  3:172979/131‑1 (MQ=255)
                                            aaCCCCCGCCACGGCCGCCCCGGGGTGGCGGGCGGTGGGGGGGCTGGGGGGGGGCGGGGCGCTGTTGGTGCTGGGGATGAGCCTGCGGGTGGGGCGGGGCACCCTGCTGGATGCCACGGTTGCCCGGCGTTGGtt                                                                                <  8:24100/135‑1 (MQ=255)
                                               tCCGGGCATGGCCGCCCCGGGGATTGGGGCGGCAGGGGGGCTGGGGGTGGTCGCGGCGATGTTGCTGCGGGGGCAGATCGTGGGGGTTGGGGGGGGAATCCAGCTGGATGCCACGGTTGCCCGGCGTTGGTTCgg                                                                             <  6:509350/135‑1 (MQ=255)
                                                 ccGGCAAGGCCGCCCTGGCGATGGTGGCGGCTGCGGGGCCGGGGGGGTGCGCGGCGCTGTGGCGGCGGGGGCTGAGCCGGCGGGTGGGGCGGGGCATCCAGATGGATGCCACGGTTGCCCGGCGTTGGTTCGGTa                                                                           <  7:158608/135‑1 (MQ=255)
                                                 ccGGCAAGGCCGCCCAGGCTATGGTGGCGGCAGCGGGGCTGGGTGGGGGCGCGGCGCTGTGGCGGCGGGGGCAGATCGTGCTGGTTGGGCGTGCCAACCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTa                                                                           <  2:25349/135‑1 (MQ=255)
                                                    gCAGGGCCGCCAGGGCTATTCTGGCGGCTGGGGGGCCGGGTGGGGGCGAGGCGCTGTGGGTGCGGGGGCTGATCCTGGGGTTGGGGGGGGGCACCCAGCAGCATGCCACGGTTGCCCGGCGTTGGTTCGGTATCg                                                                        <  2:25694/135‑1 (MQ=255)
                                                      ccggCCGCCCCGGCGATTCGGGCGGCCGGGGGGCGGGGTGGGGGCGAGGCGCAGTGGCTGCGGGGGCAGACCCTGCTGGAGGGGCGGGGCACCCCGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCt                                                                      <  1:348602/133‑1 (MQ=255)
                                                      agggCCGCCCTGGCGATTGTGGCGGCAGGGGGGCCGGGGGGGGGCGTGGCGCTGTGGGTGGGGGGGCAGATCCTGCTGGAGGGGCGTGGCACCCCGCCGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCt                                                                      <  8:267180/133‑1 (MQ=255)
                                                      acggCCGCCCCGGGGATGCGGGGGGCAGCGGGGCGGGGGGGGGGCGCGGCGCAGTGGCTGCGGGGGCAGAGCCTGGTGGAGGGGCGGGGCACCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCt                                                                      <  2:58981/133‑1 (MQ=255)
                                                      acggCCGCCCAGGCGAGTCTGGCGGCAGCGGGGCAGGGTGGGGGCGTGGCGCAGTTGCTGCTGGGGCAGATCCTGCGGTATGGGCGTGCCATCCAGCTGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCt                                                                      <  6:210464/133‑1 (MQ=255)
                                                              cctggCGATGCTGGCGGCAGGGGGGCGGGGTGGGGCCCCGGCGCTGTGGGGGGGGGGCCAGAGCCTGCGGGTGGGGCGTGGCACCCAGCTGCATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAggg                                                              <  6:470112/132‑1 (MQ=255)
                                                              cccggCGATGCTGGCGGCCGGGGGGCTGGGGGGGGGCGCGGCGCAGTGGCGGCGGGGGCAGAGCCGGCTGGAGGGGCGGGGCAGCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAggg                                                              <  4:73160/132‑1 (MQ=255)
                                                              ccAGGCTATTCGGGCGGCAGCGGGGCGGGGGGGGGGCGAGGCGCAGTGGCGGCGGGGCCAGATCCTGGTGGAGGGGGGGGCCACCCCCCCGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAggg                                                              <  2:309022/135‑1 (MQ=255)
                                                                   tctttgtGGCGGCAGGGGGGCTGGGGGTGCTCGAGGCGCTGGTGCTGCGGGGGCAGATCCTGCTGTAGGGGCGTGGCATCCCGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCa                                                         <  8:272572/129‑1 (MQ=255)
                                                                      ttCGGGCGGCAGCGGGGCAGGGTGTGGCCGAGGCGCAGTGGCTGCGGGGGCAGAGCCTGCTGGAGGGGCGGGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTg                                                      <  4:474094/135‑1 (MQ=255)
                                                                           gcggcgGGGGGGCCGGGGGTGGGCGGGGCGCTGTGGCGGGGGGGGCTGTTCCTGGTGGCGGGGGGGGGAACCCCGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATg                                                 <  5:98088/135‑1 (MQ=255)
                                                                           gcggcAGGGGGGCAGGGTGTGGGCGAGGCGCATTTGCGCCTGGGGAAGATCGTGCTGGATGGGCGTGGCATACAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATg                                                 <  2:383768/135‑1 (MQ=255)
                                                                                      gctgggTGGGGTCGTGGCGATGTTGGTGCGGGGGATGAGCCGGCGGGTGGGGGGGGGAATCCAGATGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATGGCGCACAACAt                                      <  6:265881/132‑1 (MQ=255)
                                                                                      gCAGGGTGTGGCCGCGGCGCTTTGGGGGTGGGGGAAGAGCCTGGGGGTGGGGCGGGGCAGCCTGCTGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATGGCGCACAACAt                                      <  8:175826/135‑1 (MQ=255)
                                                                                           gtgtGGGCGCGGCGCAGTGGCTGCTGGGGCAGATCCTGCTGGAGGGGCGTGCCACCCAGCCGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATGGCGCACAACATGCTCg                                 <  7:281448/135‑1 (MQ=255)
                                                                                                 gcgcggCGCTGTTGCGGGGGGGGCAGATCCGGGGGGTGGGGGGGGCCTTCCAGCTGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATGGCGCACAACATGCTCGGGCGtt                           <  7:246977/131‑1 (MQ=255)
                                                                                                 ccgcggCGCAGGTGCTGCGGGGGCAGAGCGGGGGGGAGGGGCGGGGCAGCCCGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATGGCGCACAACATGCTCGGGCGtt                           <  7:124615/131‑1 (MQ=255)
                                                                                                      gcgcTGTTGCTGCGGGGGCAGATCGTGCTGGATGGGCGTGGCAGCCAGCTGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATGGCGCACAACATGCTCGGGCGTTGCCTg                      <  2:22110/135‑1 (MQ=255)
                                                                                                       cgcAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGc                                                                                                                        <  3:78682/36‑1 (MQ=255)
                                                                                                       cgcAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGc                                                                                                                        >  4:78682/1‑36 (MQ=255)
                                                                                                           gctgctgcGGGGGCATCTCCGGCTGGAGGGGCGTGCCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATGGCGCACAACATGCTCGGGCGTTGCCTGGAGCa                 <  6:443432/133‑1 (MQ=255)
                                                                                                             tgctgcTGGGGCAGATCCTGCTGGATGGGCGTGGGCTACCGCCGGGTTCCCCGGGTGGCCCGCGGGGGGTGGGGGTCGCGGCCCCGGGTGGCCGGGGCATGGCGCCCCACATTCGGGGGCGTTGCCTGGGGGATg               >  2:404730/1‑135 (MQ=255)
                                                                                                               cggcgggggCCGAGCCTGCGGGAGGGGCGGGGCCGCCAGCTGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATGGCGCACAACATGCTCGGGCGTTGCCTGGAGCATgg              <  7:17207/129‑1 (MQ=255)
                                                                                                                    ggggCAGATCCTGCGGGGTGGGCGTGCCATCCAGTTGCATGCCATGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATGGCGCACAACATGCTCGGGCGTTGCCTGGAGCATGGCTgggg        <  2:383234/135‑1 (MQ=255)
                                                                                                                          gttCCTGCGGTTTGGGGGGGGCACCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATGGCGCACAACATGCTCGGGCGTTGCCTGGAGCATGGCTGGGGTGGCGa  <  7:296069/133‑1 (MQ=255)
                                                                                                                          gagccTGGTGGAGGGGCGTGGCAGCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATGGCGCACAACATGCTCGGGCGTTGCCTGGAGCATGGCTGGGGTGGCGa  <  1:98657/132‑1 (MQ=255)
                                                                                                                                      |                                                                                                                          
GAGGAAGTGGTGGATGTGGCCGGCTTGCAGGCCATGCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATGGCGCACAACATGCTCGGGCGTTGCCTGGAGCATGGCTGGGGTGGCGA  >  ACB122/994095‑994351

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