Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 1,610,902 T→C intergenic (‑42/+322) cspI ← / ← ydfP cold shock protein/conserved protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6061,610,9020TC96.0% 406.8 / 0.1 99intergenic (‑42/+322)cspI/ydfPcold shock protein/conserved protein
Reads supporting (aligned to +/- strand):  ref base T (3/1);  new base C (36/59);  total (39/60)
Fisher's exact test for biased strand distribution p-value = 2.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTTCACTAAACCAGTCATTTTGTTAGACATCATTATTACCTTTTGAAGAAGTTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAAT  >  REL606/1610830‑1610974
                                                                        |                                                                        
aTTTCACTAAACCAGTCATTTTGTTAGACATCATTATTACCTTTTGAAGAAGTTAGCCCTTGGGCAGAATGGCCCg                                                                       <  1:1349590/76‑1 (MQ=255)
 tttCACTAAACCAGTCATTTTGTTAGACATCATTATTACCTTTTGAAGAAGTTAGCCCTTGGGCAGAATGGCCCGa                                                                      >  2:1714382/1‑76 (MQ=255)
   tCACTAAACCAGTCATTTTGTTAGACATCATTATTACCTTTTGAAGAAGTTAGCCCTTGGGCAGAATGGCCCGaaa                                                                    <  2:1621754/76‑1 (MQ=255)
   tCACTAAACCAGTCATTTTGTTAGACATCATTATTACCTTTTGAAGAAGTTAGCCCTTGGGCAGAATGGCCCGaaa                                                                    >  2:257714/1‑76 (MQ=255)
    cACTAAACCAGTCATTTTGTTAGACATCATTATTACCTTTTGAAGAAGTTAGCCCTTGGGCAGAATGGCCCGaaaa                                                                   <  1:2112612/76‑1 (MQ=255)
    cACTAAACCAGTCATTTTGTTAGACATCATTATTACCTTTTGAAGAAGTTAGCCCTTGGGCAGAATGGCCCGaaaa                                                                   <  2:891446/76‑1 (MQ=255)
    cACTAAACCAGTCATTTTGTTAGACATCATTATTACCTTTTGAAGAAGTTAGCCCTTGGGCAGAATGGCCCGaaaa                                                                   >  1:708618/1‑76 (MQ=255)
     aCTAAACCAGTCATTTTGTTAGACATCATTATTACCTTTTGAAGAAGTTAGCCCTTGGGCAGAATGGCCCGaaaaa                                                                  <  1:1654262/76‑1 (MQ=255)
       tAAACCAGTCATTTTGTTAGACATCATTATTACCTTTTGAAGAAGTTAGCCCTTGGGCAGAATGGCCCGaaaaaaa                                                                >  1:1031359/1‑76 (MQ=255)
        aaaCCAGTCATTTTGTTAGACATCATTATTACCTTTTGAAGAAGTTAGCTCTTGGGCAGAATGGCCCGaaaaaaaa                                                               >  1:738136/1‑76 (MQ=255)
        aaaCCAGTCATTTTGTTAGACATCATTATTACCTTTTGAAGAAGTTAGCCCTTGGGCAGAATGGCCCGaaaaaaaa                                                               >  2:285789/1‑76 (MQ=255)
        aaaCCAGTCATTTTGTTAGACATCATTATTACCTTTTGAAGAAGTTAGCCCTTGGGCAGAATGGCCCGaaaaaaaa                                                               >  2:1279119/1‑76 (MQ=255)
         aaCCAGTCATTTTGTTAGACATCATTATTACCTTTTGAAGAAGTTAGCCCTTGGGCAGAATGGCCCGAAAAAAAat                                                              <  1:1010469/76‑1 (MQ=255)
         aaCCAGTCATTTTGTTAGACATCATTATTACCTTTTGAAGAAGTTAGCCCTTGGGCAGAATGGCCCGAAAAAAAat                                                              <  1:488368/76‑1 (MQ=255)
          aCCAGTCATTTTGTTAGACATCATTATTACCTTTTGAAGAAGTTAGCCCTTGGGCAGAATGGCCCGAAAAAAAata                                                             >  2:1514923/1‑76 (MQ=255)
          aCCAGTCATTTTGTTAGACATCATTATTACCTTTTGAAGAAGTTAGCCCTTGGGCAGAATGGCCCGAAAAAAAata                                                             <  1:1835157/76‑1 (MQ=255)
           ccAGTCATTTTGTTAGACATCATTATTACCTTTTGAAGAAGTTAGCCCTTGGGCAGAATGGCCCGAAAAAAAatat                                                            <  1:54488/76‑1 (MQ=255)
            cAGTCATTTTGTTAGACATCATTATTACCTTTTGAAGAAGTTAGCCCTTGGGCAGAATGGCCCGAAAAAAAATATc                                                           <  1:1534549/76‑1 (MQ=255)
            cAGTCATTTTGTTAGACATCATTATTACCTTTTGAAGAAGTTAGCCCTTGGGCAGAATGGCCCGAAAAAAAATATc                                                           <  1:916460/76‑1 (MQ=255)
            cAGTCATTTTGTTAGACATCATTATTACCTTTTGAAGAAGTTAGCCCTTGGGCAGAATGGCCCGAAAAAAAATATc                                                           <  1:424121/76‑1 (MQ=255)
            cAGTCATTTTGTTAGACATCATTATTACCTTTTGAAGAAGTTAGCCCTTGGGCAGAATGGCCCGAAAAAAAATATc                                                           >  1:1816237/1‑76 (MQ=255)
             aGTCATTTTGTTAGACATCATTATTACCTTTTGAAGAAGTTAGCCCTTGGGCAGAATGGCCCGAAAAAAAATATCa                                                          <  1:1601206/76‑1 (MQ=255)
             aGTCATTTTGTTAGACATCATTATTACCTTTTGAAGAAGTTAGCCCTTGGGCAGAATGGCCCGAAAAAAAATATCa                                                          >  1:1961733/1‑76 (MQ=255)
             aGTCATTTTGTTAGACATCATTATTACCTTTTGAAGAAGTTAGCCCTTGGGCAGAATGGCCCGAAAAAAAATATCa                                                          <  2:1518460/76‑1 (MQ=255)
                cATTTTGTTAGACATCATTATTACCTTTTGAAGAAGTTAGCCCTTGGGCAGAATGGCCCGAAAAAAAATATCagag                                                       <  2:881357/76‑1 (MQ=255)
                cATTTTGTTAGACATCATTATTACCTTTTGAAGAAGTTAGCCCTTGGGCAGAATGGCCCGAAAAAAAATATCagag                                                       >  1:1914078/1‑76 (MQ=255)
                cATTTTGTTAGACATCATTATTACCTTTTGAAGAAGTTAGCCCTTGGGCAGAATGGCCCGAAAAAAAATATCagag                                                       <  1:1114406/76‑1 (MQ=255)
                cATTTTGTTAGACATCATTATTACCTTTTGAAGAAGTTAGCCCTTGGGCAGAATGGCCCGAAAAAAAATATCagag                                                       <  2:1510676/76‑1 (MQ=255)
                cATTTTGTTAGACATCATTATTACCTTTTGAAGAAGTTAGCCCTTGCGCAGAATGGCCCGAAAAAAAATATCagag                                                       >  2:1256727/1‑76 (MQ=255)
                  ttttGTTAGACATCATTATTACCTTTTGAAGAAGTTAGCCCTTGGGCAGAATGGCCCGAAAAAAAATATCagagag                                                     <  2:1193051/76‑1 (MQ=255)
                   tttGTTAGACATCATTATTACCTTTTGAAGAAGTTAGCCCTTGGGCAGAATGGCCCGAAAAAAAATATCAgagaga                                                    <  2:1350995/76‑1 (MQ=255)
                   tttGTTAGACATCATTATTACCTTTTGAAGAAGTTAGCCCTTGGGCAGAATGGCCCGAAAAAAAATATCAgagaga                                                    >  1:545297/1‑76 (MQ=255)
                   tttGTTAGACATCATTATTACCTTTTGAAGAAGTTAGCCCTTGGGCAGAATGGCCCGAAAAAAAATATCAgagaga                                                    <  1:957301/76‑1 (MQ=255)
                   tttGTTAGACATCATTATTACCTTTTGAAGAAGTTAGCCCTTGGGCAGAATGGCCCGAAAAAAAATATCAgagaga                                                    >  1:1088560/1‑76 (MQ=255)
                   tttGTTAGACATCATTATTACCTTTTGAAGAAGTTAGCCCTTGGGCAGAATGGCCCGAAAAAAAATATCAgagaga                                                    >  1:1416700/1‑76 (MQ=255)
                   tttGTTAGACATCATTATTACCTTTTGAAGAAGTTAGCCCTTGGGCAGAATGGCCCGAAAAAAAATATCAgagaga                                                    >  1:1898543/1‑76 (MQ=255)
                    ttGTTAGACATCATTATTACCTTTTGAAGAAGTTAGCCCTTGGGCAGAATGGCCCGAAAAAAAATATCAGAGAGaa                                                   >  2:751862/1‑76 (MQ=255)
                    ttGTTAGACATCATTATTACCTTTTGAAGAAGTTAGCCCTTGGGCAGAATGGCCCGAAAAAAAATATCAGAGAGaa                                                   <  2:391867/76‑1 (MQ=255)
                       ttAGACATCATTATTACCTTTTGAAGAAGTTAGCCCTTGGGCAGAATGGCCCGAAAAAAAATATCAGAGAGaaaaa                                                >  1:1049753/1‑76 (MQ=255)
                        tAGACATCATTATTACCTTTTGAAGAAGTTAGCCCTTGGGCAGAATGGCCCGAAAAAAAATATCAGAGAGAAAAAc                                               >  1:1115397/1‑76 (MQ=255)
                        tAGACATCATTATTACCTTTTGAAGAAGTTAGCCCTTGGGCAGAATGGCCCGAAAAAAAATATCAGAGAGAAAAAc                                               >  1:480209/1‑76 (MQ=255)
                         aGACATCATTATTACCTTTTGAAGAAGTTAGCCCTTGGGCAGAATGGCCCGAAAAAAAATATCAGAGAGAAAAAcc                                              <  1:740022/76‑1 (MQ=255)
                          gACATCATTATTACCTTTTGAAGAAGTTAGCCCTTGGGCAGAATGGCCCGAAAAAAAATATCAGAGAGAAAAACca                                             >  2:1214602/1‑76 (MQ=255)
                           aCATCATTATTACCTTTTGAAGAAGTTAGCCCTTGGGCAGAATGGCCCGAAAAAAAATATCAGAGAGAAAAACcaa                                            <  1:1227129/76‑1 (MQ=255)
                           aCATCATTATTACCTTTTGAAGAAGTTAGCCCTTGGGCAGAATGGCCCGAAAAAAAATATCAGAGAGAAAAACcaa                                            <  2:1955662/76‑1 (MQ=255)
                           aCATCATTATTACCTTTTGAAGAAGTTAGCCCTTGGGCAGAATGGCCCGAAAAAAAATATCAGAGAGAAAAACcaa                                            <  1:1227119/76‑1 (MQ=255)
                            catcatTATTACCTTTTGAAGAAGTTAGCCCTTGGGCAGAATGGCCCGAAAAAAAATATCAGAGAGAAAAACcaac                                           <  1:720648/76‑1 (MQ=255)
                               catTATTACCTTTTGAAGAAGTTAGCCCTTGGGCAGAATGGCCCGAAAAAAAATATCAGAGAGAAAAACCAACAAg                                        <  1:702553/76‑1 (MQ=255)
                                  tattaCCTTTTGAAGAAGTTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGCGAAAAACCAACAAGGaa                                     >  2:417609/1‑76 (MQ=255)
                                  tattaCCTTTTGAAGAAGTTAGCCCTTGGGCAGAATGGCCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGaa                                     <  1:773434/76‑1 (MQ=255)
                                  tattaCCTTTTGAAGAAGTTAGCCCTTGGGCAGAATGGCCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGaa                                     <  1:1751360/76‑1 (MQ=255)
                                   attaCCTTTTGAAGAAGTTAGCCCTTGGGCAGAATGGCCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGaaa                                    <  2:187538/76‑1 (MQ=255)
                                    ttaCCTTTTGAAGAAGTTAGCCCTTGGGCAGAATGGCCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAAt                                   <  2:1410382/76‑1 (MQ=255)
                                      aCCTTTTGAAGAAGTTAGCCCTTGGGCAGAATGGCCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAAtct                                 >  1:1191984/1‑76 (MQ=255)
                                      aCCTTTTGAAGAAGTTAGCCCTTGGGCAGAATGGCCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAAtct                                 >  2:25654/1‑76 (MQ=255)
                                       ccTTTTGAAGAAGTTAGCCCTTGGGCAGAATGGCCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAAtctc                                <  1:140920/76‑1 (MQ=255)
                                        cTTTTGAAGAAGTTAGCCCTTGGGCAGAATGGCCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCa                               <  2:1168578/76‑1 (MQ=255)
                                        cTTTTGAAGAAGTTAGCCCTTGGGCAGAATGGCCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCa                               <  1:616749/76‑1 (MQ=255)
                                        cTTTTGAAGAAGTTAGCCCTTGGGCAGAATGGCCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCa                               <  2:1727342/76‑1 (MQ=255)
                                         ttttGAAGAAGTTAGCCCTTGGGCAGAATGGCCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCaa                              <  2:766259/76‑1 (MQ=255)
                                         ttttGAAGAAGTTAGCCCTTGGGCAGAATGGCCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCaa                              <  2:1400797/76‑1 (MQ=255)
                                         ttttGAAGAAGTTAGCCCTTGGGCAGAATGGCCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCaa                              <  2:741631/76‑1 (MQ=255)
                                          tttGAAGAAGTTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAAGCTCAag                             >  2:714048/1‑76 (MQ=255)
                                          tttGAAGAAGTTAGCCCTTGGGCAGAATGGCCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAag                             >  1:1286531/1‑76 (MQ=255)
                                              aagaagTTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGt                         <  1:1701272/76‑1 (MQ=255)
                                              aagaagTTAGCCCTTGGGCAGAATGGCCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGt                         >  1:749001/1‑76 (MQ=255)
                                              aagaagTTAGCCCTTGGGCAGAATGGCCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGt                         <  2:1775780/76‑1 (MQ=255)
                                               agaagTTAGCCCTTGGGCAGAATGGCCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTa                        <  1:574371/76‑1 (MQ=255)
                                               agaagTTAGCCCTTGGGCAGAATGGCCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTa                        <  2:1554925/76‑1 (MQ=255)
                                                 aagTTAGCCCTTGGGCAGAATGGCCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACa                      <  1:497848/76‑1 (MQ=255)
                                                 aagTTAGCCCTTGGGCAGAATGGCCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACa                      <  1:1180041/76‑1 (MQ=255)
                                                 aagTTAGCCCTTGGGCAGAATGGCCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACa                      <  1:208460/76‑1 (MQ=255)
                                                 aagTTAGCCCTTGGGCAGAATGGCCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACa                      <  1:1432604/76‑1 (MQ=255)
                                                 aagTTAGCCCTTGGGCAGAATGGCCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACa                      <  2:919661/76‑1 (MQ=255)
                                                   gTTAGCCCTTGGGCAGAATGGCCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGaaaaaa                    >  1:822545/1‑76 (MQ=255)
                                                   gTTAGCCCTTGGGCAGAATGGCCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACaaa                    <  1:958462/76‑1 (MQ=255)
                                                    ttAGCCCTTGGGCAGAATGGCCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAaat                   >  2:1004796/1‑76 (MQ=255)
                                                    ttAGCCCTTGGGCAGAATGGCCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAaat                   <  1:1872323/76‑1 (MQ=255)
                                                     tAGCCCTTGGGCAGAATGGCCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAaata                  >  2:605043/1‑76 (MQ=255)
                                                      aGCCCTTGGGCAGAATGGCCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAaataa                 >  2:1775363/1‑76 (MQ=255)
                                                      aGCCCTTGGGCAGAATGGCCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAaataa                 >  2:1344609/1‑76 (MQ=255)
                                                       gCCCTTGGGCAGAATGGCCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAaataat                <  2:907364/76‑1 (MQ=255)
                                                       gCCCTTGGGCAGAATGGCCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAaataat                <  2:326697/76‑1 (MQ=255)
                                                       gCCCTTGGGCAGAATGGCCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAaataat                >  1:1342747/1‑76 (MQ=255)
                                                       gCCCTTGGGCAGAATGGCCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAaataat                >  1:874193/1‑76 (MQ=255)
                                                       gCCCTTGGGCAGAATGGCCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAaataat                >  1:713979/1‑76 (MQ=255)
                                                        cccTTGGGCAGAATGGCCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAataata               <  2:484130/76‑1 (MQ=255)
                                                         ccTTGGGCAGAATGGCCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAAtaataa              <  2:2059054/76‑1 (MQ=255)
                                                            tGGGCAGAATGGCCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAAt           >  2:926070/1‑76 (MQ=255)
                                                            tGGGCAGAATGGCCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAAt           >  1:944218/1‑76 (MQ=255)
                                                            tGGGCAGAATGGCCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAAt           <  2:1232633/76‑1 (MQ=255)
                                                             gggCAGAATGGCCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAAtt          <  1:2002189/76‑1 (MQ=255)
                                                              agCAGAATGGTCCGAAAAAAAAATCCGAGAGAAAAACCAACAAGGAAATCTCAAGGGGTACAAATAATAAAATtat        >  2:351751/2‑76 (MQ=255)
                                                               gCAGAATGGCCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATtat        <  2:471058/76‑1 (MQ=255)
                                                                 aGAATGGCCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATaa      >  1:326453/1‑76 (MQ=255)
                                                                 aGAATGGCCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATaa      <  2:770872/76‑1 (MQ=255)
                                                                 aGAATGGCCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATaa      <  1:1210085/76‑1 (MQ=255)
                                                                   aaTGGCCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACa    <  1:2098984/76‑1 (MQ=255)
                                                                     tGGCCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAAt  <  2:1901687/76‑1 (MQ=255)
                                                                        |                                                                        
ATTTCACTAAACCAGTCATTTTGTTAGACATCATTATTACCTTTTGAAGAAGTTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAAT  >  REL606/1610830‑1610974

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