Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL606691,8520TG36.9% 272.2 / 26.6 160intergenic (+17/‑33)ybfM/ybfNpredicted outer membrane porin/predicted lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (55/46);  new base G (59/0);  total (114/46)
Fisher's exact test for biased strand distribution p-value = 4.59e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 238 total aligned reads displayed.

TCTTCCAGGATGAGCGTGACGTAAAATTTATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGTTTTGTGCCTGCCGCACGTTTGAGGAATTTGCTATGAAAAAACTGATTCTCATCGCCATGATGGCATCGGGGCTGG  >  REL606/691777‑691927
                                                                           |                                                                            
tCTTCCAGGATGAGCGTGACGTAAAATTTATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGt                                                                              >  1:7401313/1‑76 (MQ=255)
tCTTCCAGGATGAGCGTGACGTAAAATTTATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGt                                                                              >  1:415080/1‑76 (MQ=255)
 cTTCCAGGATGAGCGTGACGTAAAATTTATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGtt                                                                             <  1:1108183/76‑1 (MQ=255)
 cTTCCAGGATGAGCGTGACGTAAAATTTATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGtt                                                                             <  1:751981/76‑1 (MQ=255)
  ttCCAGGATGAGCGTGACGTAAAATTTATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGgtt                                                                            >  2:4962807/1‑76 (MQ=255)
   tCCAGGATGAGCGTGACGTAAAATTTATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGtttt                                                                           >  2:3685805/1‑76 (MQ=255)
    ccAGGATGAGCGTGACGTAAAATTTATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGTTTtg                                                                          <  1:2132854/76‑1 (MQ=255)
    ccAGGATGAGCGTGACGTAAAATTTATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGGTTtg                                                                          >  1:4932255/1‑76 (MQ=255)
     cAGGATGAGCGTGACGTAAAATTTATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGGTTtgg                                                                         >  1:4276870/1‑75 (MQ=255)
      aGGATGAGCGTGACGTAAAATTTATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGTTTtgtg                                                                        >  2:7717692/1‑76 (MQ=255)
       ggATGAGCGTGACGTAAAATTTATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGTTTTGtgc                                                                       <  1:3303656/76‑1 (MQ=255)
         aTGAGCGTGACGTAAAATTTATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGTTTtgtgtcc                                                                     >  1:2426217/1‑73 (MQ=255)
         aTGAGCGTGACGTAAAATTTATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGTTTTGtgcct                                                                     <  1:2211527/76‑1 (MQ=255)
         aTGAGCGTGACGTAAAATTTATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGTTTTGtgcct                                                                     >  1:1659094/1‑76 (MQ=255)
         aTGAGCGTGACGTAAAATTTATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGGTTTGtgcct                                                                     >  2:1216268/1‑76 (MQ=255)
          tGAGCGTGACGTAAAATTTATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGTTTTGtgcctg                                                                    <  2:2464246/76‑1 (MQ=255)
           gAGCGTGACGTAAAATTTATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGGTTTGtgcctgc                                                                  >  1:1441686/1‑76 (MQ=255)
             gggTGACATAAATTTTATGGTAATCGCACCATTCACCATCTTCTGCTGCCCGACGCGGCAGATTTTGTCCCTGCCg                                                                <  2:6147529/74‑1 (MQ=255)
             gCGTGACGTAAAATTTATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGTTTTGTGCCTGCCg                                                                <  2:4806412/76‑1 (MQ=255)
             gCGTGACGTAAAATTTATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGTTTTGTGCCTGCCg                                                                <  2:2079157/76‑1 (MQ=255)
               gggACGTAAAATTTATGGTAATCGCCCCATTCACCATCTTCTGATGCCCGACGCGGCAGGTTTTGTCCCTGCCGCa                                                              <  2:5900840/74‑1 (MQ=255)
               gTGACGTAAAATTTATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGGTTTGTGCCTGCCGcc                                                              >  1:2802548/1‑75 (MQ=255)
                tGACGTAAAATTTATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGTTTTGTGCCTGCCGCAc                                                             <  2:7041939/76‑1 (MQ=255)
                tGACGTAAAATTTATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGTTTTGTGCCTGCCGCAc                                                             >  2:1436962/1‑76 (MQ=255)
                tGACGTAAAATTTATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGTTTTGTGCCTGCCGCAc                                                             <  1:5638259/76‑1 (MQ=255)
                tGACGTAAAATTTATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGTTTTGTGCCTGCCGCAc                                                             <  1:5742669/76‑1 (MQ=255)
                 gACGTAAAATTTATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGTTTTGTGCCTGCCGCACg                                                            >  1:3737320/1‑76 (MQ=255)
                 gACGTAAAATTTATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGTTTTGTGCCTGCCGCACg                                                            >  1:1773503/1‑76 (MQ=255)
                    gTAAAATTTATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGGTTTGTGCCTTCCGCACGttt                                                         >  2:3338724/1‑76 (MQ=255)
                     tAAATTTTATGGTAATCGCCCCATTCCCCATCTTCTGACGCCGGCCGCGGCGGGTTTTGTCCCTGCCGCACGTTTg                                                        <  2:4182978/76‑1 (MQ=255)
                     tAAAATTTATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGTTTTGTGCCTGCCGCACGTTTg                                                        <  2:7334489/76‑1 (MQ=255)
                     aaaatTTTATGGTAATCCCCCCATCCACCATCTTCTGATCCCCGACGCGGCGGGTTTTGTCCCTGCCGCACGTTTg                                                        <  2:75678/75‑1 (MQ=255)
                      aaaaTTTATGTTAATCGCACCAATCTCCATCTTCTGATGCCCGACGCGAAGGGTTTGGTGCCTGCCGCACGTTTGa                                                       <  2:7541227/76‑1 (MQ=255)
                      aaaaTTTATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGGTTTGTGCCTGCCGCCCGTTTGa                                                       >  1:878321/1‑76 (MQ=255)
                      aaaaTTTATGGTAAGCGCCCCATTCACCATCTTCTGATGCCCGACGCGGCAGGTTTTGTGCCTGCCGCACGTTTGa                                                       <  2:3474329/76‑1 (MQ=255)
                       aaaTTTATGTTAATCGCCCCATTCACCATCTTCTGATGCCCGACGCGGCAGGTTTTGTGCCTGCCGCACGTTTGAg                                                      <  1:3456970/76‑1 (MQ=255)
                       aaaTTTATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGGTTTTTGGCTGCCGCACGTTTgtg                                                      >  1:4386359/1‑74 (MQ=255)
                          tttATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGTTTTGTGCCTGCCGCACGTTTGAGGaa                                                   <  2:1123273/76‑1 (MQ=255)
                           ttATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGTTTTGTGCCTGCCGCCCGGTTTAGGAAt                                                  >  2:5339665/1‑76 (MQ=255)
                           ttATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGGTTTGTGCCTGCCGCACGTTTTAGGAAt                                                  >  2:7429657/1‑76 (MQ=255)
                            tATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGGTTTTTGCCTGCCCCCCGTTTTAGGAAtt                                                 >  2:1470454/1‑76 (MQ=255)
                            tATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGGTTTGTGCCTGCCGCACGTTTGAGGAAtt                                                 >  2:6370891/1‑76 (MQ=255)
                             aTGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGTTTTGTGCCTGCCGCACGTTTGAGGAAttt                                                <  2:3498305/76‑1 (MQ=255)
                             aTGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGGTTTGTGCCTGCCGCACGGTTTAGGAAttt                                                >  2:2526343/1‑76 (MQ=255)
                              tGGTAACCGCACCTTCCCCCATCTTCTGATCCCCGACGCGGCAGGTTTTGTCCCTGCCGCACGTTTGAGGAATTTg                                               <  2:4080838/76‑1 (MQ=255)
                               ggTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGTTTTGTGCCTGCCGCACGTTTGAGGAATTTGc                                              <  2:7072880/76‑1 (MQ=255)
                               ggTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGGTTTGCGTCTGGCGCCCGGTTGACGAATTTGc                                              >  1:6415452/1‑76 (MQ=255)
                                 tAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGTTTTGTGCCTGCCCCACGTTTGAGGAATTTGTTa                                            >  1:3632154/1‑76 (MQ=255)
                                 tAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGGTTTTTGCCTGCCGCCCGTTTGAGGGATTTGCTa                                            >  1:7235102/1‑76 (MQ=255)
                                 tAACCGCCCCATTCCCCATCTTCTGATGCCCGACGGGGCAGGTTTTGTGCCTGCCGCACGTTTGAGGAATTTGCTa                                            <  1:1522836/76‑1 (MQ=255)
                                  aaTCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGTTTTGTGCCTGCCGCACGTTTGAGGGATTTGCTAt                                           >  2:6838707/1‑76 (MQ=255)
                                  aaTCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGGTTTGTGCCTGCCGCACGTTTGAGGAATTTGCTAt                                           >  1:3240131/1‑76 (MQ=255)
                                  aaTCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGGTTTGTGCCTGCCGCACGCTTGAGGAATTTTCTAt                                           >  2:5513985/1‑76 (MQ=255)
                                  aaTCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGGTTTGAGCCTGCCGCACGTTTTAGCAATTTGCtct                                           >  2:1766199/1‑74 (MQ=255)
                                  aaTCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGATTTTTTTCTTGCGGACCTTTGAGGAAATTGCTAt                                           >  2:7676297/1‑76 (MQ=255)
                                     cGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGTTTTGTTCCCGCCCCGCGTTTGACGAATTTGCTATGaa                                        >  1:6209393/1‑76 (MQ=255)
                                     cGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGTTTTGTGCCTGCCGCACGTTTGAGGAATTTGCTTTGaa                                        >  1:2826395/1‑76 (MQ=255)
                                     cGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGTTTTGTGCCTGCCCCACGGTTGAGGAATTTGCTATGaa                                        >  2:5522205/1‑76 (MQ=255)
                                       cACCATTCACCATCTTCTGATGCCCGACGCGGCAGGGTTTGTGCCTGCCGCACGTTTGAGGAATTTGCTATGaaaa                                      >  1:7116106/1‑76 (MQ=255)
                                        agcATTCACCATATTCTGATGCCCGACGCGGCAGGGTTTTTGCCTGTCGCACGGTTGAGGCATTTCCTCTGaaaaa                                     >  2:3916868/3‑76 (MQ=255)
                                               cccATCTTCTGATGCCCGACGCGGCAGGTTTTGTCCCTGCCGCACGTTTGAGGAATTTGCTATGAAAAAACTGAtt                              <  2:5408581/75‑1 (MQ=255)
                                                 cATCTTCTGATGCCCGACGCGGCAGGTTTTGTGCCTGCCGCACGTTTTAGGAATTTGCTATGAAAAAACTGATtct                            >  2:4290148/1‑76 (MQ=255)
                                                  aTCTTCTGATGCCCGACGCGGCAGGGTTTTTGCCTGCCGCACCTTTGAGGAATTTGCCATGAAAAAAATGATtttc                           >  2:4984581/1‑76 (MQ=255)
                                                  aTCTTCTGATGCCCGACGCGGCAGGGTTTGTGGCTGCCGCACGGTTGAGGAATTTGCTATGAAAAAAATGATtctc                           >  2:4864649/1‑76 (MQ=255)
                                                  aTCTTCTGATGCCCGACGCGGCAGGGTTTGTGCCTGCCGCACGGTTGAGGAATTTGCTATGAAAAAACTGATtctc                           >  2:5087201/1‑76 (MQ=255)
                                                   tcttctGATGCCCGACGCGGCAGGTTTTGTGCCTGCCGCACGTTTGAGGAAATTGCTATGAACAAACTGATTCTCa                          >  2:3167582/1‑76 (MQ=255)
                                                     ttctGATGCCCGACGCGGCAGGTTTTGTGCCTGCCGCACGTTTGAGGAATTTGCTATGAAAAAACTGATTCTCATc                        >  2:7225710/1‑76 (MQ=255)
                                                     ttctGATGCCCGACGCGGCAGGTTTTGTGCCTGCCGCACGTTTGAGGAATTTGCTATGAAAAAACTGATTCTCATc                        <  2:7465927/76‑1 (MQ=255)
                                                      tctGATGCCCGACGCGGCAGGGTTTGTGCCTGCCGCACGTTTGAGGAATTTGCTATGAAAAAACTGATTCTCATCg                       >  1:4919865/1‑76 (MQ=255)
                                                        tGATGCCCGACGCGGCAGGTTTTGTGCCTGCCGCACGTTTGAGGAATTTGCTATGAAAAAACTGATTCTCATCGcc                     >  2:470803/1‑76 (MQ=255)
                                                        tGATCCCCGCCGCGGCAGGTTTTGTGCCTGCCGCACGTTTGAGGAATTTGCTATGAAAAAACTGATTCTCATCGcc                     <  1:3013578/76‑1 (MQ=255)
                                                         gATGCCCGACGCGGCAGGTTTTGTGCCTGCCGCACGTTTGAGGAATTTGCTATGAAAAAACTGATTCTCATCGCCa                    >  2:1800223/1‑76 (MQ=255)
                                                          atccccGCCGGGGCGGTTTTTGTCCCTGCCGCACGTTTGAGGAATTTGCTATGAAAAAACTGATTCTCATCGCCat                   <  1:5362674/73‑1 (MQ=255)
                                                          aTGCCCGACGCGGCAGGTTTTGTGCCTGCCGCACGTTTGAGGAATTTGCTATGAAAAAACTGATTCTCATCGCCat                   >  1:1038368/1‑76 (MQ=255)
                                                          aTGCCCGACGCGGCAGGGTTTTTGCCTGCCGCACGTTTTGGGAAATTGCTATTAAAAAACTGATTTTCATGGCCat                   >  2:7250307/1‑76 (MQ=255)
                                                           tGCCCGCCGCGGCAGTTTTTGTGCCTGCCGCACGTTTGAGGAATTTGCTATGAAAAAACTGATTCTCATCGCCatg                  <  2:681954/76‑1 (MQ=255)
                                                            gCCGGACGCGGCAGGTTTTGTGCCTGCCGCACGTTTGAGGAATTTGCTATGAAAAAACTGATTCTCATCGCCatga                 <  2:3001263/76‑1 (MQ=255)
                                                            gCCCGACGCGGCAGGTTTTGTGCCTGCCGCACGTTTGAGGAATTTGCTATGAAAAAACTGATTCTCATCGCCatga                 >  1:6911400/1‑76 (MQ=255)
                                                            gCCCGACGCGGCAGGGTTTGTGCCTGCCGCACGTTTTAGGAATTTGCTATGAAAAAACTGATTCTCATCGCCatga                 >  2:6428173/1‑76 (MQ=255)
                                                              ccGACGCGGCAGGTTTTGTGCCTGCCGCACGTTTGAGGAATTTGCTATGAAAAAACTGATTCTCATCGCCatgatg               >  1:7641538/1‑76 (MQ=255)
                                                               cGACGCGGCAGGGTTTGTGCCTGCCGCACGTTTGAGGAATTTGCTATGAAAAAACTGATTCTCATCGCCATGATgg              >  2:6344212/1‑76 (MQ=255)
                                                                gcggcgGCAGTTTTTGTGCCTGCCGCACGTTTGAGGAATTTGCTATGAAAAAACTGATTCTCATCGCCATGATGGc             <  2:6467078/73‑1 (MQ=255)
                                                                gACGCGGCAGGTTTTGTGCCTGCCGCACGTTTGAGGAATTTGCTATGAAAAAACTGATTCTCATCGCCATGATGGc             <  1:1137543/76‑1 (MQ=255)
                                                                gACGCGGCAGGTTTTGTGCCTGCCGCACGTTTGAGGAATTTGCTATGAAAAAACTGATTCTCATCGCCATGATGGc             >  2:3989399/1‑76 (MQ=255)
                                                                 aCGCGGCGGTTTTGGTCCCTGCCGCACGTTTGAGGAATTTGCTATGAAAAAACTGATTCTCATCGCCATGATGGCa            <  2:4069785/76‑1 (MQ=255)
                                                                 aCGCGGCAGGTTTTGTGCCTGCCGCACGTTTGAGGAATTTGCTATGAAAAAACTGATTCTCATCGCCATGATGGCa            >  2:2643020/1‑76 (MQ=255)
                                                                   gcgGCAGGTTTTGTCCCTGCCGCACGTTTGAGGAATTTGCTATGAAAAAACTGATTCTCATCGCCATGATGGCATc          <  2:1648504/76‑1 (MQ=255)
                                                                   gcgGCAGGGTTTTTGGCTGCCCCACGTTTTTGGGAATTTCTATGAAAAAACTGATTTTCATCGCCATGATGGCATc          >  2:5314760/1‑76 (MQ=255)
                                                                    cgGCAGGTTTTGTGCCTGCCGCACGTTTGAGGAATTTGCTATGAAAAAACTGATTCTCATCGCCATGATGGCATCg         >  2:1026364/1‑76 (MQ=255)
                                                                    cgGCAGGGTTTGTGCCTGGCGCGCGTTTTGAGAAATTGCTATGAAAAAACAGATTCTCATCGCCATGATGGCGTCg         >  2:1161989/1‑76 (MQ=255)
                                                                     ggCAGGTTTTGTGCCTGCCGCACGTTTGAGGAATTTGCTATGAAAAAACTGATTCTCATCGCCATGATGGCATAgg        <  2:5738699/76‑1 (MQ=255)
                                                                       cAGGTTTTGTGCCTGCCGCACGTTTGAGGAATTTGCTATGAAAAAACTGATTCTCATCGCCATGATGGCATCgggg      <  2:1236664/76‑1 (MQ=255)
                                                                        aGGTTTTGTGCCTGCCGCACGTTTGAGGAATTTGCTATGAAAAAACTGATTCTCATCGCCATGATGGCATCGGGGc     >  2:3835613/1‑76 (MQ=255)
                                                                        aGGTTTTGTGCCTGCCGCACGTTTGAGGAATTTGCTATGAAAAAACTGATTCTCATCGCCATGATGGCATCGGGGc     >  2:7379112/1‑76 (MQ=255)
                                                                        aGGTTTTGTGCCTGCCGCACGTTTGAGGAATTTGCTATGAAAAAACTGATTCTCATCGCCATGATGGCATCGGGGc     >  1:6591793/1‑76 (MQ=255)
                                                                          gTTTTGTGCCTGCCGCACGTTTGAGGAATTTGCTATGAAAAAACTGATTCTCATCGCCATGATGGCATCGGGGCTg   >  2:5628216/1‑76 (MQ=255)
                                                                          gTTTTGTGCCCGCCCGCACGTTTGAGGAATTTGCTATGAAAAAACTGATTCTCATCGCCATGATGGCATCGGGGCt    <  1:6017762/76‑1 (MQ=255)
                                                                           ttttGTGCCTGCCGCACGTTTGAGGAATTTGCTATGAAAAAACTGATTCTCATCGCCATGATGGCATCGGGGCTgg  >  1:1940753/1‑76 (MQ=255)
                                                                           ttttGTGCCTGCCGCACGTTTGAGGAATTTGCTATGAAAAAACTGATTCTCATCGCCATGATGGCATCGGGGCTgg  >  2:2309054/1‑76 (MQ=255)
                                                                           ttttGTGCCTGCCGCACGTTTGAGGAATTTGCTATGAAAAAACTGATTCTCATCGCCATGATGGCATCGGGGCTgg  >  2:7685949/1‑76 (MQ=255)
                                                                           |                                                                            
TCTTCCAGGATGAGCGTGACGTAAAATTTATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGTTTTGTGCCTGCCGCACGTTTGAGGAATTTGCTATGAAAAAACTGATTCTCATCGCCATGATGGCATCGGGGCTGG  >  REL606/691777‑691927

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