Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1223,063,5210TC32.3% 38.1 / 12.6 40intergenic (+32/‑322)adh_2/pqqD_2Aldehyde dehydrogenase/PqqA binding protein
Reads supporting (aligned to +/- strand):  ref base T (8/13);  new base C (2/8);  total (18/22)
Fisher's exact test for biased strand distribution p-value = 4.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGGTGTTGGCCGTGAGACCCACAAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGC‑TTGCCCGCTCCCACA‑G‑GTTCACCGCTGGTTCCAGAACCGGTG‑AAGAACCTGTGGGGGCGGGCGAG‑CCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGT  >  ACB122/3063390‑3063654
                                                                                                                                    |                                                                                                                                         
cGGTGTTGGCCGTGAGACCCACAAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGC‑TTGCCCCCAccccc                                                                                                                                      >  3:372939/1‑135 (MQ=255)
 ggTGTTGGCCGTGAGACCCACAAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGC‑TTGCCCGCACCcacc                                                                                                                                     >  1:6156/1‑136 (MQ=255)
         ccGTGAGACCCACAAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGC‑TTGCCCGCTCCCACA‑G‑GTTCAcc                                                                                                                           <  8:96263/137‑1 (MQ=255)
         ccGTGAGACCCACAAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGC‑TTGCCCGCTCCCACA‑G‑GTTCAcc                                                                                                                           <  6:213577/137‑1 (MQ=255)
         ccGTGAGACCCACAAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGC‑TTGCCCGCTCCCACA‑G‑GTTCAcc                                                                                                                           <  2:333775/137‑1 (MQ=255)
               gaCCCACAAGATGATGCTTGACCATTATCAGCAGACCAACAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGC‑TTGCCCGCTCCCACA‑G‑GTTCACCGCTGGt                                                                                                                     <  5:330606/137‑1 (MQ=255)
                aCCCACAAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGC‑TTGCCCGCTCCCACA‑G‑GTTCACCGCTgggt                                                                                                                    >  2:28318/1‑135 (MQ=255)
                    acaAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGC‑TTGCCCGCTCCCACA‑G‑GTTCACCGCTGGTTCCAg                                                                                                                <  7:96124/137‑1 (MQ=255)
                     caAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGC‑TTGCCCGCTCCCACA‑G‑GTTCACCGCTGGTTCCAGa                                                                                                               <  4:33650/137‑1 (MQ=255)
                     caAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGC‑TTGCCCGCTCCCACA‑G‑GTTCACCGCTGGTTCCAGa                                                                                                               <  6:19000/137‑1 (MQ=255)
                      aaGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGC‑TTGCCCGCTCCCACA‑G‑GTTACCCGCTGGGTGCAgcc                                                                                                              >  6:339027/1‑135 (MQ=18)
                           gatgCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGC‑TTGCCCGCTCCCACA‑G‑GTTCACCGCTGGTTCCAGAACCGGt                                                                                                         <  4:111141/137‑1 (MQ=255)
                               cTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCCGGC‑TTTCCCCCACCCCCA‑A‑GTTTACCCCGAGTTCCAGAACCCGGg‑aaa                                                                                                    >  5:340508/1‑136 (MQ=18)
                                     cATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGC‑TTGCCCGCTCCCACA‑C‑GTTAACCGCTGGTTCCAGAACCGGTG‑AAGAACCtg                                                                                              >  8:93030/1‑137 (MQ=255)
                                          tCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGC‑GTGCCCCCCCCCCCC‑C‑GGTTACCCCGGGGTGCAGGACCCGTG‑GAGAAACCGTgggg                                                                                         >  3:50809/1‑137 (MQ=14)
                                           cagcagACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGC‑TTGCCCGCGCCCACT‑G‑GTTCACCGCTGGGTCCAGAACCGGTG‑AAGAACCTGTggggg                                                                                        >  2:100887/1‑137 (MQ=255)
                                                    aaGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGC‑TTGCCCGCTCCCACA‑G‑GTTCAcc                                                                                                                           <  8:54069/94‑1 (MQ=255)
                                                    aaGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCCCGGGC‑TTGCCCGCTCCCACA‑G‑GTTCAcc                                                                                                                           >  7:54069/1‑94 (MQ=255)
                                                              tgGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGC‑TTTCCCCCGCCCACACG‑GTTAACC‑CGGGGTGGCCAACCCGGG‑TAGAACCGGTTGGGGCGGGCGCG‑CCCCCcaaaag                                                                    >  1:5025/1‑135 (MQ=12)
                                                               ggTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGC‑TTTCCCCCACCCACA‑C‑GTTGAACGCAGGGTGCAGAAACCGGG‑AAGAACCAGTGGGGGCGGGCGCG‑CCCGCAAAAcgg                                                                   >  5:149958/1‑137 (MQ=21)
                                                                gTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTTGCCCCAACCCCCA‑‑‑GTTAAACCCTGGTTGCAGAAACCGTG‑AAAAAACCCTTGGGGCCGGCGCC‑CCCCCAAAACGGc                                                                  >  7:30004/1‑137 (MQ=9)
                                                                gTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCATTGCCC‑CCCCCCCA‑GCGTTACCCCTTGGTTGATCAACCGG‑G‑GAGAACCGGTTGGGGGGGGCGCGCCCCCCAAAAggg                                                                   >  5:365113/1‑137 (MQ=9)
                                                                 tGAGCTACCACATCAACCCGCTGGGCTTCTTCTAACCCGCGTCGCGGCCTTCGCGGGC‑TGGCCCGCTCCCACC‑G‑GTGCACCGCGGGTTCCAGAACGGGTG‑AAGAACCTGTGGGGGCGGGCGAG‑CCCGCAAAAGGGcc                                                                 <  7:375922/137‑1 (MQ=255)
                                                                  tagatACCACATCAACCCGCGGGGCTCATTCAAAACCGCGGCGGGGCCTTCGCGGGC‑TTCCCCGCCCCCACC‑G‑TTCCCCCGCTGGTCCCACAACCGGTC‑AACAACCTGGGGGGGGGGGCGAG‑CCCGCAAAAGGGCCg                                                                <  4:222884/133‑1 (MQ=255)
                                                                  gAGCAACGCCAAAAAGCCGCTGGGCTTCTTATAACCCGCGGCGCGGCCTCCGCGGGC‑TGGCCCGCCCCCACC‑G‑GTTCACCGCGGGTTCCACAACCGGTC‑AATAACCTGGGGGGGCGGGCGAG‑CCCGCAAAAGGGCCg                                                                <  2:476/137‑1 (MQ=255)
                                                                    gCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGC‑TTTCCCCCGCCCACA‑C‑GGTGACCGCGGGGTTCAGCACCCGGG‑AAAAAACAGTGGGGGCGGGGGCG‑CCCGCCAAAAGGCCCGc                                                              >  6:2131/1‑137 (MQ=16)
                                                                       atgACATCAATCCGCGGGGCTGCTACAAACCCGCGGCGCGGCCTTCGCGGGC‑TGGCCCGCCCCCACA‑G‑GTCCCCCGCTGGTCCCCCACCCGGTT‑AACAACCTGTGGGGGCGGGCGAG‑CCCGCAAAAGGGCCGGcaca                                                           <  1:231784/135‑1 (MQ=255)
                                                                           aaTCAATCCGCTGGGCTTCTTATAACCCGCGGCGCGGCCTCCGCGGGC‑TTGCCCGCTCCCACA‑G‑GTCCACCGCTGGTACCACACCCGGGC‑AAGAACCTGTGGGGGCGGGCGAG‑CCCGCAAAAGGGCCGGCAcatgca                                                       <  7:174196/136‑1 (MQ=25)
                                                                            atcacaccGCGGGGCTGCTTATAACCCGCGCCGCGGCCTCCGCGGGC‑TGGCCCGCCCCCCCC‑G‑GTGCACCGCGGGTCCCACACCCGGTC‑AACAACCTGTGGGGGCGGGCGAG‑CCCGCAAAAGGGCCGGCAcatgcat                                                      <  4:211646/131‑1 (MQ=255)
                                                                            agaaccccGCGGGGCTTCTACAAACCCGGGCCGGGGCTTGGGCGGGC‑TGCCCCGCCCCCACC‑G‑GGCCCCGGCGGGTCCCACACCCGGTA‑AACAACGTGGGGGGGGGGGCGAG‑CCCGCAAAAGGGCCGGCAcatgcat                                                      <  4:50809/131‑1 (MQ=255)
                                                                                 cccGGTGGGGTTATTATAACCCGCGGCGCGGCCTTCGCGGGC‑TGCCCCGCGCCCCCC‑G‑GTCCACCGCTGGTACAACACCCGGTT‑AATAACGTGGGGGGGGGGGCGAG‑CCCGCAAAAGGGCCGGCACATGCATGGCCt                                                 <  4:35321/136‑1 (MQ=255)
                                                                                  ccGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGC‑TTTCCCCCTCCCACA‑C‑GTTGACCGCTGGTTCCAGAACCGGGG‑AAGCACCAGTGGGGGCGGGCGCG‑CCCGCCAAATGGCCCGCCCATGCCTGGccga                                                >  4:371003/1‑135 (MQ=255)
                                                                                         gCTTCTTCTAAACGGCGGGGCGGCGTTCGCGGGC‑TTGCCCGGCCCCCCA‑G‑GTTCACCGCTGTTTCCACAACCGGTT‑AACAACGTGGGGGGGCGGGCGAG‑CCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTAc                                         <  2:86166/137‑1 (MQ=255)
                                                                                                           cgcgGCCTTCGCGGGC‑TTGCCCCCTCCCACA‑G‑GGTGACCGCTGGGTGCAGAACCGGGG‑AAGAACCGGTGGGGGCGGGCGCG‑CCCGCAAAACGGCCCGCCCAAGGCTGGCCGACCCGGACCACCGCCCCCCCCCAGcg                       >  2:51891/1‑136 (MQ=255)
                                                                                                               cccTGCGGGGGC‑GTGCCCGCCCCCCCT‑G‑GTTCACCGGTGGTACCAGAACGGGTT‑AAGAACCGGTGGGGGCGGGCGAG‑CCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGtt                   <  5:293031/136‑1 (MQ=255)
                                                                                                                    cgggggC‑TGGCCCGCTCCCACA‑G‑GTGCCCGGCTGTTTCCAGACCCGGTG‑AACAACCTGTGGGGGCGGGCGAG‑CCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCaaa              <  7:238786/134‑1 (MQ=255)
                                                                                                                     ccgGGC‑TGGCCCCCCCCCACC‑G‑GTCCCCCGCTGGTCCCA‑AACCGGGGAAACAACGGGGGGGGGCGGGCGAG‑CCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAg             <  3:153244/136‑1 (MQ=255)
                                                                                                                          ggttGCCGGCTCCCACA‑G‑GTTCACGGGGGGTTCCAGACCCGGTG‑AAGAACCTGTGGGGGCGGGCGAG‑CCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGc         <  6:86256/135‑1 (MQ=255)
                                                                                                                              gCCCCCTCCCACA‑C‑GTTGAACCCTGGGTGCAGAACCCGGG‑AAGAACCTGTGGGGGCGGGCGCG‑CCCGCAAAAAGGCCCGCCCAATCCTGGCCGACCCGGACCAACCCCCCCCCCCAGCGCGTTGCAAAGTAGCaata     >  6:293031/1‑134 (MQ=255)
                                                                                                                                 cGCTCCCCCA‑G‑GTTGACCCCGGGGTCCAGAACCGGTG‑AAGAACCTGTGGGGGCGGGCGCG‑CCCGCAAAAGGGCCGGCCCAAGCATGGCCTACCCGGACCAACGCCCCCCGCAACCGCGTTCCAAAGTAGCATTCGGt  >  8:206007/1‑137 (MQ=255)
                                                                                                                                    |                                                                                                                                         
CGGTGTTGGCCGTGAGACCCACAAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGC‑TTGCCCGCTCCCACA‑G‑GTTCACCGCTGGTTCCAGAACCGGTG‑AAGAACCTGTGGGGGCGGGCGAG‑CCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGT  >  ACB122/3063390‑3063654

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