Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL606691,8520TG38.2% 187.7 / 18.9 110intergenic (+17/‑33)ybfM/ybfNpredicted outer membrane porin/predicted lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (33/35);  new base G (42/0);  total (75/35)
Fisher's exact test for biased strand distribution p-value = 4.93e-10
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 162 total aligned reads displayed.

TCTTCCAGGATGAGCGTGACGTAAAATTTATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGTTTTGTGCCTGCCGCACGTTTGAGGAATTTGCTATGAAAAAACTGATTCTCATCGCCATGATGGCATCGGGGCTGG  >  REL606/691777‑691927
                                                                           |                                                                            
tCTTCCAGGATGAGCGTGACGTAAAATTTATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGt                                                                              >  2:2195090/1‑76 (MQ=255)
 cTTCCAGGATGAGCGTGACGTAAAATTTATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGtt                                                                             <  1:4300228/76‑1 (MQ=255)
 cTTCCAGGATGAGCGTGACGTAAAATTTATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGtt                                                                             <  1:1415334/76‑1 (MQ=255)
  ttCCAGGATGAGCGTGACGTAAAATTTATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGttt                                                                            >  2:4858898/1‑76 (MQ=255)
   tCCAGGATGAGCGTGACGTAAAATTTATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGcttt                                                                           >  2:3858490/1‑76 (MQ=255)
     cAGGATGAGCGTGACGTAAAATTTATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGTTTtgt                                                                         >  1:3896922/1‑76 (MQ=255)
     cAGGATGAGCGTGACGTAAAATTTATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGGTTtgt                                                                         >  1:119485/1‑76 (MQ=255)
      aGGATGAGCGTGACGTAAAATTTATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGGTTtgtg                                                                        >  1:970912/1‑76 (MQ=255)
      aGGATGAGCGTGACGTAAAATTTATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGGTTtgtg                                                                        >  2:550099/1‑76 (MQ=255)
       ggATGAGCGTGACGTAAAATTTATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGGTTTGtgc                                                                       >  2:1192930/1‑76 (MQ=255)
         aTGAGCGTGACGTAAACGTTACGCGAATCTGACCATTCACCTACTTCTGATGCCCGACAGGGCAATTTTTGtgcct                                                                     <  2:4240255/76‑1 (MQ=255)
         aTGAGCGTGACGTAAAATTTATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGTTTTGtgcct                                                                     <  2:1948839/76‑1 (MQ=255)
         aTGAGCGTGACGTAAAATTTATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGTTTTGtgcct                                                                     >  1:867211/1‑76 (MQ=255)
         aTGAGCGTGACGTAAAATTTATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGGTTTGtggct                                                                     >  2:1825723/1‑76 (MQ=255)
          tGAGCGTGACGTAAAATTTATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGGTTTGtggctg                                                                    >  1:4218342/1‑76 (MQ=255)
             gCGTGCCGTAAAATTTATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGTTTTGTGCCTGCCg                                                                <  2:1404451/76‑1 (MQ=255)
               gTGACGTAAAATTGATGGTAATCGCACCTTTCACCTTTTTCTGATGCCCGCCGCGGCGGGTTTTGTCCCTGCCGCa                                                              <  2:1146658/76‑1 (MQ=255)
                tGACGTAAATTTTATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCGGGTTTTGTGCCTGCCGCAc                                                             <  1:2973914/76‑1 (MQ=255)
                tGACGTAAAATTTATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGTTTTGTGCCTGCCGCAc                                                             <  1:2413952/76‑1 (MQ=255)
                tGACGTAAAATTTATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGGTTTTTGCCTGACGccc                                                             >  1:2405404/1‑74 (MQ=255)
                 gACGTAAAATTTATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGGTTTGTGCCTGCCGCACg                                                            >  1:671087/1‑76 (MQ=255)
                  aCGTAAAATTTATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGGTTTGTGCCTGCCGCACGt                                                           >  2:2023387/1‑76 (MQ=255)
                   cgaaaatTTTAGGGAAACCGCCCCTTTCCCCATCTCCTGATCCCCGACGGGGCGGTTTTTGTCCCTGCCGCACGtt                                                          <  2:3014867/73‑1 (MQ=255)
                   cGTAAAATTTATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGTTTTGTGCCTGGCGCCCGtt                                                          >  1:1973518/1‑76 (MQ=255)
                   cGTAAAATTTATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGTTTTGTGCCTGCCGCACGtt                                                          <  2:2376267/76‑1 (MQ=255)
                   cGTAAAATTTATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGGTTTTTGCCTGCCGCAccgt                                                          >  1:717055/1‑73 (MQ=255)
                    gTAAATTTTATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGTTTTGTCCCTGCCGCACGttt                                                         <  2:4334281/76‑1 (MQ=255)
                    gTAAAATTTATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGTTTTGTGCCTGCCGCACGttt                                                         <  2:2413568/76‑1 (MQ=255)
                    gTAAAATTTATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGGTTTGTGCCTGCCGCCCGttt                                                         >  2:3740241/1‑76 (MQ=255)
                      aaaaTTTTTTGTAATCTGTCCATTCACTTTTTTCATAGGCCCGACCAGGCAGGTTTGGTTCTTGCCGCACGTTTGa                                                       <  2:1890911/76‑1 (MQ=255)
                      aaaaTTTATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGTTTTGTGCCTGCCGCACGTTTGa                                                       <  2:2304158/76‑1 (MQ=255)
                      aaaaTTTATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGTTTTGTGCCTGCCGCACGGTTgg                                                       >  2:4818544/1‑75 (MQ=255)
                      aaaaTTTATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGTTTTGTCCCTGCCGCACGTTTGa                                                       <  2:3117699/76‑1 (MQ=255)
                      aaaaTTTATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGTTTTGTCCCTGCCGCACGTTTGa                                                       <  2:2077368/76‑1 (MQ=255)
                      aaaaTTTATGGAATGCGCACATATCACCATCTTCTGACGCCCGCCGCGGCTGTTTTTGTGCCTGCCGCACGTTTGa                                                       <  2:4448110/76‑1 (MQ=255)
                       aaaTTTATGGTAAGCGCGCCATGCCCCTTCTCCTGATGCCCGACCCGCAAGGTTTTGTGCCTGCCGCACGTTTGAg                                                      <  2:2529318/76‑1 (MQ=255)
                         aTTTATGGTAATCGCCCCATTCACCATCTTCTGATGCCCGACGCGGCAGGTTTTGTGCCTGCCGCACGTTTGAGGa                                                    <  2:799519/76‑1 (MQ=255)
                         aTTTATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGTTTTGTGCCTGCCGCACGTTTGAGGa                                                    >  1:3367549/1‑76 (MQ=255)
                         aTTTATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGGTTTGTGCCTGCCGGCCCGTTTggga                                                     >  1:1736075/1‑73 (MQ=255)
                         aTTTATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGGTTTGTGCCTGCCGCACGTTTGAGGa                                                    >  1:3500752/1‑76 (MQ=255)
                         aTTTATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGGTTTGTGCCTGCCGCACGGTTGAGGa                                                    >  2:76766/1‑76 (MQ=255)
                           ttATGTTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGTTTTGTGCCTGCCGCACGTTTGAGGAAt                                                  <  2:4673115/76‑1 (MQ=255)
                           ttATGGTACTCTCACCATTCCCCATCTTCTGATGCCCGCCGGGTCCGGTTTTGGGCCCGCCGCACGTTTGAGGAAt                                                  <  2:4512486/76‑1 (MQ=255)
                           ttATGGTAATCGCCCCATTCACCATCTTCTGAGGCCCGACGCGGCAGGTTTTGTGCCTGCCGCACGTTTGAGGAAt                                                  <  1:4760068/76‑1 (MQ=255)
                            tATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGTTTTGTGCCTGCCGCACGTTTGAGGAAtt                                                 <  2:1240607/76‑1 (MQ=255)
                            tATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGGTTTGTGCCTGCCGCACGTTTGAGGAAtt                                                 >  1:975312/1‑76 (MQ=255)
                             aTGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGGTTTGTGCCTGCCCCACGTTTGAGGGAttt                                                >  2:781814/1‑76 (MQ=255)
                                gTAATCGCCCCATTCACCATCTTCTGATGCCCGACGCGGCAGGTTTTGTGCCTGCCGCACGTTTGAGGAATTTGCt                                             <  2:1547441/76‑1 (MQ=255)
                                gTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGTTTTGTGCCTGCCGCACGTTTGAGGAATTTGCt                                             >  2:5060623/1‑76 (MQ=255)
                                gTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGTTTTGTGCCTGCCGCACGTTTGAGGAATTTGCt                                             >  1:4968695/1‑76 (MQ=255)
                                gTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGTTTTGTGCCTGCCGCACGGTTGAGGAATTTGCt                                             >  2:3460219/1‑76 (MQ=255)
                                gTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGGTTTGTGCCTGCCGCACGTTTGAGGGATTTGCt                                             >  2:774129/1‑76 (MQ=255)
                                gTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGGTTTGTGCCTGCCGCACGTTTGAGGAATTTGCt                                             >  2:2352408/1‑76 (MQ=255)
                                 tAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGGTTTGTGCCTGCCGCACGTTTGAGGAATTTGCTa                                            >  2:4476337/1‑76 (MQ=255)
                                  aaTCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGGTTTGTGCCTGCCGCACGTTTGAGGGATTTGCTAt                                           >  2:617575/1‑76 (MQ=255)
                                  aaTCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGGTTTGTGCCTGCCGCACGTTTGAGGGATTTGCTAt                                           >  2:1746693/1‑76 (MQ=255)
                                   aTCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGTTTTGTGCCTGCCGCCCGTTTGAGGAATTTGCTATg                                          >  1:2288501/1‑76 (MQ=255)
                                   aTCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGTTTTGTGCCTGCCGCACGTTTGAGGGATTTTCTATg                                          >  2:2898396/1‑76 (MQ=255)
                                   aTCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGGTTTGTGCCTGCCGCACGTTTGAGGAATTTGCCATg                                          >  2:827056/1‑76 (MQ=255)
                                    tCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGGTTTGGTCCTGCCGCCCGGTTGAGGGATTTGCTATGa                                         >  1:791936/1‑76 (MQ=255)
                                      gCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGGTTTGTGCCTGCCGCACGTTTGAGGAATTTGCTATGaaa                                       >  2:4383561/1‑76 (MQ=255)
                                       cACCATTCACCATCTTCTGATGCCCGACGCGGCAGGTTTTGTGCCTGCCGCACGTTTGAGGAATTTGCTATGaaaa                                      <  2:3586632/76‑1 (MQ=255)
                                       cACCATTCACCATCTTCTGATGCCCGACGCGGCAGGGTTTGTGCCTGCCGCCCGTTTGAGGGATTTGCTATGaaaa                                      >  2:2712366/1‑76 (MQ=255)
                                        aCCATTCACCATCTTCTGATGCCCGACGCGGCAGGTTTTGTGCCTGCCGCACGTTTGAGGGATTTGCTATGaaaaa                                     >  1:1010360/1‑76 (MQ=255)
                                         cctttCACCATCTTCTGATGCCCGACGCGGCAGGTTTTGTGCCTGCCGCACGTTTGAGGAATTTGCTATGaaaaaa                                    <  2:4652393/73‑1 (MQ=255)
                                          ctttCCCCTTCTTCGGACCCCCGCCGGGGGTGTTTTTGTCCCTGCCGCACGTTTGAGGAATTTGCTATGAAAAAAc                                   <  1:1730788/74‑1 (MQ=255)
                                            ttCACCATCTTCTGATGCCCGACGCGGCAGGGTTTGTGCCTGCCGCACGTTTGAGGAATTTGCTATGAAAAAACTg                                 >  2:1481553/1‑76 (MQ=255)
                                              cACCATCTTCTGATGCCCGACGCGGCAGGTTTTGTGCCTGCCGCACGTTTGAGGAATTTGCTATGAAAAAAATGAt                               >  2:3601396/1‑76 (MQ=255)
                                               aCCATCTTCTGATGCCCGACGCGGCAGGGTTTGTGCCTGCCGCACGTTTGAGGAATTTGCTATGAAAAAAATGAtt                              >  2:3307469/1‑76 (MQ=255)
                                                ccATCTTCTGATGCCCGACGCGGCAGGTTTTGTGCCTGCCGCACGTTTGAGGAATTTGCTATGAAAAAACTGATtc                             >  1:1945358/1‑76 (MQ=255)
                                                 cATCTTCTGATGCCCGACGCGGCAGGGTTTGTGCCTGCCGCCCGGTTGAGGAATTTGCTATGAAAAAACTGAtttt                            >  1:4335640/1‑74 (MQ=255)
                                                    cttctGATGCACGACGCGGCAGGTTTTGTGCCTGCCGCCCCTTTGAGGAATTTGCCATGAACAAAATGAATCTAAt                         >  1:3858462/1‑76 (MQ=255)
                                                     ttctGATGCCCGACGCGGCAGGTTTTGTGCCTGCCGCACGTTTGAGGAATTTGCTATGAAAAAACTGATTCTCATc                        <  2:2037101/76‑1 (MQ=255)
                                                     ttctGATGCCCGACGCGGCAGGGTTTGTGCCTGCCCCCCGGTTTAGGGATTTGCTATGAAAAAACTGGTTTTCATc                        >  2:4875332/1‑76 (MQ=255)
                                                      tctGATGCCCGACGCGGCAGGTTTTGTGGCTGCCGCACGGTTGAGGGATTTGCTATGAAAAAAATGATTCTCAtgg                       >  2:690177/1‑74 (MQ=255)
                                                      tctGATGCCCGACGCGGCAGGGTTTGGGGCCGCCGCACGGTTGAGGAATTTGCTATGAAAAAACTGATTCTCATCg                       >  1:1066298/1‑76 (MQ=255)
                                                       ctGATGCCCGACGCGGCAGGTTTTGTGCCTGCCGCACGTTTGAGGAATTTGCTATGAAAAAACTGATTCTCATCGc                      >  1:3031785/1‑76 (MQ=255)
                                                       ctGANGCCCGACGCGGCAGGTTTTGTGCCTGCCGCACGTTTGAGGAATTTGCTATGAAAAAACTGATTCTCATCGc                      >  2:3797640/1‑76 (MQ=255)
                                                           tccccGACGCGGCAGGTTTTGTGCCTGCCGCACGTTTGAGGAATTTGCTATGAAAAAACTGATTCTCATCGCCatg                  <  2:2370070/74‑1 (MQ=255)
                                                           tGCCCGACGCGGCAGGTTTTGTGCCTGCCGCACGTTTGAGGAATTTGCTATGAAAAAACTGATTCTCATCGCCatg                  >  2:4968799/1‑76 (MQ=255)
                                                             cccGACGCGGCAGGTTTTGTCCCTGCCGCACGTTTGAGGAATTTGCTATGAAAAAACTGATTCTCATCGCCatgat                <  1:4876165/76‑1 (MQ=255)
                                                              ccgccgCGGGAGTTTTTGTCCCCGCCGAACGTTTGAGGAATTTGCTATGAAAAAACTGATTCTCATCGCCatgatg               <  1:827056/76‑1 (MQ=255)
                                                              ccGACGCGGCAGGTTTTGTGCCTGCCGCACCTTTGAGGAATTTGCTATGAAAAAACTGATTGTCATCGCCatgatg               >  2:4392395/1‑76 (MQ=255)
                                                               cgccgcgGCAGTTTTTGTGCCGGCCTCACGTTTGAGGAATTTGCTATGAAAAAACTGATTCTCATCGCCATGATgg              <  2:3524317/73‑1 (MQ=255)
                                                               cGACGCGGCAGGTTTTGTGCCTGCCGCACGTTTGAGGAATTTGCTATGAAAAAACTGATTCTCATCGCCATGATgg              <  2:4442387/76‑1 (MQ=255)
                                                               cGACGCGGCAGGGTTTGTGCCTGCCGCACGTTTGAGGAATTTGCTATGAAAAAACTGATTCTCATCGCCATGATgg              >  2:1178650/1‑76 (MQ=255)
                                                               cGACGCGGCAGGGTTTGTGCCTGCCGCACGTTTGAGGAATTTGCTATGAAAAAACTGATTCTAATCGCCATGATgg              >  2:1032699/1‑76 (MQ=255)
                                                                gACGCGGCAGGTTTTGGGCCTGCCGCACGGTTGAGGAATTTGCTCTGAAAAAACTGATTCTCATCGCCATGATGGc             >  2:2723988/1‑76 (MQ=255)
                                                                gACGCGGCAGGGTTTGTGCCTGCCGCACGTTTGAGGAATTTGCTATGAAAAAACTGATTCTCATCGCCATGATGGc             >  2:4060019/1‑76 (MQ=255)
                                                                  cgcgGCAGGTTTTGTGCCTGCCGCACGTTTGAGGAATTTGCTATGAAAAAACTGATTCTCATCGCCATGATGGCAt           <  1:2515546/76‑1 (MQ=255)
                                                                    cgGCAGGTTTTGTGCCTGCCGCACGTTTGAGGAATTTGCTATGAAAAAACTGATTCTCATCGCCATGATGGCATCg         >  1:1311029/1‑76 (MQ=255)
                                                                     ggCAGGTTTTGTGCCTGCCGCACGTTTGAGGAATTTGCTATGAAAAAACTGATTCTCATCGCCATGATGGCATCgg        >  2:110882/1‑76 (MQ=255)
                                                                      gcgggTTTTGTGCCTGCCGCACGTTTGAGGAATTTGCTATGAAAAAACTGATTCTCATCGCCATGATGGCATCggg       <  2:2110465/73‑1 (MQ=255)
                                                                      gcgggTTTGGTGCCTGCCGCACGTTTGAGGAATTTGCTATGAAAAAACTGATTCTCATCGCCATGATGGCATCggg       <  2:1395498/73‑1 (MQ=255)
                                                                       cAGGTTTTGTGCCTGCCGCACGTTTGAGGAATTTGCTATGAAAAAACTGATTCTCATCGCCATGATGGCATCgggg      >  2:2925137/1‑76 (MQ=255)
                                                                        aGGTTTTGTGCCTGCCGCACGTTTGAGGAATTTGCTATGAAAAAACTGATTCTCATCGCCATGATGGCATCGGGGc     <  2:1034224/76‑1 (MQ=255)
                                                                        aGGTTTTGTGCCTGCCGCACGTTTGAGGAATTTGCTATGAAAAAACTGATTCTCATCGCCATGATGGCATCGGGGc     >  2:620768/1‑76 (MQ=255)
                                                                          gTTTTGTGCCTGCCGCACGTTTGAGGAATTTGCTATGAAAAAACTGATTCTCATCGCCATGATGGCATCGGGGCTg   >  1:806666/1‑76 (MQ=255)
                                                                           ttttGTGCCTGCCGCACGTTTGAGGAATTTGTTATGAAAAAACTGATTCTCATCGCCATGATGGCATCGGGGCTgg  >  1:3710597/1‑76 (MQ=255)
                                                                           tttgGTCCCTGCCGCACGTTTGAGGAATTTGCTATGAAAAAACTGATTCTCATCGCCATGATGGCATCGGGGCTgg  <  2:1329173/76‑1 (MQ=255)
                                                                           |                                                                            
TCTTCCAGGATGAGCGTGACGTAAAATTTATGGTAATCGCACCATTCACCATCTTCTGATGCCCGACGCGGCAGGTTTTGTGCCTGCCGCACGTTTGAGGAATTTGCTATGAAAAAACTGATTCTCATCGCCATGATGGCATCGGGGCTGG  >  REL606/691777‑691927

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