Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL60660,5740TG28.8% 127.4 / 10.0 52Y244S (TAC→TCC) imporganic solvent tolerance protein precursor
Reads supporting (aligned to +/- strand):  ref base T (2/35);  new base G (15/0);  total (17/35)
Fisher's exact test for biased strand distribution p-value = 3.03e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAA  >  REL606/60525‑60624
                                                 |                                                  
cgTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTa                                                   <  1:2394086/51‑1 (MQ=255)
 gTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAg                                                  <  1:2043964/51‑1 (MQ=255)
  tGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGa                                                 <  1:3372717/51‑1 (MQ=255)
   gATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGaa                                                >  1:4115778/1‑51 (MQ=255)
     tggtggCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACt                                              <  1:2862770/51‑1 (MQ=255)
     tggtggCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAACt                                              >  1:5262151/1‑51 (MQ=255)
      ggtggCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTc                                             <  1:3705265/51‑1 (MQ=255)
        tggCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCaa                                           <  1:4496855/51‑1 (MQ=255)
        tggCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAAATCaa                                           >  1:4619550/1‑51 (MQ=255)
         ggCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCaaa                                          <  1:3603023/51‑1 (MQ=255)
         ggCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCaaa                                          <  1:5417895/51‑1 (MQ=255)
         ggCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCaaa                                          <  1:2932757/51‑1 (MQ=255)
          gCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAACTCaaa                                          >  1:956808/1‑50 (MQ=255)
           cATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAagt                                        <  1:6332557/51‑1 (MQ=255)
           cATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAagt                                        <  1:5030349/51‑1 (MQ=255)
           cATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAagt                                        <  1:4923006/51‑1 (MQ=255)
            aTCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAagta                                       <  1:1509713/51‑1 (MQ=255)
            aTCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAagta                                       <  1:6776533/51‑1 (MQ=255)
             tCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAAATCaaa                                          >  1:3396212/1‑47 (MQ=255)
              ccATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAagtagt                                     <  1:6505784/51‑1 (MQ=255)
              ccATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAACTCaaa                                          >  1:5661693/1‑46 (MQ=255)
               cATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGtt                                    <  1:3635103/51‑1 (MQ=255)
               cATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGtt                                    <  1:2872718/51‑1 (MQ=255)
                  atTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTtggt                                 <  1:5562149/51‑1 (MQ=255)
                   tttGGCGCGATGTTCCAGTAATATGGCAGGGAGAAATCaaa                                          >  1:5585170/1‑41 (MQ=38)
                       gcgcGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTtggtggtgg                            <  1:2022789/51‑1 (MQ=255)
                        cgcgATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTTggtggtggt                           <  1:6357801/51‑1 (MQ=255)
                        cgcgATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTTggtggtggt                           <  1:3492088/51‑1 (MQ=255)
                        cgcgATGTTCCAGTAATATGGCAGGGAGAACTCAAAGTAGTtggggg                               >  1:2041163/1‑47 (MQ=255)
                         gcgATGTTCCAGTAATATGGCAGGGAGAACTCAAAGTAGTtggtggggg                            >  1:5555698/1‑49 (MQ=39)
                         gcgATGTTCCAGTAATATGGCAGGGAGAAATCAAAGTAGTtggtggggg                            >  1:6029401/1‑49 (MQ=38)
                         gcgATGTTCCAGTAATATGGCAGGGAGAAATCAAAATAGTtggtggggg                            >  1:5610398/1‑49 (MQ=255)
                          cgATGTTCCAGTAATATGGCAGGGAGAACTCaaa                                          >  1:4491505/1‑34 (MQ=255)
                          cgATGTTCCAGTAATATGGCAGGGAGAAATCaaa                                          >  1:965425/1‑34 (MQ=37)
                           gATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTa                        <  1:1378084/51‑1 (MQ=255)
                           gATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTa                        <  1:1391493/51‑1 (MQ=255)
                           gATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTa                        <  1:837210/51‑1 (MQ=255)
                            aTGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTAc                       <  1:6523665/51‑1 (MQ=255)
                             tGTTCCAGTAATATGGCAGGGAGAACTCaaa                                          >  1:6967430/1‑31 (MQ=255)
                              gTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACtt                     <  1:2219346/51‑1 (MQ=255)
                               ttCCAGTAATATGGCAGGGAGAACTCAAAGTAGTtggtggggg                            >  1:3841973/1‑43 (MQ=38)
                                 ccAGTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGc                  <  1:4120052/51‑1 (MQ=255)
                                 ccAGTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGc                  <  1:5679701/51‑1 (MQ=255)
                                       atatGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCgg            <  1:308700/51‑1 (MQ=255)
                                         atGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGa          <  1:2599683/51‑1 (MQ=255)
                                         atGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGa          <  1:7282140/51‑1 (MQ=255)
                                          tGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGAt         <  1:72198/51‑1 (MQ=255)
                                           ggCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATc        <  1:2191725/51‑1 (MQ=255)
                                             cAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATCaa      <  1:34187/51‑1 (MQ=255)
                                              aGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAg     <  1:4133036/51‑1 (MQ=255)
                                                gTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGaa   >  1:3891494/1‑51 (MQ=255)
                                                gTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGaa   >  1:6280803/1‑51 (MQ=255)
                                                 tAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGaaa  <  1:829727/51‑1 (MQ=255)
                                                 tAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGaaa  >  1:3197095/1‑51 (MQ=255)
                                                 |                                                  
CGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAA  >  REL606/60525‑60624

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: