Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 3,095,121 G→T 100% C66* (TGC→TGA ygiQ ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,095,1210GT100.0% 350.4 / NA 102C66* (TGC→TGAygiQhypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (38/64);  total (38/64)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 102 total aligned reads displayed.

TGGGCGATGATCCCGACGCGAAAGCCCTGCGCTTCCAGCATACGACCGCAAATCGCCATCCCGAAGCTTGGGTGATCGACATACGCGTCGCCAGT  >  REL606/3095074‑3095168
                                               |                                               
tGGGCGATGATCCCGACGCGAAAGCCCTGCGCTTCCAGCATACGACCTCaa                                              <  1:3715072/51‑1 (MQ=255)
tGGGCGATGATCCCGACGCGAAAGCCCTGCGCTTCCAGCATACGACCTCaa                                              >  1:6443582/1‑51 (MQ=255)
 gggCGATGATCCCGACGCGAAAGCCCTGCGCTTCCAGCATACGACCTCaaa                                             <  1:5025619/51‑1 (MQ=255)
 gggCGATGATCCCGACGCGAAAGCCCTGCGCTTCCAGCATACGACCTCaaa                                             <  1:10255261/51‑1 (MQ=255)
   gCGATGATCCCGACGCGAAAGCCCTGCGCTTCCAGCATACGACCTCaaa                                             <  1:8792566/49‑1 (MQ=255)
    cGATGATCCCGACGCGAAAGCCCTGCGCTTCCAGCATACGACCTCAAATCg                                          <  1:8214288/51‑1 (MQ=255)
    cGATGATCCCGACGCGAAAGCCCTGCGCTTCCAGCATACGACCTCAAATCg                                          <  1:2340063/51‑1 (MQ=255)
    cGATGATCCCGACGCGAAAGCCCTGCGCTTCCAGCATACGACCTCAAATCg                                          <  1:4903687/51‑1 (MQ=255)
    cGATGATCCCGACGCGAAAGCCCTGCGCTTCCAGCATACGACCTCAAATCg                                          >  1:2184115/1‑51 (MQ=255)
    cGATGATCCCGACGCGAAAGCCCTGCGCTTCCAGCATACGACCTCAAATCg                                          >  1:2306712/1‑51 (MQ=255)
     gatgatCCCGACGCGAAAGCCCTGCGCTTCCAGCATACGACCTCAAATCGc                                         <  1:8042382/51‑1 (MQ=255)
     gatgatCCCGACGCGAAAGCCCTGCGCTTCCAGCATACGACCTCAAATCGc                                         >  1:559394/1‑51 (MQ=255)
     gatgatCCCGACGCGAAAGCCCTGCGCTTCCAGCATACGACCTCAAATCGc                                         >  1:9347932/1‑51 (MQ=255)
     gatgatCCCGACGCGAAAGCCCTGCGCTTCCAGCATACGACCTCAAATCGc                                         <  1:4103616/51‑1 (MQ=255)
      atgatCCCGACGCGAAAGCCCTGCGCTTCCAGCATACGACCTCAAATCGcc                                        >  1:10841958/1‑51 (MQ=255)
      atgatCCCGACGCGAAAGCCCTGCGCTTCCAGCATACGACCTCAAATCGcc                                        <  1:5972626/51‑1 (MQ=255)
       tgatCCCGACGCGAAAGCCCTGCGCTTCCAGCATACGACCTCAAAt                                            <  1:5904908/46‑1 (MQ=255)
         atCCCGACGCGAAAGCCCTGCGCTTCCAGCATACGACCTCAAATCGCCATc                                     <  1:5981655/51‑1 (MQ=255)
         atCCCGACGCGAAAGCCCTGCGCTTCCAGCATACGACCTCAAATCGCCATc                                     >  1:3020804/1‑51 (MQ=255)
         atCCCGACGCGAAAGCCCTGCGCTTCCAGCATACGACCTCAAATCGCCATc                                     >  1:5835920/1‑51 (MQ=255)
         atCCCGACGCGAAAGCAATGCGCTTCCAGCATACGACCTCAAATCGCCATc                                     <  1:6962029/51‑1 (MQ=255)
          tCCCGACGCGAAAGCCCTGCGCTTCCAGCATACGACCTCAAATCGCCATcc                                    <  1:5183838/51‑1 (MQ=255)
           cccGACGCGAAAGCCCTGCGCTTCCAGCATACGACCTCAAATCGCCATccc                                   >  1:10800079/1‑51 (MQ=255)
           cccGACGCGAAAGCCCTGCGCTTCCAGCATACGACCTCAAATCGCCATccc                                   >  1:2770113/1‑51 (MQ=255)
           cccGACGCGAAAGCCCTGCGCTTCCAGCATACGACCTCAAATCGCCATccc                                   >  1:2727801/1‑51 (MQ=255)
            ccGACGCGAAAGCCCTGCGCTTCCAGCATACGACCTCAAATCGCCATCCCg                                  <  1:7445835/51‑1 (MQ=255)
            ccGACGCGAAAGCCCTGCGCTTCCAGCATACGACCTCAAATCGCCATCCCg                                  <  1:4941114/51‑1 (MQ=255)
            ccGACGCGAAAGCCCTGCGCTTCCAGCATACGACCTCAAATCGCCATCCCg                                  >  1:6514792/1‑51 (MQ=255)
                cgcgAAAGCCCTGCGCTTCCAGCATACGACCTCAAATCGCCATCCCGAAGc                              <  1:4323716/51‑1 (MQ=255)
                 gcgAAAGCCCTGCGCTTCCAGCATACGACCTCAAATCGCCATCCCGAAGCt                             <  1:10138175/51‑1 (MQ=255)
                 gcgAAAGCCCTGCGCTTCCAGCATACGACCTCAAATCGCCATCCCGAAGCt                             >  1:7002732/1‑51 (MQ=255)
                 gcgAAAGCCCTGCGCTTCCAGCATACGACCTCAAATCGCCATCCCGAAGCt                             <  1:5096259/51‑1 (MQ=255)
                 gcgAAAGCCCTGCGCTTCCAGCATACGACCTCAAATCGCCATCCCGAAGCt                             <  1:1995917/51‑1 (MQ=255)
                  cgAAAGCCCTGCGCTTCCAGCATACGACCTCAAATCGCCATCCCGAAGCtt                            >  1:4668787/1‑51 (MQ=255)
                   gAAAGCCCTGCGCTTCCAGCATACGACCTCAAATCGCCATCCCGAAGCTTg                           >  1:9846171/1‑51 (MQ=255)
                   gAAAGCCCTGCGCTTCCAGCATACGACCTCAAATCGCCATCCCGAAGCTTg                           >  1:4202921/1‑51 (MQ=255)
                   gAAAGCCCTGCGCTTCCAGCATACGACCTCAAATCGCCATCCCGAAGCTTg                           <  1:4794306/51‑1 (MQ=255)
                   gAAAGCCCTGCGCTTCCAGCATACGACCTCAAATCGCCATCCCGAAGCTTg                           >  1:7496836/1‑51 (MQ=255)
                   gAAAGCCCTGCGCTTCCAGCATACGACCTCAAATCGCCATCCCGAAGCTTg                           >  1:1011903/1‑51 (MQ=255)
                    aaaGCCCTGCGCTTCCAGCATACGACCTCAAATCGCCATCCCGAAGCTTgg                          <  1:5953623/51‑1 (MQ=255)
                      aGCCCTGCGCTTCCAGCATACGACCTCAAATCGCCATCCCGAAGCTTGGGt                        <  1:312320/51‑1 (MQ=255)
                       gCCCTGCGCTTCCAGCATACGACCTCAAATCGCCATCCCGAAGCTTGGGTg                       <  1:9565962/51‑1 (MQ=255)
                        cccTGCGCTTCCAGCATACGACCTCAAATCGCCATCCCGAAGCTTGGGTGa                      >  1:8627165/1‑51 (MQ=255)
                        cccTGCGCTTCCAGCATACGACCTCAAATCGCCATCCCGAAGCTTGGGTGa                      <  1:3526182/51‑1 (MQ=255)
                        cccTGCGCTTCCAGCATACGACCTCAAATCGCCATCCCGAAGCTTGGGTGa                      <  1:3830019/51‑1 (MQ=255)
                         ccTGCGCTTCCAGCATACGACCTCAAATCGCCATCCCGAAGCTTGGGTGAt                     >  1:10516891/1‑51 (MQ=255)
                         ccTGCGCTTCCAGCATACGACCTCAAATCGCCATCCCGAAGCTTGGGTGAt                     >  1:10286492/1‑51 (MQ=255)
                          cTGCGCTTCCAGCATACGACCTCAAATCGCCATCCCGAAGCTTGGGTGATc                    <  1:10109707/51‑1 (MQ=255)
                          cTGCGCTTCCAGCATACGACCTCAAATCGCCATCCCGAAGCTTGGGTGATc                    <  1:150925/51‑1 (MQ=255)
                          cTGCGCTTCCAGCATACGACCTCAAATCGCCATCCCGAAGCTTGGGTGATc                    <  1:5189876/51‑1 (MQ=255)
                          cTGCGCTTCCAGCATACGACCTCAAATCGCCATCCCGAAGCTTGGGTGATc                    >  1:9410463/1‑51 (MQ=255)
                           tGCGCTTCCAGCATACGACCTCAAATCGCCATCCCGAAGCTTGGGTGATCg                   <  1:4984403/51‑1 (MQ=255)
                            gcgcTTCCAGCATACGACCTCAAATCGCCATCCCGAAGCTTGGGTGATCGa                  <  1:8131035/51‑1 (MQ=255)
                            gcgcTTCCAGCATACGACCTCAAATCGCCATCCCGAAGCTTGGGTGATCGa                  >  1:2337120/1‑51 (MQ=255)
                             cgcTTCCAGCATACGACCTCAAATCGCCATCCCGAAGCTTGGGTGATCGAc                 <  1:6646621/51‑1 (MQ=255)
                               cTTCCAGCATACGACCTCAAATCGCCATCCCGAAGCTTGGGTGATCGACAt               <  1:8996374/51‑1 (MQ=255)
                               cTTCCAGCATACGACCTCAAATCGCCATCCCGAAGCTTGGGTGATCGACAt               >  1:187800/1‑51 (MQ=255)
                               cTTCCAGCATACGACCTCAAATCGCCATCCCGAAGCTTGGGTGATCGACAt               <  1:10159460/51‑1 (MQ=255)
                                ttCCAGCATACGACCTCAAATCGCCATCCCGAAGCTTGGGTGATCGACATa              >  1:3740269/1‑51 (MQ=255)
                                ttCCAGCATACGACCTCAAATCGCCATCCCGAAGCTTGGGTGATCGACATa              <  1:7219265/51‑1 (MQ=255)
                                ttCCAGCATACGACCTCAAATCGCCATCCCGAAGCTTGGGTGATCGACATa              >  1:6947594/1‑51 (MQ=255)
                                ttCCAGCATACGACCTCAAATCGCCATCCCGAAGCTTGGGTGATCGACATa              <  1:5384922/51‑1 (MQ=255)
                                 tCCAGCATACGACCTCAAATCGCCATCCCGAAGCTTGGGTGATCGACATAc             <  1:2515036/51‑1 (MQ=255)
                                  ccggCATACGACCTCAAATCGCCATCCCGAAGCt                             <  1:6766474/31‑1 (MQ=25)
                                  ccAGCATACGACCTCAAATCGCCATCCCGAAGCTTGGGTGATCGACATAcg            <  1:824460/51‑1 (MQ=255)
                                  ccAGCATACGACCTCAAATCGCCATCCCGAAGCTTGGGTGATCGACATAcg            >  1:1221625/1‑51 (MQ=255)
                                   cggCATACGACCTCAAATCGCCATCCCGAAGCTTGGGTGATCGACa                <  1:2532320/44‑1 (MQ=39)
                                   cggCATACGACCTCAAATCGCCATCCCGAAGCTTGGGTGATCGACATAc             <  1:4469720/47‑1 (MQ=39)
                                   cAGCATACGACCTCAAATCGCCATCCCGAAGCtt                            <  1:7075898/34‑1 (MQ=37)
                                   cAGCATACGACCTCAAATCGCCATCCCGAAGCTTgg                          <  1:8011663/36‑1 (MQ=38)
                                   cAGCATACGACCTCAAATCGCCATCCCGAAGCTTg                           <  1:7486063/35‑1 (MQ=255)
                                   cAGCATACGACCTCAAATCGCCATCCCGAAGCTTGGGTGa                      <  1:7213621/40‑1 (MQ=39)
                                   cAGCATACGACCTCAAATCGCCATCCCGAAGCTTGGGTGa                      <  1:9343153/40‑1 (MQ=255)
                                   cAGCATACGACCTCAAATCGCCATCCCGAAGCTTGGGTGATc                    <  1:10557333/42‑1 (MQ=255)
                                   cAGCATACGACCTCAAATCGCCATCCCGAAGCTTGGGTGATCGACATAcgc           >  1:1595804/1‑51 (MQ=255)
                                    aGCATACGACCTCAAATCGCCATCCCGAAGCTTGGGTGATCGACATAtgcg          >  1:4474007/1‑51 (MQ=255)
                                    aGCATACGACCTCAAATCGCCATCCCGAAGCTTGGGTGATCGACATAcgcg          <  1:4408933/51‑1 (MQ=255)
                                     gCATACGACCTCAAATCGCCATCCCGAAGCTTGGGTGATCGACATACGCGt         >  1:3102222/1‑51 (MQ=255)
                                     gCATACGACCTCAAATCGCCATCCCGAAGCTTGGGTGATCGACATACGCGt         <  1:3060560/51‑1 (MQ=255)
                                     gCATACGACCTCAAATCGCCATCCCGAAGCTTGGGTGATCGACATACGCGt         <  1:5920797/51‑1 (MQ=255)
                                     gCATACGACCTCAAATCGCCATCCCGAAGCTTGGGTGATCGACATACGCGt         <  1:2704821/51‑1 (MQ=255)
                                     gCATACGACCTCAAATCGCCATCCCGAAGCTTGGGTGATCGACATACGCGt         <  1:5560342/51‑1 (MQ=255)
                                      cATACGACCTCAAATCGCCATCCCGAAGCTTGGGTGATCGACATACGCGTc        <  1:7526138/51‑1 (MQ=255)
                                      cATACGACCTCAAATCGCCATCCCGAAGCTTGGGTGATCGACATACGCGTc        <  1:3233936/51‑1 (MQ=255)
                                      cATACGACCTCAAATCGCCATCCCGAAGCTTGGGTGATCGACATACGCGTc        <  1:4090541/51‑1 (MQ=255)
                                      cATACGACCTCAAATCGCCATCCCGAAGCTTGGGTGATCGACATACGCGTc        <  1:1564561/51‑1 (MQ=255)
                                       aTACGACCTCAAATCGCCATCCCGAAGCTTGGGTGATCGACATACGCGTCg       <  1:4154768/51‑1 (MQ=255)
                                        tACGACCTCAAATCGCCATCCCGAAGCTTGGGTGATCGACATACGCGTc        >  1:7873355/1‑49 (MQ=255)
                                        tACGACCTCAAATCGCCATCCCGAAGCTTGGGTGATCGACATACGCGTCGc      <  1:5707930/51‑1 (MQ=255)
                                        tACGACCTCAAATCGCCATCCCGAAGCTTGGGTGATCGACATACGCGTCGc      >  1:806912/1‑51 (MQ=255)
                                         aCGACCTCAAATCGCCATCCCGAAGCTTGGGTGATCGACATACGCGTCGcc     >  1:2343863/1‑51 (MQ=255)
                                         aCGACCTCAAATCGCCATCCCGAAGCTTGGGTGATCGACATACGCGTCGcc     >  1:8939874/1‑51 (MQ=255)
                                         aCGACCTCAAATCGCCATCCCGAAGCTTGGGTGATCGACATACGCGTCGcc     >  1:946397/1‑51 (MQ=255)
                                          cGACCTCAAATCGCCATCCCGAAGCTTGGGTGATCGACATACGCGTCGCCa    <  1:10153443/51‑1 (MQ=255)
                                          cGACCTCAAATCGCCATCCCGAAGCTTGGGTGATCGACATACGCGTCGCCa    <  1:4636978/51‑1 (MQ=255)
                                          cGACCTCAAATCGCCATCCCGAAGCTTGGGTGATCGACATACGCGTCGCCa    <  1:1234509/51‑1 (MQ=255)
                                           gACCTCAAATCGCCATCCCGAAGCTTGGGTGATCGACATACGCGTCGCCAg   <  1:6675512/51‑1 (MQ=255)
                                           gACCTCAAATCGCCATCCCGAAGCTTGGGTGATCGACATACGCGTCGCCAg   >  1:4206742/1‑51 (MQ=255)
                                            aCCTCAAATCGCCATCCCGAAGCTTGGGTGATCGACATACGCGTCGCCAGt  >  1:3951551/1‑51 (MQ=255)
                                            aCCTCAAATCGCCATCCCGAAGCTTGGGTGATCGACATACGCGTCGCCAGt  >  1:6209732/1‑51 (MQ=255)
                                               |                                               
TGGGCGATGATCCCGACGCGAAAGCCCTGCGCTTCCAGCATACGACCGCAAATCGCCATCCCGAAGCTTGGGTGATCGACATACGCGTCGCCAGT  >  REL606/3095074‑3095168

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