Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 588,490:1 +AA 100% intergenic (‑217/‑52) ybdK ← / → insJ‑1 gamma‑glutamyl:cysteine ligase/IS150 hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL606588,4901.A100.0% 196.5 / NA 69intergenic (‑217/‑52)ybdK/insJ‑1gamma‑glutamyl:cysteine ligase/IS150 hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base A (32/37);  total (32/37)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
*REL606588,4902.A100.0% 195.5 / NA 69intergenic (‑217/‑52)ybdK/insJ‑1gamma‑glutamyl:cysteine ligase/IS150 hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base A (32/37);  total (32/37)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GGCGAAGAATTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTC‑‑AAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAAT  >  REL606/588440‑588531
                                                   ||                                         
ggCGAAGAATTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTc                                             >  1:8568759/1‑51 (MQ=255)
ggCGAAGAATTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTc                                             >  1:4044461/1‑51 (MQ=255)
ggCGAAGAATTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTc                                             <  1:7593767/51‑1 (MQ=255)
 gCGAAGAATTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTC‑‑a                                          <  1:3462012/51‑1 (MQ=255)
  cGAAGAATTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTC‑‑aa                                         >  1:4903700/1‑51 (MQ=255)
   gaagaaTTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTC‑‑aaa                                        >  1:2751171/1‑50 (MQ=255)
   gaagaaTTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTC‑‑aaa                                        <  1:5798168/51‑2 (MQ=255)
     agaaTTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTC‑‑aaaag                                      >  1:7852914/1‑48 (MQ=255)
     agaaTTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTC‑‑aaaag                                      <  1:458578/51‑4 (MQ=255)
     agaaTTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTC‑‑aaaa                                       <  1:9457220/50‑3 (MQ=255)
      gaaTTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTC‑‑aaaagt                                     >  1:5335688/1‑47 (MQ=255)
      gaaTTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTC‑‑aaaag                                      >  1:4239772/1‑47 (MQ=255)
       aaTTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTC‑‑aaaagtt                                    >  1:9041050/1‑46 (MQ=255)
       aaTTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTC‑‑aaaagtt                                    <  1:7917491/51‑6 (MQ=255)
       aaTTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTC‑‑aaaagtt                                    <  1:5764141/51‑6 (MQ=255)
            tGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACt                                 <  1:6482558/51‑1 (MQ=38)
             gACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTg                                <  1:5593376/51‑1 (MQ=38)
              aCGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGc                               <  1:3108/51‑1 (MQ=38)
              aCGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGc                               <  1:5427084/51‑1 (MQ=38)
              aCGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGc                               >  1:946775/1‑51 (MQ=38)
              aCGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGc                               >  1:4135153/1‑51 (MQ=38)
               cGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCa                              <  1:1931148/51‑1 (MQ=38)
               cGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCa                              <  1:10609400/51‑1 (MQ=38)
               cGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCa                              >  1:8810932/1‑51 (MQ=38)
               cGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCa                              >  1:1267536/1‑51 (MQ=38)
                gATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCAc                             <  1:3056368/51‑1 (MQ=38)
                gATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCAc                             <  1:6610207/51‑1 (MQ=38)
                gATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCAc                             <  1:3203610/51‑1 (MQ=38)
                  ttGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCAccc                           >  1:9954815/1‑51 (MQ=38)
                  ttGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCAccc                           <  1:10870336/51‑1 (MQ=38)
                   tGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCa                          >  1:6270721/1‑51 (MQ=38)
                    gAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCAt                         >  1:10406524/1‑51 (MQ=38)
                    gAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCAt                         <  1:6548699/51‑1 (MQ=38)
                    gAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCAt                         <  1:10640069/51‑1 (MQ=38)
                    gAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCAt                         <  1:4506448/51‑1 (MQ=38)
                     aGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCAtt                        <  1:9893766/51‑1 (MQ=38)
                     aGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCAtt                        >  1:4246516/1‑51 (MQ=38)
                      gTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCAttt                       <  1:9170642/51‑1 (MQ=37)
                      gTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCAttt                       <  1:3616509/51‑1 (MQ=37)
                       tCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCAtttt                      >  1:3348436/1‑51 (MQ=37)
                       tCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCAtttt                      <  1:4161240/51‑1 (MQ=37)
                       tCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCAttt                       <  1:2913422/50‑1 (MQ=37)
                         aTGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCAcc                            >  1:10486858/1‑43 (MQ=25)
                         aTGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTttgt                    >  1:7551526/1‑51 (MQ=25)
                          tGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTttgtt                   <  1:10269523/51‑1 (MQ=25)
                          tGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTttgtt                   <  1:5084315/51‑1 (MQ=25)
                           gCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTttgttg                  <  1:10622870/51‑1 (MQ=25)
                           gCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTttgttg                  >  1:3824507/1‑51 (MQ=25)
                           gCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTttgttg                  >  1:1152440/1‑51 (MQ=25)
                           gCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTttgttg                  <  1:7742483/51‑1 (MQ=25)
                             aGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGa                <  1:5053658/51‑1 (MQ=21)
                             aGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGa                >  1:6127233/1‑51 (MQ=21)
                             aGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGa                <  1:2340265/51‑1 (MQ=21)
                              gAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGAc               <  1:8281706/51‑1 (MQ=21)
                              gAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGAc               >  1:9845396/1‑51 (MQ=21)
                              gAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGAc               >  1:9532317/1‑51 (MQ=21)
                              gAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGAc               >  1:3295407/1‑51 (MQ=21)
                               aaaaaaaCGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACg              <  1:5490599/51‑1 (MQ=21)
                               aaaaaaaCGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACg              >  1:10938373/1‑51 (MQ=21)
                               aaaaaaaCGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACg              <  1:3814326/51‑1 (MQ=21)
                               aaaaaaaCGGGTTCAGATTCAAAAGTTGTACTGCTCCCATTTTGTTGGACg              <  1:8945973/51‑1 (MQ=18)
                                aaaaaaCTGGTTCAGCTTCAAAAGTTGTACTGCGCCCATTTTGTTGGACGa             <  1:10249759/51‑1 (MQ=18)
                                aaaaaaCGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTttgtt                   >  1:9521821/1‑45 (MQ=21)
                                aaaaaaCGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGa             <  1:4143697/51‑1 (MQ=18)
                                  aaaaCGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATg           >  1:5386907/1‑51 (MQ=18)
                                   aaaCGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGa          <  1:507884/51‑1 (MQ=18)
                                    aaCGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGaa         >  1:783393/1‑51 (MQ=17)
                                    aaCGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGaa         <  1:8868692/51‑1 (MQ=17)
                                    aaCGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGaa         >  1:8430086/1‑51 (MQ=17)
                                     aCGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGaaa        >  1:640268/1‑51 (MQ=17)
                                      cGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAAt       >  1:10687610/1‑51 (MQ=17)
                                      cGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAAt       >  1:3714578/1‑51 (MQ=17)
                                      cGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAAt       <  1:8220364/51‑1 (MQ=17)
                                      cGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAAt       >  1:120574/1‑51 (MQ=17)
                                       gggTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATg      >  1:1139627/1‑51 (MQ=14)
                                        ggTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATgg     <  1:2732124/51‑1 (MQ=14)
                                        ggTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATgg     >  1:594425/1‑51 (MQ=14)
                                         gTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGa    >  1:3115700/1‑51 (MQ=14)
                                          ttCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGaa   >  1:7034265/1‑51 (MQ=14)
                                          ttCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGaa   >  1:2695409/1‑51 (MQ=14)
                                          ttCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGaa   <  1:8445595/51‑1 (MQ=14)
                                          ttCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGaa   >  1:4812376/1‑51 (MQ=14)
                                           tCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAAt  <  1:10754860/51‑1 (MQ=11)
                                           tCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAAt  <  1:9846471/51‑1 (MQ=11)
                                                   ||                                         
GGCGAAGAATTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTC‑‑AAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAAT  >  REL606/588440‑588531

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