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

ACGTTGAGTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCGGATA  >  REL606/154423‑154571
                                                                          |                                                                          
aCGTTGAGTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTg                                                                           <  1:1207637/76‑1 (MQ=255)
 cGTTGAGTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGc                                                                          <  2:450652/76‑1 (MQ=255)
   ttGAGTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAg                                                                        <  1:1431776/76‑1 (MQ=255)
     gAGTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGaa                                                                      <  1:1081708/76‑1 (MQ=255)
      aGTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGaaa                                                                     <  1:1119503/76‑1 (MQ=255)
       gTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGaaaa                                                                    <  1:132024/76‑1 (MQ=255)
         gggATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAAtt                                                                  <  2:1136104/76‑1 (MQ=255)
          ggATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAAttt                                                                 >  1:1291508/1‑76 (MQ=255)
            aTGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTAc                                                               <  1:718075/76‑1 (MQ=255)
                atatatAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTATCCGc                                                           >  2:976468/1‑76 (MQ=255)
                  atataGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAg                                                         <  2:1200982/76‑1 (MQ=255)
                    ataGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTg                                                       <  1:1163401/76‑1 (MQ=255)
                          ccTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTg                                                 <  2:725452/76‑1 (MQ=255)
                          ccTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTg                                                 <  2:738638/76‑1 (MQ=255)
                          ccTCGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTg                                                 <  2:1216702/76‑1 (MQ=255)
                           cTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCCGGGACATTGc                                                >  2:1663450/1‑76 (MQ=255)
                             ggATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGTGACATTGCtt                                              >  1:186712/1‑76 (MQ=255)
                              gATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTa                                             <  2:809639/76‑1 (MQ=255)
                               aTCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGTGACATTGCtttc                                            >  2:1483334/1‑74 (MQ=255)
                                tCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACg                                           <  2:436679/76‑1 (MQ=255)
                                tCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGCGACATTGCTTTCg                                           >  2:1636975/1‑76 (MQ=255)
                                     cAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGCGACATTGCCTACCGtttt                                      >  1:1628897/1‑76 (MQ=255)
                                        tCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTtttttt                                   >  1:1399540/1‑76 (MQ=255)
                                         cTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGTAGGAAAATTACCCGCAGAAACATTTCTTTCGGTTTTAttt                                  >  2:533817/1‑76 (MQ=255)
                                          tGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTAtttt                                 <  1:181212/76‑1 (MQ=255)
                                             tagttaTAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCCGTGACATTGCTTACGGTTTTAttttttt                              >  2:77823/1‑76 (MQ=255)
                                              agttaTAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTc                             <  2:308514/76‑1 (MQ=255)
                                                ttaTAGGTTACTGTAAATGTGGCAGGGGCCGAAAATTTACCCGCAGTGACATTTCTTTCGGTTTTATTTTTTTcta                           >  2:1590199/1‑74 (MQ=255)
                                                ttaTAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTTACCCCAGGGACATTGCTTTCGGGTTTTTTTTTTTtga                           >  2:1105615/1‑76 (MQ=255)
                                                ttaTAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCCGGGAGATTGCTTTCGGGTtttttttttttta                           >  1:942523/1‑73 (MQ=255)
                                                 tataGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTcgac                          <  1:251423/76‑1 (MQ=255)
                                                 tataGGTTACTGTAAATGTGGCAGGGGGAGGAAATTTACCCGCAGTGACATTGCTTTCGGTTTTTTTTTTTTcgac                          >  1:1186636/1‑76 (MQ=255)
                                                 tataGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGTGACATTGCTTTCGGTTTTATTTTTTTtgac                          >  1:1048124/1‑76 (MQ=255)
                                                  ataGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTcgacg                         <  2:954655/76‑1 (MQ=255)
                                                   taGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTAACCGCAGTTGCATTGCTTACGGGTTTTTTTTTTTtgacgc                        >  2:396214/1‑75 (MQ=255)
                                                   taGGTTACTGTAAATGTGGCAGGTGCAGAAAAATTACCCCCAGGGAAATTGCGTTCGGGTTTATTTTTTTtgtcga                        >  2:824573/1‑76 (MQ=255)
                                                   taGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGTGACATTGCTTTCGGTTTTTTTTTTTTtgacga                        >  2:683204/1‑76 (MQ=255)
                                                    aGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCgacgac                       <  2:869592/76‑1 (MQ=255)
                                                    aGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCgacgac                       <  1:1632653/76‑1 (MQ=255)
                                                    aGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCgacgac                       <  2:387738/76‑1 (MQ=255)
                                                       ttACTGTAAATGTGGCAGGGGCAGAAAATTTAACCGCAGTGACATTGCTTACGGTTTTTTTTTTTTCGACGACCAg                    >  2:246574/1‑76 (MQ=255)
                                                        tACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGt                   <  2:1298317/76‑1 (MQ=255)
                                                        tACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTTTTTTTTTTGACGACCAGt                   >  1:1267445/1‑76 (MQ=255)
                                                         aCTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTc                  >  1:1447904/1‑76 (MQ=255)
                                                         aCTGTAAATGTGGCAGGGGCAGAAAATTTACCCCCAGTGACATTGCTTTCGGGTTTTTTTTTTTCGACGACCCGTc                  >  2:900588/1‑76 (MQ=255)
                                                          cTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCg                 <  2:247776/76‑1 (MQ=255)
                                                          cTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCg                 <  2:664733/76‑1 (MQ=255)
                                                           tGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGa                <  1:1611288/76‑1 (MQ=255)
                                                           tGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGa                <  1:854799/76‑1 (MQ=255)
                                                             tAAATGTGGCAGGTGGAGAAAATTTACCCGCAGTGACATTTCTTACGGTTTTATTTTTTTTTCCGACCCGTCGAGc              >  2:1326575/1‑76 (MQ=255)
                                                             tAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGc              <  2:1504392/76‑1 (MQ=255)
                                                             tAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGGTTTATTTTTTTTGACGACCAGTCGAGc              >  2:1447546/1‑76 (MQ=255)
                                                              aaaTGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCa             <  2:1558362/76‑1 (MQ=255)
                                                                aTGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACt           >  2:244077/1‑76 (MQ=255)
                                                                  gtgGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCa         <  2:1144300/76‑1 (MQ=255)
                                                                   tgGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCAt        <  1:1234628/76‑1 (MQ=255)
                                                                     gCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCg      <  1:202760/76‑1 (MQ=255)
                                                                     gCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCg      >  2:230326/1‑76 (MQ=255)
                                                                       aGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCGGa    >  1:399027/1‑76 (MQ=255)
                                                                       aGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCGGa    >  2:595148/1‑76 (MQ=255)
                                                                         gTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCGGata  <  1:474992/76‑1 (MQ=255)
                                                                          |                                                                          
ACGTTGAGTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCGGATA  >  REL606/154423‑154571

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