Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1432,113,8950TG35.5% 36.6 / 11.5 31R663R (CGT→CGGCJ019_02007hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (3/17);  new base G (10/1);  total (13/18)
Fisher's exact test for biased strand distribution p-value = 6.17e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AAGAAGGACGGGGCTCGGTCGTCGCTGAGGGCGATGGCAGGCGTTGCACCTATACCGCAGGCAGTAAGGTGGATGCCAACGTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGTTCAATCTATGTGCTGGTACGCCAGCGGCGGCCTGAACTGGTGATCGACGTCAAGCAGTTGCCCGATGGCGACCTGCAGTTGACAGGTATGTTCGATGGGGTAGGCCAGCCAGACACCCAGTGGGCATTGCCAAT  >  ACB143/2113759‑2114029
                                                                                                                                        |                                                                                                                                      
aagaagGACGGGGCTCGGTCGTCGCTGAGGGCGATGGCAGGCGTTGCACCTATACCGCAGGCAGTAAGGTGGATGCCAACGTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGt                                                                                                                                        <  1:201840/137‑1 (MQ=255)
            gCTCGGTCGTCGCTGAGGGCGATGGCAGGCGTTGCACCTATACCGCAGGCAGTAAGGTGGATGCCAACGTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGTTCAATCTAtgtg                                                                                                                            >  4:210553/1‑137 (MQ=255)
                                 aTGGCAGGCGTTGCACCTATACCGCAGGCAGTAAGGTGGATGCCAACGTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGTTCAATCTATGTGCTGGTACGCCAGCGGCGGCCt                                                                                                       <  7:146935/137‑1 (MQ=255)
                                         cGTTGCACCTATACCGCAGGCAGTAAGGTGGATGCCAACGTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGTTCAATCTATGTGCTGGTACGCCAGCGGCGGCCTGAACTGGt                                                                                               <  2:107295/137‑1 (MQ=255)
                                         cGTTGCACCTATACCGCAGGCAGTAAGGTGGATGCCAACGTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGTTCAATCTATGTGCTGGTACGCCAGCGGCGGCCTGAACTGGt                                                                                               <  6:304996/137‑1 (MQ=255)
                                         cGTTGCACCTATACCGCAGGCAGAAACGTTGATGCCAACGTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGTTCAATCTATGTGCTGGTACGCCAGCGGCGGCCTGAACTGGt                                                                                               <  2:300309/137‑1 (MQ=255)
                                          gTTGCACCTATACCGCAGGCAGTAAGGTGGATGCCAACGTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGGTGTATATATTTTGTGGTTAGCCAGCGGCGGGCGGGACTGGTg                                                                                              >  3:222043/1‑137 (MQ=255)
                                                   ataCCGCAGGCAGTAAGGTGGATGCCAACGTTGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGTTCAATCTATGTGCTGGTACGCCAGCGGCGGCCTGAACTGGTGATCGACGTc                                                                                     <  5:42396/137‑1 (MQ=255)
                                                      ccGCAGGCAGTAAGGTGGATGCCAACGTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGGTGAATATATGTGCTGGTACGCCAGCGGCGGCCTGCTCTGGTGATGGACGTCAAg                                                                                  >  6:102839/1‑137 (MQ=255)
                                                           ggcagTAAGGTGGATGCCAACGTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGGGGAATATATTTTGTGGTACCGCCGCGGGGGGCGGGACTGGGGTTGGGCGTCAAGGAGtt                                                                             >  5:260441/1‑137 (MQ=255)
                                                            gcagTAAGGTGGATGCCAACGTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGGTGGATCTATGTGCTGGTCAGCCACCGGCGGCCGGAGCTGGGGTTGGACGGCCAGCCGttt                                                                            >  7:143907/1‑136 (MQ=255)
                                                               gTAAGGTGGATGCCAACGTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGGTCAATCTATGTGCTGGTGCGCCCGCGGGGGCCGGAACAGGTGATCGCCGTCAAGCAGTTGccc                                                                         >  6:297613/1‑137 (MQ=255)
                                                                        aTGCCAACGTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGTTCAATCTATGTGCTGGTACGCCAGCGGCGGCCTGAACTGGTGATCGACGTCAAGCAGTTGCCCGATGGCGAc                                                                <  8:111928/137‑1 (MQ=255)
                                                                        aTGCCAACGTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGTTCAATCTATGTGCTGGTACGCCAGCGGCGGCCTGAACTGGTGATCGACGTCAAGCAGTTGCCCAATGGCGAc                                                                <  3:210553/137‑1 (MQ=255)
                                                                           ccAACGTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGGTGAATCTAATTGCTGGTACGCCACCGGGGGGCGGCACTGGGGATCGGCGGCAACCAGTTGCCCGAGGGCGACCTg                                                             >  6:321920/1‑137 (MQ=255)
                                                                             aaCGTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGGTGAATCTATGTGGTTGTTAGCCAGCGGCGGCCTGAACTGGTGATCGACGTCAAGCAGTTGGCCGCTGGCGCGCTGcc                                                           >  4:348016/1‑136 (MQ=255)
                                                                               ggTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGTTCAATCTATGTGCTGGTACGCCAGCGGCGGCCTGAACTGGTGATCGACGTCAAGCAGTTGCCCGATGGCGACCTGCAGt                                                         <  3:295844/136‑1 (MQ=255)
                                                                                    ccTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGTTCAATCTATGTGCTGGTACGCCAGCGGCGGCCTGAACTGGTGATCGACGTCAAGCAGTTGCCCGATGGCGACCTGCAGTTGACa                                                    <  6:2996/137‑1 (MQ=255)
                                                                                    ccTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGTTCAATCTATGTGCTGGTACGCCAGCGGCGGCCTGAACTGGTGATCGACGTCAAGCAGTTGCCCGATGGCGACCTGCAGTTGACa                                                    <  3:180954/137‑1 (MQ=255)
                                                                                      tACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGTTCAATCTATGTGCTGGTACGCCAGCGGCGGCCTGAACTGGTGATCGACGTCAAGCAGTTGCCCGATGGCGACCTGCAGTTGACAgg                                                  <  3:41703/137‑1 (MQ=255)
                                                                                           gCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGGTGAATCTATTTGCTGGTACGCCAGCGGCGGCCAGAACATGTGTTCGACCTCAACCACTTGCCCGATGGCGAGCTGCAGTTGACAGGttatt                                             >  4:181159/1‑133 (MQ=255)
                                                                                                      gcaggGTACAGAACCGGGAGACCAGCGAGTCGCGTTCAATCTATGTGCTGGTACGCCAGCGGCGGCCTGAACTGGTGATCGACGTCAAGCAGTTGCCCGATGGCGACCTGCAGTTGACAGGTATGTTCGATGGGGTa                                  <  1:163009/136‑1 (MQ=255)
                                                                                                          ggTACAGAACCGGGAGACCAGCGAGTCGCGTTCAATCTATGTGCTGGTACGCCAGCGGCGGCCTGAACTGGTGATCGACGTCAAGCAGTTGCCCGATGGCGACCTGCAGTTGACATGTATGTTAGATGGGGTAGgcc                              <  7:73125/137‑1 (MQ=255)
                                                                                                              cAGAACCGGGAGACCAGCGAGTCGCGTTCAATCTATGTGCTGGTACGCCAGCGGCGGCCTGAACTGGTGATCGACGTCAAGCAGTTGCCCGATGGCGACCTGCAGTTGACAGGTATGTTCGATGGGGTAGgccagcc                          <  8:49672/137‑1 (MQ=255)
                                                                                                               aGAACCGGGAGACCAGCGAGTCGCGGTCAATCTATGTGCTGGTACGCCAGCGGCGGCCTGAACTGGTGATCGACGTCAAGCAGTTGCCCGATGGCGACCTGCAGTTGACAGGTATGTTCGATGGGGTAGgccagcca                         <  8:309249/137‑1 (MQ=255)
                                                                                                                  aCCGGGAGACCAGCGAGTCGCGGTGAATCTATGTTGTGGTTCGCCACCGGCGGCCGGAACTGGTGATCGACGTCAAGCAATTGCCCGAGGGCGACCTGCACTTGAAAGGTATGTTCGATTGGGTAGGCCCGCCAGac                      >  8:160528/1‑137 (MQ=255)
                                                                                                                  aCCGGGAGACCAGCGAGTCGCGGTGAATCTATGTGCTGGTTCGCCCGCCGCGGCCCGGACTGGGGAGAGACGTCAAGCCGTTGGCCGATGGGGCCCGGCAGTTGTCCGGTATTTTGGGTGGGGTAGGCCCGCCCGac                      >  8:204607/1‑137 (MQ=255)
                                                                                                                   ccGGGAGACCAGCGAGTCGCGTGCAATCTATGTGCTGGTACGCCAGCGGCGGCCTGAACTGGTGATCGACGTCAAGCAGTTGCCCGATGGCGACCTGCAGTTGACAGGTATGTTCGATGGGGTAGGCCAGCCAGaca                     <  7:274309/137‑1 (MQ=255)
                                                                                                                       gagaCCAGCGAGTCGCGTTGAATCTATGTGCTGGTTAGCCAGCGGCGGCCGGAACTGGTGATAGACGTCAAGCAGTTGCCCGATGGCGACCTGCAGTTGACAGGTATGTTAGATTGGGTAGGCCAGCCAGCCACCCa                 >  6:330575/1‑137 (MQ=255)
                                                                                                                       gagaCCAGCGAGTCGCGTTCAATCTATGTGCTGGTACGCCAGCGGCGGCCTGAACTGGTGATCGACGTCAAGCAGTTGCCCGATGGCGACCTGCAGTTGACAGGTATGTTCGATGGGGTAGGCCAGCCAGACACCCa                 >  7:263660/1‑137 (MQ=255)
                                                                                                                             aGCGAGTCGCGTTCAATCTATGTGCTGGTACGCCAGCGGCGGCCTGAACTGGTGATCGACGTCAAGCAGTTGCCCGATGGCGACCCGCAGTTGACAGGTATGTTCGATGGGGTAGGCCAGCCAGACACCCAGTGGGc           <  1:88185/137‑1 (MQ=255)
                                                                                                                                      cgTTCAATCTATGTGCTGGTACGCCAGCGGCGGCCTGAACTGGTGATCGACGTCAAGCAGTTGCCCGATGGCGACCTGCAGTTGACAGGTATGTTCGATGGGGTAGGCCAGCCAGACACCCAGTGGGCATTGCCAAt  <  4:56185/137‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                      
AAGAAGGACGGGGCTCGGTCGTCGCTGAGGGCGATGGCAGGCGTTGCACCTATACCGCAGGCAGTAAGGTGGATGCCAACGTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGTTCAATCTATGTGCTGGTACGCCAGCGGCGGCCTGAACTGGTGATCGACGTCAAGCAGTTGCCCGATGGCGACCTGCAGTTGACAGGTATGTTCGATGGGGTAGGCCAGCCAGACACCCAGTGGGCATTGCCAAT  >  ACB143/2113759‑2114029

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: