Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL606159,6540AC35.3% 517.7 / 44.4 258G27G (GGT→GGGyadNpredicted fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  ref base A (101/66);  new base C (0/91);  total (101/157)
Fisher's exact test for biased strand distribution p-value = 6.09e-27
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 335 total aligned reads displayed.

ACCATCTTTATTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAACCTAATTTTTTA  >  REL606/159579‑159729
                                                                           |                                                                           
aCCATCTTTATTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGacca                                                                             >  2:1929715/1‑76 (MQ=255)
aCCATCTTTATTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGacca                                                                             <  2:3504095/76‑1 (MQ=255)
aCCATCTTTATTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGacca                                                                             >  1:3467115/1‑76 (MQ=255)
 ccATCTTTATTCCCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGaccac                                                                            <  1:7028457/76‑1 (MQ=255)
   aTCTTTATTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGAccacca                                                                          >  1:4514414/1‑76 (MQ=255)
   aTCTTTATTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGAccacca                                                                          >  2:6658792/1‑76 (MQ=255)
   aTCTTTATTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGAccacca                                                                          <  2:2988422/76‑1 (MQ=255)
     cTTTATTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCAccatc                                                                        >  2:7581142/1‑76 (MQ=255)
      tttATTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCAccatcc                                                                       >  1:1541996/1‑76 (MQ=255)
       ttattaCCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACAAcaatcca                                                                      >  1:1429682/1‑76 (MQ=255)
          ttaCCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCatat                                                                   >  1:1888467/1‑76 (MQ=255)
          ttaCCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCatat                                                                   <  2:3900358/76‑1 (MQ=255)
           tcccgccgTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCCCCATCCATATc                                                                  <  2:4893323/74‑1 (MQ=255)
           taCCGCCGTCAACGCGGGTTTCACAGGTGGTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATc                                                                  >  2:6410994/1‑76 (MQ=255)
            aCCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATca                                                                 >  2:1536804/1‑76 (MQ=255)
            aCCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATca                                                                 >  2:6288909/1‑76 (MQ=255)
             ccgccgTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCCCCATCCATATcag                                                                <  2:233295/76‑1 (MQ=255)
              cgccgTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATcagc                                                               >  1:164271/1‑76 (MQ=255)
              cgccgTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATcagc                                                               >  2:6936785/1‑76 (MQ=255)
              cgccgTCAACGCGGGTTTCACAGGTGTTATAAACAACCAGACCACTGATATTTAACTGACCACCATCCATATcagc                                                               <  2:7196466/76‑1 (MQ=255)
                ccgTCAAGGCGGGTTTCACAGGTGTTATCAACAACCAGCCCACTGATATTTAACTGACCCCCATCCATATcagcag                                                             <  2:877050/76‑1 (MQ=255)
                ccgTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCCCCATCCATATcagcag                                                             <  2:5896887/76‑1 (MQ=255)
                ccgTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATcagcag                                                             >  1:2227095/1‑76 (MQ=255)
                  gTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGCCCACTGATATTTAACTGACCCCCATCCATATCAGCAGaa                                                           <  2:1055413/76‑1 (MQ=255)
                  gTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGaa                                                           >  1:868692/1‑76 (MQ=255)
                  gTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGaa                                                           <  1:5349550/76‑1 (MQ=255)
                   tCAACGCGGTTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAg                                                          <  2:4645960/76‑1 (MQ=255)
                   tCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAg                                                          <  2:426691/76‑1 (MQ=255)
                   tCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAg                                                          <  2:5471273/76‑1 (MQ=255)
                   tCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAg                                                          >  2:5799759/1‑76 (MQ=255)
                    cAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGc                                                         >  2:2615658/1‑76 (MQ=255)
                     aaCGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCa                                                        >  1:7623928/1‑76 (MQ=255)
                      aCGCGGTTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAg                                                       <  1:4398756/76‑1 (MQ=255)
                      aCGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCCCCATCCATATCAGCAGAAGCAg                                                       <  2:7766900/76‑1 (MQ=255)
                         cgGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTAc                                                    >  2:5846453/1‑76 (MQ=255)
                           gttttCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCt                                                  <  1:2915785/74‑1 (MQ=255)
                           ggTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCt                                                  <  1:5959966/76‑1 (MQ=255)
                           ggTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCt                                                  <  1:2048579/76‑1 (MQ=255)
                            gTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCTTATCAGCAGAAGCAGTACCTg                                                 <  1:6188221/76‑1 (MQ=255)
                               tCACAGGTGTTATCAACAACCAGACCACTGATTTTTAACTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCaa                                              <  1:5203440/76‑1 (MQ=255)
                                cacaGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAAc                                             >  2:5740940/1‑76 (MQ=255)
                                 acaGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAAcc                                            <  2:5249852/76‑1 (MQ=255)
                                   aGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCAt                                          >  2:1710789/1‑76 (MQ=255)
                                    ggTGTTATCAACAACCAGCCCACTGATATTTAACTGCCCACCATCCATATCAGCAGAAGCAGTACCTGCAACCAtt                                         <  1:6587452/76‑1 (MQ=255)
                                     gtgtTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTg                                        >  1:6930442/1‑76 (MQ=255)
                                       tttATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGcc                                      >  2:7885220/2‑76 (MQ=255)
                                       gtTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGcc                                      >  1:2603290/1‑76 (MQ=255)
                                        ttATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCa                                     >  1:3867953/1‑76 (MQ=255)
                                        ttATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCa                                     >  1:6017598/1‑76 (MQ=255)
                                         tataAACAACCAGACCACTGATATTTAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAg                                    <  2:4094051/76‑1 (MQ=255)
                                             aacaccCACACCACTGATTTTTAACTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATa                                <  2:5034023/76‑1 (MQ=255)
                                              acaacCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATaa                               <  2:708269/76‑1 (MQ=255)
                                               caacCAGCCCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAg                              <  1:1017441/76‑1 (MQ=255)
                                                acccAGCCCACTGATATTTACCGGCCCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGa                             <  1:5899696/74‑1 (MQ=255)
                                                aacCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGa                             >  2:1804616/1‑76 (MQ=255)
                                                  ccAGACCCCTAATATTAACCTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGAcc                           <  1:3323172/76‑1 (MQ=255)
                                                  ccAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGAcc                           >  2:3579443/1‑76 (MQ=255)
                                                   cAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCa                          >  1:6818687/1‑76 (MQ=255)
                                                    aGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAg                         <  2:7426655/76‑1 (MQ=255)
                                                    aGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAg                         >  1:4132325/1‑76 (MQ=255)
                                                    aGACCACTGATATTTAACGGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAg                         <  1:4847341/76‑1 (MQ=255)
                                                       ccACTGATTTTTAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACa                      <  2:657083/76‑1 (MQ=255)
                                                       ccACTGATATTTAACTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACa                      <  1:562145/76‑1 (MQ=255)
                                                       ccACTGATATTAAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACa                      <  2:7528705/76‑1 (MQ=255)
                                                        cACTGATATTTAACTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAg                     <  2:1348651/76‑1 (MQ=255)
                                                        cACTGATATTTAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAg                     <  1:7499809/76‑1 (MQ=255)
                                                        cACTGATATTTAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAg                     <  1:7374297/76‑1 (MQ=255)
                                                        cACTGATATTAAACGGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAg                     <  1:4018636/76‑1 (MQ=255)
                                                        cACTAATATTTAACTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAg                     <  1:72130/76‑1 (MQ=255)
                                                         aCTGATATTTAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAgg                    <  2:3696097/76‑1 (MQ=255)
                                                          cTGATTTTTAACTAACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAggg                   <  2:2425168/76‑1 (MQ=255)
                                                          cTGATATTTAACTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAggg                   <  1:1149842/76‑1 (MQ=255)
                                                          cTGATATTTAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAggg                   <  1:3599516/76‑1 (MQ=255)
                                                          cTGATATTTAACGGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAggg                   <  1:5627502/76‑1 (MQ=255)
                                                          cTGATATTAACCTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAggg                   <  2:2080232/76‑1 (MQ=255)
                                                            gATTTTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCa                 <  2:4815715/76‑1 (MQ=255)
                                                             atatTTACCCGCCCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCaa                <  1:4994943/76‑1 (MQ=255)
                                                             atatTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCaa                >  1:1057349/1‑76 (MQ=255)
                                                             atatTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCaa                >  1:5764932/1‑76 (MQ=255)
                                                             atatTAACCTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCaa                <  2:4797064/76‑1 (MQ=255)
                                                              tatTTACCGGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCaaa               <  1:2141799/76‑1 (MQ=255)
                                                              tatTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCaaa               >  1:7950235/1‑76 (MQ=255)
                                                              tatTAAACTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCaaa               <  2:4422220/76‑1 (MQ=255)
                                                              tatTAAACTGCCCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCaaa               <  2:5625106/76‑1 (MQ=255)
                                                              tatTAAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCaaa               <  1:624989/76‑1 (MQ=255)
                                                               atTTAACTGCCCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCaaaa              <  2:3287199/76‑1 (MQ=255)
                                                                tttAACTCCCCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAAc             <  2:1093125/76‑1 (MQ=255)
                                                                 ttAACTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAAcc            <  1:654091/76‑1 (MQ=255)
                                                                 ttAACTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAAcc            <  1:2757826/76‑1 (MQ=255)
                                                                 ttAACTGCCCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATaa                               <  1:2430981/57‑1 (MQ=255)
                                                                  taccTGCCCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAACCt           <  1:6281837/73‑1 (MQ=255)
                                                                  tAACTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAACCt           <  1:5432164/76‑1 (MQ=255)
                                                                  tAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAACCt           <  2:5983378/76‑1 (MQ=255)
                                                                   aaCTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAACCTa          <  1:963172/76‑1 (MQ=255)
                                                                   aaCTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAACCTa          <  1:7739363/76‑1 (MQ=255)
                                                                    aCTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAACCTaa         <  1:3119386/76‑1 (MQ=255)
                                                                    aCTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAACCTaa         >  1:3908110/1‑76 (MQ=255)
                                                                       cacccccATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAACCTAAttt      <  2:214565/75‑1 (MQ=255)
                                                                        accaccATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAACCTAAtttt     <  2:7710699/76‑1 (MQ=255)
                                                                           accaTCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAACCTAATTTTTTa  >  2:805390/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: