Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL606154,4970TG31.8% 134.9 / 20.4 66A191A (GCA→GCC) yadLpredicted fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  ref base T (8/37);  new base G (21/0);  total (29/37)
Fisher's exact test for biased strand distribution p-value = 4.82e-11
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:730973/76‑1 (MQ=255)
       tGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAAt                                                                    >  1:134223/1‑76 (MQ=255)
       tGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGNAAATGTGGCAGGTGCAGAAAAt                                                                    <  2:233254/76‑1 (MQ=255)
         ggATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAAttt                                                                  >  1:565314/1‑76 (MQ=255)
          gATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTa                                                                 <  2:1651474/76‑1 (MQ=255)
              cATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTAcccc                                                             >  1:1359558/1‑75 (MQ=255)
              cATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCg                                                             >  1:1594277/1‑76 (MQ=255)
               atatatAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGc                                                            <  2:1632035/76‑1 (MQ=255)
                 atataGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGGAGAAAAATTACCCCCAg                                                          >  2:858164/1‑76 (MQ=255)
                  tataGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGt                                                         <  1:34725/76‑1 (MQ=255)
                   ataGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGNTGCAGAAAATTTACCCGCAGTg                                                        <  2:230752/76‑1 (MQ=255)
                   ataGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTg                                                        <  2:1635738/76‑1 (MQ=255)
                    taGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGa                                                       <  1:492097/76‑1 (MQ=255)
                    taGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGNGAGAAAATTTTCCCGCAGTGa                                                       >  2:273563/1‑76 (MQ=255)
                     aGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGTGAc                                                      >  1:1031808/1‑76 (MQ=255)
                       ttCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCATGTGCAGAAAATTTACCCGCAGTGACAt                                                    >  1:807022/1‑76 (MQ=255)
                        tCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACAtt                                                   <  2:1545679/76‑1 (MQ=255)
                        tCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACAtt                                                   <  1:1305820/76‑1 (MQ=255)
                         ccTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTg                                                  >  2:644159/1‑76 (MQ=255)
                         ccTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGTGACATTg                                                  >  1:947123/1‑76 (MQ=255)
                          cTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGc                                                 <  1:1017049/76‑1 (MQ=255)
                          cTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGc                                                 <  1:86173/76‑1 (MQ=255)
                          cTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCCCAGTGACATTGc                                                 >  2:386505/1‑76 (MQ=255)
                            ggATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCtt                                               <  1:626559/76‑1 (MQ=255)
                            ggATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCtt                                               <  2:803051/76‑1 (MQ=255)
                            ggATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCtt                                               <  2:201963/76‑1 (MQ=255)
                              aTCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTAc                                             <  2:1022179/76‑1 (MQ=255)
                              aTCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGTGACATTGCTTAc                                             >  2:1645865/1‑76 (MQ=255)
                                 tGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGt                                          <  1:38839/76‑1 (MQ=255)
                                  gTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCCGTGACATTGCTTTCgggt                                         >  2:654240/1‑74 (MQ=255)
                                  gTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGTGACATTGCTTACgggt                                         >  2:1486888/1‑74 (MQ=255)
                                   tCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGTGACATTGCTTACGGttt                                        >  2:1283490/1‑76 (MQ=255)
                                   tCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGTGACATTGCTTACGGgtt                                        >  1:1212441/1‑76 (MQ=255)
                                   tCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGTGACATTGCTTACGGgtt                                        >  1:1249823/1‑76 (MQ=255)
                                    cAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGtttt                                       <  1:1239212/76‑1 (MQ=255)
                                    cAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCCGTGACATTGCTTTCGGtttt                                       >  1:744179/1‑76 (MQ=255)
                                       tCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTAtt                                    <  2:1005344/76‑1 (MQ=255)
                                       tCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTtttttt                                    >  1:66777/1‑76 (MQ=255)
                                         tGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTtttttttt                                  >  1:837397/1‑76 (MQ=255)
                                           ttagttaTAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCCGGGACATTGCTTTCGGTtttttttttt                                >  1:1111599/1‑76 (MQ=255)
                                           ttagttaTAGGTTACTGTAAATGTGGCAGGGGCAGAAAAATTACCCGCAGGGGCAATGGTTTCGGTtttttttttt                                >  2:326628/1‑76 (MQ=255)
                                            tagttaTAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACNTTGCTTACGGTTTTAttttttt                               <  2:229362/76‑1 (MQ=255)
                                            tagttaTAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTAttttttt                               <  2:1295807/76‑1 (MQ=255)
                                            tagttaTAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGTGACATTGCTTTCGGtttttttttttt                               >  1:588798/1‑76 (MQ=255)
                                            tagttaTAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGTGACATTGCTTACGGtttttttttttt                               >  2:1562105/1‑76 (MQ=255)
                                             agttaTAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTc                              <  1:1428996/76‑1 (MQ=255)
                                             agttaTAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGTGACATTGCTTACGGGTttttttttttt                              >  1:1118700/1‑75 (MQ=255)
                                             agttaTAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCCCAGGGGCATTGCTTTCGGGTttttttttttt                              >  1:979107/1‑75 (MQ=255)
                                              gttaTAGGTTACTGTAAATGTGGCAGGTGGAGAAAATTTACCCGCAGGGGCATTGGTTTCGGGTtttttttttttc                             >  2:350255/1‑74 (MQ=255)
                                              gttaTAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTcg                             <  1:717738/76‑1 (MQ=255)
                                              gttaTAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTcg                             <  1:1240598/76‑1 (MQ=255)
                                               ttaTAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTcga                            <  1:1693350/76‑1 (MQ=255)
                                                 ataGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTcgacg                          <  2:1194472/76‑1 (MQ=255)
                                                 ataGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTcgacg                          <  2:1408885/76‑1 (MQ=255)
                                                 ataGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTcgacg                          <  2:1354337/76‑1 (MQ=255)
                                                 ataGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTcgacg                          <  2:1202273/76‑1 (MQ=255)
                                                 ataGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTcgacg                          <  1:401865/76‑1 (MQ=255)
                                                 ataGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTcgaag                          <  1:654000/76‑3 (MQ=255)
                                                   aGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCgacgac                        <  2:1302702/76‑1 (MQ=255)
                                                     gTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCCCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCa                      >  2:1498721/1‑76 (MQ=255)
                                                          tGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGa                 >  2:1386057/1‑76 (MQ=255)
                                                               aTGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACt            >  2:1593268/1‑76 (MQ=255)
                                                                 gtgGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCa          <  1:214802/76‑1 (MQ=255)
                                                                  tgGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCAt         >  1:700927/1‑76 (MQ=255)
                                                                   ggCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATc        >  2:336275/1‑76 (MQ=255)
                                                                    gCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCg       <  2:611589/76‑1 (MQ=255)
                                                                    gCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGCTTTATTTTTTTCGACGACCAGTCGAGCACTCATCg       <  1:1316089/76‑1 (MQ=255)
                                                                      aGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCGGa     <  2:1333070/76‑1 (MQ=255)
                                                                      aGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCGGa     >  2:1133553/1‑76 (MQ=255)
                                                                        gTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCGGata   >  2:1498094/1‑76 (MQ=255)
                                                                        gTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCGGata   <  1:1577466/76‑1 (MQ=255)
                                                                         tGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCGGataa  <  2:1203984/76‑1 (MQ=255)
                                                                         |                                                                           
CGTTGAGTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCGGATAA  >  REL606/154424‑154572

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