Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL606159,6540AC36.1% 350.9 / 27.8 169G27G (GGT→GGGyadNpredicted fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  ref base A (76/32);  new base C (0/61);  total (76/93)
Fisher's exact test for biased strand distribution p-value = 1.54e-22
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 215 total aligned reads displayed.

ACCATCTTTATTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAACCTAATTTTTTA  >  REL606/159579‑159729
                                                                           |                                                                           
aCCATCTTTATTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGacca                                                                             <  1:1086466/76‑1 (MQ=255)
  cATCTTTATTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGaccacc                                                                           <  1:5162883/76‑1 (MQ=255)
   aTCTTTATTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGAccacca                                                                          >  2:3826050/1‑76 (MQ=255)
     cTTTATTCCCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCAccatc                                                                        <  1:4205223/76‑1 (MQ=255)
     cTTTATTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCAccatc                                                                        >  2:2825630/1‑76 (MQ=255)
     cTTTATTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCAccatc                                                                        >  2:4706149/1‑76 (MQ=255)
     cTTTATTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCAccatc                                                                        <  1:1209256/76‑1 (MQ=255)
       ttattaCCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCAccatcca                                                                      <  2:858750/76‑1 (MQ=255)
       ttattaCCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACAAccatcca                                                                      >  2:1265522/1‑76 (MQ=255)
        tattaCCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCAccatccat                                                                     <  1:3820235/76‑1 (MQ=255)
        tattaCCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCAccatccat                                                                     <  1:3684325/76‑1 (MQ=255)
        tattaCCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCAccatccat                                                                     >  1:2172968/1‑76 (MQ=255)
         attaCCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCata                                                                    >  2:3521858/1‑76 (MQ=255)
          ttcccgccgTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCatat                                                                   <  2:333126/73‑1 (MQ=255)
           tcccgccgTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCCCCATCCATATc                                                                  <  2:5189957/74‑1 (MQ=255)
           tcccgccgTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATc                                                                  <  1:1306112/74‑1 (MQ=255)
            cccgccgTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCCCCATCCATATca                                                                 <  2:4686262/75‑1 (MQ=255)
            aCCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATca                                                                 >  2:2106005/1‑76 (MQ=255)
             ccgccgTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCCCCATCCATATcag                                                                <  2:4363002/76‑1 (MQ=255)
             ccgccgTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCCCCATCCATATcag                                                                <  2:4055585/76‑1 (MQ=255)
               gccgTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATcagca                                                              >  1:4394030/1‑76 (MQ=255)
                ccgTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCCCCATCCATATcagcag                                                             <  2:1100251/76‑1 (MQ=255)
                ccgTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCCCCATCCATATcagaag                                                             <  2:4088938/76‑1 (MQ=255)
                 cgTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGa                                                            <  2:2055657/76‑1 (MQ=255)
                  gTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGaa                                                           >  1:2535713/1‑76 (MQ=255)
                  gTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGaa                                                           >  2:520305/1‑76 (MQ=255)
                   tCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAg                                                          >  2:4452148/1‑76 (MQ=255)
                   tCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAg                                                          >  1:2888236/1‑76 (MQ=255)
                    cAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCCCCATCCATATCAGCAGAAGc                                                         <  2:4891731/76‑1 (MQ=255)
                    cAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCCCCATCCATATCAGCAGAAGc                                                         <  1:876530/76‑1 (MQ=255)
                    aaaCGCGGGTTTCACAGGTGTTATCAACAACCAACCCCCGGATATTAAACTCCCCACCATCCATATCAGCAGAAGc                                                         <  2:1883399/75‑1 (MQ=255)
                       cgcgGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGt                                                      <  1:3760399/76‑1 (MQ=255)
                         cgGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTAc                                                    >  2:3421543/1‑76 (MQ=255)
                           gttttCACAGGTTTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCt                                                  <  2:677324/74‑1 (MQ=255)
                           ggTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCt                                                  <  1:2031878/76‑1 (MQ=255)
                             tttCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGc                                                >  2:1508638/1‑76 (MQ=255)
                                cacaGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAAc                                             >  1:2612091/1‑76 (MQ=255)
                                  caGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCa                                           >  1:4549737/1‑76 (MQ=255)
                                  caGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCa                                           >  1:841619/1‑76 (MQ=255)
                                   aGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCAt                                          <  1:2106533/76‑1 (MQ=255)
                                   aGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCAt                                          <  1:4813187/76‑1 (MQ=255)
                                     gtgtTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTg                                        >  1:1063405/1‑76 (MQ=255)
                                       tttATCAACAACCAGCCCACTGATTTTTAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGcc                                      <  2:4512484/75‑1 (MQ=255)
                                       gtTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGcc                                      >  1:4704150/1‑76 (MQ=255)
                                       gtTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGcc                                      <  1:3749539/76‑1 (MQ=255)
                                         tATCAACACCCGGACCACTGAATTTTAACGGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAg                                    <  2:1041769/76‑1 (MQ=255)
                                         tATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAg                                    >  2:5178717/1‑76 (MQ=255)
                                          aTCAACACCCAGACCCCTGATATTTAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGc                                   <  1:3912799/76‑1 (MQ=255)
                                           tCAACACCCAGCCCACTGATATTTAACTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCa                                  <  2:3022393/76‑1 (MQ=255)
                                           ccaacaaCCAGACCACTGATATTAAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCa                                  <  2:5112725/75‑1 (MQ=255)
                                            caacaaCCAGCCCACTGTTATTTAACTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCAt                                 <  2:2324853/76‑1 (MQ=255)
                                            caacaaCCAGACCACTGATATTTAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCAt                                 <  2:3058153/76‑1 (MQ=255)
                                            caacaaCCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCAt                                 >  2:2768215/1‑76 (MQ=255)
                                              acaccCAGACCACTGATATTTAACTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATaa                               <  1:529589/76‑1 (MQ=255)
                                              acaacCAGACCACTGATATTTAACTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATaa                               <  1:3442413/76‑1 (MQ=255)
                                              acaacCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATaa                               >  2:4292038/1‑76 (MQ=255)
                                               caacCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAg                              >  2:3108760/1‑76 (MQ=255)
                                               caacCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAg                              >  2:1451008/1‑76 (MQ=255)
                                               caacCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAg                              >  2:1370689/1‑76 (MQ=255)
                                               aacccAGCCCACTGATATTTAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAg                              <  1:161604/73‑1 (MQ=255)
                                               aaacCAGACCACTGATATTTAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAg                              <  1:4315060/75‑1 (MQ=255)
                                                aacCAGACCACTGATATTTAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGa                             <  1:1018951/76‑1 (MQ=255)
                                                aacCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGa                             >  2:3858486/1‑76 (MQ=255)
                                                  ccAGCCCACTGATATTTAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGAcc                           <  1:3588079/76‑1 (MQ=255)
                                                     gcccACTGATATTAAACTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGa                        <  1:4432150/74‑1 (MQ=255)
                                                      cccACTTATATTTAACCGCCCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGAc                       <  1:4786612/75‑1 (MQ=255)
                                                      aCCACTAATTTTTAACTGCCCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGAc                       <  2:476334/76‑1 (MQ=255)
                                                       ccACTGATATTTAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACa                      <  2:4930703/76‑1 (MQ=255)
                                                        cACTGATTTTTAGCTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAg                     <  1:4014140/76‑1 (MQ=255)
                                                        cACTGATATTTAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAg                     <  2:3220516/76‑1 (MQ=255)
                                                        cACTGATATTAAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAg                     <  2:1854086/76‑1 (MQ=255)
                                                         aCTGATATTTACCTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAgg                    <  1:1788276/76‑1 (MQ=255)
                                                         aCTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAgg                    >  2:3052089/1‑76 (MQ=255)
                                                         aCTGATATTTAACGACCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAgg                    <  2:2696459/76‑1 (MQ=255)
                                                          cTGATATTTAACTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAggg                   <  1:4077442/76‑1 (MQ=255)
                                                           tGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGc                  >  2:3197527/1‑76 (MQ=255)
                                                            gATATTTAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCa                 <  1:2478587/76‑1 (MQ=255)
                                                            gATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCa                 >  2:4563043/1‑76 (MQ=255)
                                                            gATATTTAACGGACCGCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCa                 <  2:1702643/76‑1 (MQ=255)
                                                            gATATTAACCTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCa                 <  1:3003344/76‑1 (MQ=255)
                                                             atatTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCaa                >  1:3510480/1‑76 (MQ=255)
                                                             atatTTAACTCCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCaa                <  1:22002/76‑1 (MQ=255)
                                                             atatTAAACTGCCCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCaa                <  2:250053/76‑1 (MQ=255)
                                                              tttttAACTACCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCaaa               <  1:1324989/74‑1 (MQ=255)
                                                               ttttAACGGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCaaaa              <  1:3197746/75‑1 (MQ=255)
                                                               atTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCaaaa              >  1:1976856/1‑76 (MQ=255)
                                                                ttaaaCTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAAc             <  1:590664/73‑1 (MQ=255)
                                                                 ttAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAAcc            >  2:2137906/1‑76 (MQ=255)
                                                                 ttAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAAcc            >  1:4611117/1‑76 (MQ=255)
                                                                  tAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAACCt           <  2:2457916/76‑1 (MQ=255)
                                                                   accTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAACCTa          <  1:3137317/74‑1 (MQ=255)
                                                                   aaCTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAACCTa          <  1:510173/76‑1 (MQ=255)
                                                                    ccTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAACCTaa         <  2:239422/75‑1 (MQ=255)
                                                                     cTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAACCTAAt        <  2:4394030/76‑1 (MQ=255)
                                                                     cTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAACCTAAt        >  1:5065596/1‑76 (MQ=255)
                                                                     cTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAACCTAAt        >  2:561815/1‑76 (MQ=255)
                                                                      tGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAACCTAAtt       >  2:3671645/1‑76 (MQ=255)
                                                                       gACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAACCTAAttt      >  1:1893014/1‑76 (MQ=255)
                                                                       gACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAACCTAAttt      <  2:1941458/76‑1 (MQ=255)
                                                                         ccaccaTCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAACCTAAttttt    >  2:3137347/1‑76 (MQ=255)
                                                                           |                                                                           
ACCATCTTTATTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAACCTAATTTTTTA  >  REL606/159579‑159729

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