Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6061,176,0100TG28.8% 138.1 / 13.2 59A37A (GCA→GCCfhuEferric‑rhodotorulic acid outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base T (4/38);  new base G (17/0);  total (21/38)
Fisher's exact test for biased strand distribution p-value = 2.16e-11
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:930145/76‑1 (MQ=255)
  cGCTGTAATCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTgc                                                                           <  2:1073961/76‑1 (MQ=255)
  cGCTGTAATCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTgc                                                                           <  2:2598/76‑1 (MQ=255)
  cGCTATAATCATTTTCGCCATCATCTGTAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTgc                                                                           <  1:344111/76‑1 (MQ=255)
      ttAATCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCGCCGTTTCTTCAGTGGCTGGTgcagca                                                                       <  1:394993/75‑1 (MQ=255)
      gTAATCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTgcagca                                                                       <  1:2859/76‑1 (MQ=255)
       tAATCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCaa                                                                      <  2:1473608/76‑1 (MQ=255)
       tAATCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCaa                                                                      >  2:851645/1‑76 (MQ=255)
           cATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGc                                                                  <  2:464799/76‑1 (MQ=255)
           cATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGc                                                                  <  2:38362/76‑1 (MQ=255)
            aTTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGGAGCAAAAACg                                                                 >  2:1612795/1‑76 (MQ=255)
              tttCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGc                                                               <  1:1519554/76‑1 (MQ=255)
                 cGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGGAGCAAAAGAGGgaga                                                            >  2:1295649/1‑76 (MQ=255)
                 cGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGa                                                            >  2:1207466/1‑76 (MQ=255)
                 cGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGa                                                            >  1:366017/1‑76 (MQ=255)
                   ccATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAg                                                          <  2:1311014/76‑1 (MQ=255)
                   ccATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGAAg                                                          >  1:1640406/1‑76 (MQ=255)
                    catcatCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGGAGCAAAAGCGGCAGAAgg                                                         >  2:641043/1‑76 (MQ=255)
                     atcatcTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGt                                                        <  1:388537/76‑1 (MQ=255)
                     atcatcTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGGAGCAAAAGCGGCCGAAggg                                                        >  2:1129707/1‑75 (MQ=255)
                      tcatcTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTa                                                       <  1:1421623/76‑1 (MQ=255)
                      tcatcTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTa                                                       <  1:1146997/76‑1 (MQ=255)
                       catcTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGAAGGGaa                                                      >  2:124466/1‑76 (MQ=255)
                        atcTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGGACAAGCGGGAGAGGGGAAt                                                     >  1:857950/1‑76 (MQ=255)
                        atcTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGGAGCAAAAGCGGCAGAAGGGAAt                                                     >  1:799033/1‑76 (MQ=255)
                        atcTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGAAGGGAAt                                                     >  1:857958/1‑76 (MQ=255)
                         tcTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATa                                                    <  1:1139380/76‑1 (MQ=255)
                         tcTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGAAGGGAATa                                                    >  2:1486655/1‑76 (MQ=255)
                          cTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGAAGGGAATAg                                                   >  2:1475372/1‑76 (MQ=255)
                           tGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCCGCAAAAACGGCAGAAAGGAATaag                                                  >  1:658298/1‑74 (MQ=255)
                           tGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGAAGGGAATAGt                                                  >  1:16431/1‑76 (MQ=255)
                            ggAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGGCGCAAAAGCGGCCGAAGGGAATAATg                                                 >  1:742849/1‑76 (MQ=255)
                               gCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCa                                              <  1:1422606/76‑1 (MQ=255)
                                 tgtgGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGt                                            <  1:8232/76‑1 (MQ=255)
                                 tgtgGCTGAACCCTCAACAATCACCATTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGt                                            <  1:1017610/76‑1 (MQ=255)
                                     gCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTa                                        <  2:206007/76‑1 (MQ=255)
                                      cTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTAt                                       <  1:280984/76‑1 (MQ=255)
                                       tGAACCCTCAACAATCACCGTTTCTTCAGTGGTTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATg                                      <  1:1612777/76‑1 (MQ=255)
                                        gAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGc                                     <  1:1038314/76‑1 (MQ=255)
                                        gAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGc                                     <  1:816823/76‑1 (MQ=255)
                                        gAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGc                                     <  1:632557/76‑1 (MQ=255)
                                        gAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGAAGGGAATAGTGCCAGTGCTATGc                                     >  2:94170/1‑76 (MQ=255)
                                         aaCCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCa                                    <  1:1062365/76‑1 (MQ=255)
                                         aaCCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGAAGGGAATAGTGCCAGTGCTATGCa                                    >  2:80131/1‑76 (MQ=255)
                                          aCCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCaa                                   >  2:1521506/1‑76 (MQ=255)
                                            ccTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGAAGGGAATAGTGCCAGTGCTATGCAAcc                                 >  2:917975/1‑76 (MQ=255)
                                             cTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCg                                <  1:1431391/76‑1 (MQ=255)
                                             cTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCg                                <  1:298398/76‑1 (MQ=255)
                                               caacaaTCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGGAATAGTGCCAGTGCTATGCAACCGGc                              >  1:1149492/1‑76 (MQ=255)
                                               caacaaTCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGAAGGGAATAGTGCCCGTGCTATGCAACCGGc                              >  2:410459/1‑76 (MQ=255)
                                                 acaaTCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGGAGAAGGGAATAGTGCCAGTGCTATGCAACCGGCaa                            >  2:4320/1‑76 (MQ=255)
                                                  caaTCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAag                           <  2:934166/76‑1 (MQ=255)
                                                   aaTCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAagt                          <  2:1185009/76‑1 (MQ=255)
                                                             tcttcAGTGGCTGGGGCAGCAAAAGCGGCAGAAGGGAATAGTGCCAGTGCTATGCAAACGGCAAGTAGTGACgggt                >  2:565977/1‑74 (MQ=255)
                                                              cttcAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGttt               <  1:1559388/76‑1 (MQ=255)
                                                                tcAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTgg             <  1:1490722/76‑1 (MQ=255)
                                                                 cAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGt            <  1:510838/76‑1 (MQ=255)
                                                                   gTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGa          <  2:485050/76‑1 (MQ=255)
                                                                    tGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGAt         <  2:1612029/76‑1 (MQ=255)
                                                                    tGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGAt         >  2:776208/1‑76 (MQ=255)
                                                                     ggCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGATg        <  2:1434701/76‑1 (MQ=255)
                                                                       cTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGATGGc      >  2:9270/1‑76 (MQ=255)
                                                                         ggTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGATGGCtt    <  1:208431/76‑1 (MQ=255)
                                                                          gTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGATGGCTTg   <  1:1298932/76‑1 (MQ=255)
                                                                          gTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGATGGCTTg   <  2:253193/76‑1 (MQ=255)
                                                                          gTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGATGGCTTg   <  1:1523079/76‑1 (MQ=255)
                                                                           tGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGATGGCTTGa  >  1:944041/1‑76 (MQ=255)
                                                                           |                                                                           
TACGCTGTAATCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGATGGCTTGA  >  REL606/1175935‑1176085

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