Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6064,195,8740AC28.4% 153.8 / 12.5 67T318P (ACC→CCC) aceBmalate synthase
Reads supporting (aligned to +/- strand):  ref base A (35/13);  new base C (0/19);  total (35/32)
Fisher's exact test for biased strand distribution p-value = 1.44e-08
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:1167002/76‑1 (MQ=255)
  cgTCCTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAAcc                                                                          >  1:411449/1‑76 (MQ=255)
  cgTCCTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAAcc                                                                          >  1:1244035/1‑76 (MQ=255)
   gTCCTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCt                                                                         <  1:1345132/76‑1 (MQ=255)
     ccTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGc                                                                       >  1:505662/1‑76 (MQ=255)
      cTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGcc                                                                      >  2:1568895/1‑76 (MQ=255)
      cTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGcc                                                                      >  1:1113311/1‑76 (MQ=255)
      cTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGcc                                                                      <  1:1687738/76‑1 (MQ=255)
        gCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCCCGCCTGTTGATTAAAACCTGCCAt                                                                    >  2:834205/1‑76 (MQ=255)
        gCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCAt                                                                    >  1:4581/1‑76 (MQ=255)
        gCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCAt                                                                    >  2:286872/1‑76 (MQ=255)
         ccAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTTATTAAAACCTGCCATa                                                                   <  2:1378556/76‑1 (MQ=255)
         ccAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATa                                                                   <  1:459872/76‑1 (MQ=255)
         ccAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATa                                                                   >  1:817903/1‑76 (MQ=255)
          cagacagaCAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATaa                                                                  >  1:1162012/1‑76 (MQ=255)
          cagacagaCAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATaa                                                                  >  1:847564/1‑76 (MQ=255)
             acagacagGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAAcg                                                               <  2:1071380/76‑1 (MQ=255)
             acagacagGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAAcg                                                               <  1:1683752/76‑1 (MQ=255)
              cagacagGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAAcgc                                                              <  1:1345071/76‑1 (MQ=255)
               agacagGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAAcgcg                                                             >  1:146102/1‑76 (MQ=255)
                gacagGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCCTCAACGCgg                                                            >  1:1533882/1‑76 (MQ=255)
                gacagGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCgg                                                            >  1:705405/1‑76 (MQ=255)
                 acagGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGt                                                           >  2:1681217/1‑76 (MQ=255)
                       aGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAACCCTGCCATAAACGCGGTGCtttt                                                     <  1:121602/76‑1 (MQ=255)
                        gTGACGATGGAAAACCCATCCCTGAATGCTTCCTCCCGCCTGTTGATAAAAACCTGCCATAAACGCGGTGCTTTTg                                                    <  1:1047666/76‑1 (MQ=255)
                           aCGATGGATAAACCATTCCTGAATGCTTACTCCCGCCTGTTGATTAAACCCTGCCATAAACGCGGTGCTTTTGCGa                                                 <  1:892576/76‑1 (MQ=255)
                           aCGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGa                                                 >  1:351002/1‑76 (MQ=255)
                            cGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGAt                                                >  1:972058/1‑76 (MQ=255)
                              aTGGATACCCCATCCCTTAATTCTTACTCACGCCTGTTAATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATgg                                              <  2:668617/76‑1 (MQ=255)
                                ggATAATCCATTCTGGATTGTTTACTCCCGTTTGTGGTTTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGggc                                            <  2:1190439/76‑1 (MQ=255)
                                ggATAAACCATTCCTAAATGCTTACTCCCGCCTGTTGATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATGggc                                            <  1:1244221/76‑1 (MQ=255)
                                 gATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATGggcn                                           <  1:1474850/76‑2 (MQ=255)
                                  aTAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGggcgg                                          >  2:187534/1‑76 (MQ=255)
                                  aTAAACCATTCCTGAATGCTTACTCACGCCTGTTGATAAAAACCTGCCATAAACGCGGTGCTTTTGCGATGggcgg                                          <  2:650662/76‑1 (MQ=255)
                                    aatccATTCCTGATTGCTTACCCCCCCCTTTTGATAAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCa                                        <  1:846019/73‑1 (MQ=255)
                                    aaaCCATTCCTGAATGCTTACTCACGCCTGTTGATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCa                                        <  2:1642823/76‑1 (MQ=255)
                                      cccATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATg                                      <  1:1569453/75‑1 (MQ=255)
                                      aCCATTCCTGAATGCTTCCCCCCGCTTGTTGATAAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATg                                      <  2:790577/76‑1 (MQ=255)
                                      aCCATTCCTGAATGCTTACCCACGCCTGTTGATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATg                                      <  2:1480125/76‑1 (MQ=255)
                                         aTTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATggcg                                   >  2:1130179/1‑76 (MQ=255)
                                            ccTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGgcggcg                                >  1:1118502/1‑76 (MQ=255)
                                             cTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGt                               >  1:642796/1‑76 (MQ=255)
                                              tGAATGCTTACTCACGCCTGTTGATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGtt                              <  1:1012805/76‑1 (MQ=255)
                                                aaTGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTa                            >  1:551756/1‑76 (MQ=255)
                                                aaTGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTa                            >  2:245153/1‑76 (MQ=255)
                                                 aTGCTTCCTCCCGCCTGTTGATAAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTAt                           <  1:979984/76‑1 (MQ=255)
                                                 aTGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTAt                           >  1:80396/1‑76 (MQ=255)
                                                   gCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTc                         >  1:1612259/1‑76 (MQ=255)
                                                    cTTCCTCCCGCCGGTTAATAAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTcc                        <  1:1148108/76‑1 (MQ=255)
                                                    cTTACTCACGCCTTTTGATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTcc                        <  1:1561521/76‑1 (MQ=255)
                                                    cTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTcc                        >  1:1111890/1‑76 (MQ=255)
                                                     ttACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCg                       >  2:681988/1‑76 (MQ=255)
                                                     ttACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCg                       >  1:1434645/1‑76 (MQ=255)
                                                     ttACCCACGCCTGTTAATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCg                       <  1:593019/76‑1 (MQ=255)
                                                      tACTCACGCCTGTTGATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGa                      <  1:2047/76‑1 (MQ=255)
                                                      tACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGa                      >  1:1133340/1‑76 (MQ=255)
                                                       ccTCACGCCTGTTGATAAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAg                     <  1:677842/75‑1 (MQ=255)
                                                        cTCCCGCCTGTTGATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGc                    <  2:1627181/76‑1 (MQ=255)
                                                        cTCACGCCTGTTGATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGc                    <  1:153530/76‑1 (MQ=255)
                                                         ccACGCCTGTTGATAAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCa                   <  1:289269/75‑1 (MQ=255)
                                                          cccGCCTTTTGATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCaa                  <  1:1221044/74‑1 (MQ=255)
                                                          cACGCCTGTTGATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCaa                  <  2:1073452/76‑1 (MQ=255)
                                                           ccGCCTTTTGATAAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCaaa                 <  1:334546/75‑1 (MQ=255)
                                                            cGCCTGTTGATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAg                <  2:1423959/76‑1 (MQ=255)
                                                              ccTGTTAATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGAt              <  1:1038300/76‑1 (MQ=255)
                                                                tGTTGATTAAGCCCTGCCATAAACTCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGa            <  1:120475/76‑1 (MQ=255)
                                                                tGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGa            >  1:1247836/1‑76 (MQ=255)
                                                                 gTTGATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGaa           <  1:1659688/76‑1 (MQ=255)
                                                                   tGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGAaga         >  2:206174/1‑76 (MQ=255)
                                                                   tGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGAaga         >  1:698791/1‑76 (MQ=255)
                                                                    gATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGAagag        <  2:218579/76‑1 (MQ=255)
                                                                    gATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGAagag        >  2:27627/1‑76 (MQ=255)
                                                                    gATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGAagag        >  1:1219953/1‑76 (MQ=255)
                                                                     ataaaacCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGAAGAGc       <  1:380754/73‑1 (MQ=255)
                                                                       tAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGAAGAGcac     >  1:1151306/1‑76 (MQ=255)
                                                                          aaCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGAAGAGCACAAt  >  1:1624192/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 < 29 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: