Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6062,197,8670TG42.2% 54.4 / 24.4 45A47A (GCA→GCC) cirAferric iron‑catecholate outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base T (3/23);  new base G (19/0);  total (22/23)
Fisher's exact test for biased strand distribution p-value = 6.32e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TTGAGGACATCCTTCAGATTCTGTACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGCTAACACA  >  REL606/2197793‑2197942
                                                                          |                                                                           
ttGAGGACATCCTTCAGATTCTGTACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTg                                                                            <  1:619926/76‑1 (MQ=255)
 tGAGGACATCCTTCAGATTCTGTACCGGTTTTCGCTGCACGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGc                                                                           <  2:580083/76‑1 (MQ=255)
  gAGGACATCCTTCAGATTCTGTACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCa                                                                          <  2:1138661/76‑1 (MQ=255)
      aCATCCTTCAGATTCTGTACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCtt                                                                      <  2:240432/76‑1 (MQ=255)
       cATCCTTCAGATTCTGTACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCttt                                                                     <  1:435781/76‑1 (MQ=255)
       cATCCTTCAGATTCTGTACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGGGGATCttt                                                                     >  1:302196/1‑76 (MQ=255)
        aTCCTTCAGATTCTGTACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTa                                                                    <  2:1887509/76‑1 (MQ=255)
        aTCCTTCAGATTCTGTACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTa                                                                    <  1:840462/76‑1 (MQ=255)
         tCCTTCAGATTCTGTACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTaa                                                                   <  1:1227060/76‑1 (MQ=255)
         tCCTTCAGATTCTGTACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGGGGATCtttta                                                                   >  1:1010991/1‑74 (MQ=255)
             tCAGATTCTGTACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGGGGATCTTTTAGAtt                                                               >  1:1023845/1‑76 (MQ=255)
              cAGATTCTGTACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAAGAttt                                                              <  2:1343902/76‑1 (MQ=255)
               tgcttCTGTACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAAGAtttt                                                             <  2:987967/73‑1 (MQ=255)
                 aTTCTGTACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGt                                                           <  2:206886/76‑1 (MQ=255)
                    cTGTACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTcc                                                        <  1:1734738/76‑1 (MQ=255)
                    cTGTACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTcc                                                        <  2:1768019/76‑1 (MQ=255)
                    cTGTACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGGGGGTCTTTTAGATtttttttc                                                        >  1:1175013/1‑74 (MQ=255)
                    cTGTACCGGTTTTCGCTGCAGGTCTTCCTGAGTAATGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTcc                                                        <  1:702798/76‑1 (MQ=255)
                     tGTACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCa                                                       <  2:628232/76‑1 (MQ=255)
                      gTACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGGGCCTCTTTTAGATTTTTTTcccc                                                      >  1:42919/1‑74 (MQ=255)
                        aCCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACgg                                                    <  2:1066014/76‑1 (MQ=255)
                        aCCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGGGGCTCTTTTAAATTTTTTTTCAccg                                                    >  1:807795/1‑74 (MQ=255)
                          cGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGaa                                                  <  1:1606878/76‑1 (MQ=255)
                          cGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGaa                                                  <  1:1606877/76‑1 (MQ=255)
                          cGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGGGCGTCTTTTAGATTTTGTTTCACggga                                                  >  1:263424/1‑74 (MQ=255)
                           ggTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGGGGATCTTTAAGATTTTGTTTCACgggaa                                                 >  1:147862/1‑73 (MQ=255)
                               ttCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGGGCATCTTTTAGATTTTTTTTCACGGGAGATgg                                             >  2:898809/1‑75 (MQ=255)
                               ttCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGGGCATCTTTTAGATTTTGTTTCACCGAAGaagg                                             >  1:1239946/1‑73 (MQ=255)
                                    tGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGa                                        <  2:725253/76‑1 (MQ=255)
                                     gCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGGGCATCTTTTAGATTTTTTTTCACGGGAGATGCAGTggc                                       >  2:1267242/1‑74 (MQ=255)
                                     gCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGGGCATCTTTTAGATTTTGTTCCACGGAAGATGCAGTGAc                                       >  2:556367/1‑76 (MQ=255)
                                     gCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGGGCATCTTTAAGATTTTTTTCCACGGGAGATGCCGTGaa                                       >  2:1737359/1‑75 (MQ=255)
                                      cAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACa                                      <  1:772582/76‑1 (MQ=255)
                                         gTCTTCCTGGGTAATGACGCTGATACTGGCAGGGGCATCTTTTAGATTTTTTTTCCCGGAAGATGCAGGGAAAAcc                                   >  2:1041911/1‑76 (MQ=255)
                                             tCCTGGGTAATGACGCTGATACTGGCAGGGGGGGCTTTTAGATTTTTTTTCCCGGGAAAATCAGGGACAACCAACg                               >  1:1114462/1‑76 (MQ=255)
                                               cTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTTAGATTTTGTTTCACGGAAGATGCAGGGACAACCATCggt                             >  2:1502495/1‑74 (MQ=255)
                                               cTGGGTAATGACGCTGATACTGGCAGGGGGCTCTTTTAAATTTTTTTTCACGGAAGAAGCAGTGACAACCCACGtt                             >  1:608958/1‑76 (MQ=255)
                                               cTGGGTAATGACGCTGATACTGGCAGGGGGATCTTTTAGATTTTTTTTCACCGGAGATGCAGGGGCAACCATCGtt                             >  1:1079881/1‑76 (MQ=255)
                                                tGGGTAATGACGCTGATACTGGCAGGGGCCTCTTTTAGATTTTTTTTCACGGGAGAGGGAGGGAGAAACACCGttt                            >  1:355182/1‑76 (MQ=255)
                                                tGGGTAATGACGCTGATACTGGCAGGGGCATCTTTTAGATTTTTTTTCCCGGAAGATGCAGTGACAACCATCGttt                            >  1:263818/1‑76 (MQ=255)
                                                tGGGTAATGACGCTGATACTGGCAGGGGCATCTTTTAGATTTTGTTCCCCGGAAGATGCAGTGACAAACATCGttt                            >  1:474711/1‑76 (MQ=255)
                                                 gggTAATGACGCTGATACTGGCAGGGGCTTCTTTAAGATTTTGTTTCACGGAAGATGGCGGGACAACCATCGTTTc                           >  1:82464/1‑76 (MQ=255)
                                                 gggTAATGACGCTGATACTGGCAGGGGCATCTTTTAGATTTTTTTCCACGGAAGATGCAGTGACAACCACCGTTTc                           >  2:1808797/1‑76 (MQ=255)
                                                 gggTAATGACGCTGATACTGGCAGGGGCATCTTTTAGATTTTTTTCCACGGAAGATGCAGTGACAACCAACGTTTc                           >  1:1579608/1‑76 (MQ=255)
                                                      aTGACGCTGATACTGGCAGGGGCATCTTTTAGATTTTTTTCCCCGGAAGATGCAGTGAAAAACAACGTTTCGCcat                      >  2:168587/1‑76 (MQ=255)
                                                      aTGACGCTGATACTGGCAGGGGCATCTTTTAGATTTTGTTTCACGGAAGATGGAGGGACAACCAACGTTTTGccct                      >  2:1244265/1‑74 (MQ=255)
                                                          cGCTGATACTGGCAGGGGGATCTTTTAGATTTTTTTTCACGGGAGATGCAGGGACAAACATCGTTTTGCcctcatc                  >  1:1851295/1‑76 (MQ=255)
                                                           gCTGATACTGGCAGGGGCATCTTTTAGATTTTGTTCCACGGGAGATGGAGGGAAAAAACTCCTTTCGCcatcatta                 >  2:1444289/1‑74 (MQ=255)
                                                                tACTGGCAGGGGCATCTTTTAGATTTTTTTTCACGGGAGATGCAGTGACAACCATCGTTTTGCCCTCATCCTCCAc            >  2:1784566/1‑76 (MQ=255)
                                                                  cTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCg          >  1:95172/1‑76 (MQ=255)
                                                                     gCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGCTa       <  1:1820840/76‑1 (MQ=255)
                                                                      cAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGCTaa      <  2:1057977/76‑1 (MQ=255)
                                                                        ggTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGCTAaca    <  1:142214/76‑1 (MQ=255)
                                                                         gTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGCTAacac   >  1:1693787/1‑76 (MQ=255)
                                                                          tGCATCTTTAAGATTTTGTTCNACGGAAGATGCAGTGANNACCATCGTTTCGCCATCATCATCGACCGCTAAcaca  <  1:768737/76‑1 (MQ=255)
                                                                          |                                                                           
TTGAGGACATCCTTCAGATTCTGTACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGCTAACACA  >  REL606/2197793‑2197942

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