Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6061,747,9030GT31.6% 145.9 / 67.2 76A8S (GCA→TCA) ydiLhypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (25/27);  new base T (24/0);  total (49/27)
Fisher's exact test for biased strand distribution p-value = 1.80e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.09e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AGGCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGAAA  >  REL606/1747829‑1747975
                                                                          |                                                                        
agGCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGc                                                                         >  2:704079/1‑76 (MQ=255)
 ggCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCa                                                                        >  1:1052724/1‑76 (MQ=255)
  gCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAc                                                                       >  1:305700/1‑76 (MQ=255)
   cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcg                                                                      >  1:1657916/1‑75 (MQ=255)
   cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcc                                                                      >  1:113502/1‑75 (MQ=255)
   cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcc                                                                      >  1:1238788/1‑75 (MQ=255)
   cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcc                                                                      >  1:578410/1‑75 (MQ=255)
   cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcc                                                                      >  1:1772686/1‑75 (MQ=255)
   cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAca                                                                      >  1:585312/1‑75 (MQ=255)
   cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACt                                                                      >  1:98110/1‑76 (MQ=255)
   cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACt                                                                      >  1:1204946/1‑76 (MQ=255)
   cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACt                                                                      >  2:428248/1‑76 (MQ=255)
    tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcca                                                                     >  1:1289643/1‑74 (MQ=255)
    tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcca                                                                     >  1:273253/1‑74 (MQ=255)
    tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACtt                                                                     >  1:886533/1‑75 (MQ=255)
    tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACtt                                                                     >  1:1491679/1‑75 (MQ=255)
    tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACta                                                                     >  1:268549/1‑75 (MQ=255)
    tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTc                                                                     >  2:646987/1‑76 (MQ=255)
     gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcctt                                                                    >  1:1761889/1‑73 (MQ=255)
     gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAccga                                                                    >  1:456781/1‑73 (MQ=255)
     gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAccag                                                                    >  1:180439/1‑73 (MQ=255)
     gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACtta                                                                    >  1:1023341/1‑74 (MQ=255)
     gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACtta                                                                    >  1:1620024/1‑74 (MQ=255)
     gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACtta                                                                    >  1:134211/1‑74 (MQ=255)
     gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACtta                                                                    >  1:1164874/1‑74 (MQ=255)
     gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACGcc                                                                    >  1:148115/1‑76 (MQ=255)
     gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTcc                                                                    <  1:1508334/76‑1 (MQ=255)
      aTAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCg                                                                   >  2:1476789/1‑76 (MQ=255)
      aTAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCg                                                                   <  1:362052/76‑1 (MQ=255)
       tAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGc                                                                  >  2:205679/1‑76 (MQ=255)
       tAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGc                                                                  >  2:1780982/1‑76 (MQ=255)
        aaTGTTATTTATATTTTATTGTTTTTTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACGccaat                                                                 >  1:706834/1‑73 (MQ=255)
        aaTGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGcc                                                                 <  1:1059681/76‑1 (MQ=255)
        aaTGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGcc                                                                 >  1:551286/1‑76 (MQ=255)
         aTGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCa                                                                >  2:1396870/1‑76 (MQ=255)
         aTGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCa                                                                >  2:1382219/1‑76 (MQ=255)
         aTGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCa                                                                <  1:718625/76‑1 (MQ=255)
          tGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCat                                                               >  1:747443/1‑76 (MQ=255)
          tGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCat                                                               >  1:747416/1‑76 (MQ=255)
          tGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCat                                                               >  1:399418/1‑76 (MQ=255)
            ttatttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCatat                                                             >  1:433766/1‑76 (MQ=255)
             tatttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCatatc                                                            <  1:1612474/76‑2 (MQ=255)
             tatttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATAtt                                                            <  1:507907/76‑1 (MQ=255)
             tatttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATAtt                                                            <  2:412230/76‑1 (MQ=255)
              atttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATAttt                                                           <  2:1669681/76‑1 (MQ=255)
               tttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATAtttt                                                          <  1:97708/76‑1 (MQ=255)
               tttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATAtttt                                                          >  2:617890/1‑76 (MQ=255)
                   tattttatttTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGcc                                                      <  1:908374/76‑1 (MQ=255)
                       ttATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGa                                                  >  1:1790256/1‑76 (MQ=255)
                        tATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGAc                                                 <  1:1757325/76‑1 (MQ=255)
                        tATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGAc                                                 >  2:1197064/1‑76 (MQ=255)
                        tATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGAc                                                 >  1:943130/1‑76 (MQ=255)
                           tttttATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTAt                                              <  1:676733/74‑1 (MQ=255)
                            gTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGCCTAtt                                             >  2:1510380/1‑76 (MQ=255)
                             tttATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtg                                            <  1:526044/76‑1 (MQ=255)
                               tttttCATGGAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtgat                                          <  1:730102/74‑1 (MQ=255)
                                aTTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtgatg                                         <  2:1118266/76‑1 (MQ=255)
                                  ttCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGaa                                       >  1:1036214/1‑76 (MQ=255)
                                  ttCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGaa                                       <  1:1555484/76‑1 (MQ=255)
                                  ttCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACGCCGCCATATTTTTGCCATGACTATTGATGaa                                       <  2:1061194/76‑1 (MQ=255)
                                    cATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATg                                     <  1:1784553/76‑1 (MQ=255)
                                      tGCAGCCATTACACGATGAACGCTTATGACCCCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATgcg                                   <  1:905210/76‑1 (MQ=255)
                                      tGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATgcg                                   <  1:290206/76‑1 (MQ=255)
                                       gCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATgcgc                                  >  1:1422235/1‑76 (MQ=255)
                                              ttACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTgg                           >  2:411744/1‑76 (MQ=255)
                                                acacGATGAACGCTTATGACCTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGAt                         <  1:1613087/76‑1 (MQ=255)
                                                acacAATAAACGCTAATCACCCCCAGGCCCTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGAt                         <  1:1391528/76‑1 (MQ=255)
                                                  acGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTg                       <  2:967774/76‑1 (MQ=255)
                                                  acGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTg                       <  1:1732172/76‑1 (MQ=255)
                                                     aTGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGccc                    <  1:454889/76‑1 (MQ=255)
                                                      tGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCa                   <  1:101985/76‑1 (MQ=255)
                                                       gAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAg                  <  2:1228441/76‑1 (MQ=255)
                                                       gAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAg                  <  2:1037064/76‑1 (MQ=255)
                                                             ttATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGt            <  2:67510/76‑1 (MQ=255)
                                                              tATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTg           >  1:1633385/1‑76 (MQ=255)
                                                                 gAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACa        <  2:166131/76‑1 (MQ=255)
                                                                 gAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACa        >  1:1713933/1‑76 (MQ=255)
                                                                   cccccAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGt      <  1:99592/73‑1 (MQ=255)
                                                                   accccAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGt      <  1:1499840/73‑1 (MQ=255)
                                                                   aCTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGt      >  2:774456/1‑76 (MQ=255)
                                                                    cTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTg     <  1:226958/76‑1 (MQ=255)
                                                                      ccAGGCCCTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGaa   <  1:1662714/76‑1 (MQ=255)
                                                                      ccAGCCCCTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGaa   <  1:1578051/76‑1 (MQ=255)
                                                                       cAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGaaa  >  2:719783/1‑76 (MQ=255)
                                                                          |                                                                        
AGGCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGAAA  >  REL606/1747829‑1747975

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