Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL606154,4970TG23.1% 131.4 / 10.6 52A191A (GCA→GCCyadLpredicted fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  ref base T (4/36);  new base G (12/0);  total (16/36)
Fisher's exact test for biased strand distribution p-value = 8.82e-09
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                                                                           <  1:450515/76‑1 (MQ=255)
cGTTGAGTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGc                                                                           <  1:1059398/76‑1 (MQ=255)
     aGTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGaaa                                                                      <  1:1060389/76‑1 (MQ=255)
      gTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGaaaa                                                                     <  2:1243261/76‑1 (MQ=255)
       tGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAAt                                                                    >  1:1706943/1‑76 (MQ=255)
          gATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGGAGAAAAtttt                                                                 >  2:938532/1‑75 (MQ=255)
          gATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGGAGAAAAtata                                                                 >  2:62832/1‑76 (MQ=255)
            tGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTAcc                                                               <  1:1141344/76‑1 (MQ=255)
               atatatAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGGAGAAAAATTAccccc                                                            >  2:1267647/1‑74 (MQ=255)
                   ataGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTg                                                        <  1:1421098/76‑1 (MQ=255)
                    taGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGa                                                       <  2:1737618/76‑1 (MQ=255)
                    taGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGa                                                       <  2:393950/76‑1 (MQ=255)
                    taGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCAGCAGTGa                                                       <  2:576917/76‑1 (MQ=255)
                       ttCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTAACCGCAGTGACAt                                                    >  1:1879161/1‑76 (MQ=255)
                        tCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACAtt                                                   <  2:1128272/76‑1 (MQ=255)
                        tCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACAtt                                                   <  1:417751/76‑1 (MQ=255)
                        tCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACAtt                                                   <  2:640787/76‑1 (MQ=255)
                         ccTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTg                                                  <  1:381696/76‑1 (MQ=255)
                          cTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGTGACATTGc                                                 >  1:1728284/1‑76 (MQ=255)
                           tGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTTCCCGCAGTGACATTGCt                                                >  2:1578575/1‑76 (MQ=255)
                           tGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAAATTAACCGCAGGGACATTGCt                                                >  1:1471363/1‑76 (MQ=255)
                              aTCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTAc                                             <  1:152447/76‑1 (MQ=255)
                              aTCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGGGACATTGCtttc                                             >  2:1220217/1‑74 (MQ=255)
                                cTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACgg                                           <  1:773089/76‑1 (MQ=255)
                                cTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAAATTACCCGCCGGGACATTGCTTTCgg                                           >  1:988142/1‑76 (MQ=255)
                                  gTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGtt                                         <  2:525943/76‑1 (MQ=255)
                                    cAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTTCGGtttt                                       <  2:829472/76‑1 (MQ=255)
                                    cAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGtttt                                       <  1:1229122/76‑1 (MQ=255)
                                    cAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGtttt                                       <  2:1336529/76‑1 (MQ=255)
                                    cAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGtttt                                       <  2:1230110/76‑1 (MQ=255)
                                     aGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTa                                      <  2:295081/76‑1 (MQ=255)
                                       tCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAAATTACCCGCCGGGGCATTTCTTACGGTtttttt                                    >  1:1823571/1‑76 (MQ=255)
                                            tctttaTAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTAttttttt                               <  1:928633/73‑1 (MQ=255)
                                            tagttaTAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTAttttttt                               <  2:462162/76‑1 (MQ=255)
                                              gttaTAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGTGACATTGCTTTCGGTTTTTTTTTTTTct                             >  2:218841/1‑75 (MQ=255)
                                              gttaTAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGTGAAATTGCTTTCGGGTtttttttttttt                             >  1:571543/1‑74 (MQ=255)
                                               ttaTAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTcga                            <  2:1850021/76‑1 (MQ=255)
                                               ttaTAGGTTACTGTAAATGTGGCAGGGGGAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTTTTTTTTTaga                            >  2:928633/1‑76 (MQ=255)
                                                   aGGTTACTGTAAATGTGGCAGGGGGAGAAAATTTAACCGCAGAGAAATTGCTTACGGTTTTTTTTTTTTCgaccac                        >  1:281243/1‑76 (MQ=255)
                                                       tACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGTGAACTTGCTTACGGTTTTATTTTTTTTGACGACCAGt                    >  1:387412/1‑76 (MQ=255)
                                                         cTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGTCGACCAGTCg                  <  2:962656/76‑1 (MQ=255)
                                                         cTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCg                  <  2:567209/76‑1 (MQ=255)
                                                           gTAAATGTGGCAGGGGCAGCAAATTTACCCGGAGTGACACTGCTTAAGGGTTTATTTTTTTTAACGGCCAGACGAg                >  2:1345040/1‑76 (MQ=255)
                                                            tAAATGTGGCAGGGGCAGAAAATTTACCCGCAGTGACATTGCTTACGGGTTTTTTTTTTTTGACGACCAGTCGAGc               >  2:708108/1‑76 (MQ=255)
                                                              aaTGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCAc             <  2:139296/76‑1 (MQ=255)
                                                                tgtgGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTc           <  2:924545/76‑1 (MQ=255)
                                                                tgtgGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTc           <  1:785653/76‑1 (MQ=255)
                                                                tgtgGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTc           <  2:1711567/76‑1 (MQ=255)
                                                                  tgGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCAt         <  1:1405721/76‑1 (MQ=255)
                                                                  gggCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCAt         <  1:245705/75‑1 (MQ=255)
                                                                   ggCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATc        >  2:1849291/1‑76 (MQ=255)
                                                                     cAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCgg      >  1:328935/1‑76 (MQ=255)
                                                                     cAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCgg      >  2:1083343/1‑76 (MQ=255)
                                                                       ggTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCGGat    <  1:1621729/76‑1 (MQ=255)
                                                                       ggTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCGGat    <  1:1707034/76‑1 (MQ=255)
                                                                       ggTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCGGat    <  1:132200/76‑1 (MQ=255)
                                                                        gTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCGGata   <  1:1404004/76‑1 (MQ=255)
                                                                        gTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCGGata   >  1:324353/1‑76 (MQ=255)
                                                                         tGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCGGataa  >  1:825032/1‑76 (MQ=255)
                                                                         |                                                                           
CGTTGAGTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCGGATAA  >  REL606/154424‑154572

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