Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6061,176,0100TG27.1% 138.1 / 11.4 59A37A (GCA→GCC) fhuEferric‑rhodotorulic acid outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base T (1/42);  new base G (16/0);  total (17/42)
Fisher's exact test for biased strand distribution p-value = 1.55e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TACGCTGTAATCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGATGGCTTGA  >  REL606/1175935‑1176085
                                                                           |                                                                           
tACGCTGTAATCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGt                                                                             <  1:1373721/76‑1 (MQ=255)
  cGCTGTAATCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTgc                                                                           <  2:1657277/76‑1 (MQ=255)
  cGCTGTAATCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTgc                                                                           <  2:790105/76‑1 (MQ=255)
  cGCTGTAATCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTgc                                                                           <  2:1369971/76‑1 (MQ=255)
     tGTAATCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTgcagc                                                                        <  1:1060223/76‑1 (MQ=255)
      gTAATCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGgcggca                                                                       >  2:788250/1‑76 (MQ=255)
      gTAATCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGgcagca                                                                       >  2:1473324/1‑76 (MQ=255)
        aaTCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGGAGCaaa                                                                     >  1:781811/1‑76 (MQ=255)
        aaTCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCaaa                                                                     >  1:47175/1‑76 (MQ=255)
         aTCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGGAGCaaaa                                                                    >  2:1113417/1‑76 (MQ=255)
             ttttCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGGGGCAAAAGCgg                                                                >  2:504467/1‑76 (MQ=255)
             ttttCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCGGCCAAAACgg                                                                >  2:504464/1‑76 (MQ=255)
             ttttCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCCAAAACgg                                                                >  2:1030262/1‑76 (MQ=255)
              tttCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGc                                                               <  1:806206/76‑1 (MQ=255)
              tttCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGc                                                               <  2:1091650/76‑1 (MQ=255)
              tttCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGc                                                               <  1:383128/76‑1 (MQ=255)
              tttCGCCATCATCTGGAGCTGNGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGc                                                               <  2:294317/76‑1 (MQ=255)
                tCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAg                                                             <  1:1263217/76‑1 (MQ=255)
                  gCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGGAGCAAAAGCGGCAGaa                                                           >  1:1643878/1‑76 (MQ=255)
                     atcatcTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGt                                                        <  2:61850/76‑1 (MQ=255)
                        atcTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAAt                                                     <  2:278763/76‑1 (MQ=255)
                         tcTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATa                                                    <  2:837189/76‑1 (MQ=255)
                            ggAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGGAGCAAAAGCGGCAGAAGGGAATAGTg                                                 >  1:1810740/1‑76 (MQ=255)
                            ggAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGAAGGGAATAGTg                                                 >  1:485746/1‑76 (MQ=255)
                            ggAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAACGGCAGAAGGGAATAGTg                                                 >  2:1811738/1‑76 (MQ=255)
                              aGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGcc                                               <  1:1300115/76‑1 (MQ=255)
                                cTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAg                                             <  2:1711599/76‑1 (MQ=255)
                                cTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAg                                             <  1:115433/76‑1 (MQ=255)
                                 tgtgGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGt                                            <  1:920804/76‑1 (MQ=255)
                                  gtgGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGGAGCAAAAACGGGAGAAGGGAATAGGGGCAggg                                           >  2:1857617/1‑74 (MQ=255)
                                  gtgGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGAAGGGAATAGTGCCAGTg                                           >  1:131757/1‑76 (MQ=255)
                                    ggCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGAAGGGAATAGTGCCAGGGCt                                         >  1:383041/1‑76 (MQ=255)
                                     gCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTa                                        <  2:1015324/76‑1 (MQ=255)
                                     gCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAACGGGAGGAGGGAAAAGGGCGAGGGCtt                                        >  1:353794/1‑75 (MQ=255)
                                      cTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTAt                                       <  1:133161/76‑1 (MQ=255)
                                        gAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGAAGGGAATAGTGCCAGTGCTATGc                                     >  2:467740/1‑76 (MQ=255)
                                        gAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGAAGGGAATAGTGCCAGTGCTATGc                                     >  2:1860922/1‑76 (MQ=255)
                                        gAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGAAGGGAATAGTGCCAGTGCTATGc                                     >  2:1587158/1‑76 (MQ=255)
                                           cccTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAAc                                  <  1:146505/76‑1 (MQ=255)
                                            ccTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAAcc                                 <  1:1009828/76‑1 (MQ=255)
                                            ccTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAAcc                                 <  2:368442/76‑1 (MQ=255)
                                            ccTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAAcc                                 <  1:1819793/76‑1 (MQ=255)
                                             cTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCg                                <  2:1006919/76‑1 (MQ=255)
                                              tCAACAATCACCGTTTCTTCAGTGGCTGGTGGCGCAAAAGCGGCAGAAGGGAATAGTGCCAGTGCTATGCAAACgg                               >  2:15887/1‑76 (MQ=255)
                                                aacaaTCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCa                             <  2:782545/76‑1 (MQ=255)
                                                aacaaTCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCa                             <  2:1737670/76‑1 (MQ=255)
                                                 acaaTCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCaa                            <  1:1790652/76‑1 (MQ=255)
                                                   aaTCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGGAATAGTGCCAGGGGTATGCAACCGGCCagg                          >  2:1037060/1‑75 (MQ=255)
                                                       aCCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTg                      <  2:13680/76‑1 (MQ=255)
                                                         cGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGAc                    <  1:990971/76‑1 (MQ=255)
                                                         cGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGAc                    <  1:833007/76‑1 (MQ=255)
                                                              cttcAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGttt               <  1:708780/76‑1 (MQ=255)
                                                              cttcAGTGGCTGGGGCAGCAAAAGCGGCAGAAGGGAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGgtt               >  2:1723082/1‑76 (MQ=255)
                                                              cttcAGTGGCTGGGGCAGCAAAAGCGGCAGAAGGGAATAGTGCCAGCGGTATGCAACCGGCAAGGAGTGAAGGgtt               >  2:335139/1‑76 (MQ=255)
                                                                 cAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGt            <  1:460684/76‑1 (MQ=255)
                                                                  aGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTg           <  2:1241595/76‑1 (MQ=255)
                                                                  aGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTg           <  1:223871/76‑1 (MQ=255)
                                                                  aGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTg           <  2:271259/76‑1 (MQ=255)
                                                                  aGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTg           <  1:441724/76‑1 (MQ=255)
                                                                   gTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGa          <  1:931030/76‑1 (MQ=255)
                                                                     ggCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGATg        <  2:1671064/76‑1 (MQ=255)
                                                                      gCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGATgg       <  2:1405439/76‑1 (MQ=255)
                                                                      gCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGATgg       <  1:530064/76‑1 (MQ=255)
                                                                       cTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGATGGc      <  1:178773/76‑1 (MQ=255)
                                                                        tGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGATGGCt     <  1:1256491/76‑1 (MQ=255)
                                                                         ggTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGATGGCtt    <  1:1051488/76‑1 (MQ=255)
                                                                         ggTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGATGGCtt    <  1:829262/76‑1 (MQ=255)
                                                                           tGCAGCAAAAGCGGCAGAAGGCAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGATGGCTTGa  >  2:1735926/1‑76 (MQ=255)
                                                                           |                                                                           
TACGCTGTAATCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGATGGCTTGA  >  REL606/1175935‑1176085

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: ‑