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

GCTGTAATCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGATGGCTT  >  REL606/1175938‑1176083
                                                                        |                                                                         
gCTGTAATCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTgca                                                                        <  1:1504504/76‑1 (MQ=255)
  tGTAATCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGggagc                                                                      >  1:535100/1‑76 (MQ=255)
   gTAATCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGggagca                                                                     >  2:1519998/1‑76 (MQ=255)
   gTAATCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGgcagca                                                                     >  2:1115190/1‑76 (MQ=255)
    tAATCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCaa                                                                    <  1:1340763/76‑1 (MQ=255)
    tAATCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCaa                                                                    <  2:1676889/76‑1 (MQ=255)
     aaTCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCaaa                                                                   <  1:1526518/76‑1 (MQ=255)
      aTCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCaaaa                                                                  <  1:1502034/76‑1 (MQ=255)
      aTCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCaaaa                                                                  <  2:350535/76‑1 (MQ=255)
           tttCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGc                                                             <  2:385703/76‑1 (MQ=255)
           tttCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGc                                                             >  1:777359/1‑76 (MQ=255)
              cGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTgcag                                                                       <  1:2529/63‑1 (MQ=255)
              cGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTgcag                                                                       >  2:2529/1‑63 (MQ=255)
              cGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGGAGCAAAAGCGGgaga                                                          >  1:29952/1‑76 (MQ=255)
               gCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGaa                                                         >  2:1330246/1‑76 (MQ=255)
                 catcatCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCCGCAAAAGCGGCAGAAgg                                                       >  1:738786/1‑76 (MQ=255)
                   tcatcTGGAGCTGTGGCTGAGCCCTCAACAATCACCGTTTCTGCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTa                                                     <  2:341345/76‑1 (MQ=255)
                    catcTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTaa                                                    <  1:1672519/76‑1 (MQ=255)
                    catcTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGGAGAAGGGaa                                                    >  2:1375551/1‑76 (MQ=255)
                    catcTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGAAGGGaa                                                    >  2:324538/1‑76 (MQ=255)
                     atcTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAAt                                                   <  1:858556/76‑1 (MQ=255)
                     atcTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGGAGCAAAAACGGCAGGAGGGAAt                                                   >  1:593596/1‑76 (MQ=255)
                      tcTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATa                                                  <  2:535099/76‑1 (MQ=255)
                      tcTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATa                                                  <  2:894428/76‑1 (MQ=255)
                       cTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAg                                                 <  1:1706751/76‑1 (MQ=255)
                        tGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCGGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGt                                                <  2:1274659/76‑1 (MQ=255)
                        tGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGt                                                <  1:633425/76‑1 (MQ=255)
                        tGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGt                                                <  2:709830/76‑1 (MQ=255)
                         ggAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGGAGCAAAAACGGCAGAAGGGAATAGTg                                               >  1:1595850/1‑76 (MQ=255)
                          gAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGc                                              <  1:1298643/76‑1 (MQ=255)
                            gCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCa                                            <  2:571143/76‑1 (MQ=255)
                            gCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAACGGCAGAAGGGAATAGTGCCa                                            >  1:1314094/1‑76 (MQ=255)
                             cTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAg                                           <  2:821496/76‑1 (MQ=255)
                              tgtgGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGt                                          <  1:1186970/76‑1 (MQ=255)
                              tgtgGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGt                                          <  2:875334/76‑1 (MQ=255)
                              tgtgGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGt                                          <  1:769003/76‑1 (MQ=255)
                              tgtgGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGGAGGAAAAGAGGCAGAAGGGAAAAGTGCCAgg                                          >  1:347161/1‑75 (MQ=255)
                               gtgGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGGAGGGAATAGTGCCAGTg                                         >  2:275262/1‑76 (MQ=255)
                                 ggCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCt                                       <  2:1639152/76‑1 (MQ=255)
                                  gCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTa                                      <  2:1046857/76‑1 (MQ=255)
                                   cTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTAt                                     <  2:1650981/76‑1 (MQ=255)
                                   cTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTAt                                     <  1:406207/76‑1 (MQ=255)
                                    tGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATg                                    <  1:798037/76‑1 (MQ=255)
                                    tGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGGAGCAAAAGCGGGCGAAGGGAATAGTGCCAGTGCTTTg                                    >  2:298233/1‑76 (MQ=255)
                                    tGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCCAAAGCGGCCGAAGGTAATAGTGCCCGTGCTATg                                    >  2:1237853/1‑76 (MQ=255)
                                    tGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATg                                    >  1:1419292/1‑76 (MQ=255)
                                    tGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGAAGGGAATAGTGCCAGTGCTATg                                    >  2:132521/1‑76 (MQ=255)
                                     gAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGc                                   <  1:663568/76‑1 (MQ=255)
                                     gAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGGAATAGTGCCAGTGCTATGc                                   >  2:1647206/1‑76 (MQ=255)
                                     gAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGAAGGGAATAGTGCCAGTGCTATGc                                   >  2:1271662/1‑76 (MQ=255)
                                      aaCCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCa                                  <  2:146507/76‑1 (MQ=255)
                                      aaCCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCa                                  <  1:983280/76‑1 (MQ=255)
                                      aaCCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCCGGAGGGAATAGTGGCAGGGCTATGCa                                  >  1:893698/1‑76 (MQ=255)
                                      aaCCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAACGGCAGAAGGGAATAGTGCCAGTGCTATGCa                                  >  1:634286/1‑76 (MQ=255)
                                          cTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGAAGGGAATAGGGGCAGTGCTATGCAACCg                              >  2:1024408/1‑76 (MQ=255)
                                            caacaaTCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGGAATAGTGCCAGGGCTATGCAACCGGc                            >  1:770472/1‑76 (MQ=255)
                                             aacaaTCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCa                           <  1:1108664/76‑1 (MQ=255)
                                              acaaTCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGAAGGGAATAGTGCCAGTGCTATGCAACCGGCaa                          >  1:861826/1‑76 (MQ=255)
                                               caaTCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAag                         >  2:1641473/1‑76 (MQ=255)
                                               caaTCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGAAGGGAATAGGGCCAGTGCTATGCAACCGGCAag                         >  2:614804/1‑76 (MQ=255)
                                                aaTCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAagt                        <  2:775165/76‑1 (MQ=255)
                                                 aTCACCGTTCCTTCAGTGGCTGGTGCAGCAAAAGCGACAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAagta                       <  2:1489258/76‑1 (MQ=255)
                                                  tCACCGTTTCTTCAGTGGCTGGGGGCGAGAAAGCGGCCGAAGGGAATAGGGCCCGTGCTATGCAACCGGCAagtag                      >  2:313557/1‑76 (MQ=255)
                                                     ccGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGa                   <  1:444067/76‑1 (MQ=255)
                                                          tcttcAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGtt              <  2:29952/76‑1 (MQ=255)
                                                           cttcAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGttt             <  2:593595/76‑1 (MQ=255)
                                                            ttcAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTg            <  2:152906/76‑1 (MQ=255)
                                                              cAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGt          <  1:667932/76‑1 (MQ=255)
                                                               aGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTg         <  2:1204434/76‑1 (MQ=255)
                                                                  ggCTGGTGTAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGATg      >  1:924318/1‑76 (MQ=255)
                                                                    cTGGTGCAGCGAGAGCGGCAGAAGGTAATAGTGAAAGGGCTATGGGACAGGCGAGTAGGGACGGTTGGGTGATGGc    <  2:1376270/76‑1 (MQ=255)
                                                                     tGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGATGGCt   >  2:843922/1‑76 (MQ=255)
                                                                      ggTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGATGGCtt  <  2:379909/76‑1 (MQ=255)
                                                                      ggTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGATGGCtt  <  2:133184/76‑1 (MQ=255)
                                                                        |                                                                         
GCTGTAATCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGATGGCTT  >  REL606/1175938‑1176083

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