Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6061,747,9030GT46.5% 322.4 / inf 389A8S (GCA→TCA) ydiLhypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (165/43);  new base T (181/0);  total (346/43)
Fisher's exact test for biased strand distribution p-value = 2.29e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.40e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 425 total aligned reads displayed.

CAGGCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGAAAGC  >  REL606/1747828‑1747977
                                                                           |                                                                          
  ggCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCa                                                                          >  1:5767604/1‑76 (MQ=255)
  ggCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCa                                                                          >  1:5560203/1‑76 (MQ=255)
  ggCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCa                                                                          >  1:4948783/1‑76 (MQ=255)
  ggCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCa                                                                          >  1:1672187/1‑76 (MQ=255)
  ggCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCa                                                                          >  1:4018645/1‑76 (MQ=255)
   gCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAc                                                                         >  1:1070115/1‑76 (MQ=255)
   gCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAc                                                                         >  1:7904691/1‑76 (MQ=255)
   gCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAc                                                                         >  1:4740590/1‑76 (MQ=255)
   gCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAc                                                                         >  1:4312044/1‑76 (MQ=255)
   gCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAc                                                                         >  1:3802725/1‑76 (MQ=255)
    cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcc                                                                        >  1:3343013/1‑75 (MQ=255)
    cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcc                                                                        >  1:3540480/1‑75 (MQ=255)
    cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAca                                                                        >  1:4989249/1‑75 (MQ=255)
    cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAca                                                                        >  1:4131393/1‑75 (MQ=255)
    cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACt                                                                        >  1:1281935/1‑76 (MQ=255)
    cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACt                                                                        >  1:7454380/1‑76 (MQ=255)
    cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACt                                                                        >  1:1136720/1‑76 (MQ=255)
    cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACt                                                                        >  1:6269786/1‑76 (MQ=255)
    cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACt                                                                        >  1:691212/1‑76 (MQ=255)
    cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACt                                                                        >  1:1789050/1‑76 (MQ=255)
    cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACt                                                                        >  1:371851/1‑76 (MQ=255)
    cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACt                                                                        >  2:1796369/1‑76 (MQ=255)
    cTGATAATGTTATTTATATTTTATTGTTAATTCATTCAGGCATTACACGGTGGAAGCTTTTGGACTCCCGGCACtg                                                                       >  1:3228438/1‑75 (MQ=255)
     tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcct                                                                       >  1:4698864/1‑74 (MQ=255)
     tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcct                                                                       >  1:2988073/1‑74 (MQ=255)
     tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcca                                                                       >  1:1873794/1‑74 (MQ=255)
     tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcat                                                                       >  1:5649773/1‑74 (MQ=255)
     tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcat                                                                       >  1:2202385/1‑74 (MQ=255)
     tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACtt                                                                       >  1:6830228/1‑75 (MQ=255)
     tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACtt                                                                       >  1:6977220/1‑75 (MQ=255)
     tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACtt                                                                       >  1:7065859/1‑75 (MQ=255)
     tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACta                                                                       >  1:6079721/1‑75 (MQ=255)
     tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTc                                                                       <  1:4587696/76‑1 (MQ=255)
      gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACtta                                                                      >  1:6681913/1‑74 (MQ=255)
      gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACtta                                                                      >  1:8133382/1‑74 (MQ=255)
      gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACtta                                                                      >  1:5838494/1‑74 (MQ=255)
      gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACtta                                                                      >  1:7345446/1‑74 (MQ=255)
      gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACtta                                                                      >  1:118713/1‑74 (MQ=255)
      gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACtta                                                                      >  1:5521098/1‑74 (MQ=255)
      gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACtta                                                                      >  1:7799325/1‑74 (MQ=255)
      gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACtta                                                                      >  1:19366/1‑74 (MQ=255)
      gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACtta                                                                      >  1:2138261/1‑74 (MQ=255)
       aTAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACttag                                                                     >  1:2043838/1‑73 (MQ=255)
       aTAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACttag                                                                     >  1:4396716/1‑73 (MQ=255)
       aTAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACttag                                                                     >  1:4884245/1‑73 (MQ=255)
       aTAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACttag                                                                     >  1:1745359/1‑73 (MQ=255)
       aTAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCg                                                                     <  2:2003394/76‑1 (MQ=255)
        tAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGc                                                                    >  2:4409659/1‑76 (MQ=255)
        tAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGc                                                                    >  2:7297963/1‑76 (MQ=255)
         aaTGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACGccaat                                                                   >  1:5288386/1‑73 (MQ=255)
           tGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCat                                                                 >  1:5188453/1‑76 (MQ=255)
           tGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCat                                                                 >  1:1537194/1‑76 (MQ=255)
             ttatttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACTTAGGCatat                                                               >  1:8266350/1‑76 (MQ=255)
             ttatttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCatat                                                               >  1:2265542/1‑76 (MQ=255)
              tatttatATTTTTTTGTTTTTTTAATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATAtt                                                              <  1:7645129/76‑1 (MQ=255)
              tatttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACTTAGGCatatc                                                              >  1:3931212/1‑75 (MQ=255)
              tatttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATAtt                                                              >  1:2771409/1‑76 (MQ=255)
               atttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATAttt                                                             >  2:3867599/1‑76 (MQ=255)
                 ttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATAttttt                                                           >  2:6483254/1‑76 (MQ=255)
                      ttttATTGTTTATTCCACGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCAt                                                      <  2:602839/76‑1 (MQ=255)
                       tttattttttatttCATCCAGCCATTACACGATGAACGCTTATGAACCCCAGGCCCTCCGCCATATTTTTGCCATg                                                     <  1:1015115/76‑1 (MQ=255)
                        ttATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGa                                                    >  1:2583106/1‑76 (MQ=255)
                           ttGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTa                                                 >  1:5380670/1‑76 (MQ=255)
                              tttATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtg                                              >  2:274718/1‑76 (MQ=255)
                              tttATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtg                                              >  1:3394285/1‑76 (MQ=255)
                              tttATTTCATCCAGCCATTACACGATGAACGCTTATGACCCCCAGGCACCCCGCCATATTTTTGCCATGACTATtg                                              <  1:7984549/76‑1 (MQ=255)
                              tttATTCCATCCAGCCATTACACGATAAACGTTTAGGACCCCCAGGCCCTCCGCCATATTTTTGCCATGACTATtg                                              <  1:7298685/76‑1 (MQ=255)
                               ttATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtga                                             >  1:2839284/1‑76 (MQ=255)
                                tATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtgat                                            >  1:1465186/1‑76 (MQ=255)
                                 aTTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtgatg                                           >  2:4681232/1‑76 (MQ=255)
                                 aTTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtgatg                                           >  1:5910961/1‑76 (MQ=255)
                                 aTTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtgatg                                           >  1:6005190/1‑76 (MQ=255)
                                   ttCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGaa                                         >  2:3344439/1‑76 (MQ=255)
                                    tCAGGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAAt                                        <  2:5436367/76‑1 (MQ=255)
                                     cATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATg                                       <  2:1268457/76‑1 (MQ=255)
                                     cATCCAGCCATTACACGATGAACGCTTATGAACCCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATg                                       <  1:4455482/76‑1 (MQ=255)
                                      aTGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATgc                                      >  2:1117237/1‑76 (MQ=255)
                                      aTGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATgc                                      >  1:2652638/1‑76 (MQ=255)
                                      aTGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATgc                                      >  1:6460598/1‑76 (MQ=255)
                                       tccAGCCATTACACGATGAACGCTTATGACCCCCAGGCCCTCCGCCATATTTTTGCCATGACTATTGATGAATgcg                                     <  1:6190544/74‑1 (MQ=255)
                                       tGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATgcg                                     >  1:6600371/1‑76 (MQ=255)
                                        gCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATgcgc                                    >  2:7565995/1‑76 (MQ=255)
                                        gCAGACATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATgcgc                                    >  1:4280835/1‑76 (MQ=255)
                                         cAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGcgcg                                   >  2:5875678/1‑76 (MQ=255)
                                              aTTACACGATGAACGCTTATGACCTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTg                              <  2:5547526/76‑1 (MQ=255)
                                              aTTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTg                              >  1:3682727/1‑76 (MQ=255)
                                                 acacGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGAt                           >  1:738949/1‑76 (MQ=255)
                                                  cacGATGAACTCTTATGAACCCCAGCCCCTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGAtt                          <  1:4170826/76‑1 (MQ=255)
                                                      ataaaCGTTTATGACCCCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGccc                      <  2:5189595/73‑1 (MQ=255)
                                                      aTGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGccc                      >  2:366175/1‑76 (MQ=255)
                                                       tGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCa                     <  2:6922718/76‑1 (MQ=255)
                                                         aaCGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGa                   >  2:4012744/1‑76 (MQ=255)
                                                         aaCGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGa                   >  1:7151918/1‑76 (MQ=255)
                                                           cGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGAcc                 <  2:2262790/76‑1 (MQ=255)
                                                            gtttATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCg                <  2:4280835/74‑1 (MQ=255)
                                                             cTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCgg               >  1:6505905/1‑76 (MQ=255)
                                                              ttATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGt              >  2:6665663/1‑76 (MQ=255)
                                                                       ccAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGaa     >  2:7801113/1‑76 (MQ=255)
                                                                        cAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTCACAGTGaaa    >  2:4935795/1‑76 (MQ=255)
                                                                         aGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGAAAg   <  2:161621/76‑1 (MQ=255)
                                                                           |                                                                          
CAGGCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGAAAGC  >  REL606/1747828‑1747977

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: