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

ACGCTGTAATCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGATGGCTTGA  >  REL606/1175936‑1176085
                                                                          |                                                                           
aCGCTGTAATCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTg                                                                            <  2:895623/76‑1 (MQ=255)
 cGCTGTAATCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTgc                                                                           <  2:1791190/76‑1 (MQ=255)
 cGCTGTAATCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTgc                                                                           <  2:1749483/76‑1 (MQ=255)
  gCTGTAATCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTgca                                                                          <  1:398286/76‑1 (MQ=255)
    tGTAATCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTgcagc                                                                        <  1:347929/76‑1 (MQ=255)
    tGTAATCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTgcagc                                                                        <  2:1787149/76‑1 (MQ=255)
     gTAATCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGggagca                                                                       >  2:1754089/1‑76 (MQ=255)
     gTAATCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGgcagca                                                                       >  1:548442/1‑76 (MQ=255)
       aaTCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCaaa                                                                     <  1:8367/76‑1 (MQ=255)
       aaTCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCaaa                                                                     >  1:1019692/1‑76 (MQ=255)
           aTTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGGAGCAAAAACg                                                                 >  2:479798/1‑76 (MQ=255)
                cGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGa                                                            <  1:1763864/76‑1 (MQ=255)
                 gCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGGAGCAAAAACGGCAGaa                                                           >  1:743259/1‑76 (MQ=255)
                   catcatCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGGAGAAgg                                                         >  1:1154658/1‑76 (MQ=255)
                     tcatcTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTa                                                       <  1:954631/76‑1 (MQ=255)
                     tcatcTGGAGCTGTGGCTGAACCCTCAACAATCACCATTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGAAggga                                                       >  2:1445746/1‑74 (MQ=255)
                      catcTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTaa                                                      <  1:426254/76‑1 (MQ=255)
                       atcTGGAGCTTTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAAt                                                     <  1:1625633/76‑1 (MQ=255)
                       atcTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAAt                                                     <  1:1286595/76‑1 (MQ=255)
                       atcTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAAt                                                     <  2:1818951/76‑1 (MQ=255)
                       atcTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAAt                                                     <  1:215233/76‑1 (MQ=255)
                       atcTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGAAGGGAAt                                                     >  1:797497/1‑76 (MQ=255)
                         cTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGAAGGGAATAg                                                   >  1:811437/1‑76 (MQ=255)
                          tGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGt                                                  <  2:178957/76‑1 (MQ=255)
                           ggAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTg                                                 <  1:150731/76‑1 (MQ=255)
                           ggAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCCGGAGGGAATAggg                                                 >  2:776006/1‑74 (MQ=255)
                            gAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGAAGGGAATAGTGc                                                >  2:1336420/1‑76 (MQ=255)
                              gCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCa                                              <  2:1601082/76‑1 (MQ=255)
                               cTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTTGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAg                                             <  1:117483/76‑1 (MQ=255)
                               cTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAg                                             <  1:302193/76‑1 (MQ=255)
                               cTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAg                                             <  2:737492/76‑1 (MQ=255)
                                tgtgGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGt                                            <  1:1205190/76‑1 (MQ=255)
                                 gtgGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCCAAAACGGCAGGAGGGAATAGTGGCAgcg                                           >  1:1781656/1‑74 (MQ=255)
                                 gtgGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGAAGGGAATAGTGCCACTg                                           >  2:438578/1‑76 (MQ=255)
                                    gCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTa                                        <  2:1748391/76‑1 (MQ=255)
                                     cTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTAt                                       <  1:127996/76‑1 (MQ=255)
                                      tGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGGAGGGAATAGGGCCCAGGCTATg                                      >  2:469069/1‑76 (MQ=255)
                                       gAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGAAGGGAATAGTGCCAGTGCTATGc                                     >  2:1161428/1‑76 (MQ=255)
                                       gAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAACGGCAGAAGGGAATAGTGCCAGTGCTATgg                                     >  2:1291092/1‑75 (MQ=255)
                                        aaCCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCa                                    <  1:774230/76‑1 (MQ=255)
                                        aaCCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGGAATAGTGCCAGTGCTATGCa                                    >  1:614586/1‑76 (MQ=255)
                                        aaCCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGAAGGGAATAGTGCCAGTGCTATGCa                                    >  2:1237838/1‑76 (MQ=255)
                                         aCCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGAAGGGAATAGTGGCAGTGCTATGCaa                                   >  2:459694/1‑76 (MQ=255)
                                          cccTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGGAATAGAGCCAGTGCTATGCAAc                                  >  2:116275/1‑76 (MQ=255)
                                           ccTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAAcc                                 <  1:1481837/76‑1 (MQ=255)
                                           ccTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAAcc                                 <  1:470171/76‑1 (MQ=255)
                                            cTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGAAGGGAATAGTGCCAGTGCTATGCAACCg                                >  1:797120/1‑76 (MQ=255)
                                             tCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGAAGGGAATAGTGCCAGTGCTATGCAACCgg                               >  1:1254338/1‑76 (MQ=255)
                                              caacaaTCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGc                              <  1:1449273/76‑1 (MQ=255)
                                               aacaaTCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAACGGCAGAAGGTAATAGAGCCAAAGATATGCAACCGGcc                             >  2:1596988/1‑75 (MQ=255)
                                               aacaaTCACCGTTTCTTCAGTGGCTGGGGCCGCAAAAGCGGCAGAAGGGAATAGTGCCCGGGCTATGCAACCGGCa                             >  2:1008919/1‑76 (MQ=255)
                                                      aCCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTg                      <  1:869140/76‑1 (MQ=255)
                                                       ccGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGa                     <  2:1044971/76‑1 (MQ=255)
                                                         gTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACg                   <  1:1269679/76‑1 (MQ=255)
                                                          tttCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGAAGGGAATAGTGCCAGTGCTATGCAACCCGCAAGGAGTGACgg                  >  2:624031/1‑76 (MQ=255)
                                                            tcttcAGTGGCTGGGGCAGCAAAAGCGGCAGAAGGGAATAGTGCCAGTGCTATGCAACCGGCCAGTAGTGACgggt                >  1:1732117/1‑74 (MQ=255)
                                                             cttcAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGttt               <  2:35873/76‑1 (MQ=255)
                                                               tcAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTgg             <  2:1330005/76‑1 (MQ=255)
                                                                cAGTGGCTGGGGCAGCAAAAGCGGCAGAAAGGAATAGTGGCAGGGCTATGCAACCGGCAAGGAGGGAGGGGTTggg            >  2:143968/1‑75 (MQ=255)
                                                                 aGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTg           <  1:548805/76‑1 (MQ=255)
                                                                 aGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTg           <  1:1792969/76‑1 (MQ=255)
                                                                  gTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGa          <  1:1268492/76‑1 (MQ=255)
                                                                   tGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGAt         >  1:1750247/1‑76 (MQ=255)
                                                                     gCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGATgg       >  1:1779900/1‑76 (MQ=255)
                                                                      cTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGATGGc      >  1:1023469/1‑76 (MQ=255)
                                                                       tGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGATGGCt     <  2:577869/76‑1 (MQ=255)
                                                                       tGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGATGGCt     <  2:577886/76‑1 (MQ=255)
                                                                        ggTTCAGCAAAAGCGGTAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGATGGCtt    <  2:538209/76‑1 (MQ=255)
                                                                        ggTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGATGGCtt    <  2:595768/76‑1 (MQ=255)
                                                                        ggTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGATGGCtt    <  2:1384708/76‑1 (MQ=255)
                                                                         gTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGATGGCTTg   <  2:584462/76‑1 (MQ=255)
                                                                         gTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGATGGCTTg   <  1:1223988/76‑1 (MQ=255)
                                                                          tGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGATGGCTTGa  <  2:1357368/76‑1 (MQ=255)
                                                                          |                                                                           
ACGCTGTAATCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGATGGCTTGA  >  REL606/1175936‑1176085

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