Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 588,492 +AA intergenic (‑219/‑50) 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% 118.7 / NA 42intergenic (‑217/‑52)ybdK/insJ‑1gamma‑glutamyl:cysteine ligase/IS150 hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base A (21/21);  total (21/21)
*REL606588,4902.A100.0% 118.6 / NA 42intergenic (‑217/‑52)ybdK/insJ‑1gamma‑glutamyl:cysteine ligase/IS150 hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base A (21/21);  total (21/21)

GGCGAAGAATTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTC‑‑AAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAA  >  REL606/588440‑588530
                                                   ||                                        
ggCGAAGAATTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTc                                            >  1:4044461/1‑51 (MQ=255)
 gCGAAGAATTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTC‑‑a                                         <  1:3462012/51‑1 (MQ=255)
  cGAAGAATTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTC‑‑aa                                        >  1:4903700/1‑51 (MQ=255)
   gaagaaTTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTC‑‑aaa                                       <  1:5798168/51‑2 (MQ=255)
   gaagaaTTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTC‑‑aaa                                       >  1:2751171/1‑50 (MQ=255)
     agaaTTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTC‑‑aaaag                                     <  1:458578/51‑4 (MQ=255)
      gaaTTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTC‑‑aaaagt                                    >  1:5335688/1‑47 (MQ=255)
      gaaTTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTC‑‑aaaag                                     >  1:4239772/1‑47 (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:4135153/1‑51 (MQ=38)
              aCGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGc                              <  1:5427084/51‑1 (MQ=38)
              aCGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGc                              >  1:946775/1‑51 (MQ=38)
               cGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCa                             >  1:1267536/1‑51 (MQ=38)
               cGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCa                             <  1:1931148/51‑1 (MQ=38)
                gATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCAc                            <  1:3056368/51‑1 (MQ=38)
                gATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCAc                            <  1:6610207/51‑1 (MQ=38)
                gATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCAc                            <  1:3203610/51‑1 (MQ=38)
                   tGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCa                         >  1:6270721/1‑51 (MQ=38)
                    gAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCAt                        <  1:6548699/51‑1 (MQ=38)
                    gAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCAt                        <  1:4506448/51‑1 (MQ=38)
                     aGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCAtt                       >  1:4246516/1‑51 (MQ=38)
                      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)
                          tGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTttgtt                  <  1:5084315/51‑1 (MQ=25)
                           gCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTttgttg                 >  1:1152440/1‑51 (MQ=25)
                           gCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTttgttg                 >  1:3824507/1‑51 (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:3295407/1‑51 (MQ=21)
                               aaaaaaaCGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACg             <  1:3814326/51‑1 (MQ=21)
                               aaaaaaaCGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACg             <  1:5490599/51‑1 (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)
                                     aCGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGaaa       >  1:640268/1‑51 (MQ=17)
                                      cGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAAt      >  1:120574/1‑51 (MQ=17)
                                      cGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAAt      >  1:3714578/1‑51 (MQ=17)
                                       gggTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATg     >  1:1139627/1‑51 (MQ=14)
                                        ggTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATgg    >  1:594425/1‑51 (MQ=14)
                                        ggTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATgg    <  1:2732124/51‑1 (MQ=14)
                                         gTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGa   >  1:3115700/1‑51 (MQ=14)
                                          ttCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGaa  >  1:2695409/1‑51 (MQ=14)
                                          ttCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGaa  >  1:4812376/1‑51 (MQ=14)
                                          ttCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGaa  >  1:7034265/1‑51 (MQ=14)
                                                   ||                                        
GGCGAAGAATTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTC‑‑AAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAA  >  REL606/588440‑588530

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