Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6064,393,0610TG20.9% 164.3 / 14.1 67V211G (GTG→GGG) aidBisovaleryl CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (5/48);  new base G (14/0);  total (19/48)
Fisher's exact test for biased strand distribution p-value = 1.19e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AAGCAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCTGT  >  REL606/4392988‑4393135
                                                                          |                                                                          
aaGCAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCtggtgg                                                                          >  1:551323/1‑76 (MQ=255)
aaGCAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCtggggg                                                                          >  2:1501800/1‑76 (MQ=255)
 aGCAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTggg                                                                         <  2:880864/76‑1 (MQ=255)
  gCAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTgggg                                                                        <  2:1078756/76‑1 (MQ=255)
  gCAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTgggg                                                                        <  1:1366373/76‑1 (MQ=255)
   cAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGggggggg                                                                       >  2:1065177/1‑75 (MQ=255)
    aGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGgggggggg                                                                      >  2:625804/1‑74 (MQ=255)
    aGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCa                                                                      <  1:158810/76‑1 (MQ=255)
      ggCGGTTCCGATGTTATGCGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATa                                                                    <  1:1438988/76‑1 (MQ=255)
       gCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATaa                                                                   <  1:693529/76‑1 (MQ=255)
        cGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATaaa                                                                  <  1:321790/76‑1 (MQ=255)
         ggTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAAt                                                                 <  2:783425/76‑1 (MQ=255)
         ggTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAAt                                                                 <  1:1192758/76‑1 (MQ=255)
         ggTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGGTAAAt                                                                 >  2:560912/1‑76 (MQ=255)
          gTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGGTAAATg                                                                >  1:1087047/1‑76 (MQ=255)
            tCCGATGTTGTGAGCAACACCACCCGTGCAGAGGGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGt                                                              >  2:1317578/1‑76 (MQ=255)
            tCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGt                                                              <  1:1107079/76‑1 (MQ=255)
             ccGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGtt                                                             <  2:765727/76‑1 (MQ=255)
              cGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGttt                                                            <  2:548693/76‑1 (MQ=255)
              cGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGttt                                                            <  2:1031142/76‑1 (MQ=255)
              cGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGGTAAATGGgtt                                                            >  1:813569/1‑76 (MQ=255)
              cGATGTTATGAGCAACACCACCCGTGCAGAACGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGGTAAGTGGgtt                                                            >  1:1392960/1‑76 (MQ=255)
                aTGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGATAAATGGgtttt                                                          >  1:273065/1‑76 (MQ=255)
                   ttATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGATAAATGGTTTTtttc                                                       >  1:67714/1‑76 (MQ=255)
                    tATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCg                                                      <  1:590619/76‑1 (MQ=255)
                    tATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGATAAAAGGGTTTTTTCg                                                      >  1:473479/1‑76 (MQ=255)
                            acacCCCCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCa                                              <  1:1155805/76‑1 (MQ=255)
                            acacCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGCTAAAAGGGTTNTTTTCGTTTCCgt                                               >  1:197659/1‑75 (MQ=255)
                            acacCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCa                                              <  1:1155813/76‑1 (MQ=255)
                            acacCACCCGTGCAGAGCGNCTGGAAGATGGCTCTTATCGGCTGGGGGGGGATAAATGGGTTTTTTCGGTTCCGCa                                              >  1:1433362/1‑76 (MQ=255)
                            acacCACACGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGAAGAAATGGCTTTTCTCGGTTCCGCa                                              <  1:83559/76‑1 (MQ=255)
                              accaccCGTGCAGAGCGTGTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCaaa                                            <  1:1453560/76‑1 (MQ=255)
                               ccaccCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGGTAAATGGGTTTTTTTGGGTCCGCCAAg                                           >  1:939730/1‑76 (MQ=255)
                                   ccGTGCAGAGCGTCTGGAAGATGGCTCTTATNGGCTGGTGTGGCATNAATGGNTTTTCTCGGTTCCGCAAAGCGAt                                       <  2:205209/76‑1 (MQ=255)
                                   ccGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGAt                                       <  1:455488/76‑1 (MQ=255)
                                     gTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGATAAAAGGGTTTTCTCGGTTCCGCCAAGCGAAgc                                     >  1:551141/1‑76 (MQ=255)
                                        cAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCa                                  <  1:242805/76‑1 (MQ=255)
                                         agagCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCAt                                 <  2:1154338/76‑1 (MQ=255)
                                         agagCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCAt                                 <  2:1178767/76‑1 (MQ=255)
                                         agagCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCAt                                 <  2:884419/76‑1 (MQ=255)
                                          gagCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATc                                <  1:309298/76‑1 (MQ=255)
                                           agCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCt                               <  1:160673/76‑1 (MQ=255)
                                            gCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTg                              <  2:796925/76‑1 (MQ=255)
                                             cGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTgg                             <  1:59350/76‑1 (MQ=255)
                                             cGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTgg                             <  2:695604/76‑1 (MQ=255)
                                             cGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGATAAATGGGTTTTTTCGGTTCCGCAAAGCGATGCGCATCTgg                             >  1:814421/1‑76 (MQ=255)
                                               tCTGGAAGATGGCTCTTATCGGCTGGGGGGGGATAAATGGTTTTTCTCGGTTCCGCAAAACGATGCGCATCTgggg                           >  2:968284/1‑74 (MQ=255)
                                                cTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGGGc                          >  2:1077755/1‑76 (MQ=255)
                                                cTGGAAGATGGCTCTTATCGGCTGGGGGGGCATNAATGGTTTTTCTCGGTTCCGCANAGCGTTGCCCNTCTGGTGc                          >  2:1288736/1‑76 (MQ=255)
                                                 tGGAAGATGGCTCTTATCGGCTGGTGGGGGATAAATGGGTTTTTTCGGTTCCGCAAAGCGATGCGCATCTGGGGcg                         >  2:407382/1‑75 (MQ=255)
                                                  ggAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTg                        <  1:123127/76‑1 (MQ=255)
                                                  ggAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTg                        <  1:1205052/76‑1 (MQ=255)
                                                  ggAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTg                        <  2:482081/76‑1 (MQ=255)
                                                  ggAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTATCGGTTCCGCAAAGCGATGCGCATCTGGTGCTg                        <  2:1318017/76‑1 (MQ=255)
                                                   gAAGATGGCTCTTATCGGCTGGGGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGGGCTgg                       >  2:768988/1‑76 (MQ=255)
                                                      gATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGgcgc                    <  2:188834/76‑1 (MQ=255)
                                                       aTGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCCAAGGGATGCCAATCTGGTGCTGGCGCa                   <  1:80712/76‑1 (MQ=255)
                                                       aTGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATNCGCATCTGGTGCTGGCGCa                   <  1:1426693/76‑1 (MQ=255)
                                                       aTGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCa                   <  1:269559/76‑1 (MQ=255)
                                                       aTGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCa                   <  1:911021/76‑1 (MQ=255)
                                                        tGGCTCTTATCGGCTGGTGGGGCATAAATGGGTTTTTTCGGGTCCCCAAAGCGATGCGCCTCTGGTGCTGGCGCAg                  >  2:812544/1‑76 (MQ=255)
                                                        tGGCTCTTATCGGCTGGGGGGGGATAAATGGGTTTTTTCGGGTCCGCAAAGCGATGCGCATCTGGTGCTGgcgccg                  >  2:322445/1‑74 (MQ=255)
                                                          gCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGTAGAc                <  2:1510277/76‑1 (MQ=255)
                                                          gCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGAc                <  1:1189710/76‑1 (MQ=255)
                                                           ctntTATCGGCTGGTGGGGCATAAACGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATTTGGTGCTGGCGCTGAcc               <  2:1440196/76‑1 (MQ=255)
                                                           ctctTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGAcc               <  2:1081437/76‑1 (MQ=255)
                                                             ctTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACcgc             <  1:646499/76‑1 (MQ=255)
                                                             ctTATCGGCTGGGGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGGGCTGGCGCAGACcgc             >  2:1092968/1‑76 (MQ=255)
                                                               tATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCgg           <  1:1016388/76‑1 (MQ=255)
                                                                 tCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGggt         <  1:1016896/76‑1 (MQ=255)
                                                                  cGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGggtg        <  2:698884/76‑1 (MQ=255)
                                                                   ggCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGggtgg       <  2:608124/76‑1 (MQ=255)
                                                                   ggCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGggtgg       >  1:281457/1‑76 (MQ=255)
                                                                     cTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTc     <  1:786236/76‑1 (MQ=255)
                                                                        ggtggGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCTg   <  2:1218886/76‑1 (MQ=255)
                                                                         gtggGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCTGt  >  1:1037453/1‑76 (MQ=255)
                                                                          |                                                                          
AAGCAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCTGT  >  REL606/4392988‑4393135

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: