Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6061,747,9030GT23.6% 238.7 / 53.3 89A8S (GCA→TCA) ydiLhypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (43/25);  new base T (21/0);  total (64/25)
Fisher's exact test for biased strand distribution p-value = 5.15e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 103 total aligned reads displayed.

CAGGCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGAAA  >  REL606/1747828‑1747975
                                                                           |                                                                        
cagGCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAgg                                                                          <  1:258195/76‑1 (MQ=255)
  ggCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCa                                                                        >  1:1719058/1‑76 (MQ=255)
   ncTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAc                                                                       >  1:792073/2‑76 (MQ=255)
   gCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAc                                                                       >  1:100502/1‑76 (MQ=255)
   gCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAc                                                                       >  1:58517/1‑76 (MQ=255)
   gCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAc                                                                       >  1:1253543/1‑76 (MQ=255)
   gCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAc                                                                       >  1:103564/1‑76 (MQ=255)
    cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcc                                                                      >  1:1409613/1‑75 (MQ=255)
    cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACt                                                                      >  1:1746599/1‑76 (MQ=255)
    cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACt                                                                      >  1:876259/1‑76 (MQ=255)
    cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACt                                                                      >  1:1244353/1‑76 (MQ=255)
    cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACt                                                                      <  1:1819386/76‑1 (MQ=255)
     tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcgc                                                                     >  1:165141/1‑74 (MQ=255)
     tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcca                                                                     >  1:905944/1‑74 (MQ=255)
     tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcca                                                                     >  1:1825343/1‑74 (MQ=255)
     tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTc                                                                     >  1:269354/1‑76 (MQ=255)
     tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTc                                                                     >  1:56018/1‑76 (MQ=255)
     tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTc                                                                     >  2:536519/1‑76 (MQ=255)
     tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTc                                                                     <  1:1124386/76‑1 (MQ=255)
      gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAccag                                                                    >  1:746113/1‑73 (MQ=255)
      gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACGcc                                                                    >  1:1017357/1‑76 (MQ=255)
       ataataTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCg                                                                   <  2:1802530/76‑1 (MQ=255)
       aTAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTTCAGTCACttag                                                                   >  1:291736/1‑73 (MQ=255)
       aTAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACttag                                                                   >  1:1182744/1‑73 (MQ=255)
       aTAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACGcca                                                                   >  1:289449/1‑75 (MQ=255)
       aTAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCg                                                                   <  2:1054623/76‑1 (MQ=255)
       aTAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCg                                                                   <  2:537502/76‑1 (MQ=255)
        tAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGc                                                                  >  2:777139/1‑76 (MQ=255)
         aaTGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACGccaat                                                                 >  1:9639/1‑73 (MQ=255)
         aaTGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACGccaat                                                                 >  1:1661527/1‑73 (MQ=255)
          aTGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCa                                                                >  2:113185/1‑76 (MQ=255)
           tGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCat                                                               >  1:1131948/1‑76 (MQ=255)
             ttatttttCTTTTATTGTTTATTTCATGCAGCCTTTACACGATGAACGCTTATGAACTCCAGACACTCCGGCatat                                                             >  2:1818866/1‑76 (MQ=255)
             ttatttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCatat                                                             >  2:719750/1‑76 (MQ=255)
              tatttatATTCTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATAtt                                                            >  2:824459/1‑76 (MQ=255)
               atttttATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATAttt                                                           >  2:1629945/1‑76 (MQ=255)
               atttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATAttt                                                           <  2:1468695/76‑1 (MQ=255)
                tttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACTTAGGCatatcgt                                                          >  1:1339162/1‑73 (MQ=255)
                tttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATAtttt                                                          >  2:1736431/1‑76 (MQ=255)
                tttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTACACGCACTCCGCAATAtctt                                                          >  2:1140189/1‑76 (MQ=255)
                 ttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATAttttt                                                         >  1:1209437/1‑76 (MQ=255)
                 ttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATAttttt                                                         >  2:1585333/1‑76 (MQ=255)
                  tataTTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTg                                                        <  2:1537042/76‑1 (MQ=255)
                    tatTTTATTGTTTNTTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGcc                                                      >  1:1229342/1‑76 (MQ=255)
                    tatTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGcc                                                      >  2:882821/1‑76 (MQ=255)
                    tatTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGcc                                                      >  2:280170/1‑76 (MQ=255)
                      ttttATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCAt                                                    >  2:846382/1‑76 (MQ=255)
                      ttttATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCAt                                                    >  1:1256862/1‑76 (MQ=255)
                      ttttATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCAt                                                    >  2:1362179/1‑76 (MQ=255)
                      ttttATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCAt                                                    >  2:606386/1‑76 (MQ=255)
                        ttATTTTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGa                                                  <  1:56751/76‑1 (MQ=255)
                        ttATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGa                                                  >  2:1770801/1‑76 (MQ=255)
                           ttGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTa                                               >  2:522396/1‑76 (MQ=255)
                            tGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTAt                                              >  2:1686961/1‑76 (MQ=255)
                            tGTTTATATCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTAt                                              >  1:660478/1‑76 (MQ=255)
                             gTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTAtt                                             <  1:43905/76‑1 (MQ=255)
                              tttATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCTATGACTATtg                                            <  1:1014772/76‑1 (MQ=255)
                              tttATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtg                                            >  2:1059162/1‑76 (MQ=255)
                              tttATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtg                                            <  1:1628952/76‑1 (MQ=255)
                               ttATTTCATGGAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtga                                           >  2:587484/1‑76 (MQ=255)
                                tATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtgat                                          <  1:576387/76‑1 (MQ=255)
                                 aTTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtgatg                                         >  2:404070/1‑76 (MQ=255)
                                 aTTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtgatg                                         <  2:289393/76‑1 (MQ=255)
                                 aTTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtgatg                                         >  2:200837/1‑76 (MQ=255)
                                 aTTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtgatg                                         >  1:572468/1‑76 (MQ=255)
                                   ttCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGaa                                       >  1:1795079/1‑76 (MQ=255)
                                     naTGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATg                                     >  1:886379/2‑76 (MQ=255)
                                     cATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATg                                     >  2:325091/1‑76 (MQ=255)
                                     cATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATg                                     >  1:525760/1‑76 (MQ=255)
                                       tGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATgcg                                   >  2:1046590/1‑76 (MQ=255)
                                       tGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATgcg                                   >  1:745967/1‑76 (MQ=255)
                                       cacAGCCATAACACGAAGAACGTTTATGAACTCGAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATgcg                                   <  1:1159973/74‑1 (MQ=255)
                                        gCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATgcgc                                  >  2:1610108/1‑76 (MQ=255)
                                        gCAGCCATTACACGATGAACGCTTATAAACTCCAGGCACCCCGCCATATTTTTGCCATGACTATTGATGAATgcgc                                  >  1:1341392/1‑76 (MQ=255)
                                         cAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGcgcg                                 >  2:1707943/1‑76 (MQ=255)
                                         cAGCCATTACACGATGAACGCTTAAGACCCCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGcgcg                                 <  1:1718471/76‑1 (MQ=255)
                                          aGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGa                                <  1:53853/76‑1 (MQ=255)
                                            ccATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGAcc                              >  1:962052/1‑76 (MQ=255)
                                            ccATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGAcc                              <  1:1688981/76‑1 (MQ=255)
                                            ccATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGAcc                              <  2:469218/76‑1 (MQ=255)
                                             cATTAGACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCt                             <  2:1800782/76‑1 (MQ=255)
                                              aTTACCCGAGGAGCGCTTGTGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGTGCGACCTg                            <  2:1560467/76‑1 (MQ=255)
                                               ttACACGATGAACGATTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTgg                           >  1:978722/1‑76 (MQ=255)
                                                  cacGATGAACGCTTATGACCCCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGAtt                        <  1:710478/76‑1 (MQ=255)
                                                       tGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCa                   >  1:98077/1‑76 (MQ=255)
                                                        gAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAg                  <  2:315982/76‑1 (MQ=255)
                                                               tATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTg           <  2:884838/76‑1 (MQ=255)
                                                                   accTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAg       <  1:401244/74‑1 (MQ=255)
                                                                   aaCTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAg       >  1:1210811/1‑76 (MQ=255)
                                                                   aaCTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAg       <  1:756307/76‑1 (MQ=255)
                                                                    aCTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGt      <  2:208388/76‑1 (MQ=255)
                                                                     cTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTg     >  2:1349461/1‑76 (MQ=255)
                                                                     cTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTg     <  2:684128/76‑1 (MQ=255)
                                                                     cTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGA‑‑GCGCGACCTGGATTGCCCAGACCGGTGACAGTGaa   >  2:1356185/1‑76 (MQ=255)
                                                                      tCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGa    <  2:131311/76‑1 (MQ=255)
                                                                      tCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGa    <  2:265355/76‑1 (MQ=255)
                                                                      cccAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGa    <  2:1311447/75‑1 (MQ=255)
                                                                      cccACGCCCTCCGCCAGATTTGTGCCATGACTATTGATGAGTGCGCGACCTGGATTGCCCAGGCCGGTGACAGTGa    <  2:918555/75‑1 (MQ=255)
                                                                       ccAGGCCCTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGaa   <  1:588727/76‑1 (MQ=255)
                                                                        cAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGaaa  <  1:209027/76‑1 (MQ=255)
                                                                           |                                                                        
CAGGCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGAAA  >  REL606/1747828‑1747975

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