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% 223.1 / NA 60intergenic (‑217/‑52)ybdK/insJ‑1gamma‑glutamyl:cysteine ligase/IS150 hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base A (25/35);  total (25/35)
*REL606588,4902.A100.0% 225.5 / NA 60intergenic (‑217/‑52)ybdK/insJ‑1gamma‑glutamyl:cysteine ligase/IS150 hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base A (25/35);  total (25/35)

AATTTGGCTACTTTTGCAATGTGACAAGTTATGGCAGTGCTGACTGGTGGCGAAGAATTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTC‑‑AAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATC  >  REL606/588392‑588590
                                                                                                   ||                                                                                                    
aaTTTGGCTACTTTTGCAATGTGACAAGTTATGGCAGTGCTGACTGGTGGCGAAGAATTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTC‑‑aa                                                                                                    >  2:51980/1‑101 (MQ=255)
aaTTTGGCTACTTTTGCAATGTGACAAGTTATGGCAGTGCTGACTGGTGGCGAAGAATTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTC‑‑aa                                                                                                    >  2:462077/1‑101 (MQ=255)
  tttGGCTACTTTTGCAATGTGACAAGTTATGGCAGTGCTGACTGGTGGCGAAGAATTTTCAGGATTGAGCCCTGCAGAAAAAAACGGGTTCAGCTTC‑‑aaat                                                                                                  <  2:248490/101‑3 (MQ=255)
   ttGGCTACTTTTGCAATGTGACAAGTTATGGCAGTGCTGACTGGTGGCGAAGAATTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTC‑‑aaaag                                                                                                 >  2:1707819/1‑98 (MQ=255)
   ttGGCTACTTTTGCAATGTGACAAGTTATGGCAGTGCTGACTGGTGGCGAAGAATTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTC‑‑aaaag                                                                                                 >  1:1207026/1‑98 (MQ=255)
   ttGGCTACTTTTGCAATGTGACAAGTTATGGCAGTGCTGACTGGTGGCGAAGAATTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTC‑‑aaaag                                                                                                 <  1:709283/101‑4 (MQ=255)
     ggCTACTTTTGCAATGTGACAAGTTATGGCAGTGCTGACTGGTGGCGAAGAATTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTC‑‑aaaagtt                                                                                               >  1:1144778/1‑96 (MQ=255)
      gCTACTTTTGCAATGTGACAAGTTATGGCAGTGCTGACTGGTGGCGAAGAATTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTC‑‑aaaagttg                                                                                              <  1:793799/101‑7 (MQ=255)
        tACTTTTGCAATGTGACAAGTTATGGCAGTGCTGACTGGTGGCGAAGAATTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTC‑‑aaaagttgta                                                                                            <  2:1299108/101‑9 (MQ=255)
        tACTTTTGCAATGTGACAAGTTATGGCAGTGCTGACTGGTGGCGAAGAATTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTC‑‑aaaagttgta                                                                                            <  1:1576758/101‑9 (MQ=255)
          cTTTTGCAATGTGACAAGTTATGGCAGTGCTGACTGGTGGCGAAGAATTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACt                                                                                            >  1:1071524/1‑101 (MQ=255)
          cTTTTGCAATGTGACAAGTTATGGCAGTGCTGACTGGTGGCGAAGAATTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACt                                                                                            >  2:386056/1‑101 (MQ=255)
           ttttGCAATGTGACAAGTTATGGCAGTGCTGACTGGTGGCGAAGAATTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTg                                                                                           <  1:141075/101‑1 (MQ=255)
           ttttGCAATGTGACAAGTTATGGCAGTGCTGACTGGTGGCGAAGAATTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTg                                                                                           <  2:456947/101‑1 (MQ=255)
               gCAATGTGACAAGTTATGGCAGTGCTGACTGGTGGCGAAGAATTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCAcc                                                                                       <  2:1143427/101‑1 (MQ=255)
                     tgACAAGTTATGGCAGTGCTGACTGGTGGCGAAGAATTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCAtttt                                                                                 <  2:223293/101‑1 (MQ=255)
                      gACAAGTTATGGCAGTGCTGACTGGTGGCGAAGAATTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTttg                                                                                >  2:45135/1‑101 (MQ=255)
                          aGTTATGGCAGTGCTGACTGGTGGCGAAGAATTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTgg                                                                            <  2:1325004/101‑1 (MQ=255)
                           gTTATGGCAGTGCTGACTGGTGGCGAAGAATTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGa                                                                           >  2:480836/1‑101 (MQ=255)
                            ttATGGCAGTGCTGACTGGTGGCGAAGAATTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGAc                                                                          <  2:282955/101‑1 (MQ=255)
                             tATGGCAGTGCTGACTGGTGGCGAAGAATTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACg                                                                         >  1:260713/1‑101 (MQ=255)
                              aTGGCAGTGCTGACTGGTGGCGAAGAATTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGa                                                                        >  2:1069464/1‑101 (MQ=255)
                               tGGCAGTGCTGACTGGTGGCGAAGAATTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGAt                                                                       >  1:796551/1‑101 (MQ=255)
                               tGGCAGTGCTGACTGGTGGCGAAGAATTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGAt                                                                       <  1:516656/101‑1 (MQ=255)
                                ggCAGTGCTGACTGGTGGCGAAGAATTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATg                                                                      <  2:50687/101‑1 (MQ=255)
                                    gTGCTGACTGGTGGCGAAGAATTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAAt                                                                  >  2:1608150/1‑101 (MQ=255)
                                    atGCTGACTGGTGGCGAAGAATTTTGACGATTGAGTCATGCAGAAAATAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAAt                                                                  <  1:1689166/100‑1 (MQ=38)
                                       cTGACTGGTGGCGAAGAATTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGa                                                               >  2:64776/1‑101 (MQ=255)
                                        tGACTGGTGGCGAAGAATTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGaa                                                              <  1:993538/101‑1 (MQ=255)
                                         gACTGGTGGCGAAGAATTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAAt                                                             <  1:494039/101‑1 (MQ=255)
                                         gACTGGTGGCGAAGAATTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAAt                                                             <  2:1071519/101‑1 (MQ=255)
                                         gACTGGTGGCGAAGAATTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAAt                                                             <  1:635589/101‑1 (MQ=255)
                                         gACTGGTGGCGAAGAATTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAAt                                                             <  2:593399/101‑1 (MQ=255)
                                          aCTGGTGGCGAAGAATTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATa                                                            >  1:755916/1‑101 (MQ=255)
                                          aCTGGTGGCGAAGAATTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATa                                                            <  1:110053/101‑1 (MQ=255)
                                               tggCGAAGAATTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGcccc                                                       <  1:1114907/101‑1 (MQ=255)
                                                ggCGAAGAATTTTGACGATTGAGTCATGCAGAAAAAAACGAGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCt                                                      <  1:4459/101‑1 (MQ=25)
                                                 gCGAAGAATTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTa                                                     >  2:1200857/1‑101 (MQ=255)
                                                 gCGAAGAATTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTa                                                     >  2:322804/1‑101 (MQ=255)
                                                    aagaaTTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAata                                                  <  1:265244/101‑1 (MQ=255)
                                                     agaaTTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAatat                                                 <  2:129529/101‑1 (MQ=255)
                                                     agaaTTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAatat                                                 <  1:1608156/101‑1 (MQ=255)
                                                     agaaTTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAatat                                                 <  1:706575/101‑1 (MQ=255)
                                                       aaTTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGt                                               <  2:260713/101‑1 (MQ=255)
                                                         ttttGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCa                                             <  1:1586390/101‑1 (MQ=255)
                                                          tttGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCaa                                            <  1:1707825/101‑1 (MQ=255)
                                                           ttGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCaaa                                           >  2:1491358/1‑101 (MQ=255)
                                                           ttGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCaaa                                           <  1:1561171/101‑1 (MQ=255)
                                                            tGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAg                                          <  1:537846/101‑1 (MQ=255)
                                                             gACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGc                                         >  2:752190/1‑101 (MQ=255)
                                                             caCGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTTCTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGc                                         <  1:1136370/100‑1 (MQ=21)
                                                                gATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCaa                                      >  2:526209/1‑101 (MQ=255)
                                                                gATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCaa                                      <  1:348501/101‑1 (MQ=255)
                                                                  ttGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCaaaa                                    >  1:1644016/1‑101 (MQ=255)
                                                                  ttGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCaaaa                                    <  2:340584/101‑1 (MQ=255)
                                                                  ttGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCaaaa                                    <  1:1123274/101‑1 (MQ=255)
                                                                   tGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAAt                                   <  2:750352/101‑1 (MQ=255)
                                                                      gTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATAcc                                <  2:1214981/101‑1 (MQ=255)
                                                                        cATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCt                              <  2:1004170/101‑1 (MQ=255)
                                                                         aTGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCtt                             >  1:449146/1‑101 (MQ=255)
                                                                         aTGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCtt                             >  1:1657168/1‑101 (MQ=255)
                                                                         aTGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCtt                             >  2:1141109/1‑101 (MQ=255)
                                                                           gCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCtttt                           >  2:544281/1‑101 (MQ=255)
                                                                             aGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGa                         <  1:111760/101‑1 (MQ=255)
                                                                                 aaaaaCGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAg                     >  2:495303/1‑101 (MQ=255)
                                                                                  aaaaCGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAgc                    >  2:105975/1‑101 (MQ=255)
                                                                                    aaCGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAgcgc                  <  1:1278006/101‑1 (MQ=255)
                                                                                       gggTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCtt               <  2:417472/101‑1 (MQ=255)
                                                                                          ttCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGaa            >  2:670744/1‑101 (MQ=255)
                                                                                          ttCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGaa            >  2:1247734/1‑101 (MQ=255)
                                                                                           tCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAg           >  2:762206/1‑101 (MQ=255)
                                                                                                  c‑‑AAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATc  >  1:211017/1‑101 (MQ=33)
                                                                                                   ||                                                                                                    
AATTTGGCTACTTTTGCAATGTGACAAGTTATGGCAGTGCTGACTGGTGGCGAAGAATTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTC‑‑AAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATC  >  REL606/588392‑588590

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