Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6064,195,8740AC29.4% 109.2 / 11.8 52T318P (ACC→CCC) aceBmalate synthase
Reads supporting (aligned to +/- strand):  ref base A (26/10);  new base C (1/14);  total (27/25)
Fisher's exact test for biased strand distribution p-value = 2.25e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGCGTCCTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGAAGAGCACAAT  >  REL606/4195799‑4195948
                                                                           |                                                                          
cgcgTCCTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTaaaa                                                                            <  1:630687/76‑1 (MQ=255)
  ggTCCTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAAcc                                                                          <  2:1202856/75‑1 (MQ=255)
  cgTCCTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAAcc                                                                          <  2:1202866/76‑1 (MQ=255)
   gTCCTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCt                                                                         >  2:287359/1‑76 (MQ=255)
   gTCCTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCt                                                                         <  1:556890/76‑1 (MQ=255)
   gTCCTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCt                                                                         <  1:872312/76‑1 (MQ=255)
    tCCTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTg                                                                        <  2:652482/76‑1 (MQ=255)
     ccTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGc                                                                       >  2:920841/1‑76 (MQ=255)
       tGCCAGATAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCa                                                                     <  2:420181/76‑1 (MQ=255)
       tGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCa                                                                     >  1:120258/1‑76 (MQ=255)
        gccagccaGACAGGCAGTGACGATGGATAACCCATTCCTGAATGCTTTCCCACGCCTGTTAATTAAAGCCTGCCAt                                                                    <  2:243526/76‑1 (MQ=255)
        gCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCAt                                                                    <  2:579922/76‑1 (MQ=255)
        gCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCAt                                                                    >  1:1154118/1‑76 (MQ=255)
          cagacagaCAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATaa                                                                  >  2:363104/1‑76 (MQ=255)
          cagacagaCAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATaa                                                                  <  1:153291/76‑1 (MQ=255)
           agacagacAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCGNTTGATTAAAACCTACCATaac                                                                 >  2:270463/1‑75 (MQ=255)
            gacagacaGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAACCCTGCCATAAAc                                                                <  1:989408/76‑1 (MQ=255)
                  cagGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAACCCTGCCATAAACGCGGTg                                                          <  2:597175/76‑1 (MQ=255)
                  cagGCAGTGACGATGGAAAAACCATCCCTGAATGCTTACTCACGCCTGTTAATTAAAACCTGCCATAAACGCGGTg                                                          <  2:810559/76‑1 (MQ=255)
                   agGCAGTGACGATGGATAAACCATACCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGc                                                         >  2:1059390/1‑76 (MQ=255)
                     gCAGTGACGATGGATAAACCATTCCTGAATGCTTANNNACGCCTGTTGATTAAAANCTGCCATAAACGCGGTGCtt                                                       >  1:843379/1‑76 (MQ=255)
                     gCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTTTTGATTAAAACCTGCCATAAACGCGGTGCtt                                                       <  1:1559927/76‑1 (MQ=255)
                     gCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCtt                                                       >  1:1164655/1‑76 (MQ=255)
                          gACGATGGATAAACCATTCCTGAATGCTTACCCACGCCTGTTGATTAAACCCTGCCATAAACGCGGTGCTTTTGCg                                                  <  1:617145/76‑1 (MQ=255)
                            cGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGAt                                                >  1:1200911/1‑76 (MQ=255)
                              aTGGATAAACCATTCCTGAATGCTTCCTCCCGCCTGTTGATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATgg                                              <  2:1102012/76‑1 (MQ=255)
                              aTGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAACACTGCCATACACGCGGTGCTTTTGCGACgg                                              >  2:394907/1‑76 (MQ=255)
                              aTGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATgg                                              >  2:715764/1‑76 (MQ=255)
                               tGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATggg                                             >  1:1039567/1‑76 (MQ=255)
                               tGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATggg                                             >  1:147436/1‑76 (MQ=255)
                               tGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATggg                                             >  1:1651774/1‑76 (MQ=255)
                                ggATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGggc                                            >  2:1479459/1‑76 (MQ=255)
                                ggATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGggc                                            >  1:267376/1‑76 (MQ=255)
                                  aTAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGggcgg                                          >  1:1207377/1‑76 (MQ=255)
                                  aTAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGggcgg                                          >  1:318126/1‑76 (MQ=255)
                                    aaaCCATCCCTGAATGCTTACCCCCGCCTGTTGATAAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCa                                        <  1:586976/76‑1 (MQ=255)
                                     aaCCATTCCTAAATGCTTACTCACGCCTTTTGATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCAt                                       <  1:264905/76‑1 (MQ=255)
                                        cATCCCTGAATGCTTACTCCCGCCTGTTGATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATggc                                    <  2:957404/76‑1 (MQ=255)
                                            ccTGAATGCTTACTCCCGCCTGTTGATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGgcggcg                                <  1:154307/76‑1 (MQ=255)
                                             cTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGt                               >  1:1128885/1‑76 (MQ=255)
                                              tGAATGCTTCCCCACGCCTGTTAATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGtt                              <  2:82264/76‑1 (MQ=255)
                                               gAATGCTTCCTCCCGCCTGTTGATAAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGttt                             <  1:1507467/76‑1 (MQ=255)
                                               gAATGCTTACTCCCGCCTGTTGATAAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGttt                             <  1:515969/76‑1 (MQ=255)
                                               gAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGttt                             >  1:940771/1‑76 (MQ=255)
                                                 aTGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTAt                           >  1:255470/1‑76 (MQ=255)
                                                  tGCTTCCCCCCGCCTTTTGATAAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTAtt                          <  2:1267867/76‑1 (MQ=255)
                                                     ttACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCg                       >  2:594484/1‑76 (MQ=255)
                                                         tCACGCATTTTAATAAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCa                   <  1:1369850/76‑1 (MQ=255)
                                                         ccACGCCTTTTGATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCa                   <  1:699750/75‑1 (MQ=255)
                                                          cgcGCCTGTTGATAAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCaa                  <  1:1500226/74‑1 (MQ=255)
                                                            cGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAg                >  1:1182502/1‑76 (MQ=255)
                                                              cttGTTGATAAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGAt              <  2:1049775/74‑1 (MQ=255)
                                                                tGTTGATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATgn            <  1:463590/76‑2 (MQ=255)
                                                                tGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGa            >  1:151703/1‑76 (MQ=255)
                                                                  ttGATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGAag          <  1:1386653/76‑1 (MQ=255)
                                                                  ttGATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGAag          <  1:1096885/76‑1 (MQ=255)
                                                                  ttGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGAag          >  1:10103/1‑76 (MQ=255)
                                                                  ttGATAAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGTGGCATGGCGGCGTTTATTCCGAGCAAAGATGAag          <  1:928609/76‑1 (MQ=255)
                                                                    gATAAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGAagag        <  1:842989/76‑1 (MQ=255)
                                                                       tAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGAAGAGcac     >  1:1464723/1‑76 (MQ=255)
                                                                       aaaacCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGAAGAGcac     <  1:715764/75‑1 (MQ=255)
                                                                        aaaaCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGAAGAGcaca    >  2:1552211/1‑76 (MQ=255)
                                                                        aaaaCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGAAGAGcaca    <  1:1520167/76‑1 (MQ=255)
                                                                         aaaCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGAAGAGCACaa   >  2:6277/1‑76 (MQ=255)
                                                                          aaCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGAAGAGCACAAt  >  1:1301386/1‑76 (MQ=255)
                                                                           |                                                                          
CGCGTCCTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGAAGAGCACAAT  >  REL606/4195799‑4195948

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