Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6061,747,9030GT21.9% 172.2 / 30.3 64A8S (GCA→TCA) ydiLhypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (32/18);  new base T (14/0);  total (46/18)
Fisher's exact test for biased strand distribution p-value = 6.56e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.24e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CAGGCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGAAAG  >  REL606/1747828‑1747976
                                                                           |                                                                         
cagGCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGAAGAACGCTTATGAACTCCAgg                                                                           >  2:510919/1‑76 (MQ=255)
   gCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAc                                                                        >  1:767802/1‑76 (MQ=255)
   gCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAc                                                                        >  1:402303/1‑76 (MQ=255)
   gCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCAc                                                                        >  2:28261/1‑76 (MQ=255)
    cTGATAATGTTGTTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcc                                                                       >  1:956050/1‑75 (MQ=255)
    cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcc                                                                       >  1:300574/1‑75 (MQ=255)
    cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcc                                                                       >  1:712907/1‑75 (MQ=255)
    cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAca                                                                       >  1:893537/1‑75 (MQ=255)
    cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACt                                                                       >  1:654245/1‑76 (MQ=255)
     tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACtt                                                                      >  1:413946/1‑75 (MQ=255)
     tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTc                                                                      >  2:768319/1‑76 (MQ=255)
     tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTc                                                                      >  1:103077/1‑76 (MQ=255)
     tGATAATGTTATTAATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGTACTCCAGTCACtt                                                                      >  1:78255/1‑75 (MQ=255)
      gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcctt                                                                     >  1:268863/1‑73 (MQ=255)
      gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACGcc                                                                     >  1:631187/1‑76 (MQ=255)
      gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTcc                                                                     <  1:1181581/76‑1 (MQ=255)
      gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAAATCCAGTCAcatc                                                                     >  1:348130/1‑73 (MQ=255)
       aTAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACttag                                                                    >  1:264763/1‑73 (MQ=255)
       aTAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGGACTCCg                                                                    <  1:40475/76‑1 (MQ=255)
        tAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACGccaa                                                                   >  1:735513/1‑74 (MQ=255)
        tAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGc                                                                   <  1:189981/76‑1 (MQ=255)
         aaTGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTACAGTCACGccaat                                                                  >  1:1032380/1‑73 (MQ=255)
            gTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCata                                                               <  1:290088/76‑1 (MQ=255)
            gTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCata                                                               >  2:341292/1‑76 (MQ=255)
             ttatttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCatat                                                              >  1:1168220/1‑76 (MQ=255)
             ttatttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCatat                                                              >  2:571981/1‑76 (MQ=255)
                tttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATAtttt                                                           >  1:539984/1‑76 (MQ=255)
                tttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATAtttt                                                           >  1:539973/1‑76 (MQ=255)
                 ttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATAttttt                                                          >  2:212094/1‑76 (MQ=255)
                     atTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCa                                                      >  1:328483/1‑76 (MQ=255)
                      ttttATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCAt                                                     >  2:569047/1‑76 (MQ=255)
                       tttATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATg                                                    >  2:896420/1‑76 (MQ=255)
                       tttATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTNTTGCCATg                                                    >  1:1059848/1‑76 (MQ=255)
                       tttATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTNTTGCCATg                                                    >  1:1059861/1‑76 (MQ=255)
                        ttATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGa                                                   <  1:529307/76‑1 (MQ=255)
                        ttATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGa                                                   >  1:539610/1‑76 (MQ=255)
                        ttATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGa                                                   >  2:555024/1‑76 (MQ=255)
                         tATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGAc                                                  >  2:1076595/1‑76 (MQ=255)
                           ttGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTa                                                >  2:969395/1‑76 (MQ=255)
                            tGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTAt                                               <  1:280191/76‑1 (MQ=255)
                             gTTTTTTTCATGCAGCCATTACACGATGAACGCTTATAACCTCCAGGCACTCCGCCATATTTTTGCCATGACTAtt                                              <  1:297432/76‑1 (MQ=255)
                             gTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTAtt                                              <  2:626286/76‑1 (MQ=255)
                             gTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTAtt                                              <  2:1119241/76‑1 (MQ=255)
                              tttATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtg                                             >  2:115771/1‑76 (MQ=255)
                               ttATTTCATGCAGCCATTGCATGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtga                                            <  1:1140302/76‑1 (MQ=255)
                               ttATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtga                                            >  2:1112779/1‑76 (MQ=255)
                               ttATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtga                                            >  1:442197/1‑76 (MQ=255)
                               ctgtttCATGCAGNCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtga                                            <  1:162995/73‑1 (MQ=255)
                               ctctttCATGCAGNCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtga                                            <  1:162985/73‑1 (MQ=255)
                                 aTTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATGTTATTGCCATGACTATtgatg                                          >  1:1133363/1‑76 (MQ=255)
                                 aTTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTCCCATGACTATtgatg                                          >  2:40576/1‑76 (MQ=255)
                                   ttCATCCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGaa                                        <  1:689689/76‑1 (MQ=255)
                                     cATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATg                                      >  1:572112/1‑76 (MQ=255)
                                      aTGCAGCCATTACACGATGACCGTTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATgc                                     <  1:743289/76‑1 (MQ=255)
                                      aTGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATgc                                     >  2:267114/1‑76 (MQ=255)
                                         cAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGcgcg                                  >  1:977134/1‑76 (MQ=255)
                                         cAGCCATTACACGATGAACGCTTATGAACTCCAGGCACCCCGCCATATTTTTGCCATGACTATTGATGAATGcgcg                                  >  1:1073678/1‑76 (MQ=255)
                                          aGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGa                                 <  2:59758/76‑1 (MQ=255)
                                            ccATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGAcc                               >  1:491444/1‑76 (MQ=255)
                                             cATTACACGACGAACGCTTATGAACTCCAGGCACTCCGCCGTATTTTTGCCATGACTATTGATGAATGCGCGGCCt                              <  2:597609/76‑1 (MQ=255)
                                                tACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGa                           <  2:54208/76‑1 (MQ=255)
                                                 acacGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGAt                          >  1:975320/1‑76 (MQ=255)
                                                      aTGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGccc                     <  2:210222/76‑1 (MQ=255)
                                                           agCTTATGTACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATCGCCCAGAcc                >  1:39759/2‑76 (MQ=255)
                                                             cTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCgg              >  1:381148/1‑76 (MQ=255)
                                                             cTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCgg              >  2:98205/1‑76 (MQ=255)
                                                              ttATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGt             <  2:1168220/76‑1 (MQ=255)
                                                                        cAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGaaa   <  2:312208/76‑1 (MQ=255)
                                                                         aGGCACTCNGCCATATTTTNGCCATGACTATTGATGAATGCGNGACCTGGATTGCCCAGACCGGTGACAGTGAAAg  <  2:1007050/76‑1 (MQ=255)
                                                                           |                                                                         
CAGGCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGAAAG  >  REL606/1747828‑1747976

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