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.

CGGTGTTGGCCGTGAGACCCACAAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGT  >  ACB122/3063390‑3063654
                                                                                                                                    |                                                                                                                                         
cGGTGTTGGCCGTGAGACCCACAAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCCCAccccc                                                                                                                                      >  3:372939/1‑135 (MQ=255)
 ggTGTTGGCCGTGAGACCCACAAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCACCcacc                                                                                                                                     >  1:6156/1‑136 (MQ=255)
         ccGTGAGACCCACAAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCAcc                                                                                                                           <  8:96263/137‑1 (MQ=255)
         ccGTGAGACCCACAAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCAcc                                                                                                                           <  6:213577/137‑1 (MQ=255)
         ccGTGAGACCCACAAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCAcc                                                                                                                           <  2:333775/137‑1 (MQ=255)
               gaCCCACAAGATGATGCTTGACCATTATCAGCAGACCAACAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGt                                                                                                                     <  5:330606/137‑1 (MQ=255)
                aCCCACAAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTgggt                                                                                                                    >  2:28318/1‑135 (MQ=255)
                    acaAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAg                                                                                                                <  7:96124/137‑1 (MQ=255)
                     caAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGa                                                                                                               <  4:33650/137‑1 (MQ=255)
                     caAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGa                                                                                                               <  6:19000/137‑1 (MQ=255)
                      aaGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTACCCGCTGGGTGCAgcc                                                                                                              >  6:339027/1‑135 (MQ=18)
                           gatgCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGt                                                                                                         <  4:111141/137‑1 (MQ=255)
                               cTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCCGGCTTTCCCCCACCCCCAAGTTTACCCCGAGTTCCAGAACCCGGgaaa                                                                                                    >  5:340508/1‑136 (MQ=18)
                                     cATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACACGTTAACCGCTGGTTCCAGAACCGGTGAAGAACCtg                                                                                              >  8:93030/1‑137 (MQ=255)
                                          tCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCGTGCCCCCCCCCCCCCGGTTACCCCGGGGTGCAGGACCCGTGGAGAAACCGTgggg                                                                                         >  3:50809/1‑137 (MQ=14)
                                           cagcagACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCGCCCACTGGTTCACCGCTGGGTCCAGAACCGGTGAAGAACCTGTggggg                                                                                        >  2:100887/1‑137 (MQ=255)
                                                    aaGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCAcc                                                                                                                           <  8:54069/94‑1 (MQ=255)
                                                    aaGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCCCGGGCTTGCCCGCTCCCACAGGTTCAcc                                                                                                                           >  7:54069/1‑94 (MQ=255)
                                                              tgGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTTCCCCCGCCCACACGGTTAACCCGGGGTGGCCAACCCGGGTAGAACCGGTTGGGGCGGGCGCGCCCCCcaaaag                                                                    >  1:5025/1‑135 (MQ=12)
                                                               ggTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTTCCCCCACCCACACGTTGAACGCAGGGTGCAGAAACCGGGAAGAACCAGTGGGGGCGGGCGCGCCCGCAAAAcgg                                                                   >  5:149958/1‑137 (MQ=21)
                                                                gTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTTGCCCCAACCCCCAGTTAAACCCTGGTTGCAGAAACCGTGAAAAAACCCTTGGGGCCGGCGCCCCCCCAAAACGGc                                                                  >  7:30004/1‑137 (MQ=9)
                                                                gTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCATTGCCCCCCCCCCAGCGTTACCCCTTGGTTGATCAACCGGGGAGAACCGGTTGGGGGGGGCGCGCCCCCCAAAAggg                                                                   >  5:365113/1‑137 (MQ=9)
                                                                 tGAGCTACCACATCAACCCGCTGGGCTTCTTCTAACCCGCGTCGCGGCCTTCGCGGGCTGGCCCGCTCCCACCGGTGCACCGCGGGTTCCAGAACGGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGcc                                                                 <  7:375922/137‑1 (MQ=255)
                                                                  tagatACCACATCAACCCGCGGGGCTCATTCAAAACCGCGGCGGGGCCTTCGCGGGCTTCCCCGCCCCCACCGTTCCCCCGCTGGTCCCACAACCGGTCAACAACCTGGGGGGGGGGGCGAGCCCGCAAAAGGGCCg                                                                <  4:222884/133‑1 (MQ=255)
                                                                  gAGCAACGCCAAAAAGCCGCTGGGCTTCTTATAACCCGCGGCGCGGCCTCCGCGGGCTGGCCCGCCCCCACCGGTTCACCGCGGGTTCCACAACCGGTCAATAACCTGGGGGGGCGGGCGAGCCCGCAAAAGGGCCg                                                                <  2:476/137‑1 (MQ=255)
                                                                    gCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTTCCCCCGCCCACACGGTGACCGCGGGGTTCAGCACCCGGGAAAAAACAGTGGGGGCGGGGGCGCCCGCCAAAAGGCCCGc                                                              >  6:2131/1‑137 (MQ=16)
                                                                       atgACATCAATCCGCGGGGCTGCTACAAACCCGCGGCGCGGCCTTCGCGGGCTGGCCCGCCCCCACAGGTCCCCCGCTGGTCCCCCACCCGGTTAACAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGcaca                                                           <  1:231784/135‑1 (MQ=255)
                                                                           aaTCAATCCGCTGGGCTTCTTATAACCCGCGGCGCGGCCTCCGCGGGCTTGCCCGCTCCCACAGGTCCACCGCTGGTACCACACCCGGGCAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCAcatgca                                                       <  7:174196/136‑1 (MQ=25)
                                                                            atcacaccGCGGGGCTGCTTATAACCCGCGCCGCGGCCTCCGCGGGCTGGCCCGCCCCCCCCGGTGCACCGCGGGTCCCACACCCGGTCAACAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCAcatgcat                                                      <  4:211646/131‑1 (MQ=255)
                                                                            agaaccccGCGGGGCTTCTACAAACCCGGGCCGGGGCTTGGGCGGGCTGCCCCGCCCCCACCGGGCCCCGGCGGGTCCCACACCCGGTAAACAACGTGGGGGGGGGGGCGAGCCCGCAAAAGGGCCGGCAcatgcat                                                      <  4:50809/131‑1 (MQ=255)
                                                                                 cccGGTGGGGTTATTATAACCCGCGGCGCGGCCTTCGCGGGCTGCCCCGCGCCCCCCGGTCCACCGCTGGTACAACACCCGGTTAATAACGTGGGGGGGGGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCt                                                 <  4:35321/136‑1 (MQ=255)
                                                                                  ccGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTTCCCCCTCCCACACGTTGACCGCTGGTTCCAGAACCGGGGAAGCACCAGTGGGGGCGGGCGCGCCCGCCAAATGGCCCGCCCATGCCTGGccga                                                >  4:371003/1‑135 (MQ=255)
                                                                                         gCTTCTTCTAAACGGCGGGGCGGCGTTCGCGGGCTTGCCCGGCCCCCCAGGTTCACCGCTGTTTCCACAACCGGTTAACAACGTGGGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTAc                                         <  2:86166/137‑1 (MQ=255)
                                                                                                           cgcgGCCTTCGCGGGCTTGCCCCCTCCCACAGGGTGACCGCTGGGTGCAGAACCGGGGAAGAACCGGTGGGGGCGGGCGCGCCCGCAAAACGGCCCGCCCAAGGCTGGCCGACCCGGACCACCGCCCCCCCCCAGcg                       >  2:51891/1‑136 (MQ=255)
                                                                                                               cccTGCGGGGGCGTGCCCGCCCCCCCTGGTTCACCGGTGGTACCAGAACGGGTTAAGAACCGGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGtt                   <  5:293031/136‑1 (MQ=255)
                                                                                                                    cgggggCTGGCCCGCTCCCACAGGTGCCCGGCTGTTTCCAGACCCGGTGAACAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCaaa              <  7:238786/134‑1 (MQ=255)
                                                                                                                     ccgGGCTGGCCCCCCCCCACCGGTCCCCCGCTGGTCCCAAACCGGGGAAACAACGGGGGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAg             <  3:153244/136‑1 (MQ=255)
                                                                                                                          ggttGCCGGCTCCCACAGGTTCACGGGGGGTTCCAGACCCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGc         <  6:86256/135‑1 (MQ=255)
                                                                                                                              gCCCCCTCCCACACGTTGAACCCTGGGTGCAGAACCCGGGAAGAACCTGTGGGGGCGGGCGCGCCCGCAAAAAGGCCCGCCCAATCCTGGCCGACCCGGACCAACCCCCCCCCCCAGCGCGTTGCAAAGTAGCaata     >  6:293031/1‑134 (MQ=255)
                                                                                                                                 cGCTCCCCCAGGTTGACCCCGGGGTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGCGCCCGCAAAAGGGCCGGCCCAAGCATGGCCTACCCGGACCAACGCCCCCCGCAACCGCGTTCCAAAGTAGCATTCGGt  >  8:206007/1‑137 (MQ=255)
                                                                                                                                    |                                                                                                                                         
CGGTGTTGGCCGTGAGACCCACAAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGT  >  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: