Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL606154,4970TG26.2% 143.2 / 10.8 61A191A (GCA→GCCyadLpredicted fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  ref base T (9/36);  new base G (16/0);  total (25/36)
Fisher's exact test for biased strand distribution p-value = 1.01e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TACGTTGAGTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCGGATAA  >  REL606/154422‑154572
                                                                           |                                                                           
tACGTTGAGTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGt                                                                             <  1:68930/76‑1 (MQ=255)
   gTTGAGTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCa                                                                          >  2:873581/1‑76 (MQ=255)
     tGAGTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGgaga                                                                        >  2:662915/1‑76 (MQ=255)
       aGTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGaaa                                                                      <  2:1824048/76‑1 (MQ=255)
         tGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAAt                                                                    <  2:959247/76‑1 (MQ=255)
            gATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAANTGTGGCAGGTGCAGAAAATTTa                                                                 <  2:258100/76‑1 (MQ=255)
            gATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTa                                                                 >  1:1639606/1‑76 (MQ=255)
            gATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAAATTa                                                                 >  1:132656/1‑76 (MQ=255)
              tGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTAcc                                                               <  1:469013/76‑1 (MQ=255)
              tGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTAcc                                                               <  1:1359726/76‑1 (MQ=255)
               gCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAAATTAccc                                                              >  1:1222519/1‑76 (MQ=255)
                cATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCg                                                             <  1:115487/76‑1 (MQ=255)
                cATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCg                                                             <  2:1804872/76‑1 (MQ=255)
                 atatatAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGc                                                            <  2:1541118/76‑1 (MQ=255)
                 atatatAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAAATTAccccc                                                            >  1:214120/1‑74 (MQ=255)
                    tataGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAAATTACCCGCAgg                                                         >  1:135972/1‑75 (MQ=255)
                     ataGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTg                                                        >  1:126156/1‑76 (MQ=255)
                     ataGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTg                                                        <  2:1827006/76‑1 (MQ=255)
                       aGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGGAGAAAATTTACCCGCCGCGAc                                                      >  2:1185001/1‑76 (MQ=255)
                        gTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACa                                                     <  1:628212/76‑1 (MQ=255)
                          tCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACAtt                                                   <  1:36720/76‑1 (MQ=255)
                           ccTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGTGACATTg                                                  >  1:239332/1‑76 (MQ=255)
                            cTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGGGGAAAATTTACCCCCAGCGACATTGc                                                 >  2:516904/1‑76 (MQ=255)
                             tGGATCTGTCAGTCTGTTAGTTATAGGTTACTGAAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTgat                                                <  1:1637098/76‑3 (MQ=255)
                               gATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTa                                              <  2:871469/76‑1 (MQ=255)
                                 tCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACg                                            <  1:1529458/76‑1 (MQ=255)
                                  cTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCGGGGGCATTGCTTATgg                                           >  2:1908261/1‑76 (MQ=255)
                                    gTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGGGGAAAAATTAACCCCGGAGACATTGCTTTCGGtt                                         >  2:823513/1‑76 (MQ=255)
                                    gTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGGGACATTGCTTTCgggt                                         >  1:1371447/1‑74 (MQ=255)
                                     tCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGttt                                        <  1:1243520/76‑1 (MQ=255)
                                     tCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGTGACATTGCTTACGGttt                                        >  1:1767582/1‑76 (MQ=255)
                                      cAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGtttt                                       <  2:964050/76‑1 (MQ=255)
                                      cAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTAACCGCAGTGACATTGCTTACGGtttt                                       >  1:15112/1‑76 (MQ=255)
                                        gTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTAt                                     <  2:1303978/76‑1 (MQ=255)
                                         tCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTAtt                                    <  2:289679/76‑1 (MQ=255)
                                         tCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGGGACATTGCTTACGGTtttttt                                    >  1:1882869/1‑76 (MQ=255)
                                          cTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTAttt                                   <  1:29009/76‑1 (MQ=255)
                                           tGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTAtttt                                  <  2:1738865/76‑1 (MQ=255)
                                           tGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTAtttt                                  <  2:1240561/76‑1 (MQ=255)
                                           ngttagttaTAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCCCAGTGACATTGCTTTCGGGtttttttt                                  >  1:31254/2‑76 (MQ=255)
                                            gttagttaTAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTAttttt                                 <  1:1500483/76‑1 (MQ=255)
                                                gttaTAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTcg                             <  1:1313527/76‑1 (MQ=255)
                                                gttaTAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTTTTTTTTTcg                             >  1:1648681/1‑76 (MQ=255)
                                                    taGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTcgacga                         <  2:574209/76‑1 (MQ=255)
                                                    taGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTcgacga                         <  2:244788/76‑1 (MQ=255)
                                                    taGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCCCAGGGACATTGCTTTCGGTTTTTTTTTTTTttacga                         >  1:163407/1‑76 (MQ=255)
                                                    taGGTTACTGTAAATGTGGCAGGGGCAGAAAAATTACCCGCAGTGACATTGCTTTCGGGTTTTTTTTTTTtgacga                         >  1:1288129/1‑76 (MQ=255)
                                                       gTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCa                      <  2:1130479/76‑1 (MQ=255)
                                                       gTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCa                      <  1:980945/76‑1 (MQ=255)
                                                        ttACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAg                     <  1:515364/76‑1 (MQ=255)
                                                          aCTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTTTTTTTTTTGACGACCAGTc                   >  1:24374/1‑76 (MQ=255)
                                                           cTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCg                  <  2:1732721/76‑1 (MQ=255)
                                                             gTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAg                <  2:1232545/76‑1 (MQ=255)
                                                              tAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGc               <  1:829908/76‑1 (MQ=255)
                                                               aaaTGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCa              <  2:917192/76‑1 (MQ=255)
                                                                aaTGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCAc             <  1:1152032/76‑1 (MQ=255)
                                                                 aTGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACt            <  2:1551034/76‑1 (MQ=255)
                                                                 aTGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACt            <  2:1526129/76‑1 (MQ=255)
                                                                  tgtgGCAGGGGCAGAAAATTTACCCCCCGGGACATTTCTTTCGGGTTTTTTTTTTTTTACGACCCGTCGAGCACTc           >  2:403298/1‑76 (MQ=255)
                                                                    tgGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCCCTCAt         >  2:527737/1‑76 (MQ=255)
                                                                    tgGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCAt         >  1:578643/1‑76 (MQ=255)
                                                                     ggCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATc        >  1:181438/1‑76 (MQ=255)
                                                                     ggCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATc        >  1:181455/1‑76 (MQ=255)
                                                                      gCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCg       <  2:779232/76‑1 (MQ=255)
                                                                      gCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCg       >  1:1829511/1‑76 (MQ=255)
                                                                       cAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCgg      >  2:227638/1‑76 (MQ=255)
                                                                           tGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCGGataa  <  2:1626476/76‑1 (MQ=255)
                                                                           tGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCGGataa  <  2:1560182/76‑1 (MQ=255)
                                                                           |                                                                           
TACGTTGAGTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCGGATAA  >  REL606/154422‑154572

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