Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6062,197,8670TG30.9% 87.4 / 15.7 42A47A (GCA→GCC) cirAferric iron‑catecholate outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base T (4/25);  new base G (13/0);  total (17/25)
Fisher's exact test for biased strand distribution p-value = 9.33e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TTTGAGGACATCCTTCAGATTCTGTACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGCTAACAC  >  REL606/2197792‑2197941
                                                                           |                                                                          
tttGAGGACAGCCTTCAGATTCTGTACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGt                                                                            <  1:643686/76‑1 (MQ=255)
 ttGAGGACATCCTTCAGATTCTGTACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTg                                                                           <  2:412485/76‑1 (MQ=255)
   gAGGACATCCTTCAGATTCTGTACCGGTTTTCGCTGCAGTTCTTCCTGGGTAATGACGCTGATACTGGCAGGGGct                                                                         >  1:1488264/1‑75 (MQ=255)
    aGGACATCCTTCAGATTCTGTACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGGGCAt                                                                        >  2:535067/1‑76 (MQ=255)
      gACATCCTTCAGATTCTGTACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCt                                                                      <  2:1224763/76‑1 (MQ=255)
        cATCCTTCAGATTCTGTACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGGGTCttt                                                                    >  1:1042739/1‑76 (MQ=255)
         aTCCTTCAGATTCTGTACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTa                                                                   <  1:675867/76‑1 (MQ=255)
           ccTTCAGATTCTGTACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGGGGATCTTTAAg                                                                 >  2:263824/1‑76 (MQ=255)
            cTTCAGATTCTGTACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGGGGATCTTTTAGa                                                                >  1:1157328/1‑76 (MQ=255)
                aGATTCTGTACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGGGCCTCTTTTAGAtttt                                                            >  2:1463624/1‑76 (MQ=255)
                  aTTCTGTACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGt                                                          <  1:2020321/76‑1 (MQ=255)
                  aTTCTGTACCGGTTTTCGCTGCAGGGCTTCCTGGGTAATGACGCTGATACTGGCAGGCGCCTCTttttgtttttgt                                                          >  1:362195/1‑76 (MQ=255)
                     cTGTACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGGGCATCTTTAAAATTTTGtttc                                                       >  1:1782818/1‑74 (MQ=255)
                      tGTACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGGGCATCTTTTAGATTTTGTTCCa                                                      >  2:1771276/1‑76 (MQ=255)
                       gTACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCAc                                                     <  1:1179953/76‑1 (MQ=255)
                       gTACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGGGGATCTTTTAGATTTTTTTCCAc                                                     >  1:313437/1‑76 (MQ=255)
                         aCCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGGGGGTCTTTTAGATTTTTTTTCACgg                                                   >  1:668659/1‑76 (MQ=255)
                         aCCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGGGGATCTTTTAGATTTTTTTTCAccg                                                   >  1:371968/1‑74 (MQ=255)
                          ccGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGGGGCTCTTTTAGATTTTTTTCCACGGa                                                  >  2:717205/1‑76 (MQ=255)
                           cGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGaa                                                 <  1:2013519/76‑1 (MQ=255)
                            ggTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAg                                                <  1:44550/76‑1 (MQ=255)
                             gTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGGGGATCTTTTTGATTTTTTTTCACGGgaga                                               >  1:794132/1‑76 (MQ=255)
                             gTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGGGCATCTTTAAAAATTTTTTTCACGCAAGa                                               >  1:1608871/1‑76 (MQ=255)
                              ttttCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGAt                                              <  2:1945153/76‑1 (MQ=255)
                              ttttCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGGGCATCTTTAAGATTTTTTTCCCCGGAAGAt                                              >  2:1412792/1‑76 (MQ=255)
                                 tCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGGGCATCTTTTAGATTTTGTTTCACGGAAGATGCa                                           >  2:543822/1‑76 (MQ=255)
                                   gCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGt                                         <  1:696031/76‑1 (MQ=255)
                                     tGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGGGGATCTTTAAGATTTTGTTTCACGGAAGATGCAGTGa                                       >  1:278399/1‑76 (MQ=255)
                                      gCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGGGGATTTTTTAGATTTTGTTTCACGGGAGAAGGAGGGAc                                      >  2:966009/1‑76 (MQ=255)
                                       cAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACa                                     <  1:1925464/76‑1 (MQ=255)
                                       cAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGGGCATCTTTTAGATTTTGTTCCACGGAAGACGCAGTgtga                                     >  1:1193002/1‑73 (MQ=255)
                                        aGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACaa                                    <  1:461443/76‑1 (MQ=255)
                                         ggTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAAc                                   <  1:1931144/76‑1 (MQ=255)
                                          gTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAAcc                                  <  1:643253/76‑1 (MQ=255)
                                          gTCTTCCTGGGTAATGACGCTGATACTGGCAGGGGCATCTTTTAGATTTTTTTTCACCGAAGATGCAGTGACaaac                                  >  1:1311095/1‑74 (MQ=255)
                                            cTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCAt                                <  2:1562854/76‑1 (MQ=255)
                                             ttCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATc                               <  1:645579/76‑1 (MQ=255)
                                              tCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGGAGACGCAGTGACAACCATCg                              >  1:1516815/1‑76 (MQ=255)
                                                cTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACACCCATCGtt                            <  2:1903304/76‑1 (MQ=255)
                                                cTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGtt                            <  1:539753/76‑1 (MQ=255)
                                                cTGGGTAATGACGCTGATACTGGCAGGGGCATCTTTTAGATTTTGTTTCACGGAAGATGCAGTGACAACCATCGtt                            >  1:1397158/1‑76 (MQ=255)
                                                 tGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGttt                           <  2:1247927/76‑1 (MQ=255)
                                                    gTAATGACGCTGATACTGGCAGGGGCATCTTTTAGAATTTTTTCCCCGGAAGATGCAGTGACAACCATCGTTTCGc                        >  2:1914399/1‑76 (MQ=255)
                                                         gACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCcatca                   <  1:443580/76‑1 (MQ=255)
                                                            gCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCcatcatca                <  1:615562/76‑1 (MQ=255)
                                                              tGATACTGGCAGGGGCATCTTTTAGATTTTTTTTCACGGAAGATGCAGTGACAACCATCGTTTTGCCctcatcatc              >  2:1389289/1‑76 (MQ=255)
                                                                gtACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGa            <  1:1423444/75‑1 (MQ=255)
                                                                 tACTGGCAGGGGCATCTTTAAGATTTTTTTTCACGGAAGATGGAGTGACAACCAACGTTTTGCCCTCATCATCGAc           >  2:1312275/1‑76 (MQ=255)
                                                                  aCTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGAcc          <  1:290532/76‑1 (MQ=255)
                                                                    tGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGc        >  1:963144/1‑76 (MQ=255)
                                                                      gCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACATCCATCGTTTCGCCATCATCATCGACCGCTa      <  2:2027809/76‑1 (MQ=255)
                                                                      gCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGCTa      <  2:876378/76‑1 (MQ=255)
                                                                       cAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGCTaa     <  1:758431/76‑1 (MQ=255)
                                                                       cAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGCTaa     >  2:1010937/1‑76 (MQ=255)
                                                                        aGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGCTAac    >  1:960806/1‑76 (MQ=255)
                                                                          gTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGCTAacac  <  1:538540/76‑1 (MQ=255)
                                                                           |                                                                          
TTTGAGGACATCCTTCAGATTCTGTACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGCTAACAC  >  REL606/2197792‑2197941

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