Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL606159,6540AC35.5% 537.8 / 46.0 262G27G (GGT→GGGyadNpredicted fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  ref base A (105/64);  new base C (0/93);  total (105/157)
Fisher's exact test for biased strand distribution p-value = 1.71e-28
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 340 total aligned reads displayed.

ACCATCTTTATTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAACCTAATTTTTTA  >  REL606/159579‑159729
                                                                           |                                                                           
aCCATCTTTATTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGacca                                                                             <  1:998829/76‑1 (MQ=255)
 ccATCTTTATTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGaccac                                                                            <  2:7679289/76‑1 (MQ=255)
   aTCTTTATTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGAccacca                                                                          >  2:8053799/1‑76 (MQ=255)
   aTCTTTATTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGAccacca                                                                          >  2:4914910/1‑76 (MQ=255)
    tCTTTATTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCAccat                                                                         <  1:4101673/76‑1 (MQ=255)
    tCTTTATTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCAccat                                                                         <  2:8709483/76‑1 (MQ=255)
     cTTTATTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCAccatc                                                                        >  2:1222301/1‑76 (MQ=255)
     cTTTATTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCAccatc                                                                        >  2:114982/1‑76 (MQ=255)
       ttattaCCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCAccatcca                                                                      >  1:8938351/1‑76 (MQ=255)
        tattcCCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCAccatccat                                                                     <  2:6394708/76‑1 (MQ=255)
        tattaCCGCCGTCAACGCGGGTTTCACAGGTTTTATCAACAACCAGACCACTGATATTTAACTGACCAccatccat                                                                     <  1:5770448/76‑1 (MQ=255)
        tattaCCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCAccatccat                                                                     >  1:4550874/1‑76 (MQ=255)
         tttaCCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCCCCATCCata                                                                    <  2:2698669/75‑1 (MQ=255)
         attaCCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCata                                                                    >  2:3094891/1‑76 (MQ=255)
          ttaCCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGAGCACCATCCatat                                                                   <  2:3012536/76‑1 (MQ=255)
          ttaCCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCatat                                                                   <  2:8859940/76‑1 (MQ=255)
          ttaCCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCatat                                                                   >  1:1045414/1‑76 (MQ=255)
           taCCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATc                                                                  >  2:4335427/1‑76 (MQ=255)
              cgccgTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCCCCATCCATATcagc                                                               <  2:5600110/76‑1 (MQ=255)
               gccgTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATcagca                                                              >  1:6258146/1‑76 (MQ=255)
                ccgTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCCCCATCCATATcagcag                                                             <  1:7150017/76‑1 (MQ=255)
                ccgTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATcagcag                                                             >  2:8597111/1‑76 (MQ=255)
                   tCAACGCGGGTTTCACAGGTTTTATCAACAACCAGCCCACTGATATTTAACTGCCCCCCATCCATATCAGCAGAAg                                                          <  2:748993/76‑1 (MQ=255)
                   tCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCCCCATCCATATCAGCAGAAg                                                          <  2:8438818/76‑1 (MQ=255)
                      aCGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAg                                                       <  1:693444/76‑1 (MQ=255)
                       cgcgGTTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCCCCATCCATATCAGCAGAAGCAGt                                                      <  2:8179613/76‑1 (MQ=255)
                       cgcgGGTTTCCAAGGTGTTACCAACATCCAGACCGGTGATTTTTAGCTACCCCCCATCCATATCAGCAGAAGCAGt                                                      <  1:6472007/76‑1 (MQ=255)
                       cgcgGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCCCCATCCATATCAGCAGAAGCAGt                                                      <  1:1669872/76‑1 (MQ=255)
                        gcgGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCCCCATCCATATCAGCAGAAGCAGTa                                                     <  2:7254661/76‑1 (MQ=255)
                        gcgGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTa                                                     >  2:3241212/1‑76 (MQ=255)
                         cgGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTAc                                                    >  2:1983220/1‑76 (MQ=255)
                         cgGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTAc                                                    >  1:5697215/1‑76 (MQ=255)
                            gTTTCACAGGTTTTATCAACAACCAGACCACTGATATTTAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTg                                                 <  2:2640862/76‑1 (MQ=255)
                                cacaGGTTTTAACAAGACCAAGCCCCCTGATATTTAACTGCCCACCATCCATATCAGCAGAAGCAGTACCTGCAAc                                             <  2:1433823/76‑1 (MQ=255)
                                 acaGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAAcc                                            >  2:3630659/1‑76 (MQ=255)
                                 acaGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAAcc                                            >  1:2191187/1‑76 (MQ=255)
                                  caGGTGTTATCAAAAACCAGCCCACTGATATTTACCTGTCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCa                                           <  2:7084707/76‑1 (MQ=255)
                                   gggTTTTATCAACAACCAGACCACTGATATTTAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCAt                                          <  2:5249224/75‑1 (MQ=255)
                                   aGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCAt                                          >  2:2277367/1‑76 (MQ=255)
                                   aGGAGTTATCAACAACCAGACCACGGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCAt                                          <  1:6676275/76‑1 (MQ=255)
                                    ggTGTTATCAACACCCAGACCCCTAATTTTTAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCAtt                                         <  2:412481/76‑1 (MQ=255)
                                     ttgtAATCAACAACCAGACCACTGATATTAAACTCACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTg                                        <  2:246997/75‑1 (MQ=255)
                                     gtgtTATCAACAACCAGCCCACTGATATTTAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTg                                        <  1:2579606/76‑1 (MQ=255)
                                        ttATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCa                                     >  2:7432042/1‑76 (MQ=255)
                                         taacaacacCCAGCCCACTGATATTTAACTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAg                                    <  2:5456333/73‑1 (MQ=255)
                                         tATCAACACCCAGCCCACTGATATTTAACTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAg                                    <  2:422761/76‑1 (MQ=255)
                                         tATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAg                                    >  2:7952481/1‑76 (MQ=255)
                                         tATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAg                                    >  1:3954264/1‑76 (MQ=255)
                                         aaTCAAAACCCGGCCCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAg                                    <  1:5904055/75‑1 (MQ=255)
                                          aTCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGc                                   >  1:266157/1‑76 (MQ=255)
                                          aTCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGc                                   >  2:4742422/1‑76 (MQ=255)
                                            caacacCCAGACCACTGATATTTAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCAt                                 <  1:2103128/76‑1 (MQ=255)
                                             aacaacCAGACCACTGATATTTAACTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATa                                <  1:8083115/76‑1 (MQ=255)
                                             aacaacCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATa                                <  1:1506361/76‑1 (MQ=255)
                                              acaccCAGCCCACTGATATTTACCTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATaa                               <  1:3989665/76‑1 (MQ=255)
                                              acaacCAGCCCACTGATATTTAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATaa                               <  1:1288181/76‑1 (MQ=255)
                                               cacccAGCCCACTGATATTTAACTGCCCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAg                              <  1:3328812/73‑1 (MQ=255)
                                               caacCAGACCACTGATTTTTAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAg                              <  1:2331555/76‑1 (MQ=255)
                                               caacCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAg                              >  2:3441858/1‑76 (MQ=255)
                                                aacCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGa                             >  1:1764311/1‑76 (MQ=255)
                                                 acCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGAc                            >  1:7637432/1‑76 (MQ=255)
                                                  ccAGCCCACTGATATTTAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGAcc                           <  1:5009498/76‑1 (MQ=255)
                                                   cAGCCCACTAATTTTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCa                          <  2:4612453/76‑1 (MQ=255)
                                                   cAGACCCCTTATATTTAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCa                          <  2:162778/76‑1 (MQ=255)
                                                    agcccACTGATATTTAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAg                         <  1:1222122/73‑1 (MQ=255)
                                                    agcccACTGATATTTAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAg                         <  1:8652524/73‑1 (MQ=255)
                                                    aGACCACTGATATTTAACTGCCCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAg                         <  1:7421837/76‑1 (MQ=255)
                                                    aGACCACTGATATTTAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAg                         <  2:599379/76‑1 (MQ=255)
                                                      cccACTGATATTTAACTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGAc                       <  2:3750501/75‑1 (MQ=255)
                                                        cACTGATATTTAACTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAg                     <  1:4956317/76‑1 (MQ=255)
                                                         cttGATTTTTAACTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAgg                    <  1:6850821/74‑1 (MQ=255)
                                                         aCTGATTTTTAACTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAgg                    <  2:447342/76‑1 (MQ=255)
                                                           tGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGc                  <  2:1755247/76‑1 (MQ=255)
                                                            tatatTAAACTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCa                 <  1:8434557/75‑1 (MQ=255)
                                                            gATATTTAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCa                 <  1:8740718/76‑1 (MQ=255)
                                                            gATATTTAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCa                 <  1:2858962/76‑1 (MQ=255)
                                                             atatTTAACTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCaa                <  1:3110775/76‑1 (MQ=255)
                                                             atatTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCaa                <  2:1804899/76‑1 (MQ=255)
                                                             atatTTAACGGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCaa                <  1:4754545/76‑1 (MQ=255)
                                                             aattttAACTGCCCACCACCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCaa                <  1:6147138/73‑1 (MQ=255)
                                                              tatTTAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCaaa               <  1:4194304/76‑1 (MQ=255)
                                                              tatTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCaaa               >  2:5342537/1‑76 (MQ=255)
                                                              tatTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCaaa               >  1:4313737/1‑76 (MQ=255)
                                                               ttttAACGGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCaaaa              <  1:5161818/75‑1 (MQ=255)
                                                               atTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCaaaa              >  1:8044402/1‑76 (MQ=255)
                                                               atTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCaaaa              >  1:5849208/1‑76 (MQ=255)
                                                               atTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCaaaa              >  2:612548/1‑76 (MQ=255)
                                                                ttgaaCTGCCCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAAc             <  2:5838360/73‑1 (MQ=255)
                                                                 ttAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAAcc            <  2:6881617/76‑1 (MQ=255)
                                                                 ttAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAAcc            <  1:783483/76‑1 (MQ=255)
                                                                 ttAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAAcc            >  2:6665660/1‑76 (MQ=255)
                                                                  tAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAACCt           <  1:5354101/76‑1 (MQ=255)
                                                                  tAACGGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAACCt           <  1:378906/76‑1 (MQ=255)
                                                                    aCTGCCCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAACCTaa         <  1:3379485/76‑1 (MQ=255)
                                                                     cggACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAACCTAAt        <  1:8872969/74‑1 (MQ=255)
                                                                     cTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAACCTAAt        <  2:2892182/76‑1 (MQ=255)
                                                                     cTGCCCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAACCTAAt        <  2:509139/76‑1 (MQ=255)
                                                                     cTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAACCTAAt        <  1:3561524/76‑1 (MQ=255)
                                                                           accaTCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAACCTAATTTTTTa  >  1:1988236/1‑76 (MQ=255)
                                                                           accaTCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAACCTAATTTTTTa  >  2:4023879/1‑76 (MQ=255)
                                                                           |                                                                           
ACCATCTTTATTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAACCTAATTTTTTA  >  REL606/159579‑159729

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