Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL606273,5810TG37.5% 80.9 / 17.0 56G136G (GGT→GGGyagUconserved inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (15/20);  new base G (21/0);  total (36/20)
Fisher's exact test for biased strand distribution p-value = 5.54e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TGTTTGCTGTCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACACCACCTCTGTTTGATCTTCC  >  REL606/273508‑273655
                                                                         |                                                                          
tgtTTGCTGTCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGc                                                                          <  2:1026704/76‑1 (MQ=255)
   ttGCTGTCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGGttt                                                                       >  2:1371015/1‑76 (MQ=255)
     gCTGTCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAg                                                                     <  2:976508/76‑1 (MQ=255)
     gCTGTCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAg                                                                     <  2:723418/76‑1 (MQ=255)
       tGTCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTTGcc                                                                   >  1:1088028/1‑76 (MQ=255)
        gTCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGccc                                                                  >  1:101861/1‑76 (MQ=255)
        gTCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCCTTTGccc                                                                  >  2:503088/1‑76 (MQ=255)
         tCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCa                                                                 <  1:817692/76‑1 (MQ=255)
         tCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCa                                                                 <  2:213590/76‑1 (MQ=255)
          cGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCaa                                                                <  2:559780/76‑1 (MQ=255)
          cGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTTGCCCaa                                                                >  1:680860/1‑76 (MQ=255)
          cGGTTATTGTGTGGTCGCTGAAGTATTCCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTTTCCCaa                                                                >  2:435875/1‑76 (MQ=255)
             ttATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCCACtt                                                             >  2:1034743/1‑76 (MQ=255)
             ttATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTAGCCCCAcct                                                             >  2:950181/1‑74 (MQ=255)
              tATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTAGCCCAACttt                                                            >  1:195801/1‑76 (MQ=255)
               aTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCCTTAGCCCAACtttt                                                           >  1:1527400/1‑76 (MQ=255)
                ttGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTTGCCCCACttttt                                                          >  2:221177/1‑76 (MQ=255)
                 tgtgtgGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGGTTTTGCCCAACtttttt                                                         >  2:928432/1‑75 (MQ=255)
                 tgtgtgGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTAGCCCAACtttttt                                                         >  1:1117242/1‑75 (MQ=255)
                 tgtgtgGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTAGCCCAACcttttt                                                         >  1:1046374/1‑75 (MQ=255)
                  gtgtgGTCGCTGAAGTATTTCCAAAACTTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTTGCCCCACttttttt                                                        >  1:988300/1‑74 (MQ=255)
                    gtgGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTTGCCCCACttttttttt                                                      >  2:652757/1‑75 (MQ=255)
                    gtgGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGGTTTTGCCCACCTTTTTGttt                                                      >  2:598180/1‑75 (MQ=255)
                    gtgGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTTGCCCAACTTTTTGttt                                                      >  2:706329/1‑75 (MQ=255)
                     tgGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTTGCCCCACTTTTTGttta                                                     >  1:565777/1‑74 (MQ=255)
                      ggTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCCACTTTTTGTTTat                                                    >  2:1159821/1‑76 (MQ=255)
                      ggTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTAGCCCAACTTTTTGTTCat                                                    >  1:197529/1‑76 (MQ=255)
                         cGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTCGCCCAACTTTTTGTTCATATg                                                 >  2:1401129/1‑76 (MQ=255)
                          gCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTTGCCCCACTTTTTTTTCATATGa                                                >  1:333162/1‑76 (MQ=255)
                            tGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCCACCTTTTTTTCATATGAtt                                              >  1:449873/1‑76 (MQ=255)
                             gAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTTTTCATATGAttt                                             >  2:1354357/1‑76 (MQ=255)
                             gAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGAttt                                             <  2:146551/76‑1 (MQ=255)
                             gAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTAGCCCCACTTTTTTTTTATATGAttt                                             >  1:830495/1‑76 (MQ=255)
                              aaGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGAtttc                                            <  1:984852/76‑1 (MQ=255)
                               aGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGAtttct                                           <  2:870161/76‑1 (MQ=255)
                                gTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGAtttctt                                          <  1:1467874/76‑1 (MQ=255)
                                  aTTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTTTTCATATGAtttctttc                                        >  1:96969/1‑76 (MQ=255)
                                   tttCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTTGCCCCACTTTTTTTTTTTTTGAttttttttc                                       >  2:805108/1‑74 (MQ=255)
                                      ccAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCctc                                    >  1:1406107/1‑76 (MQ=255)
                                       cAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTTTTCATATGATTTCTTTTCctct                                   >  2:190044/1‑76 (MQ=255)
                                          aaaTTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTTGCCCAACTTTTTTTTTATATGATTTTTTTTCCtctcct                                >  1:883179/1‑74 (MQ=255)
                                          aaaTTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCAt                                <  2:924101/76‑1 (MQ=255)
                                          aaaTTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCAt                                <  2:1498905/76‑1 (MQ=255)
                                           aaTTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTTGCCCAACTTTTTTTTCATATGATTTTTTTTCCTCTCCTg                               >  2:950338/1‑76 (MQ=255)
                                               aaaCTCTGGCAGGGCTTACTGGCAGGGGCTTTTGCCCAACTTTTTGTTTATATGATTTTTTTTCCTCTCATggggc                           >  2:777201/1‑74 (MQ=255)
                                               aaaCTCTGGCAGGGCTTACTGGCAGGGGCTTTAGCCCAACTTTTTTTTCATATGATTTTTTTTCCTCTCATGGGAc                           >  1:1341734/1‑76 (MQ=255)
                                                aaCTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCCACTTTTTGTTCATATGATTTCTTTTCCTCTCATGGGACt                          >  1:1273043/1‑76 (MQ=255)
                                                  ctctGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTTTTTTCCTCTCATGGGACTGa                        >  1:209532/1‑76 (MQ=255)
                                                  ctctGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGa                        <  1:324241/76‑1 (MQ=255)
                                                  ctctGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGa                        <  2:827392/76‑1 (MQ=255)
                                                   tctGGCAGGGCTTACTGGCAGGGGCTTTAGCCCCACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGac                       >  1:1484944/1‑76 (MQ=255)
                                                    ctGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTTATATGATTTTTTTTCCTCTCATGGGACTGaca                      >  2:1034615/1‑76 (MQ=255)
                                                    ctGGCAGGGCTTACTGGCAGGGGCTTTTGCCCAACTTTTTGTTCATATGATTTTTTTCCCTCTCATGGGGCTGaca                      >  1:306754/1‑76 (MQ=255)
                                                      ggCAGGGCTTACTGGCAGGTGCTTTAGCCGAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGAcacc                    <  2:447161/76‑1 (MQ=255)
                                                      ggCAGGGCTTACTGGCAGGTGCTTTAGCCCAACCTTTTTTTCATATGAATTCTTTCCCCCCCCCGGGACCGAcacc                    >  1:987324/1‑76 (MQ=255)
                                                         aGGGCTTACTGGCAGGGGCTTTTGCCCCACTTTTTTTTCTTATTATTTCTTTTCCTCTCCTGGGACTGACaccacc                 >  1:375090/1‑76 (MQ=255)
                                                          gggCTTACTGGCAGGGGCTTTAGCCCAACTTTTTTTTCATATGATTTCTTTTCCTCTCCTGGGACTGACaccaccc                >  2:291387/1‑75 (MQ=255)
                                                             cTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACACCACCTCTg             >  2:359927/1‑76 (MQ=255)
                                                              ttACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACACCACCTCTGt            <  2:1149257/76‑1 (MQ=255)
                                                               tACTGGCAGGTGGTTTTGCCCAACTTTTTTTTTATATGATTTTTTTTCCTCTCATGGGACTGACACCACCTCTGtt           >  2:914178/1‑76 (MQ=255)
                                                                 cTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACACCACCTCTGTTTg         >  1:642710/1‑76 (MQ=255)
                                                                 cTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACACCACCTCTGTTTg         <  2:404007/76‑1 (MQ=255)
                                                                  tGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACACCACCTCTGTTTGa        <  2:90387/76‑1 (MQ=255)
                                                                   ggCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACACCACCTCTGTTTGAt       <  2:1038961/76‑1 (MQ=255)
                                                                     cAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACACCACCTCTGTTTGATCt     >  2:1047989/1‑76 (MQ=255)
                                                                      aGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACACCACCTCTGTTTGATCtt    <  2:389418/76‑1 (MQ=255)
                                                                       ggTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACACCACCTCTGTTTGATCTTc   >  2:491586/1‑76 (MQ=255)
                                                                        gTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACACCACCTCTGTTTGATCTTcc  <  2:446600/76‑1 (MQ=255)
                                                                         |                                                                          
TGTTTGCTGTCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACACCACCTCTGTTTGATCTTCC  >  REL606/273508‑273655

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