Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1431,721,8780AC45.0% 18.8 / 14.0 40G95G (GGT→GGGCJ019_01617hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (19/3);  new base C (2/16);  total (21/19)
Fisher's exact test for biased strand distribution p-value = 2.84e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTCCATTTCCCTACAGCCTTGGGTAACCAAGT  >  ACB143/1721747‑1722012
                                                                                                                                   |                                                                                                                                      
gTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCCgg                                                                                                                                   >  7:226747/1‑137 (MQ=255)
gTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACCCCAGGTGATTTCCTCTTCACCGTCGGCGCTTATGCGGACCCAGCGGTCAGCAAATTCTTCCCCGTATTACGAGTACCACCCACCCgg                                                                                                                                   >  1:103870/1‑137 (MQ=255)
    aCCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTACTCGACCCAGCCACCCGgcgcg                                                                                                                               >  5:335078/1‑137 (MQ=255)
        ataccTACTGAAACTACCACGCCTGCCTTTTCATTATGTACGCCTGGTTATTTCCTCTTCCCGGGCGGCGCTGATGCGGACCCCGCGGGACGCTTTTTCTCCCCCGTATTCCTCGCCCCAGCCCCCCGGCgcgcagc                                                                                                                           <  6:244841/134‑1 (MQ=255)
             tACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCCCCCGGCGCGCAGCGCAca                                                                                                                      >  7:17127/1‑136 (MQ=255)
                          aaGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGGCGCGCAGCGCACTTCCACACCCAgcc                                                                                                         >  3:129119/1‑136 (MQ=255)
                            gCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCCGCCACCCGGCGCGCAGCGCACTTCCACACCCAGcgcg                                                                                                       >  8:278440/1‑137 (MQ=255)
                               tGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGGCGCGCAGCGCACTTCCCCACCCAGCCCGGcc                                                                                                    >  5:246188/1‑136 (MQ=255)
                                          gatgTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGcgcg                                                                                         >  7:193802/1‑137 (MQ=255)
                                              tACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGgcg                                                                                     >  6:252468/1‑137 (MQ=255)
                                                        gATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTGAGCATCTTCTTCCCCGTCTTCCTAGACCAAGCCACCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGAGGCGAGTGCGCAGGCAACGt                                                                           >  3:166419/1‑137 (MQ=255)
                                                        gATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCCCCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGt                                                                           >  2:31304/1‑137 (MQ=255)
                                                         aTTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGGCGCGCAGCGCACTTCCACACCCCGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTc                                                                          >  5:122565/1‑137 (MQ=255)
                                                         aTTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTc                                                                          >  1:184572/1‑137 (MQ=255)
                                                                 ttCACCGTCGGCGCTTATGCGGCCCCAGCGGTAAGAATCTTCTCCCCCGTCTACCTCCACCCCCCCCCCCGGCGCGCCGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGccc                                                                  <  4:304444/137‑1 (MQ=255)
                                                                         cGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCCGCCACCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCggc                                                          >  4:226814/1‑135 (MQ=255)
                                                                             gcTGTTGCGCCCCCCGCGGTCAGAATGTTCTTCCCCGTATCCCCCCACCCCCCCCCCGGGCGCGCCGCCCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGt                                                      <  2:114471/137‑1 (MQ=255)
                                                                             gcTGATGCGGACCCAGGGGTCAGCATCTTATTCCCCGTCTCCCGCGACCCCCCCCCCCGGCGCGCCGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGt                                                      <  4:265878/137‑1 (MQ=255)
                                                                             gcTGATGCGCACCCCGCGGTCAGAATCTTCTTCCCCGTCTTCCTCCCCCCACCCCCCCGGCGCGCCGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGt                                                      <  8:271456/137‑1 (MQ=255)
                                                                                 ttGCGGACCCTGCGGTCAGCATGTTCTACCCCGTCTTCCTCCACCCCGCCACCCGGCGCGCCGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCt                                                  <  6:358465/136‑1 (MQ=255)
                                                                                 aggcgcACCCAGCGGTCATCATCTTCTGCCCCGTCTTCCTCCCCCCACCCACCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCt                                                  <  6:69111/135‑1 (MQ=255)
                                                                                  ggcgcACCCCGCGGTCAGCATCTTCTGCCCGGTCTTCCTCGCCCCCGCCCCCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCtt                                                 <  5:151759/136‑1 (MQ=255)
                                                                                           cccGGTCACCATATTATCCCCCGTCTCCCTCGCCCCCCCCCCCCGGCGCCCCGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATAc                                        <  8:107704/135‑1 (MQ=255)
                                                                                              ggTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAgg                                     >  4:242147/1‑137 (MQ=255)
                                                                                                   gCATGTTCTTCCCCGTCTTCCCCCACCCCCCCCCCCGGCGCGCCGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGtt                                <  4:338734/137‑1 (MQ=255)
                                                                                                    aaTCTTCTTCCCCGTCTACCCCGACCCACCCCCCCGGCGCCCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTc                               <  3:72173/136‑1 (MQ=255)
                                                                                                    aaTCTTCTACCCCGTCTGCCGCGCCCCCCCCCCCCGGCGCGCAGCCCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTc                               <  4:156871/136‑1 (MQ=255)
                                                                                                      tcttctTCCCCGTCTTCCTCGACCCAGCCACCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTCCa                             >  8:257331/1‑137 (MQ=255)
                                                                                                      tcttatTCCCCGTCTTCCTCGACCCAGCCACCCGGCGCGCAGCGCACTTACACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACCTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTcaa                             >  6:328523/1‑135 (MQ=255)
                                                                                                           ttcCCCGTCTACCGCCACCCACCCCCCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTCCATTTcc                        <  5:52614/137‑1 (MQ=255)
                                                                                                           tgccccGGCTTCCCCGACCCCCCCCCCGGGCGCGCACCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTCCATTTcc                        <  6:224536/135‑1 (MQ=255)
                                                                                                                 gTCTTCCTCGACCCAGCCACCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTCCATTTACCTACaa                  >  5:224536/1‑136 (MQ=255)
                                                                                                                  tCTTCCTCGCCCCAGCCCCCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTCCATTTCCCTACAGc                 <  8:285856/137‑1 (MQ=255)
                                                                                                                    ttCCTCCCCCCAGCCCCCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTCCATTTCCCTACAGCCt               <  1:31304/137‑1 (MQ=255)
                                                                                                                     tCCTCGACCCAGCCACCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTCCATTTCCCTACAGCCtt              >  2:16784/1‑137 (MQ=255)
                                                                                                                       cTCGACCCAGCCACCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTCCATTTCCCTACAGCCTTgg            >  1:8002/1‑137 (MQ=255)
                                                                                                                             cccgCCCCCCGGCGCGCCGCCCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTCCATTTCCCTACAGCCTTGGGTAAcc      <  6:184758/134‑1 (MQ=255)
                                                                                                                             ccACCCCCCCGGCGCGCCGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTCCATTTCCCTACAGCCTTGGGTAAcc      <  7:271962/137‑1 (MQ=255)
                                                                                                                               aGCCACCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTCCATTTCCCTACAGCCTTGGGTAACCaa    >  8:347871/1‑137 (MQ=255)
                                                                                                                                 ccACCCGGCGCGCCGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTCCATTTCCCTACAGCCTTGGGTAACCAAGt  <  7:278440/137‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                      
GTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTCCATTTCCCTACAGCCTTGGGTAACCAAGT  >  ACB143/1721747‑1722012

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: