Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL606154,4970TG35.2% 97.8 / 21.2 54A191A (GCA→GCC) yadLpredicted fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  ref base T (6/29);  new base G (19/0);  total (25/29)
Fisher's exact test for biased strand distribution p-value = 9.64e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGTTGAGTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCGGATAA  >  REL606/154424‑154572
                                                                         |                                                                           
cGTTGAGTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGc                                                                           <  2:1383966/76‑1 (MQ=255)
  ttGAGTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAg                                                                         <  1:1232622/76‑1 (MQ=255)
  ttGAGTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAg                                                                         >  1:733071/1‑76 (MQ=255)
   tGAGTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGa                                                                        <  2:672761/76‑1 (MQ=255)
      gTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGaaaa                                                                     <  1:735757/76‑1 (MQ=255)
      gTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGaaaa                                                                     >  1:756226/1‑76 (MQ=255)
        gggATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAAtt                                                                   >  1:666919/1‑76 (MQ=255)
         ggATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAAttt                                                                  <  2:1139086/76‑1 (MQ=255)
            tGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTAcc                                                               <  1:194218/76‑1 (MQ=255)
            tGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTAcc                                                               <  1:345163/76‑1 (MQ=255)
             gCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTAccc                                                              <  1:454192/76‑1 (MQ=255)
               atatatAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGGAGAAAATTTACCCGc                                                            >  2:1467367/1‑76 (MQ=255)
                 atataGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAg                                                          <  2:525864/76‑1 (MQ=255)
                 atataGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAAATTATCCCCAg                                                          >  1:1261515/1‑76 (MQ=255)
                  tataGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGGAGAAAATTTACCCGccgg                                                         >  2:164004/1‑73 (MQ=255)
                   ataGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTg                                                        <  1:1463078/76‑1 (MQ=255)
                   ataGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTg                                                        <  1:1291596/76‑1 (MQ=255)
                      gTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACa                                                     <  1:958227/76‑1 (MQ=255)
                       ttCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACAt                                                    <  1:904483/76‑1 (MQ=255)
                        tCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACAtt                                                   <  2:506414/76‑1 (MQ=255)
                        tCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACAtt                                                   <  2:671148/76‑1 (MQ=255)
                         ccTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTg                                                  <  1:1089700/76‑1 (MQ=255)
                           tGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCt                                                <  2:745945/76‑1 (MQ=255)
                            ggATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCtt                                               <  1:456208/76‑1 (MQ=255)
                                 tGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGTGACATTGCTTTCGGt                                          >  1:176713/1‑76 (MQ=255)
                                 tGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGTGACATTGCTTTCGGt                                          >  2:970343/1‑76 (MQ=255)
                                 tGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAAATTACCCGCCGTGACATTGCTTACGGt                                          >  1:1154098/1‑76 (MQ=255)
                                   tCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGTGACATTGCTTACGGttt                                        >  2:346302/1‑76 (MQ=255)
                                       tCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGTGACATTGCTTACGGGtttttt                                    >  1:1323868/1‑76 (MQ=255)
                                         tGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTAtttt                                  <  2:961649/76‑1 (MQ=255)
                                           ttagttaTAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTTTCCGCAGGGAAATTGCTTACGGGtttttttttt                                >  2:474422/1‑76 (MQ=255)
                                           ttagttaTAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGTGACATTGCTTTCGGTtttttttttt                                >  1:1151966/1‑76 (MQ=255)
                                           ttagttaTAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTtttttttttt                                >  2:566967/1‑76 (MQ=255)
                                              gttaTATGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGTGACATTGCTTACGGGTtttttttttttg                             >  1:694953/1‑74 (MQ=255)
                                               ttaTAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGGGGCATTGCTTATGGGTtttttttttttgt                            >  2:517360/1‑73 (MQ=255)
                                                   aGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCgacgac                        <  2:419506/76‑1 (MQ=255)
                                                     gTTACTGTAAATGTGGCAGGTGCAGAAAAATTACCCGCTGTGAACATGCTTTCGGGTTTTTTTTTTTTGGCGAccc                      >  2:33786/1‑75 (MQ=255)
                                                      ttACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCCGTGACATTGCTTACGGTTTTATTTTTTTTGACGACCAg                     >  1:11767/1‑76 (MQ=255)
                                                       tACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGt                    >  2:369538/1‑76 (MQ=255)
                                                         cTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCg                  <  1:1375310/76‑1 (MQ=255)
                                                         cTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCg                  <  2:410504/76‑1 (MQ=255)
                                                         cTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCg                  <  2:488250/76‑1 (MQ=255)
                                                           gTAAATGTGGCAGGTGGAGAAAATTTACCCGCAGGGGGATTTCTTTCGGTTTTATTTTTTTCGTCGACCAGTCGAg                >  2:1159178/1‑76 (MQ=255)
                                                           gTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGGCATTGCTTTCGGTTTTTTTTTTTTTGACGCCCAGGCGAg                >  2:591725/1‑76 (MQ=255)
                                                           gTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAg                <  2:1212182/76‑1 (MQ=255)
                                                           gTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAg                >  1:1514934/1‑76 (MQ=255)
                                                             gaaTGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTTTTTTTTTCGACGACCAGTCGAGCa              >  2:367835/2‑76 (MQ=255)
                                                             aaaTGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCa              <  1:1522934/76‑1 (MQ=255)
                                                                tgtgGCAGGGGGAGAAAATTTACCCGCCGGGACATTGCTTTCGGGTTTTTTTTTTTTTGCGACCCGTCGAGCACTc           >  2:1035481/1‑76 (MQ=255)
                                                                 gtgGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCa          <  1:1333456/76‑1 (MQ=255)
                                                                  tgGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCAt         >  2:452399/1‑76 (MQ=255)
                                                                    gCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCg       <  1:486920/76‑1 (MQ=255)
                                                                     cAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCgg      <  2:8029/76‑1 (MQ=255)
                                                                        gTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCGGata   >  1:1494334/1‑76 (MQ=255)
                                                                        gTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCGGata   <  2:923952/76‑1 (MQ=255)
                                                                         tGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCGGataa  <  1:252461/76‑1 (MQ=255)
                                                                         |                                                                           
CGTTGAGTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCGGATAA  >  REL606/154424‑154572

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